Free Access
Issue
A&A
Volume 640, August 2020
Article Number A123
Number of page(s) 15
Section Stellar atmospheres
DOI https://doi.org/10.1051/0004-6361/202038511
Published online 26 August 2020

© ESO 2020

1 Introduction

The precise determination of the atmospheric parameters and chemical composition of stars plays a crucial role in many astrophysical fields. In particular, this is of primary importance for exoplanetary studies to fully understand the main observational correlations between the properties of exoplanets and the characteristics of their host star; these correlations include the giant planet-metallicity relation and the trends observed between the condensation temperature (Tc) and abundances ratios [X/Fe] (see e.g. Nissen 2015; Brewer et al. 2016; Adibekyan 2019 and references therein).

High-resolution stellar spectroscopy is one of the most powerful tools to fully characterise a star. This technique allows us to determine with great precision the physical properties of stars, for example effective temperatures (Teff), and the chemical abundances of various atomic and molecular species (see Jofré et al. 2019 for a complete review). However, in some cases the spectroscopic analysis of stars is not trivial and the results can be affected by multiple issues. In particular, young stars (< 200 Myr) have higher activity levels, which can alter the structure of the atmosphere, complicate the derivation of the stellar parameters and consequently of the metallicity, and can hamper the detectability of planets (e.g. Carleo et al. 2018, 2020).

The presence of active chromospheres or intense photospheric magnetic fields (Folsom et al. 2016) may alter the spectral line formation. Recently, Baratella et al. (2020) show that the apparent sub-solar metallicity observed for the young stars in the solar neighbourhood may be related to an over-estimation of the microturbulence velocity (ξ) parameter. This is a free fictitious parameter representing small-scale motions of matter in the photospheric layers of the star and it is introduced in 1D spectroscopic analysis to account for the difference between the observed and predicted equivalent widths (EWs), when models account only for thermal and damping broadening. Weaker lines are less affected by this parameter, which is calculated by forcing lines, usually of iron (Fe), to give the same abundance. However, higher values of ξ lead to systematic under-estimation of the elemental abundances, so that the stars belonging to young associations and open clusters (OCs; <200 Myr) could appear more metal-poor than what it is predicted by Galactic chemical evolution models (D’Orazi et al. 2011; Biazzo et al. 2011a; Spina et al. 2017). For example, James et al. (2006) analysed young stars observed in three star-forming regions (τ ~ 10 Myr) and reported extremely high ξ values for pre-main-sequence stars of up to 2.5 km s−1, which is expected for more evolved stars. Similar results have been also reported by Santos et al. (2008). Viana Almeida et al. (2009) analysed stars in 11 young associations (ages less than 100 Myr) and again found ξ values up to 2.6 km s−1. Moreover, they reported a small trend of increasing ξ at decreasing Teff. Several other authors confirmed the anomalous values of ξ found for young stars, but they reported close-to-solar values of the metallicity. For instance, D’Orazi et al. (2009) and Biazzo et al. (2011a,b) reported values up to ~2 km s−1 for stars belonging to the Orion complex. However, for some stars they also reported values of metallicity varying from − 0.15 to 0.01 dex, despite the large scatter in the ξ parameter measurements (e.g. in Biazzo et al. 2011b a star in λ Orionis has ξ = 2.1 km s−1, but [Fe/H] = 0.00 dex).

Recently, Yana Galarza et al. (2019) show that the observed EWs of iron lines vary with the activity phase of the young solar analogue HD 59967 (age ~400 Myr). In particular, they find that the line strength increases when the star is more active, producing variations of ξ and iron abundance along with the stellar cycle. Moreover, these authors demonstrate that such variations of EWs depend on the optical depth of line formation and, marginally, on the Landè gL factor, which measures the sensitivity of a spectral line to magnetic fields. Spina et al. (2020) conducted the same study as Yana Galarza et al. (2019) to a sample of 211 solar-analogue stars observed with High Accuracy Radial velocity Planet Searcher (HARPS) and find similar results. Yana Galarza et al. (2019), Baratella et al. (2020), and Spina et al. (2020) demonstrate how iron lines forming in the upper layers of the photosphere of young stars can yield higher abundances due to the possible influence of the more intense chromospheric or photospheric magnetic fields.

Higher levels of stellar activity can also affect the abundances of some elements when they are derived using high-excitation potential lines (see Schuler et al. 2015 and references therein). It has been observed that young (<200 Myr) and cool (Teff ≲ 5400 K) stars display anomalous abundances of oxygen O I triplet (χ = 9.15 eV) and sulfur (S I line at 6053 Å with χ = 7.87 eV) (Schuler et al. 2004; Teske et al. 2013; Ramírez et al. 2013). In particular, the abundances increase at decreasing Teff, reaching values of 0.8–1.0 dex over solar for the coolest stars (Teff ~ 4700 K). Similarly, for the same kind of stars and in the same Teff regime, differences between the neutral and ionised species of the same element of the order of +0.8 dex have been observed for Fe and Ti (Schuler et al. 2006; D’Orazi & Randich 2009). Such differences can produce unreliable results, in particular for the derivation of log g, which should be decreased in order to satisfy the ionisation equilibrium. These effects may be caused either by non-local thermodynamic equilibrium (NLTE) departures, for which the high-energy levels are not correctly modelled, or by the presence of unidentified blends (see also Tsantaki et al. 2019), or a combination of both. Aleo et al. (2017) argue that the large differences between Fe I and Fe I may be related to blending of Fe II lines that become more severe at decreasing Teff. These results were corroborated by Takeda & Honda (2020), who also concluded that the O I overabundance obtained from the oxygen triplet by Schuler et al. (2006) might be due to the different Teff scale and to over-estimation of the strength of the lines in coolest stars. Even though we are aware of the issues related to the spectroscopic analysis of young and cool stars and we are starting to shed light on the topic, we still lack a definitive solution. However, we can overcome these problems with strategic choices of the line list to use in the analysis, for example with a refined selection based on the EWs (Spina et al. 2020) or using new approaches (Baratella et al. 2020).

Thanks to the advent of large spectroscopic surveys, such as Gaia-ESO Survey (Gilmore et al. 2012) or GALactic Archaeology with HERMES (GALAH; De Silva et al. 2015), the number of stellar spectra has increased enormously. Along with the increasing number of available spectra, the need arose to assess the precision and accuracy of spectroscopic analysis techniques (Jofré et al. 2019). Until recently, the study of stellar spectra mainly involved the analysis of data in the optical band, covering the wavelength range from ~4000 to ~7000 Å. However, the advent of high-resolution near-infrared (NIR) spectroscopy allowed us to extend the analysis of stellar spectra at longer wavelengths as well, and to test the validity of optical and NIR analysis techniques (Marfil et al. 2020). This is particularly important in the study of young and intermediate-age stars, for which stellar activity and other effects can alter the derivation of atmospheric parameters and, specifically, the chemical composition.

For five years, the Global Architecture of Planetary Systems (GAPS) project (Covino et al. 2013; Poretti et al. 2016)searched for planets through a radial velocity (RV) technique with High Accuracy Radial velocity Planet Searcher forthe Northern emisphere (HARPS-N; Cosentino et al. 2014) at Telescopio Nazionale Galileo (TNG, Roque de los Muchachos, La Palma) around different types of stars, including the characterisation of selected planet-host stars. Recently, a new phase of the project started with the aim of exploiting the full capabilities of the GIARPS mode (Claudi et al. 2017). This means that we can study and fully characterise planetary systems by analysing GIAno-B (Oliva et al. 2006) and haRPS stellar spectra acquired simultaneously. In this context, theGAPS Young Objects (GAPS-YO) project (Carleo et al. 2020) aims to monitor and study young (< 100 Myr) and intermediate-age (<700 Myr) stars to search for and characterise hot and warm planets down to sub-Neptune mass in formation or at an early stage of their evolution.

In this first paper of a series, we present the results of spectral characterisation in terms of astrophysical parameters and elemental abundances of stars observed within the GAPS-YO project. Our analysis includes the Sun, two RV standard stars (HD 3765 and HD 15922), and seven more stars members of intermediate-age stellar clusters and moving groups. In Sect. 2 we present the data we analysed and we report some information on the selected stars. Our analysis is separated between the optical and NIR spectral ranges. In particular, in Sect. 3 we describe the new method applied to derive atmospheric parameters and elemental abundances from optical spectra. These parameters were used to derive abundances of neutral C, Na, Mg, Al, Si, Ca, Ti, Fe, and Ni from NIR spectral lines (Sect. 4). In Sect. 5 we present the resulting chemical abundances of various atomic species and discuss the scientific implications. In Sect. 6 we present our conclusions.

