The GAPS Programme at TNG

Context. Exoplanetary properties strongly depend on stellar properties: to know the planet with accuracy and precision it is necessary to know the star as accurately and precisely as possible


Introduction
Until around ten years ago many of the known extrasolar planets were revealed by the use of Doppler radial velocity (RV). Alone, the RVs only yield partial information on orbital elements of the planets and their minimum masses, and no insight is obtained about the planetary physical properties, like their true masses, radii, and mean densities. This additional information is available in the case of planetary transits. In recent years, major A&A 664, A161 (2022) of the parent stars and, in turn, of the planet properties. The connection between star and planet is perpetually interweaved, so one cannot be studied without accounting for the other. For instance, when the planet is transiting, it is possible to have hints of the internal planetary structure and the gas and/or ice and/or rock ratios from the bulk planetary density. However, the required accurate estimates of planetary mass and radius and density necessarily rely on precise determination of mass and radius of the hosting star. The derivation of stellar mass (and radius) is strongly connected to the effective temperature (T eff ), surface gravity (log g), and iron abundance ([Fe/H]) of the star, and is also dependent on the evolutionary tracks considered. Therefore, planetary properties are critically dependent on the properties of the hosting stars (see, e.g., Torres et al. 2012;Santos et al. 2013;Sousa et al. 2015;Maldonado et al. 2018Maldonado et al. , 2019. Moreover, several studies have pointed out the existence of correlations between the characteristics of the host stars and the properties and frequencies of their planetary systems, in particular for giant transiting planets, which are the focus of the present work. As a result, the evidence supporting the correlation between the stellar metallicity and occurrence rate of giant planets (e.g., Santos et al. 2004;Valenti & Fischer 2005) and the weakening of this correlation toward lower regimes of planetary mass (e.g., Sousa et al. 2008;Ghezzi et al. 2010) or for wide-orbit planets (Swastik et al. 2021), the connection between planet radius and stellar metallicity (Buchhave et al. 2014), the correlation between stellar metallicity and planetary heavy-metal content (e.g., Guillot et al. 2006), the trend between orbit eccentricity and star metallicity (Dawson & Murray-Clay 2013), the role of the abundances of α-elements (Robinson et al. 2006;Gonzalez 2009;Adibekyan et al. 2012cAdibekyan et al. , 2015 or of the Li depletion observed for massive planet-hosting stars (e.g., Delgado  is becoming more numerous thanks to the constant increase in new planet discoveries. Again, for the information on the hosting star's chemical composition it is important to separate the signatures left on the planet during its formation and migration from those due to the star. The various physical processes associated with forming giant planets (e.g., planet-disk interaction, planet-planet scattering, in situ formation, multiple generation of embryos) are thought to result in differences in atmospheric composition depending on the enrichment by chemical elements present at the formation site or accreted during migration (Voelkel et al. 2022, and references therein). In particular, some studies suggest that the planetary carbon-to-oxygen ratio and metallicity with respect to the corresponding stellar values could provide constraints on the original formation region of the planet with respect to the H 2 O, CO 2 , and CO snow lines, and on the time when planet migrated to its present orbit (e.g., Öberg et al. 2011). For instance, an enhanced C/O ratio of the planet compared to its host star was found to be produced when the planet formed far beyond the water snowline, predominantly by gas accretion, and then underwent a subsequent disk-free (high-eccentricity) inward migration (Madhusudhan et al. 2014). Therefore, C/O ratio in planet-host stars can provide key information about the protoplanetary disk regions in which the planet was formed, as abundances of the volatiles in the disk gas and solids are heavily affected by the disk radial temperature profile and, therefore, by the distance from the host star while planets accrete. On the contrary, other elemental ratios, like magnesium-to-silicon (Mg/Si) which governs the distribution of silicates in the protoplanetary disk, do not depend as strongly on the distance to the stars as the C/O ratio does (Thiabaud et al. 2015a). However, recent studies highlighted how the use of multiple elemental ratios involving elements with a high contrast in volatility (e.g., S/N, Si/N, S/O, N/O) can provide more detailed and robust constraints on the formation and migration history of giant planets than is possible with C/O alone (Turrini et al. 2021a,b). The same studies argued how planetary elemental ratios normalized to stellar abundances can provide more unequivocal indications and allow for a more straightforward comparison between different planets orbiting different host stars (Turrini et al. 2021a,b;Kolecki & Wang 2021).
With this in mind, the requirement for homogeneity and precision for the stellar parameters and elemental abundances becomes even more crucial; fundamental parameters of large samples of planet-hosting stars are often found in the literature as the result of analysis performed by different methodologies, resulting in an inhomogeneous census of stars with planets. All these arguments highlight how the homogeneous and precise determination of the fundamental properties and elemental abundances of the hosting stars is pivotal for achieving a comprehensive characterization of exoplanets. We therefore analyzed a sample of transiting planet host stars with the HARPS-N spectrograph at the Telescopio Nazionale Galileo (TNG) within the Global Architecture of Planetary Systems (GAPS; Covino et al. 2013) project. Our immediate aims are twofold: to apply an accurate procedure to derive stellar parameters, global properties, abundances of multiple elements, and kinematic properties of transiting planet host stars as homogeneously and precisely as possible using high-quality data, and to study the possible relationships between the astrophysical, kinematic, and chemical parameters of exoplanet host stars and the properties of their transiting planets, thus providing necessary information for future studies of their exoplanets with new facilities. We are aware that our procedure is based on non-automatic tools and therefore it is time consuming, but we think that this type of approach can be used as a benchmark analysis for interpreting the composition, the origin, and evolution of planets with current and future theoretical models and statistical studies. For instance, over the last decades, thanks to the successful photometric space missions (CoRoT, Kepler/K2, and TESS), a remarkable synergy has emerged between ground-based spectroscopy and asteroseismic techniques for the accurate determination of the fundamental parameters of exoplanet host stars (e.g., di Mauro et al. 2011;Chaplin et al. 2013). Similarly, the efforts made for a large sample of data applying automatic procedures to derive stellar parameters and iron abundances, like those presented by Sousa et al. (2021), are very useful for the statistical approach.
The outline of this paper is as follows. We first present in Sect. 2 the spectroscopic dataset. In Sect. 3 we describe the measurements of the stellar parameters and elemental abundances for 26 species. We then present our results and discuss the behavior of the elemental abundances of the stars with respect to their kinematic or global properties and with respect to planetary properties in Sect. 4. In Sect. 5 we draw our conclusions.

Stellar sample, observations, and data reduction
The stellar sample was selected within two GAPS subprograms aimed at searching for additional companions in known systems and at determining the Rossiter-McLaughlin effect in transiting systems (see, e.g., Bonomo et al. 2017). Within this sample we selected targets with spectral types from F5 to K7 (see Table A.1) and with rotational velocity (v sin i) known from the NASA Exoplanet Archive less than 10 km s −1 . This was done to avoid strong problems due to line blending or the presence of molecular lines, which must be dealt with via different procedures than those adopted in this work. In the end, we analyzed a total of 28 targets with 9.3 < ∼ V < ∼ 13.4 mag, of which 13 from the HATNet Exoplanet Survey (HAT-P; Bakos et al. 2004) in the north, 5 from the Wide Angle Search for Planets (WASP; Kane et al. 2004), 2 in common between HATNet and WASP (i.e., WASP11 or HAT-P-10A and HAT-P-30 or WASP51), 1 from the Kilodegree Extremely Little Telescope survey (KELT; Pepper et al. 2007), 2 from the XO Project (McCullough et al. 2005), 2 from of the Qatar Exoplanet Survey (Alsubai et al. 2013), and 1 from the Trans-Atlantic Exoplanet Survey (TrES; O'Donovan 2007). Twenty-seven stars host massive planets with masses (M p ) from ∼0.2 to ∼7.3 M Jup , and 1 star (HAT-P-26) hosts a Neptune-mass planet with ∼0.06 M Jup . All stars but one (namely, XO-2S, observed as a wide companion of the known planet host; Desidera et al. 2014) host transiting planets. Table A.1 lists the basic information on the final sample taken from the literature, together with some characteristics of the planets. Interestingly, the transits of some of our targets (e.g., HAT-P-12, HAT-P-26, WASP-43) will be observed within Cycle 1 of the Guaranteed Time Observations of the James Webb Space Telescope 1 (JWST).
Observations were performed between the end of 2012 and 2016 with the high-resolution HARPS-N at TNG spectrograph (R ∼ 115 000, λ ∼ 3900−6900 Å; Cosentino et al. 2012). Solar spectra were also performed through observations of the asteroid Vesta. Spectra reduction was performed using the HARPS-N instrument Data Reduction Software pipeline (see the radial velocity curves of the sample in Bonomo et al. 2017). The standard steps for data reduction and the appropriate crosscorrelation masks were applied to each target. The final high signal-to-noise ratios (S NR) merged spectra for each target star and for the Sun were obtained by co-adding, after the proper shift, to the rest frame by the corresponding RV. All individual spectra of the given star and of Vesta reached a S NR (per pixel at λ ∼ 6000 Å) between 100 and 300 for the targets (except for Qatar-2, the faintest one in the sample, with a mean S NR around 50) and S NR ∼ 300 for the solar spectrum.
Due to the importance of deriving accurate oxygen abundances for as many targets as possible (see, e.g., discussions in Sects. 4. 6-4.9), and since very useful oxygen lines like those of the O I triplet at ∼777 nm are not present within the HARPS-N spectral range, we decided to search for processed FEROS at ESO spectra (R ∼ 48 000, λ ∼ 3600−9000 Å; Kaufer et al. 1999) available in the ESO archive. In the end, we found FEROS spectra for the Sun and seven stars: HAT-P-30, WASP-54, HAT-P-17, HAT-P-26, HAT-P-20, Qatar 2, WASP-43. The typical S NR of these spectra was greater than ∼50 at λ ∼ 6000 Å. These spectra also allowed us to derive nitrogen abundances from high-excitation permitted lines (see Sect. 3.2). For homogeneity reasons, we verified that for the targets observed with FEROS and using the same method applied for HARPS-N spectra based on the spectral synthesis of the oxygen λ6300.3 Å line and the CN molecule (see Sect. 3.2), we derived O and N abundances very close to those obtained with HARPS-N spectra (within the uncertainties). On the other hand, we verified that for the two targets for which we could measure oxygen and nitrogen both using FEROS and HARPS-N diagnostics we obtained very similar results, within the errors.
We note that for some of the targets within the GAPS project we performed a preliminary analysis of some astrophysical parameters (effective temperature, surface gravity, and iron abundance), but it was non-homogeneous and based on previous version of tools, model atmospheres, and line lists (see Covino et al. 2013;Desidera et al. 2014;Esposito et al. 2014Esposito et al. , 2017Damasso et al. 2015a,b;Biazzo et al. 2015;Sozzetti et al. 2015;Mancini et al. 2015Mancini et al. , 2018Mancini et al. , 2022. We note here that the analysis performed in this work is absolutely new, and innovative, homogeneous, and aimed at a characterization of our sample of stars that is as broad as possible.

Stellar parameters and iron abundance
We derived stellar parameters and iron abundances using the code MOOG (Sneden 1973;version 2017) with the driver abfind, which assumes local thermodynamic equilibrium (LTE) and treating the radiative and Stark broadenings in a standard way. For collisional broadening, we used the Barklem et al. (2000) prescriptions for damping values. We used plane-parallel model atmospheres linearly interpolated from the ATLAS9 grids of Castelli & Kurucz (2003), with solar-scaled chemical composition and new opacities (ODFNEW).
Effective temperature T eff , surface gravity log g, microturbulence velocity ξ, and iron abundance [Fe/H] were measured through a method based on equivalent widths (EWs) of iron lines. We adopted the list of iron lines by Biazzo et al. (2012) and corrected the atomic parameters of the Fe II line at λ6516.08 Å with the most recent NIST values (National Institute of Standards and Technology; Kramida 2019) because they led to the most reliable solar iron abundance measurements. The iron line list was built to minimize potential correlations between the atmospheric parameters, by including lines of different strengths at a given excitation potential and by having a wide distribution of excitation potentials. In addition, we keep in our line list only the iron lines that could be reliably measured at our spectral resolution. In the end, a total of 82 Fe I+11 Fe II lines were considered. The EWs of the target stars were measured by means of a direct integration or Gaussian fitting procedure using the IRAF 2 SPLOT task. We discarded strong lines (EW > 150 mÅ) and those lines with fitting errors larger than 2σ. Each line equivalent width was controlled and measured several times and particular attention was paid to the continuum tracing. The continuum placement of each stellar line was defined looking at the continuum position of the same line in the solar spectrum and the same criteria (both for the continuum definition and the intervals selected for the integration) were adopted. The values of T eff and ξ were determined by imposing the condition that the Fe I abundance does not depend on the excitation potential and equivalent width of the lines, while log g was obtained from the Fe I/Fe II ionization equilibrium. We used an iterative procedure by changing the parameters at steps of 5 K in T eff , 0.01 km s −1 in ξ, and 0.01 dex in log g, and requiring that the slope of the Fe I abundance with respect to the excitation potential (for T eff ) or EW (for ξ) was close to zero and that the difference between the mean iron abundance obtained from the Fe I and Fe II lines was lower than 0.01 dex (for log g). We therefore derived stellar parameters with internal accuracy (at 3σ) ranging from 15 to 90 K in T eff , from 0.09 to 0.19 dex in log g, from 0.02 to 0.41 km s −1 in ξ, and from 0.07 to 0.15 dex in [Fe/H] (see Table 1). As a sanity Table 1. Final stellar parameters (effective temperature, surface gravity, microturbulence velocity, rotational velocity), and iron and lithium abundances derived in this work (for Fe I and Fe II the number of lines used to derive the abundances is in parentheses), together with macroturbulence velocity V macro and mean radial velocity V rad .

Abundance of other elements and rotational velocity
For elements other than iron, we applied spectroscopic techniques based on line EWs and spectral synthesis depending of the element. We also measured the rotational velocity through spectral synthesis.

Analysis based on equivalent widths
Once the stellar parameters and iron abundances were measured, we computed the abundances of the other elements ([X/H] 3 ) using the MOOG code (Sneden 1973;version 2017) and the drivers abfind and blends for the treatment of the lines without and with hyperfine structure (HFS). In particular, in addition to Fe we computed the abundance of 25 elements: Li, C, N, O, Na, Mg, Al, Si, S, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, Zr, Ba, La, Nd, and Eu. As was done for Fe, we measured the abundances of two ionization states also for Ti and Cr, while for the other elements only the abundance of one species (first or second ionization state) was measured (see Table A .2).
As was done for the iron lines, the EW of each spectral line was measured on the one-dimensional spectra interactively using the splot task in IRAF. The location of the local continuum was carefully selected tracing as closely as possible the same position for the spectral line of each star and the asteroid Vesta. This was done with the aim of minimizing the error on the selection of the continuum. We also excluded features affected by telluric absorption.
Following the prescriptions by Biazzo et al. (2015), we started from the line list of Na I, Mg I, Al I, Si I, Ca I, Ti I, Ti II, Cr I, Cr II, Ni I, and Zn I by Biazzo et al. (2012) complemented with additional lines and atomic parameters for Na I, Al I, Si I, Ti I, Ti II, Cr I, Ni I, and Zn I taken from Schuler et al. (2011) and Sozzetti et al. (2006). For the Mg I line at λ4730 Å and the Al I line at λ6698.67 Å we considered the NIST (Kramida 2019) and the Meléndez et al. (2014) atomic parameters. In the cases of C I, S I, Sc II, V I, Mn I, Co I, and Cu I, we considered the line lists by Kurucz (1993), Schuler et al. (2011), Kramida (2019), Johnson et al. (2006), and Scott et al. (2015), where the hyperfine structure by Johnson et al. (2006) and Kurucz (1993) was adopted for Sc, V, Mn, Cu, and Co. Solar isotopic ratios by Anders & Grevesse (1989) were considered for Cu (i.e., 69.17% for 63 Cu and 30.83% for 65 Cu). For the s-process elements Y II and Zr II (first peak), Ba II and La II (second peak), the mixed s/r-process element Nd II, and the r-process element Eu II we considered the line lists by Johnson et al. (2006), Ljung et al. (2006), Prochaska et al. (2000), Lawler et al. (2001b), andDen Hartog et al. (2003), where the HFS by Gallagher et al. (2010) and Lawler et al. (2001a) was adopted for Ba and La, respectively. Solar isotopic ratios by Anders & Grevesse (1989)  For the targets observed also with FEROS, we measured oxygen and nitrogen abundances through line EWs. The oxygen abundance was estimated using the O I triplet of permitted lines at λ7771.94, λ7774.17, λ7775.39 Å with atomic parameters by NIST (Kramida 2019), and considering the non-LTE (NLTE) corrections by Amarsi et al. (2015). To our knowledge, no NLTE correction are present for [Fe/H] > 0 and for T eff < 5000 K, and the effect is important for T eff differences also of 50 K (while it is negligible for [Fe/H] differences within 0.4 dex). This is why we considered the O I abundances measured for the Sun and the stars WASP-54, HAT-P-17, HAT-P-26, and the slightly metal-rich star HAT-P-30, while we excluded from the analysis the other (cooler) targets observed with FEROS (i.e., HAT-P-20, Qatar-2, and WASP-43). The importance of these effects has been also demonstrated for hot exoplanetary atmospheres (Borsa et al. 2021).
Thanks to the FEROS spectra, we also measured nitrogen abundance using three high-excitation permitted lines (e.g., λ7442.3, λ7468.3, λ8216.3 Å), considering the atomic parameters by Caffau et al. (2009). Stars with T eff < ∼ 5200 K are too cool to have detectable N I lines, while within the warmer targets observed with FEROS only the spectra of Vesta and HAT-P-30 had sufficient S NR to measure the nitrogen abundance. No correction for NLTE effects was applied for the N abundance of our sample because to our knowledge it is only available for the Sun (see Caffau et al. 2009).
In addition to O I, due to the relatively wide range in the stellar parameters of our targets (mainly T eff and [Fe/H]), we also applied NLTE corrections to each line abundance of C I, Na I, Mn I, and Co I, following the prescriptions given by Amarsi et al. (2019), Lind et al. (2011), Bergemann & Gehren (2008, and Bergemann et al. (2010), respectively. For the other elements, no NLTE departure was considered because corrections were negligible or not reported in the literature for the lines used in this work and for the T eff and [Fe/H] ranges of our targets. As for the iron lines, the final elemental abundance of each target was plotted as a function of stellar T eff , log g, ξ, and v sin i to look for possible trends, for example due to line blendings. When trends were found for a specific element (in particular at lower and higher T eff ), each line of that element was plotted as a function of T eff , log g, ξ, and v sin i to recognize the presence of effects due mainly to line blending or the bad quality of the spectrum for that specific line, which was then removed. In the end, no trend was found for the final abundances and stellar parameters, with the exception of Cr II and C I for which a residual trend with effective temperature remained. These two elements are discussed in the next paragraph.

Analysis based on spectral synthesis
The projected rotational velocity (v sin i) and the lithium abundance (log n(Li)) were measured by applying the spectral synthesis technique. We used the synth driver within the MOOG code (Sneden 1973;version 2017) and considered synthetic spectra obtained from the Castelli & Kurucz (2003) grids of model atmospheres at the stellar parameters (T eff , log g, ξ, [Fe/H]) derived in Sect. 3.1. We applied the same method to derive the abundances of carbon, nitrogen, and oxygen, which are also elements analyzed through line EWs. In particular, we convolved the synthetic spectra with a Gaussian profile corresponding to the resolution of HARPS-N of R ∼ 115 000, taking into account the optical limb-darkening coefficients by Claret (2019) at the stellar T eff , log g, ξ, and [Fe/H]. Moreover, the empirical relationships obtained using the asteroseismic rotational velocities by Doyle et al. (2014) from the Kepler data were used to derive the macroturbulence velocity (V macro ) for T eff > 5700 K, while empirical relationships by Brewer et al. (2016) were considered for T eff < 5700 K 4 .
For the v sin i we synthesized spectral lines in two regions around 6200 and 6700 Å, as done in Barbato et al. (2019), until the minimum of residuals between stellar and synthetic spectra were reached. The final values are listed in Table 1, where the errors take into account uncertainties in stellar parameters (T eff , log g, ξ, [Fe/H]) and spectral continuum definition (for the uncertainties in the continuum position, see next paragraph).
The forbidden [O I] line at 6300.3 Å was considered for deriving the O abundance through spectral synthesis and using the above-mentioned code and model atmospheres. Atomic parameters for the O line and the nearby (blended) Ni I line were taken from Caffau et al. (2008) and Johansson et al. (2003), respectively (see also Bertran de Lis et al. 2015, and references therein). Nickel abundances were fixed to the values we derived above. No NLTE correction were made because the λ6300.3 Å line is not affected by deviation from LTE (Caffau et al. 2008).
Then, adopting the estimated O abundance, C was measured from 12 CH and 13 CH molecular features around 4320 Å (see the example in Fig. 1). Molecular parameters were taken from Masseron et al. (2014). Atomic lines were re-adjusted on the solar spectrum. The C abundances from these CH bands were measured for all targets but two cool stars (Qatar-2 and WASP-43), mainly because of the very critical continuum placement.
Finally, the N abundance was obtained from the CN features at ∼4215 Å, adopting the measured C and O abundances derived through spectral synthesis. The fitting procedure was differential  (i.e., adopting the same spectral range and features in each star), and it was repeated until the minimum of residuals was obtained. The line lists for the 12 C 14 N, 13 C 14 N, and 12 C 15 N isotopologues were taken from Sneden et al. (2014). Given the sensitivity to T eff , the N abundance from the CN band was measured for stars in the 4600−5300 K T eff range.
Since CNO abundances were derived with both EWs and spectral synthesis, as final abundance we considered the weighted average coming from these two methods. In particular, Fig. 2 shows the C abundances derived from C I lines (open diamonds) and those derived from the CH band (filled circles). We find a good agreement between the [C/H] measured from these two diagnostics for T eff > 5000 K, with the largest differences for cooler stars, for which the atomic carbon lines provide abundances larger than twice the typical error on [C/H]. Similar results were found by Baratella et al. (2020) andDelgado Mena et al. (2021), among others. The latter authors mentioned a possible dependence of the trend with the metallicity (slight increase in [C/Fe] toward lower metallicity for [C/H] derived through atomic lines, and a flatter trend for abundances obtained from molecular lines). We do not find similar trends, most probably due to our small stellar sample. We find instead an evident trend with the Mittag chromospheric activity index log R HK derived by Claudi et al. (in prep). Similar results were obtained by Baratella et al. (2020), who also found a positive correlation between atomic [C/H] and log R HK , justified as possible unknown blends in the optical lines becoming more important in active stars than in quiet stars. On the contrary, we do not find any trend between log R HK and other abundances or ξ, thus further validating our method of deriving elemental abundances (see also discussion in Baratella et al. 2020). In the end, we decided to consider a weighted average between C abundances derived with the atomic and the CH band for T eff > 5000 K, while for T eff < 5000 K we considered only the values from the CH band as C abundances 5 .
The lithium abundance (log n(Li)) was derived through spectral synthesis of the absorption line at λ6707.8 Å, which was present in seven stars (KELT-6, HAT-P-14, WASP-38, HAT-P-30, WASP-13, HAT-P-4, HAT-P-3; see Sect. 4.5). In particular, we used the lithium line list by Reddy et al. (2002) in the vicinity (±0.5 Å) of the Li 6707.8 Å line, implemented with the VALD database (Kupka et al. 2000) for wavelengths farther from the line center. The line list by Reddy et al. (2002) considers the isotopes 6 Li and 7 Li of the λ6707.8 Å line allowing us to estimate the lithium isotopic 6 Li/ 7 Li ratio. The lithium abundance derived through the MOOG code was then corrected for the departure from LTE considering the non-LTE calculation of Lind et al. (2009). A double-check of the final results was done computing the lithium abundances also using the Li equivalent widths and converting them in abundances through the curves of growth by Lind et al. (2009). We found consistent values with both methods within the uncertainties. In Table 1 we list only the results obtained through spectral synthesis. The errors on log n(Li) were derived by adding quadratically the uncertainties due to the abundance measurements and those due to stellar parameters.
The results for all solar elemental abundances are given in Table 2 together with those given by Asplund et al. (2021) and obtained using a 3D radiative-hydrodynamic model of the solar surface convection and atmosphere, and correction for departures from LTE conditions when necessary. This table shows good agreement between our solar abundance values and the literature results, with a mean difference between values derived in this work and in the literature of 0.00 ± 0.04 dex.

Uncertainties in elemental abundances
Derived abundances are mainly affected by uncertainties in atomic parameters, stellar parameters, and measured EWs (or continuum position when spectral synthesis was applied).
Uncertainties in atomic parameters, such as the transition probability (log g f ), should cancel out since our analysis is carried out differentially with respect to the Sun. Errors due to uncertainties in stellar parameters (T eff , ξ, log g, [Fe/H]) were estimated first by assessing errors on stellar parameters themselves and then by varying each parameter separately, while keeping the other two unchanged. As shown in Table 1 and explained in Sect. 3.1, the uncertainties in stellar parameters are Notes. Z represents the atomic number and T cond the 50% condensation temperature values derived by Lodders (2003).
in the range 15-90 K for T eff (with a standard deviation σ of 25 K), 0.09-0.19 dex for log g (σ = 0.03 dex), 0.02-0.41 km s −1 for ξ (σ = 0.14 km s −1 ), and 0.07-0.15 dex for [Fe/H] (σ = 0.02 dex). Due to the small values of the standard deviations across the wide range of stellar effective temperature for the errors on T eff , log g, and [Fe/H], and due to the smaller influence of errors on ξ for almost all elements, when compared to the other sources of uncertainties, we decided to consider typical uncertainties in stellar parameters of 60 K, 0.10 dex, 0.15 km s −1 , and ∼0.08 dex in T eff , log g, ξ, and [Fe/H], respectively. We thus assumed these values as uncertainties in stellar parameters. Errors on abundances [X/H] due to uncertainties in stellar parameters are summarized in Table 3 for two different stars in our sample from which we derived all elemental abundances: one of the coolest stars (HAT-P-12 or WASP-10) and one of the warmest stars (HAT-P-30). Total errors on log n(Li) are instead reported in Table 1.
The errors due to uncertainties in EWs are well represented by the standard deviation around the mean abundance determined from all the lines. These errors are listed in Tables 1 and A.2, where the uncertainties in [X/H] were obtained by quadratically adding the error for the target and the error for the Sun. When only one line was measured, the error on [X/H] is the standard deviation of three independent EW measurements obtained taking different positions of the continuum. The number of lines employed for the abundance analysis is listed in Tables 1 and A.2 in brackets.
For the uncertainties due to the definition of continuum position when spectral synthesis was applied, random errors affecting our best-fit procedure were evaluated by changing the continuum position until the standard deviation (observed minus synthetic spectra) was two times larger than the best-fit value, where residuals of our best-fit solutions are typically smaller than 0.02. These error budgets are listed in Table A Figure 3 shows the position of our targets in the Kiel T efflog g diagram along with stellar model tracks for three different metallicities (i.e., Z = 0.5 Z , 1.5 Z , 2.5 Z ) and ages (i.e., log(Age/yr) = 8. 7, 9.3, 10.1) spanning the values of our stellar sample (see Sect. 3.1 for our iron abundance measurement and Sect. 4.4 for our age estimation using chemical clocks). The stars are color-coded with respect to [Fe/H], from approximately −0.3 dex up to +0.4 dex, as shown in the histogram. Our homogeneous procedure led to stellar parameters agreeing with theoretical tracks within the uncertainties. Some evidence that our log g is slightly greater than that derived from the tracks at a mean level of ∼0.1 dex seems to be present for stars cooler than ∼4900 K, but this is compatible with the internal error on log g and is probably due to the smaller number of Fe II lines for lower T eff . Similar findings were obtained by Brucalassi et al. (2021) and Magrini et al. (2022) for targets of the ARIEL mission. A161, page 7 of 30 A&A 664, A161 (2022) Table 3. Internal errors on abundance determination due to uncertainties in stellar parameters for one of the coolest stars (HAT-P-12 for all elements but C and N; WASP-10 for C and N) and for one of the warmest stars (HAT-P-30) in our sample.

Final parameters and comparison with previous works
HAT-P-12  Table 4. Mean difference and standard deviation between the values of the atmospheric parameters derived in the literature and those obtained in this work (n is the number of stars in common).
−25 ± 82 0.06 ± 0.07 −0.04 ± 0.11 0.03 ± 0.33 ( * ) −0.1 ± 0.7 28 Discovery papers ( ) 20 ± 60 0.06 ± 0.06 0.00 ± 0.08 ...  Table 4. In particular, we consider the comparisons with the values of the exoplanet discovery papers derived through different methods, the stellar parameters obtained by Torres et al. (2012) averaging three methods based on spectral synthesis, line equivalent widths, and crosscorrelation against a library of spectra, the values by Brewer et al. (2016) derived using spectral synthesis, and those listed in the SWEET-Cat (Stars With ExoplanETs Catalogue; Santos et al. 2013;Sousa et al. 2021). SWEET-Cat, for the first time described in Santos et al. (2013), is a continuously updated catalogue of stellar parameters for planet-hosting stars derived, whenever possible, using the same methodology. A very recent version of the catalogue, whose parameters were re-derived using better quality spectra and following the same homogeneous procedure, is available (see Sousa et al. 2021). In the new version of the catalogue, recent Gaia EDR3 parallaxes were also considered to derive accurate surface gravity of the host stars.
As a result of the comparison of our T eff , log g, and [Fe/H], and v sin i with the values given by the literature and listed in Table 4, the standard deviation of the difference is 60−90 K, 0.07−0.23 dex, 0.04−0.11 dex, and ∼0.5−0.7 km s −1 , respectively. For the microturbulence velocity, small differences are present among the comparisons with the two versions of SWEET-Cat (see Table 4), while within the discovery papers, most values of this parameter were fixed, and thus we cannot make a meaningful comparison; however, we can consider the relationship by Adibekyan et al. (2012a), which is based on the dependence of ξ on T eff and log g for stars with 4500−6500 K, 3.0−5.0 in log g, and −1.4 < [Fe/H] < 0.5 dex. We find a very good agreement between the values computed through the mentioned relationship and those derived through our MOOG analysis, with mean difference of 0.10 ± 0.18 km s −1 .
Surface gravities were also derived by Mortier et al. (2013) and Sousa et al. (2021) using alternative methods to those based on spectroscopy. In particular, considering the surface gravities derived by Mortier et al. (2013) for 11 targets in common with our sample using a method based on photometric light curves and the previous version of the SWEET-Cat, the difference with respect to our values is of 0.04 ± 0.11 dex (see the inset in the second panel of Fig. B.1), while considering the log g derived by Sousa et al. (2021) through Gaia parallaxes, the standard deviations strongly improve with respect to the spectroscopic values listed by the authors, becoming 0.09 ± 0.09 dex. Therefore, our analysis based on a careful method leads to log g results that are very close to the values derived in accurate ways (through transit light curves or Gaia parallaxes), even if we are aware that the kind of analysis performed in this work is very time-consuming and cannot be easily applied to surveys of hundreds to thousands of targets. A similar method to derive trigonometric log g is applied in Brucalassi et al. (2021) for ∼150 targets within the ARIEL reference sample (see also Tinetti et al. 2021;Danielski et al. 2021;Magrini et al. 2022). Due to the high quality of the Gaia parallax, the authors suggest adopting trigonometric log g as a viable possibility for big stellar samples for which some spectroscopic methods based on automatic tools tend to under-or overestimate the surface gravity at low and high temperatures.
Concerning the iron abundances, from Fig. B.1 we recognize a presence of most discarding values when the comparison is made with respect to the discovery papers. In particular, the targets with a difference of more than 1σ are those from which the microturbulence velocity was not derived or was fixed (i.e., XO-2S, TRES-4, Qatar-2, HAT-P-20) or for part of the WASP sample. For this subsample of WASP targets (namely, WASP-38, WASP-60, WASP-54, WASP-39), T eff and log g were derived via a method based on different strong lines; by chance, this method seems to lead to lower T eff , and therefore lower [Fe/H] than our values. We do not go deeper into detail because we are aware that these analyses date back to more than ten years ago. Here, we only note that the most recent analysis points toward findings that are more similar to our results. Mortier et al. (2013) and Brewer et al. (2016) respectively measured abundances of 12 and 14 elements other than iron. For the 11 targets in common with Mortier et al. (2013), we find a mean difference in [X/H] of −0.01 ± 0.07 dex for Al, Ca, Co, Cr, Cr II, Mg, Mn, Na, Ni, Sc II, Si, Ti, Ti II, and V, while for the five targets in common with Brewer et al. (2016)  Comparing our results with those achieved by these authors, we find mean differences of 0.03 ± 0.09, −0.04 ± 0.08, and 0.00 ± 0.07 dex with respect to Burke et al. (2007) and Teske et al. (2014Teske et al. ( , 2019, respectively. Therefore, a general agreement between our abundance values and those from the literature is obtained for the targets in common.  chronological indicator of nucleosynthesis. In Fig. 4 Adibekyan et al. 2012c). Here we cannot conclude whether the observed dispersions are intrinsic or due, at least in part, to statistical or methodological reasons (e.g., number of lines, S NR). We only highlight how for all elements the general pattern of [X/H] and [X/Fe] with the iron abundance is similar to those of the sample of dwarf stars in the Galactic disk, without clear evidence of peculiar behavior for our planet-host stars when compared with field stars. Therefore, the [X/H]-[Fe/H] trends of our sample seem to reflect the Galactic chemical evolution in the solar neighborhood.

Elemental abundance versus condensation temperature
In Fig. 5 we plot for each star the elemental abundance [X/H] as a function of the condensation temperature T cond , from the warmest target to the coolest. It was indeed reported that stars hosting high-mass planets are expected to be more enriched in refractory elements (i.e., elements with high condensation A161, page 10 of 30 temperatures) and deficient in low T cond volatile elements (e.g., Smith et al. 2001). This should happen because any accretion event, occurring very close to the star (i.e., in high-temperature environments) would add refractory elements, which condense at high T cond with respect to volatiles (see, e.g., Sozzetti et al. 2006, and references therein). Similarly, Meléndez et al. (2009) concluded that solar twins without close-in giant planets chemically resemble the Sun, with depletion of refractory elements relative to the volatiles, suggesting that the presence of such planets might prevent the formation of an Earth-sized planet. Revealing these trends requires developing a very accurate differential analyses that is precise for binary systems (e.g., Gratton et Biazzo et al. 2015;Maia et al. 2019) and members of stellar clusters (e.g., Yong et al. 2013) or through the use of a set of comparison stars (Liu et al. 2020;Tautvaišienė et al. 2022) because many observational uncertainties can be considered common-mode effects. If no comparison target is available, an unambiguous explanation for these trends is difficult to reach because they could also reflect the wide diversity of exoplanetary systems, as well as a variety of scenarios which could occur within the circumstellar disk (Spina et al. 2016b), or could be associated with the correlation of elemental abundances with the age and birthplace in the Galaxy (e.g., Adibekyan et al. 2014). González Hernández et al. (2013) have found that after removing the Galactic chemical evolution effects from a sample of A161, page 11 of 30 A&A 664, A161 (2022) main-sequence objects, stars with and without planets show similar mean abundance patterns. We therefore applied an approach similar to the latter authors and adopted Spina et al. (2016a) for the correction of Galactic chemical evolution (GCE). In Fig. 5 we show the [X/H] abundances as a function of T cond before (filled dots and dashed lines) and after (open squares and dotted lines) the removal of the GCE effects. The Spearman (ρ) and Kendall (τ) statistical significance computed with IDL 6 for the linear regression analysis before and after the GCE removal are displayed respectively in the bottom right and bottom left corners of the plots. Excluding Qatar-2 and WASP-43 for which we could not measure the abundances of volatile elements, among the other 26 targets there are 10 showing a trend after the GCE removal, with both ρ and τ lower than 0.05, and a wide range of planetary masses (from 0.06 M J to 2.6 M J ) and effective temperatures (from 4650 to 6715 K). In particular, HAT-P-12 and HAT-P-26, respectively hosting planets with M p ∼ 0.2 M J and 0.06 M J , show the most pronounced significant trend at a level of ρ ∼ τ ∼ 2−3 × 10 −5 up to 2−5 × 10 −4 and among the highest values of positive slopes, while WASP-60 (with a planet of 0.5 M J ) is the only one within this subsample of ten targets showing a negative trend (i.e., decreasing refractory-to-volatile abundance ratios). On the one hand, we note that HAT-P-26 and HAT-P-12 are also targets with kinematic properties consistent with transition thin-thick disk (see Sect. 4.3) and with high values of [α/Fe] and Mg/Si (see Sect. 4.7), which, together with the possible [X/H]-T cond trend, could indicate a pattern in the formation of the nuclei of their planets. On the other hand, we note that the binary components XO-2N and XO-2S, here analyzed as single stars, do not show a reliable correlation between elemental abundances and T cond , while precise and accurate differential analysis demonstrated that the difference in elemental abundance between these two binary components shows a trend with the condensations temperature indicating possible ingestion of material by XO-2N or depletion in XO-2S (see Teske et al. 2015;Biazzo et al. 2015). Finally, excluding Qatar-2 and WASP-43 for which we could not measure elemental abundance of volatiles, we find higher values of T cond -[X/H] slopes for cooler and older stars (see Sect. 4.4) with higher log g, regardless of planetary mass and not evident without the GCE removal. We also find tentative evidence that stars with smaller galactocentric distance and greater Galactic eccentricity (see Sect. 4.3) have steeper slopes. Similar results were also found by other authors (e.g., Adibekyan et al. 2014;Tautvaišienė et al. 2022). In other words, stellar parameters and Galactic position are crucial to establishing the stellar chemical pattern of the stars and only a strictly differential analysis (like those performed for binary stars) can remove spurious trends and help to draw definitive conclusions, which is not the case in this work.

Chemical and kinematic properties
To study possible chemokinematic peculiarities of planethosting stars, we computed the stellar Galactic space velocities. The space velocity components UVW were derived with respect to the local standard of rest (LSR), correcting for the solar motion derived by Coşkunoǧlu et al. (2011): (U , V , W ) = (−8.50, 13.38, 6.49) km s −1 . Parallaxes (π) and propers motions (µ α , µ δ ) were taken from Gaia EDR3 (Gaia Collaboration 2016, 2021), mean radial velocities ( V rad ) were obtained from the HARPS-N spectra, and ICRS coordinates at epoch=2000 were taken from the SIMBAD Astronomical Database. We considered the general outline of Johnson & Soderblom (1987) in a left-hand coordinate system (i.e., with U pointing toward the Galactic anti-center, V toward the local direction of rotation in the plane of the Galaxy, and W toward the north Galactic pole). The uncertainties were obtained considering the prescriptions by Gagné et al. (2014), thus we used the full covariance matrix taking into account the error contributions of V rad , π, µ α , and µ δ . Combining the measurement errors on parallaxes, proper motions, and radial velocities, the resulting average uncertainties in the U, V, W velocities are about 0.15 km s −1 . The values derived for each target are listed in Table A.3 and the Boettlinger diagram in the (U, V) plane is shown in the left panel of Fig. 6, where the boundary separating the young disk (YD) and the old disk (OD) stars according to Eggen (1996) is displayed with a solid line. The YD locus contains associations younger than ∼1 Gyr (see Gagné et al. 2018). Nine targets, namely HAT-P-14, TRES-4, WASP-54, WASP-13, HAT-P-3, WASP-11, HAT-P-18, WASP-10, and HAT-P-20, seem to belong to the young disk, and three others (WASP-43, HAT-P-4, WASP-38) are very close to the YD boundary. In the right panel of the same figure we show the Toomre diagram, which is a representation of the combined vertical and radial kinetic energies versus the rotational energy. Low-velocity stars, within a total velocity v tot = (U 2 LSR + V 2 LSR + W 2 LSR ) 1/2 = 50 km s −1 are, to a first approximation, mainly thin disk stars, while stars with 70 < ∼ v tot < ∼ 180 km s −1 are likely to be thick disk stars (see Bensby et al. 2014, and references therein). Eleven of our targets, namely HAT-P-14, WASP-38, HAT-P-29, WASP-54, WASP-13, HAT-P-3, WASP-11, HAT-P-18, WASP-10, HAT-P-20, and WASP-43, are very close to the Sun, with v tot < ∼ 20 km s −1 , while five exoplanet-hosting stars (XO-2N, XO-2S, Qatar-2, HAT-P-26, HAT-P-12) have v tot > 70 km s −1 , compatible with thick disk stars.
We also calculated the thick-to-thin disk probability ratios. In particular, we considered the prescriptions of Bensby et al. (2014) for the Gaussian distributions of random velocities of the different stellar populations. To get the probability D and T D for the thin and thick disk that a given star belongs to a specific population, we considered the asymmetric drift, the velocity dispersion, and the fractions of each population listed in their Table A.1. By then dividing the thick disk probability with the thin disk probability, we get the probability for the thick disk-tothin disk (T D/D) membership. Bensby et al. (2014) require that for a star to be a candidate thick disk star its probability must be at least twice that of being a thin disk star (i.e., T D/D > 2), and vice versa for a candidate thin disk star T D/D < 0.5. All our targets show T D/D < 0.5 with the exception of five stars: two targets with probability ratios between the thin and thick disks (i.e., HAT-P-21 and HAT-P-26) with T D/D ∼ 0.5; three targets with probabilities consistent with the thick disk, namely the star HAT-P-12 with a value of ∼3.7, and the XO-2 binary system with T D/D ∼ 13. All targets within or near the YD boundary show T D/D close to zero. The same occurs for targets with v tot < ∼ 20 km s −1 . Following Bensby et al. (2014), with the aim to further investigate the chemokinematic properties of our sample, we show in Fig. 7 Bensby et al. (2014) for dwarf thin and thick disk stars in the solar neighborhood. We also note that within the four stars hosting planets with masses smaller than 5 M Nep (i.e., XO-2S, HAT-P-26, HAT-P-18, HAT-P-12) all but HAT-P-18 show v tot > 20 km s −1 . Moreover, their mean distance from the Sun is smaller (by more than 60 pc) than that of stars hosting planets with higher masses. Similar results were found by Adibekyan et al. (2012b), who noted that, as expected, low-mass planets are easier to find at smaller distances due to the higher apparent magnitudes of their hosts. We therefore used a similar approach to that proposed by these authors. In particular, we calculated the maximum height a star can reach above the Galactic plane (Z max ), the current Galactic eccentricity (e G ), the peri-/apo-(R peri , R apo ) center radii of an orbit, and the galactocentric distance (R GC ) with the galpy 7 package, a python package for Galactic-dynamics calculations. We assumed the built-in model MWpotential2014 for the Milky Way's gravitational potential (see Bovy 2015). We set the distance of the Sun from the Galactic center to R 0 = 8.0 kpc (Bovy et al. 2012) and its height above the plane to z 0 = 0.025 kpc (Jurić et al. 2008), and used the parallaxes and proper motions from Gaia EDR3 to transform the celestial coordinates in galactocentric radius and height above the Galactic plane. We also used our mean radial velocities from HARPS-N to obtain the orbital parameters and considered the chemical ages as derived in Sect. 4.4. Table A.3, together with UVW velocities and T D/D probabilities, also lists the output results from the galpy package.
First, we note that all targets with T D/D > ∼ 0.5 also show |∆(R mean − R GC )| > ∼ 1, where R mean is the mean position of the stellar Galactic orbits (calculated as the average of R peri and R apo ). This could indicate that stars with R mean very different from R GC could have experienced greater migration and therefore could have a higher probability of belonging to the thin-thick disk transition or thick disk (see Magrini et al. 2022; and references therein). Then, we find that on average stars hosting planets with masses lower than 5 M Nep have higher Galactic eccentricity (by ∼0.2) and greater age (by ∼5 Gyr), while the difference in Z max of only ∼45 pc is within the average uncertainties (of around 150 pc). These low-mass planet-hosting stars also have 7 http://github.com/jobovy/galpy A161, page 13 of 30 A&A 664, A161 (2022) lower V LSR and W LSR space velocity components than stars hosting higher-mass planets. Moreover, comparing their [Mg/Fe], [Si/Fe], and [Ti/Fe] ratios, we find that the stars with M p < 5 M Nep planets are more enhanced by Mg (with a difference of 0.09 dex), by Si (difference of 0.04 dex), and Ti (difference of 0.09 dex) when compared to stars with higher-mass planets. These elements are all tracers of the rocky component of the cores of the Neptune-mass planets, thus possibly indicating that Neptunian and/or super-Earth host stars tend to belong to the thicker disk when compared with Jupiter-mass hosting stars. Even if we are cautious about these findings because they could depend on several selection effects (e.g., stellar magnetic activity), our results seem to give support to the findings by Adibekyan et al. (2012b), according to which stars hosting low-mass planets tend to belong to a thicker disk.
In Fig. 8  of 20 km s −1 in the Toomre diagram, and the T D/D ratio is very low (∼0.02), thus appearing to belong to the thin disk (see Fig. 6). HAT-P-12 shows a relatively high proper motion, it is compatible with thick disk stars in the Toomre diagram, and T D/D ∼ 3.7 (i.e., slightly higher than the thick disk threshold). In conclusion, our chemokinematic analysis for these two targets seems to be consistent with KELT-6 clearly belonging to the thin disk and HAT-P-12 between the thin and thick disk. Moreover, our findings support the conclusion by Adibekyan et al. (2012b), according to which planet-hosting metal-poor stars (like KELT-6) can have high [α/Fe] even belonging to the thin disk.
Furthermore, we note that the only target within our sample hosting a Neptunian and/or super-Earth planet, namely the solarmetallicity star HAT-P-26, shows the highest value of [α/Fe], when compared to the sample by Adibekyan et al. (2012b) with similar planetary masses and iron abundance (triangles in Fig. 8). The relatively high abundance values of α elements for this target are similar to those of their Neptunian-super-Earth host stars with iron abundance values around −0.3 dex. Its position in Fig. 8 is indeed consistent with the stars with high-α content, as defined by Adibekyan et al. (2012c). Again, this target shows high space velocities (see Fig. 6) and in the Toomre diagram, accordingly to Bensby et al. (2014), it is in the region populated by thick disk stars. Finally, its Z max is around 565 pc (one of the highest in the sample) and the T D/D ratio is ∼0.5. Therefore, the chemokinematic analysis for this target seems to be consistent with a star close to the thin-thick disk transition, also supported by its relatively old age (see Sect. 4.4).
Two other targets with α content higher than the limit defined by Adibekyan et al. (2012c) are HAT-P-22 and HAT-P-3. The kinematic position in the (U, V) plane and the Toomre diagram for the super-solar star HAT-P-3 are consistent with the thin disk, also supported by the T D/D ratio close to zero. HAT-P-22 shows v tot ∼ 50 km s −1 and a T D/D ratio of ∼0.08. Again, both components of the XO-2 binary system are slightly above the above-mentioned chemical limit and also have high Galactic space velocities. Moreover, their position in the Toomre diagram is within the locus of the thick disk stars, and the T D/D ratio is ∼13, higher than the threshold established by Bensby et al. (2014) for potential thick disk stars. From the analysis of the Galactic orbits we find an eccentric orbit (e G ∼ 0.44) with a maximum height above the Galactic plane of ∼104 pc. The relatively low Z max led in the past to conclude that the binary system is confined to the Galactic thin disk (see Burke et al. 2007;Damasso et al. 2015a), but the other chemokinematic indicators suggest that the XO-2 binary system is at least in the thin disk-thick disk transition.

Chemical and isochronal ages
We computed stellar ages using elemental abundance ratios. It was indeed demonstrated that abundance ratios of pairs of elements produced over different timescales (e.g., [Y/Mg], [Y/Al]) can be used as valuable indicators of stellar age (Nissen 2015). Their [X/Fe] ratios show opposite behaviors with respect to stellar age (e.g., [Mg/Fe] and [Y/Fe] abundances respectively decrease and increase with stellar age). Therefore, their ratio, for instance [Y/Mg], shows a steep increasing trend with stellar age (see Casali et al. 2020, and references therein). The latter authors derived relations in the form [A/B] = c + x 1 · [Fe/H] + x 2 · Age, with [A/B] generic abundance ratio used as a chemical clock (see their paper for a wide description of the method). Here we considered the multivariate linear regression parameters c, x 1 , and x 2 for all abundance ratios in common   Table 6). We derived mean ages from all the abundance ratios and rejected those values discrepant by more than 1σ from the average. Moreover, we also computed stellar ages (and masses) from isochrone fittings (and evolutionary tracks). We therefore considered the PARSEC 9 models by Bressan et al. (2012) and the PARAM interface 10 (version 1.3; da Silva et al. 2006). This code considers as input some observational parameters (effective temperature, parallax, apparent V magnitude, and iron abundance) to perform a Bayesian determination of the most likely stellar intrinsic properties, appropriately weighting all the isochrone sections that are compatible with the observational parameters. A flat distribution of ages with a range of 0.1-15 Gyr was considered as priors for the analysis. We considered as effective temperature and iron abundance those values we derived as described in Sect. 3.1. The parallax was taken from Gaia EDR3, while the V magnitude was computed from the Gaia EDR3 G magnitudes, G BP and G RP colors using the appropriate photometric relationships (Gaia Collaboration 2021) and the reddening maps by Capitanio et al. (2017).
In Fig. 9 the comparison between the ages obtained with all the abundance ratios as a function of the isochronal ages obtained from PARSEC models is shown. In the same figure the range of ages listed in Bonomo et al. (2017) for each star and derived through stellar evolutionary tracks is shown for comparison. We tried multiple combinations of abundance ratios to find the mean chemical ages by using the best abundance ratios discussed in Casali et al. (2020), and the better agreement with isochronal ages from PARSEC was found considering the two ratios [Y/Al] and [Y/Mg]. In Table 5 we list the chemical ages obtained by using all abundances and these two abundance ratios, together with the isochronal ages from the PARAM tool. 9 PAdova and TRieste Stellar Evolutionary Code. 10 http://stev.oapd.inaf.it/cgi-bin/param_1.3 In particular, we find a mean difference between the ages derived from all the available abundance ratios and those obtained through PARSEC models of 2.7 Gyr with a standard deviation of 2.4 Gyr. Six targets (Qatar-2, HAT-P-12, HAT-P-26, Qatar-1, HAT-P-3, WASP-43) fall outside the 1σ limit, while three (KELT-6, WASP-10, WASP-39) are close to it. The comparison with Bonomo et al. (2017) gives an age difference of 2.9 Gyr with a standard deviation of 2.7 Gyr. If we consider the two best abundance ratios ([Y/Al] and [Y/Mg]; see, e.g., Casali et al. (2020), and references therein) the agreement is even better than that obtained using all abundance ratios, with a mean age difference of 2.0 Gyr and a standard deviation of 2.0 Gyr. In this case the position of Qatar-2, HAT-P-12, WASP-43, and KELT-6 is closer to that of the isochronal age, with a visibly improved general agreement. The comparison with Bonomo et al. (2017) in this case gives an age difference of 2.5 Gyr with a standard deviation of 2.3 Gyr. We note here that the most discrepant target, HAT-P-26, is also one of the stars with high values of space velocities (see Sect. 4.3), with v tot > 70 km s −1 , and with the [α/Fe] and T D/D ratios compatible with the thin to thick disk transition. This justifies its chemically old origin. Other stars (HAT-P-12, Qatar-1, Qatar-2) are in a region of the log g-T eff diagram for which the age determination from both the isochrones and the chemical indicators is problematic mainly due to their relatively cool effective temperature.

Lithium abundance
We find that seven targets show lithium in their spectra, with HAT-P-3 the most uncertain case (see Fig. 10). KELT-6 shows a value of ∼1.2 dex; WASP-38, HAT-P-14, and WASP-13 have lithium abundances around 2 dex; and HAT-P-4 and HAT-P-30 have values log n(Li) ∼ 2.8−3.0 dex. Lithium in WASP-38, WASP-13, HAT-P-4, and HAT-P-30 was also detected by Mortier et al. (2013) and Enoch et al. (2011), who measured abundances very close to our values. Our targets with lithium A161, page 15 of 30 A&A 664, A161 (2022) HAT-P-3 9.6 ± 1.5 8.6 ± 1.6 4.7 ± 4.4 0.88 ± 0.03 HAT-P-4 4.8 ± 0.9 3.8 ± 0.9 3.8 ± 0.7 1.28 ± 0.05 HAT-P-12 13.8 ± 2.2 11.8 ± 3.0 6.0 ± 5.2 0.69 ± 0.02 HAT-P-14 2.4 ± 1.5 1.0 ± 0.4 1.2 ± 0.4 1.41 ± 0.03 HAT-P-15 3.6 ± 1.4 3.9 ± 0.6 2.9 ± 2.9 0.98 ± 0.04 HAT-P-17 8.7 ± 0.4 8.3 ± 0.6 7.3 ± 5.1 0.87 ± 0.04 HAT-P-18 7.4 ± 2.0 6.8 ± 2.4 5.3 ± 4.8 0.78 ± 0.03 HAT-P-20 7.5 ± 1.2 3.6 ± 3.4 5.4 ± 4.9 0.74 ± 0.02 HAT-P-21 11.2 ± 0. 9 9.4 ± 0.5 10.5 ± 1.7 0.97 ± 0.04 HAT-P-22 9.5 ± 1.1 9.7 ± 1.1 10.8 ± 3.1 0.92 ± 0.03 HAT-P-26 13.5 ± 1.4 14.2 ± 0.5 7.0 ± 5.4 0.84 ± 0.03 HAT-P-29 2.1 ± 0.4 1.2 ± 1.0 1.4 ± 1.1 1.20 ± 0.03 HAT-P-30 2.2 ± 0.8 1.9 ± 1.0 1.7 ± 1.0 1.25 ± 0.04 HAT-P-36 3.0 ± 0.8 3.3 ± 0.5 5.8 ± 4.0 1.00 ± 0.05 KELT-6 8.8 ± 1.9 5.5 ± 3.1 4.8 ± 0.8 1.14 ± 0.05 Qatar-1 12.0 ± 1.7 11.4 ± 2.5 6. are outside the region known as the lithium desert, located at ∼5900−6200 K and log n(Li) < 2 dex (see Fig. 11), in agreement with Ramírez et al. (2012) and López-Valdivia et al. (2015). By chance, the number of stars with T eff > 5900 K with lithium content higher than ∼2.0 dex (i.e., 5) exceeds the number of stars with depleted lithium (i.e., 1), in accordance with Pavlenko et al. (2018), and justified by their thinner envelopes. On the other hand, all targets with clearly detected lithium have T D/D ≤ 0.05, with the most metal-poor KELT-6 having the lower value of log n(Li). This behavior could also reflect the pattern observed for Galactic thin stars, in addition to the depletion dependences related to their different stellar parameters (see Ramírez et al. 2012). We find, as expected, that the targets with lithium show on average higher v sin i (∼5.0 km s −1 ) than stars without lithium detected (∼2.6 km s −1 ), higher T eff (∼6090 K vs. ∼5260 K), and higher stellar mass (∼1.10 M vs. ∼0.97 M ), lower log g (4.2 dex vs. 4.4 dex), lower [Fe/H] (0.08 dex vs. 0.13 dex), lower chemical age (3.9 Gyr vs. 6.7 Gyr), and lower planetary mass (1.1 M Jup vs. 1.5 M Jup ). Moreover, within the sample with detected lithium, we recognize a tendency for the star with lower Li content  to be also that with a chemically derived older age. Similar dependence on lithium abundance on stellar parameters were also found by Delgado  for solar twins observed with HARPS and by Pavlenko et al. (2018) for CHEOPS dwarf stars. The dependence of Li abundance on T eff and on stellar mass is mainly due to the fact that high T eff (and therefore higher-mass) stars have thinner envelopes, which naturally leads to higher Li, while cooler stars have deeper convective envelopes which allow for Li to get into hot enough regions for processing to occur (see, e.g., the pioneering work by Boesgaard et al. 1998). The barely noticeable effect on log g could be due to stronger lithium depletion in the atmospheres of older stars (and therefore with higher surface gravity) or enhanced mixing in stars with deeper convective envelopes; alternatively, those stars with lower log g could be slightly evolved and thus their T eff at the main sequence could have been higher; therefore, they did not destroy so much Li due to their thinner convective envelopes (see Pavlenko et al. 2018). The slightly lower mean [Fe/H] of stars with Li compared to stars without detectable lithium could be caused by deeper convective envelopes expected for high opacities or Galactic chemical evolution (see Delgado . Even if our sample is biased toward relatively slow rotating stars, we observe a difference in v sin i between targets with and without lithium compatible with depletion induced by rotation, as suggested in the rotational evolution models by Bouvier (2008) and observed by Pavlenko et al. (2018) for CHEOPS dwarf stars with planets. We also find higher planetary masses for stars without lithium detected, which seems to indicate that destruction of Li is greater when the planet is more massive. Similar results were found by Delgado , who suggested that this could be justified within a scenario where the disk is affecting the evolution of angular momentum, with a stronger effect for a more massive disk, a condition needed to form a giant planet (Bouvier 2008). Moreover, when a giant planet forms in the disk, the accretion processes are expected to be more frequent and violent and produce Li destruction because of the temperature increase at the base of the convective envelope (Baraffe et al. 2015).
In Fig. 11 we plot the lithium abundance versus T eff and stellar age derived through chemical indicators. Despite the relatively wide T eff range (∼5200−6700 K), the plot shows an evident negative trend of Li abundance versus age, in which stars with higher lithium content show younger chemically derived ages. For the chemically old (∼8 Gyr) star HAT-P-3 we cannot make conclusions about the possible age derived through the Li because of its upper limit in log n(Li). KELT-6 is placed on the left of the warmest part of M67 members with similar T eff and its age derived through [Y/Al] and [Y/Mg] elemental abundance ratios is 5.5 Gyr, consistent with the open cluster M67 (Pasquini et al. 2008). The position of WASP-13 within the T eff -log n(Li) diagram is compatible with M67 stars, giving support to the chemically derived age of ∼4.6 Gyr. HAT-P-4 and HAT-P-30 are placed close to the Hyades group, while their chemical age is ∼2−4 Gyr, with HAT-P-4 older than HAT-P-30. Finally, HAT-P-14 and WASP-38 appear to be respectively in the warm and in the cool side of the Hyades Li dip, with a chemical age of ∼1−2 Gyr.
Finally, in our fit of the lithium line we tried to include the 6 Li/ 7 Li ratio as a free parameter. This was done because the inclusion of this isotopic ratio improves the fit of the λ6707.8 Å line and also because the determination of 6 Li/ 7 Li would improve our knowledge of the stars in our sample. Standard and non-standard stellar evolution models predict in solar-type stars a destruction of 6 Li at the base of the convective envelope (Talon & Charbonnel 2005, and references therein), and hence the presence of 6   . Spectral synthesis around the Li region for the seven stars with detected lithium (with the exception of HAT-P-3 for which only an upper limit in Li abundance was derived). The best-fit NLTE log n(Li) and its margins of error superimposed on the observed spectrum are shown with solid red and dotted light blue lines, respectively. Israelian et al. 2001;Mott et al. 2017). On the other hand, no 6 Li/ 7 Li was detected in other stars with planets (Reddy et al. 2002;Ghezzi et al. 2009;Harutyunyan et al. 2018). Unfortunately, our spectra are not of sufficiently high S NR to determine 6 Li/ 7 Li ratio with high enough precision. Usually, the 6 Li/ 7 Li isotopic ratio is determined from analysis of spectra with S NR > 600 (Mott et al. 2017). Trying to derive the 6 Li/ 7 Li ratio from our MOOG analysis (see Sect.3.2), we tentatively find a 6 Li/ 7 Li value of around 0.07 for the solar spectrum, which is between the values of ∼0.05 and ∼0.08 derived for the Sun by Asplund et al. (2021) and Lang (1974), respectively. For our subsample with lithium detected, we infer from our analysis a 6 Li/ 7 Li lower than the solar value, with typical values ranging from ∼0.01−0.02 for HAT-P-3, KELT-6, and WASP-38 up to ∼0.03 for HAT-P-4, HAT-P-14, HAT-P-30, and WASP-13. Similar low values of 6 Li/ 7 Li ratios are compatible with null results, in the sense that we do not find a significant amount of 6 Li in our sample of stars with detected Li. We tried to use other line lists (like that by Meléndez et al. 2012) and obtained similar findings; therefore, we think that higher S NR is needed for more reliable 6 Li/ 7 Li measurements, as suggested by Mott et al. (2017). [O/α] shows a possible decreasing trend (with a Spearman statistical significance ρ ∼ 0.2) with increasing iron abundance, as also found for FG-type stars observed within the HARPS GTO (gray dots; Bertran de . This kind of trend is expected if O abundance scales with the iron abundance, as suggested above. Finally, hints of enhancements in sulfur are observed for the components of the XO-2 binary system (at [Fe/H] ∼ 0.3−0.4 dex).
As mentioned in the introduction, volatile elements like C, N, O, and S, can be used as proxy of the star-planet formation history (Turrini et al. 2021a). These authors have found that using the ratios of C/N, N/O, and C/O for the planets and for their hosting stars breaks the degeneracy in the formation and migration tracks of giant planets, while the ratio S/N provides an additional independent probe into the metallicity of giant planets and their accretion of solids. We computed for our stellar sample the mean elemental ratios for two different stellar metallicity regimes around the iron abundance peak of our targets, which is [Fe/H] ∼ 0.12 dex (see Table 6). As the comparison of these mean values highlights, the elemental ratios of interest for planetary studies can vary by as much as ∼10-20% between stars with different metallicity, for example with C/N going from values of around ∼2.8 for solar-metallicity or metal-poor targets up to ∼3.1 for super metal-rich targets. As such, the use of reference solar values would introduce significant biases in the interpretation of the planetary compositional data (see Turrini et al. 2021a,b). As an illustrative example, a planet with a C/N ratio of ∼3.4 or C/O ratio of ∼0.6 orbiting a star with metallicity belonging to the ≤1.3 Z group would be interpreted as possessing a solar value of this ratio when compared to the Sun, while its C/N value could actually be ∼1.   impose solar abundance values to study the bulk exoplanetary properties.

[Mg/Si] vs. [Fe/H] and Mg/Si vs. C/O
In Fig. 13 we show how the [Mg/Si] ratio depends on stellar iron abundance. We also overplot the results by Adibekyan et al. (2015) for planet-hosting stars, keeping in mind that the authors declared that they were not able to find differences between stars with and without high-mass planets. Our targets show values consistent with those by Adibekyan et al. (2015), with a few outliers that are placed at similar positions of undetected planethosting stars by the same authors (see their Fig. 1). Moreover, we note that the lowest-mass planet-hosting star in our sample, HAT-P-26, shows the highest [Mg/Si] ratio, in agreement with Adibekyan et al. (2015), who declare that low-mass planets are more prevalent around stars with high [Mg/Si]. This could be due to the fact that [Mg/Si] probably plays a very important role in the formation of low-mass planets; for example, high Mg abundances could mitigate the lower iron abundances or metallicities and make core accretion comparatively more efficient. Moreover, there could be also a dependence of the planetary structure on the Galactic chemical evolution, as for this target we recognized from chemodynamical diagnostics a possible thin to thick disk transition origin within the Galactic plane (see Sect. 4.3). A clearer picture can become evident if we also consider volatile elements. While refractory elements (as traced by the abundance ratios Mg/Si 11 and Fe/Si) condense close to the star and their abundance ratios remain constant throughout most of the disk, the same is not true for volatile elements like C/O (see Thiabaud et al. 2015b). Such elemental ratios are important because they govern the distribution and formation of chemical species in the protoplanetary disk, and hence the mineralogy of planets. In Fig. 14 Brewer et al. (2016) for FGK stars in the solar neighborhood. We also note a tendency for stars cooler than T eff < 5000 K to have mean C/O ratios smaller than those of warmer stars, with a difference of ∼0.08. No clear difference is evident for the Mg/Si ratio.
In regions of high C/O, planets form primarily from carbonates, and in regions of low C/O, the Mg/Si determines the types of silicates that dominate the compositions (e.g., Brewer et al. 2016, and references therein). This means that the C/O ratio controls the distribution of Si among carbide and oxide species. If C/O is greater than 0.8, Si exists in solid form primarily as SiC, and also graphite and TiC will be formed; for C/O ratios below 0.8, Si is present in rock-forming minerals as SiO 4− 4 or SiO 2 , serving as seeds for Mg silicates for which the exact composition will be controlled by the Mg/Si value (Bond et al. 2010;Thiabaud et al. 2015b). Moreover, Thiabaud et al. (2015a) have shown that the condensation of volatile species as a function of radial distance allows for C/O enrichment in specific parts of the protoplanetary disk of up to four times the solar values, leading to the formation of planets that can be enriched in C/O in their envelope up to three times the solar value. This is the case of HD209458b observed by Giacobbe et al. (2021) for which a scenario of planet formation beyond the water snowline and migration toward its host star through disk or disk-free migration was hypothesized (see also Brewer et al. 2017). At the same time, Mg/Si governs the distribution of silicates: for Mg/Si < 1, Mg forms orthopyroxene (MgSiO 3 ) and the excess Si is present as other silicate species such as feldspars (CaAl 2 Si 2 O 8 , NaAlSiO 8 ) or olivine (Mg 2 SiO 4 ); for Mg/Si values ranging from 1 to 2, Mg is distributed between olivine and pyroxene; for Mg/Si > 2, all available Si is consumed to form olivine with excess Mg available to bond with other minerals, mostly oxides such as MgO or MgS (Bond et al. 2010;Thiabaud et al. 2015b). The peak of the Mg/Si-C/O distribution for our targets is therefore consistent with Si which will take solid form as SiO 4− 4 and SiO 2 and Mg equally distributed between pyroxene and olivine.
Considering the whole sample of 28 targets for which Mg/Si was determined, the star hosting the planet with mass <30 M ⊕ (HAT-P-26) has the highest Mg/Si value (∼1.5) 12 , while 75% of the higher-mass companion host sample show Mg/Si values between 1.0 and 2.0, which means that Mg is equally distributed between pyroxene and olivine. For the 25% of high-mass planet hosts with Mg/Si values below 1.0, Mg and Si will form mainly orthopyroxene, whereas the remaining Si will take other forms, such as feldspars or olivine. No stars with Mg/Si > 2.0 were found. We note here that targets like HAT-P-26 and HAT-P-12, with their high values of Mg/Si, are also among those resulting chemically older in our sample. We indeed find a weak trend between Mg/Si and isochronal ages, which could mean that the determination of this chemical ratio could be favored for older (and chromospherically less active) stars (see Sect. 4.4 for the measurement of chemical and isochronal ages). Conclusions of this kind should be confirmed by statistically more significant samples.
If we consider the sample of 18 targets for which both Mg/Si and C/O ratios were measured, 100% of our planet-hosting stars have C/O values lower than 0.8, and 11% of the sample stars have C/O < 0.4. This means that Si will be present in rock-forming minerals as SiO 4− 4 and SiO 2 . In these cases, silicate mineralogy will be controlled by the Mg/Si ratio. Within these 18 targets, 15 high-mass planet hosts show Mg/Si values between 1.0 and 2.0, two targets have Mg/Si < 1.0, and the low-mass planet-hosting star HAT-P-26 shows a value of ∼1.5. This supports the finding by Suárez-Andrés et al. (2018), who claimed that low-mass planets are likely to be found in the 1.0-1.5 Mg/Si regime, although mixed with stars with high-mass planets.
We also note here that most targets in our sample show subsolar values of Mg/Si and C/O ratios. In particular, if we consider the solar values of Mg/Si = 1.17 ± 0.08, C/O = 0.57 ± 0.04, we find that 18/28 (i.e., 64%) and 11/18 (i.e., 61%) of our targets show respectively Mg/Si and C/O ratios lower than the solar values at the 1σ level. As already mentioned in Sect. 4.6, this highlights once again how the use of solar values as reference stellar abundances could introduce biases in the interpretation of the planetary compositional models (see Turrini et al. 2021a,b).

Stellar abundances versus planetary properties
We show in Fig. 15 [X/Fe] versus planetary mass 13 for some volatile elements (i.e., C, N, O, S), with the aim of looking 13 In the cases of multiple planets, such as KELT-6, XO-2S, and HAT-P-17, we considered the planet b, which is the closest planet and that for possible relations between stellar abundances and planetary properties. The masses of the planets range between 0.058 and 7.27 M Jup . In the figure, to better visualize possible trends, we represent with asterisks bins at increasing mass steps (i.e.,0.4,0.6,1.2,2.4 M Jup ) in order to have similar number of targets per bin. Due to the few targets for which we could measure N abundances, we are not able to draw any conclusion for this element. For the other elements, we see a probable decreasing trend of [C/Fe], [O/Fe], and [S/Fe] with increasing M p . For these trends we calculated the Spearman associated statistical significance (ρ), finding ρ ∼ 9 × 10 −6 for the [O/Fe]-M p relation and ρ ∼ 0.07 for sulfur and carbon. We note that for the planet XO2-Sb we plot M p sin i p (instead of M p ), but we find similar results also when excluding this target, in particular for oxygen. Again, since at least part of these trends could be due to the star hosting the lowest-mass planet, we computed the same statistical significance after excluding HAT-P-26 and we find ρ around 0.18, 6 × 10 −5 , and 0.07 for C, O, S, respectively. This means that a more significant correlation is present for the [O/Fe] ratio versus M p with respect to S and C versus M p , with an evident decreasing step toward low [O/Fe] values for M p > 0.5 M Jup . A flat tendency was found for C by Suárez-Andrés et al. (2017) and a probable increasing trend was found by Suárez-Andrés et al. (2016) for N. We note that most of the targets hosting lowmass planets and showing higher [O/Fe] are also those resulting chemically old and possibly belonging to the thin-thick disk transition. Moreover, we note that the number of stars in those papers and in this work is not statistically significant, but if our findings are confirmed through bigger samples it could imply that the formation of low-mass planets is favored at the highest values of stellar volatile elements. A possible, yet speculative, interpretation of these trends could be offered by the giant planet formation process in the framework of the pebble accretion scenario. Giant planets of Jovian or super-Jovian mass are expected to form early in the lifetime of circumstellar disks, when the high disk mass accretion rate is capable of supporting the rapid growth of their cores and sustaining their gas accretion rates (see Johansen et al. 2019;Tanaka et al. 2020 for a discussion). As a result, such massive planets can also form in disks comparatively poor in the abundant volatile elements O and C. Less massive planets are expected to form over longer timescales, in circumstellar disks characterized by lower mass accretion rates (Hartmann et al. 1998;Johansen et al. 2019). Such planets may therefore form more easily around volatile-rich stars, so that the higher abundance of ices can partly compensate the lower disk mass accretion rates supporting their growth. In this framework, the stronger trend observed between the planetary mass and the [O/Fe] ratio could be explained by the larger contribution of O to the mass fraction of heavy elements in stars and their disks (in the Sun O provides about 45% of the mass of heavy elements, while C and S only 16 and 2% respectively; Lodders 2010), and its lower volatility compared to C (e.g., Turrini et al. 2021a, and references therein). This means that any increase in the abundance of O would affect wider orbital regions and would have a larger impact on the availability of solid material of planetforming disks than equal increases in the abundance of C and S, thus resulting more effective in promoting the formation of low-mass planets. On the contrary, the possible trends discussed above cannot be explained as easily in terms of effects linked to the planetary multiplicity, for example by the more massive giant planets blocking the flux of pebbles and promoting the causing the transit (in the cases of KELT-6 and HAT-P-17). We verified that this choice, motivated by homogeneity reasons, does not change our results and conclusions.  formation of low-mass planets on outer orbits. All giant planets in the sample, with the possible exclusion of HAT-P-26, are significantly more massive than the pebble isolation masses typical of the inner regions of circumstellar disks (Johansen et al. 2019), so they should all be capable of blocking the pebble flux at some point during their migration. Furthermore, due to the lower volatility of S with respect to O (e.g., Turrini et al. 2021a,b and references therein), this process should result in lower values of S/Fe than O/Fe at higher planetary masses, a signature which is not observed in the data. Specifically, the sublimation of the O frozen as ice in the pebbles (e.g., water ice sublimating once the inward-drifting pebbles cross the water snowline) results in a smaller fraction of O being trapped in the pebbles with respect to S and Fe (which remain in solid form until closer to the star), making the trapping process less effective for O than for S and Fe.
We also analyzed Mg/Si and C/O as a function of planetary mass. In Fig. 16 we see a hint of a trend between Mg/Si and M p (with slope of ∼−0.18 and a Spearman associated statistical significance ρ ∼ 0.3, which became ∼−0.09 and 0.6 after excluding HAT-P-26) and a still less evident C/O-M p relation (slope of ∼0.16 and ρ ∼ 0.5, which became ∼0.31 and ∼0.2 if we do not consider the lowest-mass planet HAT-P-26b). Similar results are obtained after excluding XO-2S for which we know M p sin i p (and not M p ). A possible downward trend of Mg/Si-M p and a slight increasing trend of C/O-M p relationships were reported by Mishenina et al. (2021), but no conclusion was drawn because of the large scatter of the ratios at certain planetary mass ranges. Very recently, Tautvaišienė et al. (2022) find a weak negative C/O slope and a slightly more negative Mg/Si slope toward stars with high-mass planets. Here, we note that for icy and giant planets, Mg/Si (and also Fe/Si) in stars should be a direct information about the planetary composition, as no differences are expected between star and planet in terms of Mg/Si, as proposed by Thiabaud et al. (2015a). This is because giant planets are mainly formed outside the ice line, a region where all refractory material has condensed. Our findings could imply higher values of Mg/Si for lower-mass planets, like HAT-P-26b and HAT-P-12b, which indeed were found to have respectively super-solar and either solar or super-solar metallicity (see Kawashima & Min 2021). We note here that the solar-metallicity star HAT-P-26 and the metal-poor HAT-P-12 are also targets for which we find old chemical ages and the possibility of belonging to the Galactic thin-thick disk transition (see Sects. 4.4,4.3); therefore, their relatively high Mg/Si ratios could be related to their position in the Galactic disk. Regarding C/O, an indirect relation between the planet and the star should be present. This is mainly because final planetary C/O depends on the location and timescale of formation, how much of the atmosphere is accreted from gas versus solids, and how isolated the atmosphere is from mixing with core materials (see, e.g., Teske et al. 2014 andGiacobbe et al. 2021, andreferences therein). In Sect. 4.9 we try to give some possible pathways for exoplanets in common with our hosting stars for which C/O and N/O elemental ratios were computed in the literature.
In Fig. 17 we show the distribution of the iron abundance of our stellar sample in terms of the planetary orbit eccentricity, mass, radius, density, and stellar mass. We find a tendency for high-eccentricity planets to be around more metal-rich stars, as also found by Dawson & Murray-Clay (2013) and Mills et al. (2019). Even if we are aware of our not statistically significant sample, we note that the mean [Fe/H], and also the mean values of [α/H] and Mg/Si, are greater by ∼0.03−0.05 dex for planets with e > 0.1 when compared with lower-eccentricity planets. This tendency appears not to be influenced by the planetary radii because the planet radii of our sample are within a range (6−21 R E ) for which no stellar metallicity-planet radius correlation was found (see Buchhave et al. 2014). In order to assign a confidence level of our result, we used a one-side 2 × 2 Fisher's exact test 14 (Agresti 1992). Choosing the divisions at e = 0.1 for high-and low-eccentricity orbits and [Fe/H] = 0.00 dex for metal-poor and metal-rich stars, we find a p-value of 0.29 as the chance that random data would yield this trend, indicating a probability of correlation of 71%. With this in mind, we also find some evidence that our metal-rich stars have wider ranges of planetary masses and denser planets (mean ρ P = 1.68 g cm −3 ) are around stars with greater [Fe/H] (and therefore in more eccentric orbits) when compared with planets around more metal-poor stars with e < 0.1 (mean ρ P = 1.11 g cm −3 ). We note here that we are cautious about these possible trends both because the number of our targets is relatively small and also because each correlation depends on the interplay of many planetary parameters and stellar properties. The disk gas therefore becomes enriched in N with respect to O and C the farther we move from the star, while solids become increasingly enriched in C and O with respect to N. As a result, the farther from the host star the giant planets start their migration, the more Table 7. Stellar elemental ratios derived in this work (Cols. (2)-(4)), planetary elemental ratios computed by Kawashima & Min (2021) (Cols. (5)- (7)), and planetary-to-stellar elemental ratios defined by Turrini et al. (2021a) (Cols. (8)-(10)). 12 2.09 ± 0.16 0.31 ± 0.18 0.15 ± 0.16 5.33 ± 0.18 0.80 ± 0.15 0.15 ± 0.10 2.55 ± 0.24 2.58 ± 0.23 1.00 ± 0.19 WASP-10 2.00 ± 0.21 0.37 ± 0.18 0.19 ± 0.18 6.31 ± 0.32 0.82 ± 0.30 0.13 ± 0.10 3.16 ± 0.38 2.22 ± 0.35 0.68 ± 0.21 HAT-P-26 3.80 ± 0.15 0.55 ± 0.12 0.14 ± 0.12 1.79 ± 0.16 0.25 ± 0.13 0.14 ± 0.10 0.47 ± 0.22 0.45 ± 0.18 1.00 ± 0.16 WASP-39 3.73 ± 0.12 0.51 ± 0.12 0.14 ± 0.14 1.73 ± 0.16 0.26 ± 0.12 0.15 ± 0.10 0.46 ± 0.20 0.51 ± 0.17 1.07 ± 0.17 Notes. i. Planetary C/N were derived by dividing C/O by N/O; ii. When we were not able to measure one of the C, N, O stellar abundances, we considered as [X/H] a value consistent with the abundance of elements with the closer T cond (within the volatile elements C, N, and O; see Table 2 and Sect. 4.2). their C/N and N/O ratios will diverge from the stellar values due to the different accretion efficiencies of gas and solids.
The use of planetary elemental ratios normalized to the respective stellar ratios makes it easier to extract the information on the nature of giant planets and constrain whether their metallicity is dominated by the accretion of gas or of solids (see Turrini et al. 2021a,b for additional discussion). Four of the stars analyzed in the present work (HAT-P-26, WASP-10, HAT-P-12, and WASP-39) host planets for which Kawashima & Min (2021) derived metallicity, N/O, and C/O through the spectral disequilibrium retrieval models. Dividing these planetary ratios by those we derived for the hosting stars (which we label X/Y * 15 ), we can therefore gain insights into the formation pathways of these planets.
We list in Table 7 the elemental ratios as defined by Turrini et al. (2021a). Values of C/N * > C/O * > N/O * imply that the budget of heavy elements is dominated by the accretion of solids, while N/O * > C/O * > C/N * implies its accretion is mostly from the disk gas. The authors also claim that in both cases, the separation between the values of the three normalized ratios will increase with the extent of disk-driven migration experienced by the giant planet. Specifically, the farther from the star a giant planet starts its formation, the more significant the difference between its C/N * , C/O * , and N/O * . For our sample of four targets we note that large uncertainties associated with the measurements of the X/Y * ratio cannot allow us to draw a definitive conclusion, and therefore our discussion is mainly qualitative and indicative of possible planetary formation scenarios.
The pattern of the abundance ratios observed for HAT-P-26 and WASP-39 points toward C/N * lower than N/O * , with possible constant and increasing trends between C/N * and C/O * , respectively. Given their C/N * and N/O * abundance patterns, these planets likely underwent migration to get to their present orbits (both at ∼0.048 au) and accreted mostly gas along their path (Turrini et al. 2021a). The comparison of the C/N * , C/O * , and N/O * values is consistent with a scenario where both planets started forming outside the CO 2 snowline and accreted most of their gas inward of it (Turrini et al. 2021a), but where the gas was enriched in O by the evaporation of O-rich ices from the inward drifting dust (Booth & Ilee 2019). The two host stars have solar iron abundances, while their planets show super-solar and solar metallicity, respectively for HAT-P-26b and WASP-39b (see MacDonald & Madhusudhan 2019; Kawashima & Min 2021). Since the accretion of non-enriched gas would result in substellar metallicity values (Turrini et al. 2021a,b), the metallicity of these planets appears consistent with the accretion of gas enriched in heavy elements. The estimated abundance patterns tentatively favor the accretion of gas in a disk that underwent chemical reset (Pacetti et al., in prep.).
For the planets orbiting HAT-P-12 and WASP-10 we find C/N * greater than N/O * , with possible constant and decreasing trends between C/N * and C/O * , respectively. The C/N * , C/O * , and N/O * abundance patterns of the planet orbiting WASP-10 are consistent with extensive migration toward its present orbit (at ∼0.038 au) and the accretion of heavy elements being dominated by the accretion of solids (Turrini et al. 2021a). The high value of the C/O * ratio is suggestive of an additional contribution to the planetary C budget by the accretion of C-enriched gas, which would favor the accretion of gas between the CO 2 and CH 4 snowlines (Booth & Ilee 2019). The C/N * , C/O * , and N/O * abundance patterns of the planet orbiting HAT-P-12 are also consistent with extensive migration and the accretion of heavy elements being dominated by the accretion of solids. In this case, however, the higher values of the C/O * and N/O * ratios suggest that the planetary budgets of both C and N were affected by the accretion of gas enriched by the evaporation of ices from the inward drifting dust. This would point toward the planet accreting a significant fraction of its gas between the N 2 and CO 2 snowlines (Booth & Ilee 2019). In such a scenario, the planet orbiting HAT-P-12 would have started its formation farther out than its counterpart orbiting WASP-10.

Conclusions
In this paper we presented for the first time a wide and comprehensive characterization of a sample of 27 transiting planet host stars within the GAPS programme obtained through a homogeneous and accurate spectroscopic procedure, which is mainly based on high-resolution HARPS-N at TNG spectra (and a few FEROS at ESO spectra). We analyzed this sample for the first time with the aim of deriving different stellar properties, abundances of many elements, and kinematic information. Our main results can be summarized as follows: -We obtained the atmospheric parameters and the abundances of 26 elements from lithium to europium, together with their kinematic properties. The chemokinematic analysis allowed us to recognize that most of the targets appear to belong to the Galactic thin disk, and a few stars likely belong to a thin-thick disk transition. The lithium line is present in seven stars. -From the analysis of some elemental ratios, most of them including s-process elements and α-elements, we derived, for the first time in exoplanet hosts, stellar ages often consistent with those obtained through theoretical isochrones. -From the analysis of the Mg/Si ratios, we find that most of the targets show values consistent with a distribution of Mg between olivine and pyroxene, and a few show Mg forming orthopyroxene. The C/O ratio for all targets is lower than 0.8, compatible with Si present in rock-forming minerals, with a slight tendency of higher values for stars hosting lower-mass planets, even if we note that at least part of this trend could be related to the position of these targets in the Galactic disk.  iii. The solar-metallicity HAT-P-21 shows T D/D and Z max compatible with thin-thick disk transition and seems to originate from the outer Galactic disk. -Finally, since detailed knowledge of the formation of a planet requires accurate knowledge of chemical abundances of its host star, we tried to discuss the formation and migration mechanisms of those targets for which abundances of planets hosted by stars analyzed in the present work were obtained for the same elemental ratios. We suggest that the planets orbiting HAT-P-26 and WASP-39 started forming outside the CO 2 snowline, while those around HAT-P-12 and WASP-10 probably formed between the CO 2 and CH 4 showlines and between the N 2 and CO 2 showlines, respectively. We think that analyses like those performed in this work will be necessary for future studies on planetary composition that take into account host star composition, in particular for transiting planet host stars, for which more information about the system formation, migration, and evolution can be retrieved. Metallicity and Mg/Si, C/O, C/N, S/N, N/O ratios are important indicators of planet formation; therefore, future high-precision observations are essential to further explore the trend between stellar and planetary properties toward understanding the formation mechanisms of planets. For instance, some of the planets will be observed by JWST, and therefore this kind of studies are motivated by the prospects of the chemical characterization of exoplanets, and how the chemical compositions of planet host stars relate to those of their planets. Forthcoming JWST observations and other upcoming infrared spectroscopic missions (like ARIEL) will allow us to draw more robust conclusions, in particular regarding the level of precision for planetary abundances, useful to provide definitive conclusions, for example on planetary formation and evolution. Similarly, precise spectroscopic studies together with high-quality photometric data like those acquired with TESS or in the near future with PLATO will be useful to clarify evolutionary stages and to better characterize planetary systems in a global way.    (12) HAT-P-26 0.12±0.11 (2,s) ... 0.14±0.07(3) 0.04±0.04(2) 0.18±0.09(3) 0.25±0.08(2) 0.08±0.08(12) 0.08±0.06(3) 0.07±0.12(9) 0.18±0.11(3) 0.17±0.11(21) 0.20±0.07(4) 0.20±0.07 (12) HAT-P-29 0.08±0.03 (2,s) ... 0.17±0.15(s) 0.13±0.03(4) 0.13±0.11(4) 0.11±0.10(3) 0.21±0.09(11) 0.25±0. 16(3) 0.24±0.08(13) 0.18±0.10(3) 0.17±0.06(20) 0.20±0.08 (5) 0.16±0.09 (9) HAT-P-30 −0.04±0.03(2,s) 0.02±0.08(2) −0.02±0.06(3) −0.01±0.06(4) 0.01±0.09(4) −0.10±0.08(2) 0.04±0.07(11) 0.15±0.14(3) 0.09±0.11(14) 0.03±0.08(3) 0.04±0.06(19) 0.05±0.11(6) 0.04±0.08 (7) HAT-P-36 0.  (2) 0.26±0.08(11) 0.29±0.13(2) 0.19±0.08(15) 0.26±0.08 (7) 0.19±0.07 (6) WASP-10  (2) 0.12±0.09 (9) ... 0.10±0.14(6) 0.15±0.14(3) 0.14±0.13 (16) (6) (2) 0.36±0.13(1) 0.22±0.14(3) ...