2 Sample selection and spectroscopic data

In this work, we analysed high-resolution spectra of seven young and intermediate-age stars observed in the GAPS-YO project. We selected spectra with high signal-to-noise ratios (S/N), low rotational velocities (v sin i < 15 km s−1) to avoid line blending, and spectral type F-G-K. We excluded from the analysis stars with spectral types later than K to avoid problems with the molecular bands. The selected targets are as follows: TYC 1991-1235-1, HIP 61205, TYC 1989-0049-1, and TYC 1989-147-1, which belong to the Coma Berenices OC (Mermilliod et al. 2008), with an age of ~600 Myr; HD 167389 and HD 59747, which are part of the Ursa Major moving group (Montes et al. 2001), with an age of ~500 Myr; and HD 70573 of Hercules Lyra moving group (López-Santiago et al. 2006), with an age of ~200 Myr. We also analysed, for validation, the spectra of two old stars observed as RV standard stars, HD 3765 and HD 159222. We reported some information on the selected targets in Table 1.

The spectra were acquired with HARPS-N and GIANO-B spectrographs placed at the 3.6 m INAF-Telescopio Nazionale Galileo (TNG) in La Palma. The HARPS-N spectrograph is the northern counterpart of HARPS at the La Silla Observatory (Chile), mounted at the Nasmyth-B focus of the TNG. With a resolving power R ~115 000 and large wavelength coverage in the optical range (0.38–0.69 μm), it allows us to obtain very precise (less than 1 m s−1) RV measurements, thanks to an accurate control system that minimises pressure and temperature variations and prevents spectral drifts due to environmental conditions. The GIANO-B spectrograph is a high-resolution (R ~ 45 000–50 000) NIR spectrograph covering the wavelength range from 0.95 to 2.45 μm that is placed at Nasmyth-B focus of the TNG. The configuration of the two spectrographs allows us to observe the stars simultaneously in the optical and NIR wavelengths in the GIARPS mode. We analysed GIARPS spectra of the Sun, HD 3765, HD 159222, TYC 1991-1235-1, HIP 61205, and HD 167389. We also analysed the optical spectra of the remaining four additional targets.

HARPS-N data are reduced with the standard Data Reduction Software (DRS). Since the spectra were collected by the GAPS-YO collaboration to obtain time series for RV monitoring, the available HARPS-N data for each target were then combined to obtain a co-added spectrum with SN > 100 (Malavolta et al. 2016). The NIR data reduction was performed with the pipeline GOFIO (Rainer et al. 2018; Harutyunyan et al. 2018), while the telluric correction was performed following the method described in Carleo et al. (2016). We verified that co-adding HARPS-N spectra did not introduce any systematic errors, thanks to the high stability of the instrument over several months. This was not the case for GIANO-B spectra, for which we decided to consider the highest S/N observation (SN > 70) for each star. Since we performed a differential analysis with respect to the Sun, we derived our solar abundance scale by analysing the HARPS-N and GIANO-B spectra of Ganymede, which have a SN = 145 at 607 nm and SN = 180 at 1500 nm, respectively.

Table 1

Some basic information of both standard stars and members in young associations analysed in this work.

3 Optical analysis

For the analysis of the HARPS-N optical spectra, we employed the same approach as in Baratella et al. (2020), which exploits the use of Ti lines to derive the atmospheric parameters. On average, Ti lines form deeper in the photospherethan Fe lines, so they are less affected by the chromosphere, which is more active in young stars. In this way, wecan overcome the issues affecting the analysis of young stars, which have been already presented in Sect. 1. Briefly, the new spectroscopic method is based on the use of Ti and Fe lines to derive Teff by imposing the excitation equilibrium and the use of Ti lines only to derive log g and ξ by imposing the ionisation equilibrium and by removing the trend between the single line abundances and the reduced equivalent width (REW1), respectively.

For the analysis, we used the LTE code MOOG2 (version 2017, Sneden 1973; Sobeck et al. 2011). We estimated the abundances of C I, Na I, Mg I, Al I, Si I, Ca I, Ti I, Ti II, Cr I, Cr II, Fe I, Fe II, Ni I, and Zn I using the EW method by running the abfind driver. We adopted the same line list used in D’Orazi et al. (2020) that includes 86 Fe I lines, 17 Fe II lines, 57 Ti I lines, 22 Ti II lines, and 42 more lines of different atomic species; Table A.1 provides a complete line list with the atomic data. We added two C I lines to the original line list, taking into account the atomic data from Amarsi et al. (2019), specifically lines 5380.34 and 6587.61Å. We used the Barklem prescriptions for damping values (see Barklem et al. 2000 and references therein).

We measured EWs for all lines via the software ARESv2 (Sousa et al. 2007)3, which calculates EWs through a Gaussian fitting of the line. We discarded the lines with fitting errors larger than 10% and those lines with EWs > 120 mÅ. In this way, we removed strong lines for which the Gaussian approximation is not adequate. We used 1D model atmospheres linearly interpolated from the ATLAS9 grid of Castelli & Kurucz (2003), with new opacities (ODFNEW). We estimated the input values of the Teff and of surface gravities (log g) as in Baratella et al. (2020). The Teff estimates were obtained via 2MASS photometry (Cutri et al. 2003) in the calibrated relation by Casagrande et al. (2010) that is valid for (JK) de-reddened colours. The initial values of the surface gravities (trigonometric gravities, log gtrig) were estimated using the classical equation, based on Gaia DR2 distances as calculated by Bailer-Jones et al. (2018) (see Table 1). Instead, the initial values of ξ were derived using the relation by Dutra-Ferreira et al. (2016), calibrated for dwarf stars, that is ξ(kms1)=0.998+3.16×104X0.253Y2.86×104XY+0.165Y2, \begin{equation*} \begin{split} \xi(\textrm{km\,s}^{-1}) = &\, 0.998&#x002B;3.16\times10^{-4}\,X-0.253\,Y\\ & -2.86\times10^{-4}\,X\,Y&#x002B;0.165Y^2, \end{split} \end{equation*}(1)

where X = T(JK) − 5500 (K) and Y = log g − 4.0 (dex). In the calibrated relation used to derive the Teff, the input metallicity was assumed to be solar, which was later confirmed by the chemical abundances analysis.

For the derivation of Teff we required that the slope of the trend between the Fe + Ti individual line abundances and χ is lower than its error. We adopted the same criteria for ξ, derived from the relation between Ti lines abundances and REWs. Instead, for the log g we required that the difference between Ti I and Ti II is lower than the quadratic sum of the errors on the abundances as calculated by MOOG. The uncertainties on Teff and ξ were calculated by varying each quantity until the slopes of the relative trends are larger than their errors, while for log g it was calculated by varying these parameters until the difference between neutral and ionised species is larger than the total error. The uncertainties on the abundances include the internal errors due to EWs measurements (σ1) and the contribution of the atmospheric parameters (σ2), which is calculated by varying Teff, log g, and ξ one by one by their uncertainties, and calculating the difference with the new abundances.

For the solar atmospheric parameters, we obtained Teff = 5790 ± 75 K, log g = 4.40 ± 0.05 dex and ξ = 0.93 ± 0.05 km s−1. We reportedthe solar abundances of each element in Table 2, where the uncertainties are the quadratic sum of the σ1 and σ2 contributions. As shown, our abundances are in very good agreement with the solar abundances of Asplund et al. (2009). The final values of atmospheric parameters and the derived abundances of neutral and ionised Fe and Ti for the stars in our sample are reported in Table 3. We also calculated the abundance ratios [X/Fe] as [X/Fe] = [X/H]−[Fe/H] (in particular, for the ionised species [X/Fe]II = [X/H]II−[Fe/H]II). The final abundance ratios are reported in Table 4: for star HD 70573, we could not derive the abundances of C and Al because of the relatively high v sin i. The analysis of HD 3765 was not trivial. We derived Teff = 5001 ± 75 K, but if the ionisation equilibrium was satisfied for Ti, this was not the case for Fe and Cr. In particular, we obtained a difference of +0.11 and +0.18 dex between ionised and neutral Fe and Cr species, respectively. The same issue was raised by Ramírez et al. (2007), who obtained a difference between Fe I and Fe II of +0.18 dex. Moreover, we also obtained an anomalously large value of [C/H] = +0.36 ± 0.05 ± 0.09 dex. Similar values of carbon abundances were also obtained for the other cool stars in our sample. This behaviour is discussed extensively in Sect. 5.

Table 2

Solar abundances derived in the present work from the analysis of HARPS-N and GIANO-B spectra.

4 Near-infrared analysis

The GIANO-B spectra were acquired for HD 3765, HD 159222, TYC 1991-1235-1, HIP 61205 and HD 167389. The NIR abundances were measured through spectral synthesis via the driver synth in MOOG. We measured abundances of Na I, Mg I, Al I, Si I, Ca I, Ti I, Fe I, and Ni I using the same line list as D’Orazi et al. (2020). Moreover, we added one C I line to the original line list, the 16021.7 Å line, for which atomic data were taken from Shetrone et al. (2015). Since in general the number of lines of the ionised atomic species is significantly lower in the NIR than in the optical part of the spectrum, we used the atmospheric parameters derived from the analysis in the optical part. To derive the abundances, we synthesised a region of 1000 Å and determined the best instrumental profile. Then we focussed on the line of interest and we derived the given abundance with 0.1 dex steps to find the best-fit profile that minimises the sum of the squared residuals between the synthetic and the observed spectra. In Table 2 we reported the mean values obtained from the analysis of the solar spectrum. As shown, the agreement with the optical values is extremely good, also validating the results of the atmospheric parameters we obtained from the optical analysis. The final abundances for the stars for which we analysed GIANO-B spectra are reported in Table 5. For C, Na, Al, and Mg, we measured only one line, so the uncertainties σ1 and σ2 on the abundances in the table account for the uncertainties on the fitting procedure and the sensitivity of [X/H] to changes in the atmospheric parameters, respectively. Instead, for the remaining elements for which we measured more than one line, we reported the mean values of the abundances, where σ1 is the error on the mean and σ2 related to the atmospheric parameters, respectively.

Recently, Marfil et al. (2020, hereafter M20) analysed CARMENES spectra of a sample of F-G-K stars (wavelength coverage between 5200–17100 Å and R = 95 000 and R = 80 000 in the optical and NIR channel, respectively) with the EWs method, using an extended line list that also comprises Fe I and Fe II lines in the NIR part (216 and 1 lines, respectively) to assess the impact of the NIR lines on the derived stellar parameters. We compared the log gf values of our line list and those of the authors and we find that the values are nearly the same, as shown in Fig. 1. The mean difference between our values and those in M20 is +0.02 ± 0.08 for optical and −0.13 ± 0.26 for NIR; thus we expect to obtain the same results as M20. We applied to our solar spectrum the same procedure and the same line list (the one that the M20 authors optimised for metal-rich dwarfs) as in M20. Since the wavelength coverage of CARMENES spectrograph is different than that of by HARPS-N and GIANO-B, we measured a total of 165 Fe I and Fe I lines adopted from M20 in the solar spectrum, in particular 125 lines in VIS and 40 lines in NIR. The analysis of optical+NIR spectra produced Teff = 5790 ± 50 K, log g = 4.50 ± 0.10 dex, ξ = 0.70 ± 0.10 km s−1, with log(Fe I) = 7.53 ± 0.01 ± 0.04 and log(Fe II) = 7.54 ± 0.05 ± 0.04. These values are very similar to what we obtained from the optical analysis alone and using our line list. The M20 authors, instead, find nearly the same Teff and log g, but ξ = 1.31 ± 0.09 km s−1. We derived the atmospheric parameters applying the same criteria as the code STEPAR (Tabernero et al. 2019) and we find Teff = 5750 ± 75 K, log g = 4.40 ± 0.05 dex and ξ = 0.77 ± 0.15 km s−1, confirming what we previously found. We believe that such a large discrepancy of the ξ values is mainly due to the EW measurements of lines in the NIR. The NIR may pose a challenge when it comes to measuring EWs: for example, telluric lines in emission that remain after the correction, which that the placement of the continuum and a variable S/N that is smaller at shorter wavelengths (Marfil, private communication). Thus, it is possible that we measured different EWs than M20 and these differences may be responsible for the discrepancy in ξ values.

Table 3

Input values of the atmospheric parameters and results of the spectroscopic analysis derived for the stars in our sample from the optical analysis.

Table 4

Abundance ratios obtained from the optical analysis.

Table 5

Mean values of the [X/H] ratios derived from the NIR analysis for four stars in our sample.

thumbnail Fig. 1

Comparison of log gf between ourline lists and those used in Marfil et al. (2020, M20).

5 Discussion

5.1 Stellar parameters and elemental abundances

The final values of the atmospheric parameters for all the stars analysed in our sample, as well as abundances of Fe I, Fe II, Ti I, and Ti II, are reported in Table 3. All the abundance ratios obtained from the optical analysis are reported in Table 4, where the errors were derived as described in Sect. 3. Given the relatively wide range in Teff of our targets, we applied NLTE corrections to Na and C abundances derived from the optical range, following the prescriptions given by Lind et al. (2011) and Amarsi et al. (2019), respectively. The final NIR abundances and uncertainties are reported in Table 5.

The atmospheric parameters we derived with the new approach agree well with the input estimates used in the analysis, as shown in Fig. 2. We calculated the mean difference between the initial guesses and the final spectroscopic values for each parameter. The temperatures are in excellent agreement; ΔTeff = 33 ± 64 K. For the comparison between the spectroscopic and trigonometric gravities, we find a mean difference of −0.08 ± 0.06 dex. As already noted by several authors (Sozzetti et al. 2007; Tsantaki et al. 2013; Maldonado et al. 2015), the spectroscopic gravities tend to be under-estimated with respect to the trigonometric values, especially for log g > 4.50 dex, where Teff < 5000 K. This is again a manifestation of the ionisation balance problem affecting cool dwarf stars, that is enhanced as the stellar age decreases. As an additional check, for the stars with significant difference between spectroscopic and trigonometric gravities the Gaia DR2 astrometric solutions4 (including the parallax) are all well behaved, based on the reduced unit weight error (RUWE) metric (see e.g. Lindegren et al. 2018). This further argues for the discrepancy arising because of above-mentioned limits in the spectroscopic measurements. Finally, regarding the ξ values, we find a mean difference of − 0.004 ± 0.066 km s−1.

In Figs. 3 and 4, we report our abundances as a function of Teff : the open symbols refer to the values we obtained from the optical analysis, while the red symbols represent the results from the NIR analysis. In particular, the diamond symbol is TYC 1991-1235-1, the star symbol stands for HIP 61205, and the pentagon symbol represents HD 167389. The lack of systematic trends between the derived optical abundances and Teff estimates validates our derivation of the atmospheric parameters. We calculated the Pearson correlation coefficient for all the trends in the two figures: none of these is statistically significant at p-value < 0.1 with the exception of Cr II. Despite the low [Mg/H] value obtained for TYC 1989-0049-1, equal to −0.17 ± 0.03 ± 0.04, the trend [Mg/H] versus Teff has a Pearson correlation coefficient r = 0.548, and p = 0.2. For this star only optical spectra are available, so we could not compare this low value with the NIR estimate. We only measured the abundance for the line 4730 Å for which NLTE corrections for the Sun are of the order of 0.01 dex, as calculated by Zhao et al. (2016). The 5711 Å line is strong in the spectrum of this star, has an EW of 143 mÅ, and according to our selection criteria, this line was excluded from the line list for the derivation of the abundance. According to Osorio et al. (2015), the NLTE corrections for a star such as TYC 1989-0049-1 are very small, of the order of −0.006 dex and we expect the NLTE corrections of the line 4730 Å are of the same order.

The ionization equilibrium is satisfied for Ti and also for Fe for the stars in our sample, as shown in Table 3. Interestingly, this is not true for Cr, for which we find an anti-correlation with Teff. As shown in Fig. 4, the Cr II abundances increase at decreasing Teff, especially for stars with Teff ≲ 5400 K. As already mentioned in Sect. 1, this can be explained by the over-ionisation effect. Differences between the neutral and ionised species for some atomic species, such as Fe, Ti, and Cr, have been observed in cool dwarfs with Teff ≲ 5400 K (King et al. 2000; Ramírez et al. 2007; D’Orazi & Randich 2009; Schuler et al. 2010). These differences can reach values up to 0.6–0.8 dex in stars younger than 100 Myr and consequently the value of log g should be decreased. The over-ionization effect is seen in cool dwarf stars, both in OC stars and in field stars (Bensby et al. 2014; Tsantaki et al. 2019). We can see this effect in the cool (Teff = 5001 K) standard star HD 3765 (age ~5 Gyr). While the ionisation equilibrium is satisfied for Ti, suggesting a good estimate of log g for this star, wefind a large discrepancy between Fe I and Fe II of about +0.11 dex, as already noted in Ramírez et al. (2007). A similar discrepancy is also seen for the Cr abundances, for which we obtained a difference of +0.18 dex between Cr I and Cr II, as shown in Table 4. In the analysis, we used two Cr II lines that are not blended with known contaminants, according to Lawler et al. (2017). We note that the NIR Fe abundance is in agreement with the Fe II estimate in the optical range. The reason for such observed discrepancies is still unknown: this may be due to the limitations of 1D-LTE model atmospheres, 3D effects, stellar activity, or a combination of these.

The agreement between the optical and NIR abundances for the stars is overall good, within the uncertainties. Our results corroborate the previous findings of Caffau et al. (2019), who derived chemical abundances for different species of 40 stars by analysing spectra from GIANO in its previous configuration (fibre-fed). However, we noted that for the star HIP 61205 we obtained larger discrepancies for Mg, Si and Ni between the optical and NIR abundance; the latter values are nearly solar. Such discrepancies could be related to the different number of lines used to derive the abundances: more in the optical range (16 lines) than in the NIR (1 line).

In Fig. 5, the abundance ratios [X/H] are plotted as a function of the condensation temperature TC, taken from Lodders (2003). Given the wide range in Teff covered by the stars in our sample (~1500 K), we could not perform a strictly differential analysis with respect to stars of the same association and with similar Teff (see e.g. Meléndez et al. 2009, 2014). For each trend we calculated the Pearson and Spearman correlation coefficients and we find that the trend is significant at p < 0.05 in both cases for HD 167389 alone. It has been suggested that the positive slopes observed in [X/H]-TC plots might be the result of accretion onto the star of refractory material present in the circumstellar disc or a signature of planet engulfment episodes. Instead, for the other stars, the correlation coefficients are not in agreement; thus no exhaustive conclusions can be drawn in those cases.

thumbnail Fig. 2

Comparison between the input estimates of the atmospheric parameters and the derived spectroscopic values. The dash-dotted line represents the 1:1 relation.

thumbnail Fig. 3

Abundances of Fe I and Ti I as function of Teff. The open symbols represent the optical measurements, while the red symbols indicate the NIR measurements. The diamond indicates TYC 1991-1235-1, the star symbol represents HIP 61205, and the pentagon denotes HD 167389, for which we have both optical and NIR measurements. The Pearson correlation coefficient of the trend in the top panel is r = 0.67, which is not significant at p < 0.1. For Ti the Pearson correlation coefficient is r = −0.03, which is not significant at p < 0.1.

5.2 Comparison with previous studies

Our abundance measurements are in overall fair agreement with other studies found in the literature, as shown in Table 6, where we reported the mean values for each cluster and results from different studies. Regarding the Coma Berenices cluster, we find that our measurements are in good agreement with Blanco-Cuaresma et al. (2015) for all the atomic species. Netopil et al. (2016) report a mean [Fe/H] equal to 0.00 ± 0.08, averaging over different estimates in the literature. Other studies on the chemical composition of this cluster are mainly focussed on the analysis of A-F type stars that have temperatures higher than 6000 K. Burkhart & Coupry (2000) reported a mean [Fe/H] = − 0.07 ± 0.05 dex for 1 F-type star; on the contrary, Gebran et al. (2008) found ⟨[Fe∕H]⟩ = 0.07 ± 0.09 dex for 11 F-type stars analysed through the spectral synthesis technique. The large difference between these two studies could be related to the different line lists and techniques employed and also to the different number of stars analysed. Friel & Boesgaard (1992) analysed high-resolution, high S/N spectra of 14 F-G type stars, with Teff > 5950 K, through the EW method. In particular, they analysed the spectral window 6078–7755 Å, where they measured 8 Fe I lines. These authors found a mean metallicity of −0.05 ± 0.05 dex, again in good agreement with our measurements, despite the different type of stars analysed. Regarding the individual stars, we find thatour measurement for HIP 61205 confirmed the results of Brewer et al. (2016). The authors find Teff = 5796 K, log g = 4.51 dex and [Fe/H] = −0.02, which are in excellent agreement with our estimates. Brewer et al. (2016) also derived abundances for different atomic species, in particular they find that [C/Fe] = −0.04, [Na/Fe] = −0.12, [Mg/Fe] = −0.06, [Al/Fe] = −0.12, [Si/Fe] = −0.04, [Ca/Fe] = 0.02, [Ti/Fe] = −0.02, [Cr/Fe] = 0.01, and [Ni/Fe] = −0.09. All abundance ratios agree very well with our estimates, as shown in Table 4. To our knowledge, for TYC 1991-1235-1, TYC 1989-0049-1, and TYC 1989-147-1, there are no previous studies on abundances in the literature.

For the UMa moving group, our estimate of mean [Fe/H] = −0.01 ± 0.01 is in fair agreement with the results from Soderblom & Mayor (1993), King & Schuler (2005), and Monier (2005), which reported mean values equal to −0.08 ± 0.09, −0.06 ± 0.05, and −0.05 ± 0.02, respectively. Our results also confirm what Biazzo et al. (2012) and Tabernero et al. (2017) find; these authors analysed stars similar to our sample employing the EW method. As shown in Table 6, our results agree well with the two studies. HD 167389 has been analysed by Ammler-von Eiff & Guenther (2009) and Tabernero et al. (2017). In particular, the former derived the stellar parameters and abundances of Fe and Mg through spectral synthesis: they find Teff = 5895 ± 80 K, log g = 4.37 ± 0.15 dex, ξ = 0.99 ± 0.20 km s−1, [Fe/H] = −0.02 ± 0.07, [Mg/Fe] = −0.03 ± 0.05, in excellent agreement with our estimates. Tabernero et al. (2017) analysed candidate members of the UMa group to confirm their membership through chemical tagging by employing the EW analysis method. The authors find for HD 167389 Teff = 5978 K, log g = 4.56 dex and [Fe/H] = +0.01, confirming our results. Moreover, they derived abundances for various atomic species, finding [Na/Fe] = −0.06 ± 0.01, [Mg/Fe] = −0.07 ± 0.03, [Al/Fe] = −0.05 ± 0.01, [Si/Fe] = −0.02 ± 0.01, [Ca/Fe] = 0.03 ± 0.01, [Ti/Fe] = 0.01 ± 0.01, [Cr/Fe] = −0.01 ± 0.01, and [Ni/Fe] = −0.04 ± 0.01. HD 59747 was analysed by Ammler-von Eiff & Guenther (2009), who derived the stellar parameters, [Fe/H], and [Mg/H] through spectral synthesis fitting. The authors found Teff = 5094 K, log g = 4.55 dex, [Fe/H] = −0.03, and [Mg/Fe] = −0.01, in excellent agreement with our results.

In the Her-Lyr association, we analysed only the star HD 70573, which is also the star with the highest v sin i in our sample. This star was also analysed by Gonzalez et al. (2010) in the standard way, that is using Fe (neutral and ionised) lines to derive the atmospheric parameters. These authors found Teff = 5807 ± 85 K, log g = 4.35 ± 0.08 dex, ξ = 1.80 ± 0.16 km s−1 and [Fe/H]= −0.05 ± 0.06. These results are also confirmed by Ghezzi et al. (2010), who analysed the star in the same way, finding Teff = 5884 ± 26 K, log g = 4.57 ± 0.08 dex, ξ = 1.69 ± 0.06 km s−1 and [Fe/H]= −0.04 ± 0.03. Our results confirm the values reported in the two different studies, with the exception of ξ, for which we find a lower value equal to 1.10 ± 0.10 km s−1 ; there is a difference of the order of 0.6–0.7 km s−1 between the two studies. Such a discrepancy could be due to the different line lists used, in particular to differences in the atomic data, but they also seem to confirm the results of Baratella et al. (2020), regarding the possible overestimation of the ξ parameter when using iron lines. As shown in Table 6, our results confirm those found by Brewer et al. (2016). However we note large differences for some elements, such as Na, Mg, and Ni, which could be due to differences in the line list used, in particular to differences in the atomic line parameters, and to the different spectroscopic analysis technique employed.

thumbnail Fig. 4

Individual values of [X/H] as a function of spectroscopic estimates of Teff, derived from the analysis of the optical spectra (open symbols) and from the analysis of NIR spectra (red symbols). The symbols for the three stars for which we analysed GIARPS spectra are the same as in Fig. 3. All trends have Pearson correlation coefficients that are not significant at p < 0.1, apart from Cr II (see the text for details).

thumbnail Fig. 5

Values of [X/H] as a function of the condensation temperature TC, taken from Lodders (2003).

Table 6

Mean values of the abundances ratios for each cluster and comparison with literature studies.

5.3 Carbon abundances

We measured C abundance using four different indicators: two high-excitation potential lines at 5380.337 and 6587.61 Å (atomic data and NLTE corrections from Amarsi et al. 2019); the NIR line at 16021.7 Å (atomic data from Shetrone et al. 2015); and the CH molecular band at 4300 Å, for which molecular line data come from Plez (priv. comm.). We note that the line 16021 Å suffers from a blend with Fe and Si lines in the blue wing. This blend is not significant for the Sun, TYC 1991-1235-1, and HIP 61205, but it becomes more important in HD 167389 (Teff = 6000 K), because the Si feature is a high-excitation line (χ = 7.035 eV), which strengthens at these temperatures. For this reason, we could not derive the abundance in the NIR range for this star. The abundances of C I lines in the optical part were calculated with the EW method, while the abundances for the NIR lines and for the CH feature were calculated through spectral synthesis. We applied NLTE corrections to C I optical abundances following Amarsi et al. (2019). Despite the wide range in Teff covered by the stars we analysed, NLTE corrections of optical lines are small, typically of the order of −0.01 dex.

Carbon is one of the most important elements for life on Earth and also for planetary formation models. In particular, deriving reliable estimates of the C/O ratio is crucial, since it provides clues to where the planets formed in the protoplanetary disc and possible subsequent radial migration (Brewer et al. 2017). However, we could not derive abundances of oxygen(O) from OH molecular features in the NIR, because in F-G stars molecular absorption is less important, weakening the lines (Souto et al. 2018). For the K-type star TYC 1989-0049-1 GIANO-B spectrum is not available, hampering the determination of abundances through the NIR molecular lines. Also, extremely high-resolution, high S/N spectra are required to be able to measure OH lines (Meléndez 2004). Additionally, the HARPS-N spectra does not allow us to cover the O I triplet at 7773 Å, which is ideal for solar-type stars. Despite the forest of CN lines in the solar spectrum covering the blue, red, and NIR part of it (Sneden & Lambert 1982), we did not derive N abundance from those molecular lines. First of all, for solar-type stars the best tools to derive reliable N abundance are high-excitation N I atomic lines at 7400–8720 Å (Asplund et al. 2009; Sneden et al. 2014), which is not covered by HARPS-N (and GIANO-B) spectra. Moreover, we notice that N has little impact on the molecular equilibrium.Thus, we derive only C abundance, fixing O and N abundances to solar values, which is a reasonable assumption for our sample that is comprised of intermediate-age, thin disc, main-sequence stars (e.g. Bensby et al. 2014).

The values for the Sun are reported in Table 2; for C abundance inferred from the CH features we obtained log n(C) = 8.35 ± 0.08, that is marginally lower than that obtained from the atomic lines. Such difference can be becasue the atomic data of the CH feature are not so precise. In Fig. 6 the different C abundance estimates are reported as a function of Teff. The empty symbols represent the values obtained in the optical range, the red symbols stand for the NIR measurements, and the blue symbols indicate the C values from CH molecular band. The different stars are represented by the different symbols, as described in the caption of the figure. As shown, we obtained different trends from the different lines. The increasing C I abundances at decreasing Teff for the optical measurements is noteworthy. However, this is not seen for the NIR abundances and the values derived from the CH. Since NLTEcorrections of C I optical abundances are negligible for the stars in our sample, we believe that its trend with Teff is due to over-excitation effects. Schuler et al. (2015) find similar behaviour when deriving C abundance from two high-excitation lines, which have χ similar to our adopted lines, for a star with Teff = 5406 K. The authors find that the C abundance inferred from atomic lines is +0.16 dex higher than that derived from the C2 feature. The interesting aspect is that even if the NIR C line has a high-excitation energy (χ = 9.631 eV), we do not see the same effect in the abundances. A possible explanation is that at 1.6 μm there is the H absorption minimum, so we see the deepest photospheric layers of the atmosphere where LTE is a good approximation to compute the populations of atomic levels. Schuler et al. (2015) argue that the over-excitation could be explained as a NLTE effect, specifically it could be the result of our incapacity to properly model the population of high-energy levels under LTE approximation. Thus, the C abundances obtained from the NIR line, even if it has a high-excitation energy, could be the real C abundances of the stars. Moreover, while for star TYC 1991-1235-1, Teff = 5070 K, we obtained a difference between optical and NIR abundances of +0.26 dex, for HIP 61205, Teff = 5825 K, the behaviour is reversed; in this star, the optical C abundanceis smaller than the NIR estimate by −0.22 dex. Unfortunately, we could not confirm this trend in the whole temperature range, since GIANO-B spectra are not available for all the stars. As a further test, we measured C abundance from the CH band at 4300 Å. We find that the trend with Teff is not statistically meaningful (with a p-value > 0.1) and, most importantly, the over-excitation effect observed for the C I abundances is not present. Also, especially for TYC 1991-1235-1, the C abundance from CH is in better agreement with the NIR estimate than the optical. We suggest that the values obtained from molecular features in the optical for very young stars are more reliable estimates of C abundances, as already suggested by Schuler et al. (2015). For this reason, the [C/Fe] values in Table 4 are calculated with the C abundance derived from CH lines.

thumbnail Fig. 6

Abundances of C I as a function of Teff, derived from the optical analysis (empty symbols), from the NIR line (red symbols), and from CH band at 4300 Å. The different symbols are the different estimates for the same star: the diamond represents TYC 1991-1235-1, the star HIP 61205, the pentagon HD 167389, the circle HD 70573, the triangle HD 59747, the square TYC 1989-0049-1, and finally the x-shaped symbol indicates TYC 1989-147-1.

5.4 Effects of stellar activity

The over-excitation and over-ionisation effects that we observed for the C and Cr II abundances are among the principal problems affecting the analysis of young cool stars. Such effects are more evident in young, intermediate-age (τ ≲ 800 Myr) and cool dwarf stars (Teff ≲5400 K), which are more intense at decreasing ages and temperatures. Young stars are more active and they have more intense chromospheric or photospheric magnetic fields than older stars. The main effect of local magnetic fields on spectral lines is the broadening of their profile through the Zeeman effect that causes a splitting of the spectral line into its multipletcomponents. This effect is directly proportional to the wavelength and to the value of the Landè gL factor. The latter parameter measures the sensitivity of an atomic transition to the magnetic fields, meaning the higher the gL factor, the more sensitive the line is to Zeeman splitting. The C I NIR (16021 Å) and optical (5380 Åand 6588 Å) lines have a Landè factor equal to 1.15, 1.0, and 1.0, respectively. For the Cr II lines at 4848.23 and 5237.33 Å the gL is equal to 1.25 and 1.30, respectively. According to Shchukina et al. (2015), these lines are insensitive to the presence of magnetic fields, which produce the stronger effects in lines with gL ~ 2.0 typically. Moreover, the Zeeman splitting has two main effects on the spectral line. On one hand, it produces a broadening of the profile and an increase of the EW. On the other hand, the line weakens, with a decrease of its depth; thus,the two effects compensate for each other (Reiners et al. 2013). In this case, we can exclude the Zeeman effect as a possible explanation of the observed trends in Figs. 4 and 6.

In Fig. 7, we plot the estimates of the C abundances from the two atomic lines in the optical range and from the CH molecule as a function of the activity index logRHK$\log R ^{\prime}_{\rm{HK}}$. The latter values were calculated with Yabi5 interface (Hunter et al. 2012), following the prescription of Noyes et al. (1984) and through the procedure described by Lovis et al. (2011). Yabi is a Python web application installed at IA2 in Trieste that allows authorised users to run the HARPS-N DRS pipeline on proprietary data with custom input parameters. Since we analysed co-added spectra of the stars, we calculated the mean values of the activity indexes, averaging over the spectra we used in the co-adding procedure. The logRHK$\log R ^{\prime}_{\rm{HK}}$ indexes are reported in Table 1. As shown, the C I abundances from atomic lines in the optical range have a positive correlation with the logRHK$\log R ^{\prime}_{\rm{HK}}$, with a Pearson correlation coefficient of r = 0.69 and p-value = 0.08. On the contrary, the abundance values derived from the CH features do not correlate with the activity indexes, again suggesting that these values are more reliable estimates of C abundances. We also find a similar behaviour for Cr II abundances, as shown in Fig. 8. We may envisage different, plausible scenarios to explain this peculiar trend. The effect of the chromospheric emission in active stars is observed mainly in the Lyman-α lines of H atom, the Ca II H and K lines, Mg II lines, and He lines. In particular, the photons from the Lyman-α, with an energy of 10.2 eV, can ionise Cr atoms that have a first ionisation potential of 6.77 eV. So, the population of Cr II atoms is larger than Cr I and this can qualitatively explain the increase of abundance at increasing levels of activity. In the case of C I lines, a possible explanation of the overabundance is the presence of unknown blends in the optical lines that become stronger as the Teff decreases. Moreover, these blends could be more important in active stars than in quiet stars, if a significant part of the flux comes from cool regions, such as photospheric dark spots, where the lines are stronger and/or more sensible to the temperature. Another possible explanation is that the increase of the population of the atomic levels (from which the 5380 and 6587 Å lines form) is mainly due to UV continuum photons between 1450 and 1650 Å, which increase in intensity at increasing levels of activity (Linsky et al. 2012). These photons might be responsible of the larger population of the levels from which the 5380 and 6587 Å lines are formed. Although we do not have a definitive explanation to the over-ionisation/excitation effects, the solutions proposed seem intriguing and worthy of a detailed investigation.

We also find that the ξ values we obtained using Ti lines do not seem to correlate with logRHK$\log R ^{\prime}_{\rm{HK}}$, as shown in Fig. 9. The symbols in this figure are colour coded according to the Teff. As already known, the ξ velocity increases systematically towards higher Teff and lower log g. In particular, in dwarf stars (log g ~ 4.50 dex) the ξ values are of the order of 0.70km s−1 at Teff ~ 4500 K (Steffen et al. 2013). Thus, the trend of ξ with Teff observed in Fig. 9 is expected. We calculated the Pearson correlation coefficient for the trend, that is equal to r = −0.49, with p = 0.26; thus it is not significant at p < 0.10. This result further validates our method and what was previously found by Baratella et al. (2020).

thumbnail Fig. 7

Abundances of C I derived from atomic lines in the optical range (empty symbols) and from the CH molecular features (blue symbols) as a function of the activity index logRHK$\log R^{\prime}_{\rm{HK}}$. The symbols are the same as in Fig. 6.

thumbnail Fig. 8

Abundances of Cr II as a function of activity index logRHK$\log R^{\prime}_{\rm{HK}}$. The symbols are the same as in Fig. 6.

thumbnail Fig. 9

Values of ξ parameterderived with the new approach as a function of the chromospheric activity index logRHK$\log R ^{\prime}_{\rm{HK}}$. The symbolsare colour-coded according to the Teff.

6 Conclusions

In this first paper of a series, we presented the preliminary results of an extensive analysis of optical and NIR spectra of stars observed by the GAPS-YO programme. In particular, we derived the atmospheric parameters and the chemical composition of seven target stars, the Sun, and two RV standard stars, HD 159222 and HD 3765, using a new spectroscopic approach to overcome analytical issues related to the relatively young ages of the stars.

The analysis of young and intermediate-age stars, in particular in the cool temperature regime (Teff ≲ 5400 K), is not trivial, owing to a series of effects still unexplained from a theoretical point of view. For these reasons, we applied the same methodology as in Baratella et al. (2020) for the analysis of the optical HARPS-N spectra. In general, our derived spectroscopic estimates of the atmospheric parameters are in excellent agreement with the initial guesses. The atmospheric parameters we obtained from the optical analysis were used to derive the abundances in the NIR part, through the spectral synthesis technique and using the same line list as in D’Orazi et al. (2020). We derived abundances for 11 atomic species, both α-, proton-capture and iron-peak elements. Overall, we find a good agreement between optical and NIR abundances. The lack of trends between [X/H] and Teff confirm that our analysis is reliable, with the exception of Cr II, for which instead we observed increasing abundances at decreasing temperatures. This trend confirms the previous findings of Schuler et al. (2006, 2010) about the over-ionisation effects. Our derivation of C I abundances from optical atomic lines reveals a similar effect. The two lines used have high-excitation potential and they yield higher abundances (up to almost +1.0 dex) at decreasing Teff. In the NIR, we analysed another high-excitation line, 16021 Å, but only in two stars, TYC 1991-1235-1 and HIP 61205. Despite what we obtained from the optical lines, in the NIR we do not see the same effect as in the optical. Schuler et al. (2015) find a trend similar to what we observed for two C lines with χ > 7eV and these authors suggested that the C abundance from C2 features is more reliable. In a similar way, we derived C abundances from CH molecular band at 4300 Å. At variance with what we obtained from the atomic lines, we did not observe the same trend for the new abundance determinations, in agreement with the findings of Schuler et al. (2015). We suggest that for very young and cool stars the C abundance derived from molecular lines is more reliable. The over-ionisation/excitation effects could be explained by a combination of different factors, such as the higher level of activity due to the young age of the stars and the presence of intense local chromospheric and/or photospheric magnetic fields that can alter the line profiles. Indeed, we find a positive correlation between the C abundances derived from the atomic lines in the optical range and the activity indexes logRHK$\log R^{\prime}_{\rm{HK}}$, suggesting that these effects are related to higher activity levels. This behaviour is not seen in the C estimates from CH molecular features. We also find a positive correlation between the Cr II values and logRHK$\log R^{\prime}_{\rm{HK}}$. However, as already pointed out in previous studies (Baratella et al. 2020; Spina et al. 2020), the main causes are still unknown and they may be a combination of different factors, most likely a combination of more intense chromospheric or photospheric magnetic fields. Finding a theoretical explanation to these issues is beyond the scopes of this paper, but the topic is interesting and deserves a deeper investigation.

Acknowledgements

We thank the anonymous referee for her/his very helpful comments and suggestions. The authors acknowledgesupport by INAF/WOW and INAF/FRONTIERA through the “Progetti Premiali” funding scheme of the Italian Ministry of Education, University, and Research.

Appendix A Optical line list

The line list used in the analysis of the HARPS-N optical spectra is shown in Table A.1. The source of oscillator strengths include the NIST database, Lawler et al. (2013) for Ti lines, line lists published by D’Orazi et al. (2017), and Ruffoni et al. (2014) for Fe I.

Table A.1

Line list for the HARPS-N spectra.

References

  1. Adibekyan, V. 2019, Geosciences, 9, 105 [NASA ADS] [CrossRef] [Google Scholar]
  2. Aleo, P. D., Sobotka, A. C., & Ramírez, I. 2017, ApJ, 846, 24 [CrossRef] [Google Scholar]
  3. Amarsi, A. M., Nissen, P. E., & Skúladóttir, Á. 2019, A&A, 630, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Ammler-von Eiff, M., & Guenther, E. W. 2009, A&A, 508, 677 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Mantelet, G., & Andrae, R. 2018, AJ, 156, 58 [NASA ADS] [CrossRef] [Google Scholar]
  7. Baratella, M., D’Orazi, V., Carraro, G., et al. 2020, A&A, 634, A34 [CrossRef] [EDP Sciences] [Google Scholar]
  8. Barklem, P. S., Piskunov, N., & O’Mara, B. J. 2000, A&AS, 142, 467 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
  9. Bensby, T., Feltzing, S., & Oey, M. S. 2014, A&A, 562, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Biazzo, K., Randich, S., & Palla, F. 2011a, A&A, 525, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Biazzo, K., Randich, S., Palla, F., & Briceño, C. 2011b, A&A, 530, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Biazzo, K., D’Orazi, V., Desidera, S., et al. 2012, MNRAS, 427, 2905 [NASA ADS] [CrossRef] [Google Scholar]
  13. Blanco-Cuaresma, S., Soubiran, C., Heiter, U., et al. 2015, A&A, 577, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Brewer, J. M., Fischer, D. A., Valenti, J. A., & Piskunov, N. 2016, ApJS, 225, 32 [Google Scholar]
  15. Brewer, J. M., Fischer, D. A., & Madhusudhan, N. 2017, AJ, 153, 83 [NASA ADS] [CrossRef] [Google Scholar]
  16. Burkhart, C., & Coupry, M. F. 2000, A&A, 354, 216 [NASA ADS] [Google Scholar]
  17. Caffau, E., Bonifacio, P., Oliva, E., et al. 2019, A&A, 622, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Carleo, I., Sanna, N., Gratton, R., et al. 2016, Exp. Astron., 41, 351 [NASA ADS] [CrossRef] [Google Scholar]
  19. Carleo, I., Benatti, S., Lanza, A. F., et al. 2018, A&A, 613, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Carleo, I., Malavolta, L., Lanza, A. F., et al. 2020, A&A, 638, A5 [CrossRef] [EDP Sciences] [Google Scholar]
  21. Casagrande, L., Ramírez, I., Meléndez, J., Bessell, M., & Asplund, M. 2010, A&A, 512, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Castelli, F., & Kurucz, R. L. 2003, IAU Symp., 210, A20 [Google Scholar]
  23. Claudi, R., Benatti, S., Carleo, I., et al. 2017, Eur. Phys. J. Plus, 132, 364 [CrossRef] [Google Scholar]
  24. Cosentino, R., Lovis, C., Pepe, F., et al. 2014, SPIE Conf. Ser., 9147, 91478C [Google Scholar]
  25. Covino, E., Esposito, M., Barbieri, M., et al. 2013, A&A, 554, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, 2MASS All Sky Catalog of point sources [Google Scholar]
  27. De Silva, G. M., Freeman, K. C., Bland-Hawthorn, J., et al. 2015, MNRAS, 449, 2604 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  28. D’Orazi, V., & Randich, S. 2009, A&A, 501, 553 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. D’Orazi, V., Randich, S., Flaccomio, E., et al. 2009, A&A, 501, 973 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. D’Orazi, V., Biazzo, K., & Randich, S. 2011, A&A, 526, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. D’Orazi, V., Desidera, S., Gratton, R. G., et al. 2017, A&A, 598, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. D’Orazi, V., Oliva, E., Bragaglia, A., et al. 2020, A&A, 633, A38 [CrossRef] [EDP Sciences] [Google Scholar]
  33. Dutra-Ferreira, L., Pasquini, L., Smiljanic, R., Porto de Mello, G. F., & Steffen, M. 2016, A&A, 585, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Eisenbeiss, T., Ammler-von Eiff, M., Roell, T., et al. 2013, A&A, 556, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Folsom, C. P., Petit, P., Bouvier, J., et al. 2016, MNRAS, 457, 580 [NASA ADS] [CrossRef] [Google Scholar]
  36. Friel, E. D., & Boesgaard, A. M. 1992, ApJ, 387, 170 [NASA ADS] [CrossRef] [Google Scholar]
  37. Gebran, M., Monier, R., & Richard, O. 2008, A&A, 479, 189 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Ghezzi, L., Cunha, K., Smith, V. V., et al. 2010, ApJ, 720, 1290 [NASA ADS] [CrossRef] [Google Scholar]
  39. Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25 [NASA ADS] [Google Scholar]
  40. Gonzalez, G., Carlson, M. K., & Tobin, R. W. 2010, MNRAS, 403, 1368 [NASA ADS] [CrossRef] [Google Scholar]
  41. Harutyunyan, A., Rainer, M., Hernandez, N., et al. 2018, SPIE Conf. Ser., 10706, 1070642 [Google Scholar]
  42. Hunter, A. A., Macgregor, A. B., Szabo, T. O., Wellington, C. A., & Bellgard, M. I. 2012, Source Code Biol. Med., 7, 1 [CrossRef] [Google Scholar]
  43. James, D. J., Melo, C., Santos, N. C., & Bouvier, J. 2006, A&A, 446, 971 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Jofré, P., Heiter, U., & Soubiran, C. 2019, ARA&A, 57, 571 [Google Scholar]
  45. King, J. R., & Schuler, S. C. 2005, PASP, 117, 911 [NASA ADS] [CrossRef] [Google Scholar]
  46. King, J. R., Soderblom, D. R., Fischer, D., & Jones, B. F. 2000, ApJ, 533, 944 [NASA ADS] [CrossRef] [Google Scholar]
  47. Lawler, J. E., Guzman, A., Wood, M. P., Sneden, C., & Cowan, J. J. 2013, ApJS, 205, 11 [NASA ADS] [CrossRef] [Google Scholar]
  48. Lawler, J. E., Sneden, C., Nave, G., et al. 2017, ApJS, 228, 10 [CrossRef] [Google Scholar]
  49. Lind, K., Asplund, M., Barklem, P. S., & Belyaev, A. K. 2011, A&A, 528, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Linsky, J. L., Bushinsky, R., Ayres, T., Fontenla, J., & France, K. 2012, ApJ, 745, 25 [NASA ADS] [CrossRef] [Google Scholar]
  52. Lodders, K. 2003, ApJ, 591, 1220 [NASA ADS] [CrossRef] [Google Scholar]
  53. López-Santiago, J., Montes, D., Crespo-Chacón, I., & Fernández-Figueroa, M. J. 2006, ApJ, 643, 1160 [NASA ADS] [CrossRef] [Google Scholar]
  54. Lovis, C., Dumusque, X., Santos, N. C., et al. 2011, ArXiv e-prints, [arXiv:1107.5325] [Google Scholar]
  55. Luck, R. E. 2017, AJ, 153, 21 [NASA ADS] [CrossRef] [Google Scholar]
  56. Malavolta, L., Nascimbeni, V., Piotto, G., et al. 2016, A&A, 588, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Maldonado, J., Eiroa, C., Villaver, E., Montesinos, B., & Mora, A. 2015, A&A, 579, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Marfil, E., Tabernero, H. M., Montes, D., et al. 2020, MNRAS, 492, 5470 [CrossRef] [Google Scholar]
  59. Marsden, S. C., Petit, P., Jeffers, S. V., et al. 2014, MNRAS, 444, 3517 [NASA ADS] [CrossRef] [Google Scholar]
  60. Martínez-Arnáiz, R., Maldonado, J., Montes, D., Eiroa, C., & Montesinos, B. 2010, A&A, 520, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Meléndez, J. 2004, ApJ, 615, 1042 [NASA ADS] [CrossRef] [Google Scholar]
  62. Meléndez, J., Asplund, M., Gustafsson, B., & Yong, D. 2009, ApJ, 704, L66 [NASA ADS] [CrossRef] [Google Scholar]
  63. Meléndez, J., Ramírez, I., Karakas, A. a. I., et al. 2014, ApJ, 791, 14 [NASA ADS] [CrossRef] [Google Scholar]
  64. Mermilliod, J. C., Grenon, M., & Mayor, M. 2008, A&A, 491, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Mermilliod, J. C., Mayor, M., & Udry, S. 2009, A&A, 498, 949 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Monier, R. 2005, A&A, 442, 563 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Montes, D., López-Santiago, J., Gálvez, M. C., et al. 2001, MNRAS, 328, 45 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  68. Netopil, M., Paunzen, E., Heiter, U., & Soubiran, C. 2016, A&A, 585, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Nissen, P. E. 2015, A&A, 579, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & Vaughan, A. H. 1984, ApJ, 279, 763 [NASA ADS] [CrossRef] [Google Scholar]
  71. Oliva, E., Origlia, L., Baffa, C., et al. 2006, SPIE Conf. Ser., 6269, 626919 [Google Scholar]
  72. Osorio, Y., Barklem, P. S., Lind, K., et al. 2015, A&A, 579, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Poretti, E., Boccato, C., Claudi, R., et al. 2016, Mem. Soc. Astron. It., 87, 141 [NASA ADS] [Google Scholar]
  74. Rainer, M., Harutyunyan, A., Carleo, I., et al. 2018, SPIE Conf. Ser., 10702, 1070266 [Google Scholar]
  75. Ramírez, I., Allende Prieto, C., & Lambert, D. L. 2007, A&A, 465, 271 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Ramírez, I., Allende Prieto, C., & Lambert, D. L. 2013, ApJ, 764, 78 [NASA ADS] [CrossRef] [Google Scholar]
  77. Reiners, A., Shulyak, D., Anglada-Escudé, G., et al. 2013, A&A, 552, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Ruffoni,M. P., Den Hartog, E. A., Lawler, J. E., et al. 2014, MNRAS, 441, 3127 [NASA ADS] [CrossRef] [Google Scholar]
  79. Santos, N. C., Melo, C., James, D. J., et al. 2008, A&A, 480, 889 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Schuler, S. C., King, J. R., Hobbs, L. M., & Pinsonneault, M. H. 2004, ApJ, 602, L117 [NASA ADS] [CrossRef] [Google Scholar]
  81. Schuler, S. C., Hatzes, A. P., King, J. R., Kürster, M., & The, L.-S. 2006, AJ, 131, 1057 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  82. Schuler, S. C., Plunkett, A. L., King, J. R., & Pinsonneault, M. H. 2010, PASP, 122, 766 [NASA ADS] [CrossRef] [Google Scholar]
  83. Schuler, S. C., Vaz, Z. A., Katime Santrich, O. o. J., et al. 2015, ApJ, 815, 5 [NASA ADS] [CrossRef] [Google Scholar]
  84. Shchukina, N. G., Sukhorukov, A. V., & Trujillo Bueno, J. 2015, IAU Symp., 305, 368 [Google Scholar]
  85. Shetrone, M., Bizyaev, D., Lawler, J. E., et al. 2015, ApJS, 221, 24 [Google Scholar]
  86. Silaj, J., & Landstreet, J. D. 2014, A&A, 566, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Sneden, C. A. 1973, PhD thesis, The University of Texas at Austin, Texas, USA [Google Scholar]
  88. Sneden, C., & Lambert, D. L. 1982, ApJ, 259, 381 [NASA ADS] [CrossRef] [Google Scholar]
  89. Sneden, C., Lucatello, S., Ram, R. S., Brooke, J. S. A., & Bernath, P. 2014, ApJS, 214, 26 [Google Scholar]
  90. Sobeck, J. S., Kraft, R. P., Sneden, C., et al. 2011, AJ, 141, 175 [NASA ADS] [CrossRef] [Google Scholar]
  91. Soderblom, D. R., & Mayor, M. 1993, AJ, 105, 226 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  92. Sousa, S. G., Santos, N. C., Israelian, G., Mayor, M., & Monteiro, M. J. P. F. G. 2007, A&A, 469, 783 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Souto, D., Cunha, K., Smith, V. V., et al. 2018, ApJ, 857, 14 [NASA ADS] [CrossRef] [Google Scholar]
  94. Sozzetti, A., Torres, G., Charbonneau, D., et al. 2007, ApJ, 664, 1190 [NASA ADS] [CrossRef] [Google Scholar]
  95. Spina, L., Randich, S., Magrini, L., et al. 2017, A&A, 601, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. Spina, L., Nordlander, T., Casey, A. R., et al. 2020, ApJ, 895, 52 [CrossRef] [Google Scholar]
  97. Steffen, M., Caffau, E., & Ludwig, H. G. 2013, Mem. Soc. Astron. It. Suppl., 24, 37 [Google Scholar]
  98. Tabernero, H. M., Montes, D., González Hernández, J. I., & Ammler-von Eiff, M. 2017, A&A, 597, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Tabernero, H. M., Marfil, E., Montes, D., & González Hernández, J. I. 2019, A&A, 628, A131 [CrossRef] [EDP Sciences] [Google Scholar]
  100. Takeda, Y., & Honda, S. 2020, AJ, 159, 174 [CrossRef] [Google Scholar]
  101. Teske, J. K., Cunha, K., Schuler, S. C., Griffith, C. A., & Smith, V. V. 2013, ApJ, 778, 132 [NASA ADS] [CrossRef] [Google Scholar]
  102. Tsantaki, M., Sousa, S. G., Adibekyan, V. Z., et al. 2013, A&A, 555, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Tsantaki, M., Santos, N. C., Sousa, S. G., et al. 2019, MNRAS, 485, 2772 [NASA ADS] [CrossRef] [Google Scholar]
  104. Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  105. Viana Almeida, P., Santos, N. C., Melo, C., et al. 2009, A&A, 501, 965 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  106. Yana Galarza, J., Meléndez, J., Lorenzo-Oliveira, D., et al. 2019, MNRAS, 490, L86 [NASA ADS] [CrossRef] [Google Scholar]
  107. Zhao, G., Mashonkina, L., Yan, H. L., et al. 2016, ApJ, 833, 225 [Google Scholar]

1

REW = log(EW/λ).

4

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.

All Tables

Table 1

Some basic information of both standard stars and members in young associations analysed in this work.

Table 2

Solar abundances derived in the present work from the analysis of HARPS-N and GIANO-B spectra.

Table 3

Input values of the atmospheric parameters and results of the spectroscopic analysis derived for the stars in our sample from the optical analysis.

Table 4

Abundance ratios obtained from the optical analysis.

Table 5

Mean values of the [X/H] ratios derived from the NIR analysis for four stars in our sample.

Table 6

Mean values of the abundances ratios for each cluster and comparison with literature studies.

Table A.1

Line list for the HARPS-N spectra.

All Figures

thumbnail Fig. 1

Comparison of log gf between ourline lists and those used in Marfil et al. (2020, M20).

In the text
thumbnail Fig. 2

Comparison between the input estimates of the atmospheric parameters and the derived spectroscopic values. The dash-dotted line represents the 1:1 relation.

In the text
thumbnail Fig. 3

Abundances of Fe I and Ti I as function of Teff. The open symbols represent the optical measurements, while the red symbols indicate the NIR measurements. The diamond indicates TYC 1991-1235-1, the star symbol represents HIP 61205, and the pentagon denotes HD 167389, for which we have both optical and NIR measurements. The Pearson correlation coefficient of the trend in the top panel is r = 0.67, which is not significant at p < 0.1. For Ti the Pearson correlation coefficient is r = −0.03, which is not significant at p < 0.1.

In the text
thumbnail Fig. 4

Individual values of [X/H] as a function of spectroscopic estimates of Teff, derived from the analysis of the optical spectra (open symbols) and from the analysis of NIR spectra (red symbols). The symbols for the three stars for which we analysed GIARPS spectra are the same as in Fig. 3. All trends have Pearson correlation coefficients that are not significant at p < 0.1, apart from Cr II (see the text for details).

In the text
thumbnail Fig. 5

Values of [X/H] as a function of the condensation temperature TC, taken from Lodders (2003).

In the text
thumbnail Fig. 6

Abundances of C I as a function of Teff, derived from the optical analysis (empty symbols), from the NIR line (red symbols), and from CH band at 4300 Å. The different symbols are the different estimates for the same star: the diamond represents TYC 1991-1235-1, the star HIP 61205, the pentagon HD 167389, the circle HD 70573, the triangle HD 59747, the square TYC 1989-0049-1, and finally the x-shaped symbol indicates TYC 1989-147-1.

In the text
thumbnail Fig. 7

Abundances of C I derived from atomic lines in the optical range (empty symbols) and from the CH molecular features (blue symbols) as a function of the activity index logRHK$\log R^{\prime}_{\rm{HK}}$. The symbols are the same as in Fig. 6.

In the text
thumbnail Fig. 8

Abundances of Cr II as a function of activity index logRHK$\log R^{\prime}_{\rm{HK}}$. The symbols are the same as in Fig. 6.

In the text
thumbnail Fig. 9

Values of ξ parameterderived with the new approach as a function of the chromospheric activity index logRHK$\log R ^{\prime}_{\rm{HK}}$. The symbolsare colour-coded according to the Teff.

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.