Context. The study of exoplanets at different evolutionary stages can shed light on their formation migration and evolution. The determination of exoplanet properties depends on the properties of their host stars. It is therefore important to characterise the host stars for accurate knowledge on their planets. Aims. Our final goal is to derive in a homogeneous and accurate way the stellar atmospheric parameters and elemental abundances of ten young TESS Transiting Exoplanet Survey Satellite transiting planethosting GK stars followed up with the HARPSN High Accuracy Radial velocity Planet Searcher for the Northern hemisphere at TNG Telescopio Nazionale Galileo spectrograph within the Global Architecture of Planetary Systems GAPS programme. Methods. We derived stellar kinematic properties atmospheric parameters and abundances of 18 elements. Depending on stellar parameters and chemical elements we used methods based on line equivalent widths and spectral synthesis. Lithium line measurements were used as approximate age estimations. We exploited chemical abundances and their ratios to derive information on planetary composition. Results. Elemental abundances and kinematic properties are consistent with the nearby Galactic thin disk. All targets show COlt0.8 and 1.0ltMgSilt1.5 compatible with an interpretation of silicate mantles made of a mixture of pyroxene and olivine assemblages for any rocky planets around them. The FeMg ratios with values of 0.71.0 show a propensity for the planets to have big iron cores. All stars hosting very lowmass planets show MgSi values consistent with the Earth values thus demonstrating their similar mantle composition. Hot Jupiter host stars show a lower content of OSi which could be related to the lower presence of water content. We confirm a trend found in the literature between stellar OFe and total planetary mass implying an important role of the O in shaping the mass fraction of heavy elements in stars and their disks. Conclusions. Exploring the composition of planets through the use of elemental abundances of their hosting stars holds promise for future investigations particularly with large samples. Meanwhile the detailed host star abundances provided can be employed for further studies on the composition of the planets within the current sample when their atmospheres will be exploited.
['techniques: spectroscopic', 'stars: abundances', 'stars: atmospheres', 'stars: fundamental parameters', 'planetary systems', 'Astrophysics - Solar and Stellar Astrophysics', 'Astrophysics - Earth and Planetary Astrophysics']
The GAPS Programme at TNG LXI. Atmospheric parameters and elemental abundances of TESS young exoplanet host stars
{'LXI. Atmospheric parameters and elemental abundances of TESS young exoplanet host stars /star': "- S. Filomeno 1 , 2 , 3 , K. Biazzo 1 , M. Baratella 4 , S. Benatti 5 , V. D'Orazi 2 , 6 , S. Desidera 6 , L. Mancini 2 , 7 , 8 , S. Messina 9 , D. Polychroni 10 , D. Turrini 7 , 11 , L. Cabona 6 , I. Carleo 12 , 13 , 7 , M. Damasso 7 , L. Malavolta 6 , 14 , G. Mantovan 14 , 6 , D. Nardiello 14 , 6 , G. Scandariato 9 , A. Sozzetti 7 , T. Zingales 6 , 14 , G. Andreuzzi 15 , 1 , S. Antoniucci 1 , A. Bignamini 10 , A. S. Bonomo 7 , R. Claudi 6 , 16 , R. Cosentino 15 , A. F. M. Fiorenzano 15 , S. Fonte 11 , A. Harutyunyan 15 10  \nand C. Knapic  \n- 1 INAF-Osservatorio Astronomico di Roma, Via Frascati 33, I-00040 Monte Porzio Catone (RM), Italy e-mail: simone.filomeno@inaf.it\n- 2 Dipartimento di Fisica, Università di Roma Tor Vergata, Via della Ricerca Scientifica 1, I-00133 Roma, Italy\n- 3 Dipartimento di Fisica, Sapienza Università di Roma, Piazzale Aldo Moro 5, I-00185 Roma, Italy\n- 4 ESO - European Southern Observatory, Alonso de Córdova 3107, Casilla 19, Santiago 19001, Chile\n- 5 INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, I-90134 Palermo, Italy\n- 6 INAF - Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, I-35122 Padova, Italy\n- 7 INAF - Osservatorio Astrofisico di Torino, Via Osservatorio 20, I-10025 Pino Torinese, Italy,\n- 8 Max Planck Institute for Astronomy, Königstuhl 17, D-69117, Heidelberg, Germany\n- 9 INAF - Osservatorio Astrofisico di Catania, Via S. Sofia 78, I-95123 Catania, Italy\n- 10 INAF - Osservatorio Astronomico di Trieste, Via Tiepolo 11, I-34143 Trieste, Italy\n- 11 INAF - Istituto di Astrofisica e Planetologia Spaziali, Via del Fosso del Cavaliere 100, I-00133 Roma, Italy\n- 12 Instituto de Astrofísica de Canarias (IAC), Calle Vía Láctea s / n, 38200 La Laguna, Tenerife, Spain\n- 13 Departamento de Astrofísica, Universidad de La Laguna (ULL), 38206 La Laguna, Tenerife, Spain\n- 14 Dipartimento di Fisica e Astronomia, Universit'a di Padova, Vicolo dell'Osservatorio 2, I-35122 Padova, Italy\n- 15 Fundación Galileo Galilei - INAF, Rambla José Ana Fernandez Pérez 7, E-38712 Bre˜na Baja (TF), Spain\n- 16 Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy  \nReceived ... / accepted ...", 'ABSTRACT': 'Context. The study of exoplanets at di ff erent evolutionary stages can shed light on their formation, migration, and evolution. The determination of exoplanet properties depends on the properties of their host stars. It is therefore important to characterise the host stars for accurate knowledge on their planets.  \nAims. Our final goal is to derive, in a homogeneous and accurate way, the stellar atmospheric parameters and elemental abundances of ten young TESS (Transiting Exoplanet Survey Satellite) transiting planet-hosting GK stars followed up with the HARPS-N (High Accuracy Radial velocity Planet Searcher for the Northern hemisphere) at TNG (Telescopio Nazionale Galileo ) spectrograph within the Global Architecture of Planetary Systems (GAPS) programme.  \nMethods. We derived stellar kinematic properties, atmospheric parameters, and abundances of 18 elements. Depending on stellar parameters and chemical elements, we used methods based on line equivalent widths and spectral synthesis. Lithium line measurements were used as approximate age estimations. We exploited chemical abundances and their ratios to derive information on planetary composition.  \nResults. Elemental abundances and kinematic properties are consistent with the nearby Galactic thin disk. All targets show C / O < 0.8 and 1.0 < Mg / Si < 1.5, compatible with an interpretation of silicate mantles made of a mixture of pyroxene and olivine assemblages for any rocky planets around them. The Fe / Mg ratios, with values of ∼ 0.7-1.0, show a propensity for the planets to have big (iron) cores. All stars hosting very low-mass planets show Mg / Si values consistent with the Earth values, thus demonstrating their similar mantle composition. Hot Jupiter host stars show a lower content of O / Si, which could be related to the lower presence of water content. We confirm a trend found in the literature between stellar [O / Fe] and total planetary mass, implying an important role of the O in shaping the mass fraction of heavy elements in stars and their disks.  \nConclusions. Exploring the composition of planets through the use of elemental abundances of their hosting stars holds promise for future investigations, particularly with large samples. Meanwhile, the detailed host star abundances provided can be employed for further studies on the composition of the planets within the current sample, when their atmospheres will be exploited.  \nKey words. Stars: abundances, fundamental parameters, atmospheres - Techniques: spectroscopic - Planetary systems  \n- ,', '1. Introduction': 'Precise stellar parameters and elemental abundances of planethosting stars are crucial for correctly interpreting exoplanet observations. Ground-based high-resolution spectroscopy provides the main avenue for determining the stellar parameters of exoplanet hosts. In recent years, much e ff ort by many researchers in the community has pushed this work forward. Most of them are focussed on old FGK (Sousa et al. 2021, Magrini et al. 2022, Biazzo et al. 2022, MacDougall et al. 2023, Strassmeier et al. 2023, Yun et al. 2024) and M (Maldonado et al. 2020, Marfil et al. 2021) dwarf planet-hosting stars. In this context, the characterisation of transiting exoplanet host stars is of particular interest because they represent unique laboratories to test theories of exoplanet formation and evolution with a relatively high precision and possibly trace any potential migration they underwent (Biazzo et al. 2022, Kolecki & Wang 2022).  \nPlanet formation models and theories can be constrained through the properties of the stellar atmospheres, as the composition of the host star can be used as a proxy for the composition of the initial protoplanetary disk from which the planets were born. Therefore, by comparing planetary and stellar atmospheric chemical composition, we can gain insight into where and how planets formed and possibly from where they migrated. This has, in turn, led modelling studies to focus on predicting the enrichment of volatiles, such as carbon and oxygen, in giant planet atmospheres (Öberg et al. 2011,Madhusudhan et al. 2014,Fonte et al. 2023). Although this framework provides useful constraints on the potential formation location of gas giant exoplanets, carbon and oxygen measurements alone are not enough to determine from where a given gas giant planet originates. The degeneracy in the information provided by, for example, the C / O ratio can be broken by characterising the abundances of both volatile elements (such as sulphur, in addition to carbon and oxygen) and refractory elements (such as magnesium and silicon), which condense at higher temperatures. This can provide stronger constraints on the formation region and extent of planet migration. Therefore, the ratio of refractory content relative to volatile content in the atmosphere can give us direct insights into the nature of the solids that were incorporated into the planetary envelope (Turrini et al. 2021, 2022; Pacetti et al. 2022). This is particularly interesting for young stars because, compared with older counterparts that possibly experienced migration throughout the Galactic disk, their chemical content should have preserved much of the original chemical composition of the cloud from which they formed (Spitoni et al. 2023, Magrini et al. 2023). Therefore, the chemistry of young exoplanet-hosting stars, and thus of their planets, is strongly connected to the time and place of its birth, unless the gas from which they formed has undergone peculiar enrichment.  \nIn this paper, we conduct a homogeneous and accurate stellar parameter and elemental abundance analysis of several atomic species, including volatiles and refractories (such as carbon, oxygen, sulphur, magnesium, silicon, titanium, and iron), which are processed in planetary formation and migration. The analysis was performed as in previous works (Biazzo et al. 2022) in the  \ncontext of the Italian Global Architecture of Planetary Systems (GAPS) programme (Covino et al. 2013) and within the subprogramme devoted to relatively young transiting planet host stars (Baratella et al. 2020a, Carleo et al. 2020). Inside the whole sample of targets in this sub-programme, we studied a subset of ten relatively young dwarf stars for which no similar extensive analysis has been previously done involving all elements we measured. Our final aim is to provide accurate stellar parameters and abundances of multiple elements using high-quality spectroscopic data, to possibly link stellar elemental abundances to the composition of their planets, and to investigate possible relationships between stellar abundances and the properties of their transiting companions. This will represent a step forward in providing necessary information for future studies on the atmospheric abundance of relatively young planets with facilities such as JWST (James Webb Space Telescope) (Gardner et al. 2023) and the next-coming Ariel (Atmospheric Remote-Sensing Infrared Exoplanet Large-survey) telescope (Tinetti et al. 2021).  \nThe paper is organised as follows. In Sects. 2 and 3 we describe the spectroscopic dataset and the measurements of stellar parameters and elemental abundances. In Sect. 4 we present and discuss our results in the context of chemical evolution, kinematics, and planetary abundances and properties. Finally, in Sect. 5 we draw our conclusions.', '2. Stellar sample, observations, and data reduction': 'The GAPS-Young Objects (YOs) sub-programme has the goal to study young and intermediate-aged planetary systems (with typical ages of ≤ 1Gyr) and to characterise their exoplanets down to sub-Neptune masses with transit and radial velocity analysis techniques, after accurate determinations of the host star properties. Within the GAPS-YOs sample, we selected TESS (Transiting Exoplanet Survey Satellite) Objects of Interest (TOIs; Ricker et al. 2014) hosting at least one known transiting planet and showing indications of young ages from our systematic check of their stellar properties (see Carleo et al. 2020). Then, within this sample, we selected targets with e ff ective temperatures ( T e ff ) higher than 5000 K and rotational velocities ( v sin i ) lower than 10 km / s known from the NASA Exoplanet Archive. This was done to avoid problems caused by the presence of molecular lines or by line blending, which need a di ff erent approach than that adopted in this work (Biazzo et al. 2022, and references therein). We also considered two targets not belonging to the GAPS-YOs sub-programme but with similar characteristics and observed within the GAPS collaboration (see Table 1). The final sample consists of ten targets with V magnitude between ∼ 6.9 and ∼ 12.0 mag. All these targets have one to six confirmed planets with masses ranging from ∼ 2 . 4 M E to ∼ 10 M J (rocky to giant planets) as shown in more detail in the discovery papers mentioned in Tables 1.  \nObservations were performed between the beginning of 2020 and the end of 2023 with the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N; resolving power R ∼ 115000, spectral range λ ∼ 3830-6900Å; Cosentino et al. 2012) placed at the 3.6 m Telescopio Nazionale Galileo (TNG) in La Palma, Spain. Spectra reduction was performed by using the standard HARPS-N Data Reduction Software (Pepe et al. 2002) and all the spectra collected for each target were subsequently combined to obtain a co-added spectrum with signal-to-noise ratio S / N > 100 at ∼ 6000Å.  \nSince we are interested in the oxygen abundance and one of its diagnostics is the near-IR triplet at ∼ 7770Å, that is outside  \nthe HARPS-N spectral range, we searched for spectra acquired with the Tillinghast Reflector Echelle Spectrograph (TRES; R ∼ 44000, λ ∼ 3850-9100Å; Szentgyorgyi & Furész 2007) spectrograph at the Tillinghast telescope (Mt. Hopkins, Arizona) or with the CHIRON ( R ∼ 80000, λ ∼ 4150-8800Å; Tokovinin et al. 2013) spectrograph at the Cerro Tololo Inter-American Observatory (Chile) and available in the Exoplanet Follow-up Observing Programme (ExoFOP) website. In the end, we found TRES spectra for eight stars (TOI-1430, TOI-1726, TOI-2076, TOI-4515, TOI-5082, TOI-5398, TOI-5543, TOI-1136) and CHIRON spectra for one target (TOI-179) with typical S / N of ∼ 3040 at around 6000 Å.', '3. Analysis and results': 'The analysis of the HARPS-N spectra can be divided into two main steps: the measurement of the stellar atmospheric parameters (e ff ective temperature T e ff , surface gravity log g , and microturbulence velocity ξ ), iron and titanium abundances ([Fe / H], [Ti / H]), and the measurement of other photospheric elemental abundances fixing the previously derived parameters. Since we performed a di ff erential analysis relative to the Sun, we took as reference the solar parameters and elemental abundances derived by Biazzo et al. (2022, see that paper for details).', '3.1. Stellarparametersandabundanceofironandtitanium': 'The first step of the analysis was performed following the same approach and same line list as in Baratella et al. (2020a), which exploits both Fe and Ti lines. This procedure is preferred to the standard approach (based only on the Fe lines), especially for the study of young and active stars, like those analysed in this work. The main reason for this preference is that young stars have a high chromospheric variability which a ff ects spectral lines. Since Ti lines form deeper in the photosphere than Fe lines, they are less a ff ected by the chromosphere (see also Baratella et al. 2020b for details).  \nFor the measurement of atmospheric parameters, Fe and Ti abundances we employed pyMOOGi (Adamow 2017) which is a python wrapper for the MOOG code (Sneden 1973, version 2019). We used the equivalent widths (EWs) method by running the driver abfind , which assumes a local thermal equilibrium (LTE). We then considered the 1D model atmospheres linearly interpolated from the ATLAS9 grid of Castelli & Kurucz (2003), with solar-scaled chemical composition and new opacities (ODFNEW). The EWs of the target stars were measured for each line using the ARESv2 code (Sousa et al. 2015), but substantial double checks (with the splot task in IRAF 1 ) were performed, particularly in the blue part of the spectrum because of the intrinsic di ffi culties in optimal continuum tracing. In this case, depending on the line profile, Gaussian fitting procedure or direct integration was performed. During the analysis, when the EWs were greater than 150mÅ (strong lines), the related lines were discarded. The same was done for those lines with fitting errors larger than 2 σ . Once the line EWs were measured, T e ff was derived by imposing the excitation equilibrium, therefore reducing to zero the slope between the abundances of both Fe i and Ti i lines and the excitation potential. Then, Ti i and  \nTi ii lines were used to derive log g by imposing the ionisation equilibrium (di ff erence between the mean Ti i and Ti ii abundances lower than 0.01 dex), while ξ was obtained by removing the trend between the Ti i abundances and the reduced equivalent width REW = log( EW /λ ) (Baratella et al. 2020a). We used an iterative procedure by changing parameters at steps of 5 K in T e ff , 0.01 km / s in ξ , and 0.01 dex in log g . The uncertainties on T e ff and ξ were calculated by varying each quantity until the slopes of the relative trends were larger than three times their value with the optimised parameters, while the uncertainty of log g was calculated by varying this parameter until the di ff erence between neutral and ionised Ti became larger than the highest error between the abundances of the two species. We finally derived the stellar parameters with internal accuracy (at 3 σ ) ranging from 30 to 100 K in T e ff , from 0.05 to 0.17 dex in log g , from 0.07 to 0.15 km / s in ξ , from 0.08 to 0.12 dex in [Fe / H], and from 0.07 to 0.14 dex in [Ti / H]. In Table 1 we list the final stellar parameters ( T e ff , log g , ξ ), and the Fe and Ti abundances, which we derived within the GAPS-YOs sub-programme. We note that for seven out of ten targets, the analysis of the stellar parameters was performed within GAPS-YOs for planet validation studies already published or to be submitted. Instead, for the remaining targets, the stellar parameters were derived in a homogeneous way in this work and for the first time within the GAPS project 2 . Moreover, for the whole sample, no abundance analysis of multiple elements was previously performed.', '3.2. Otherelementalabundancesthroughequivalentwidths': "Once the stellar parameters and iron together with titanium abundances were computed, we measured the abundances ([X / H] 3 ) of the following other 16 elements: lithium (Li), carbon (C), oxygen (O), sodium (Na), magnesium (Mg), aluminium (Al), silicon (Si), sulphur (S), calcium (Ca), chromium (Cr), nickel (Ni), zinc (Zn), yttrium (Y), zirconium (Zr), cerium (Ce) 4 , and neodymium (Nd), with Cr abundances measured in the first two ionisation states, like Fe and Ti (see Table 2). For all the 16 elements mentioned above, we adopted the line list in Biazzo et al. (2022) and applied the EW method. Taking into account the same grids of model atmospheres, tools, and procedures as those described in Sect. 3.1, we measured the EWs of each line and discarded the strong lines with EWs greater than 150 mÅ and those with fitting errors larger than 2 σ .  \nFor those targets observed with the TRES and CHIRON spectrographs, we measured the oxygen abundance through the O i triplet lines (7771.94Å, 7774.17 Å, 7775.39 Å) and applied the non-LTE (NLTE) corrections by Amarsi et al. (2015).  \nIn addition to O i triplet lines, we also applied NLTE correction to each line abundance of C i and Na i and taking into account the prescriptions by Amarsi et al. (2019) and Lind et al. (2011), respectively. For the other elements, the corrections were  \n/epsilon1  \nwhere log  \n(  \nX  \n) is the absolute abundance.  \nnegligible or not reported in the literature (see Biazzo et al. 2022 for details).  \nWe also measured the abundance of lithium by applying the EW procedure to the 6707.8Å line. We find that all our targets show the lithium 6707.81Å line in their spectra, with the exception of TOI-5543. Table 1 lists the values of our measurements of lithium equivalent widths and abundances. In the case of TOI-1430 and TOI-4515, we could only derive for log /epsilon1 (Li) upper limits of 0.1 dex and 0.5 dex, respectively. Then, we considered the departure from LTE considering the NLTE calculation of Lind et al. (2009). Final NLTE Li abundances, together with their uncertainties due to both EW measurements and stellar parameters added quadratically, are reported in Table 1. A doublecheck of the final results was done also computing the lithium abundances using the spectral synthesis method (see Sect. 3.3).  \nFinally, we looked for possible trends of final elemental abundances (also for Fe and Ti derived in the previous subsection) with the stellar parameters ( T e ff , log g , ξ , v sin i ), and with the activity index log R ' HK 5 listed in Table 1. The log R ' HK values were measured in this work for each target adopting the colour index B-V extracted from the tables of Pecaut & Mamajek (2013) (version 2022) by using our T e ff values. No relationship was found for our sample.", '3.3. Rotationalvelocitiesandelementalabundancesthrough spectralsynthesis': 'The projected rotational velocity ( v sin i ) was measured through the spectral synthesis technique and considering the synth driver of the same codes used in Sect. 3.1. We then created synthetic spectra from the Castelli & Kurucz (2003) grids of model atmospheres fixing the stellar parameters T e ff , log g , ξ , and [Fe / H]to those derived in Sect. 3.1. We convolved the resulting synthetic spectra with a Gaussian profile corresponding to the resolving power of HARPS-N of R ∼ 115 000, taking into account the optical limb-darkening coe ffi cients by Claret (2019) at the T e ff , log g , ξ , and [Fe / H] of our targets. As done in Biazzo et al. (2022), we applied the empirical relationships for the macroturbulence velocity ( V macro) by Doyle et al. (2014) for stellar T e ff> 5700 K, while empirical relationships by Brewer et al. (2016) were considered for T e ff< 5700 K (see that paper for details). Finally, we synthesised spectral lines in three regions around 5400, 6200, and 6700 Å, as done in, for example, Carleo et al. (2024), until the minimum of residuals between stellar and synthetic spectra were reached. The final values are listed in Table 1, where uncertainties in stellar parameters together with spectral continuum definition were taken into account for the v sin i error estimates (for the uncertainties in the continuum position, see Sect. 3.4).  \nWe also applied the same spectral synthesis method to derive the oxygen abundance using the forbidden [O i ] line at 6300.3Å. Atomic parameters for the oxygen line and the nearby (blended) Ni i line, which exhibits an isotopic structure, were taken from Ca ff au et al. (2008) and Johansson et al. (2003), respectively. Since nickel abundances were needed for the synthesis of the oxygen line profile, their values were fixed to the ones we derived above (see Sect. 3.2). Since the 6300.3Åline is not a ff ected by deviation from LTE, we did not apply any NLTE corrections (see Ca ff au et al. 2008 for details). In the end, when TRES and CHIRON spectra were available, O abundances were derived with both EWs and spectral synthesis methods, therefore  \nas final abundances, we considered the weighted average coming from these two methods.  \nFinally, the spectral synthesis technique was also adopted for the lithium abundance (log /epsilon1 (Li)) measurements. Therefore, we used the line list by Reddy et al. (2002) in the vicinity of the Li 6707.8 Å line, implemented with the VALD database (Kupka et al. 2000) for wavelengths farther from the line centre, as done in Biazzo et al. (2022). The lithium abundance derived through the same codes and model atmospheres used in Sect. 3.1 was then corrected for the departure from LTE considering the NLTE calculations of Lind et al. (2009). We found consistent values with both methods based on line EW and synthesis within the uncertainties. In Table 1 we list only the results obtained through the EW method, with the exception of TOI-1430 and TOI-4515 for which only the synthesis method allowed us to derive upper limits in log /epsilon1 (Li).', '3.4. Uncertaintiesinelementalabundances': "The errors in the abundances [X / H]are mainly due to uncertainties in atomic parameters, stellar parameters, and EW measurements or continuum position (for the spectral synthesis method).  \nUncertainties in atomic parameters should be negligible since our analysis is performed di ff erentially with respect to the Sun.  \nErrors due to stellar parameter uncertainties were estimated by varying each parameter separately by its assessed error while keeping the others unchanged. As shown in Table 1, the uncertainties in e ff ective temperature are in the range 30-100 K (with standard deviation σ = 20K), 0.05-0.17dex for log g (with σ = 0.04 dex), 0.07-0.15km / s for ξ (with σ = 0.03 km / s), 0.080.12 dex for [Fe / H] (with σ = 0.01 dex), and 0.07-0.14dex for [Ti / H] ( σ = 0.02 dex). Due to the small values of the standard deviations and the relatively small range in all stellar parameters and Fe and Ti abundances (among the others with very similar σ ), we decided to consider typical uncertainties in stellar parameters of 70 K, 0.10 dex, 0.10 km / s, 0.10 dex in T e ff , log g , ξ , and [Fe / H], respectively. Errors in elemental abundances [X / H] due to uncertainties in individual stellar parameters are listed in Table 3 for, respectively, the coolest and the warmest targets in our sample.  \nThe errors in log /epsilon1 (X) due to uncertainties in EWs were obtained by computing the standard deviation around the mean abundance determined from all the measured lines of each element. From these values, the uncertainties in [X / H] were obtained by the quadratic sum of the error for the target and the error for the Sun, as shown in Tables 1 and 2 together with the number of lines employed for the abundance analysis (in brackets). When only one line was measured (such as for the lithium element), the error we used was the standard deviation of three independent EW measurements obtained considering three different positions of the continuum.  \nThe uncertainties on the oxygen abundances derived through the spectral synthesis method are of two types: i . errors related to the best-fit determination, which were evaluated by changing the continuum position until the standard deviation was two times larger than the best-fit value; typical uncertainties, in this case, are around 0.12 dex; ii . errors related to the uncertainty in the nickel abundance measurements, as the oxygen line at 6300.3 Å is blended with a Ni line (see Sect. 3.3; typical uncertainties due to Ni are around 0.05 dex). These sources of internal errors were quadratically added to obtain the total error (see Table 2).  \nTable 1. Final stellar parameters, Fe and Ti abundances (the numbers in parenthesis indicate the number of lines used), EW Li and NLTE lithium abundance, the mean activity index log R ' HK , and the stellar rotation period P rot of our sample.  \nReferences: a Desidera et al. (2023), b This Work, c Nardiello et al., in prep., d Damasso et al. (2023), e Damasso et al. (2024, in press), f Carleo et al. (2024), g Mantovan et al. (2024a), h Malavolta et al. in prep., i Dai et al. (2023), j Newton et al. (2023), k Nardiello et al. (2022).  \nTable 2. Final elemental abundances of our stellar sample.  \nNote: The number of lines used in the analysis and the cases for which we applied the spectral synthesis method (s) are shown in parentheses.", '4.1. ElementalabundancesinthecontextoftheGalacticdisk': 'Each Galactic component (thin and thick disk, bulge, and halo) is characterised by di ff erent chemical abundance patterns when compared to the iron abundance. The observed di ff erences in these abundance patterns are the consequence of a variety of star formation histories. Therefore, possible relationships between elemental abundances and Fe abundances are generally used to trace the chemical evolution of our Galaxy (Bensby et al. 2014). In particular, the study of low-mass stars can be very useful in explaining the history of the Galactic chemical evolution. The chemistry of a star is indeed intimately connected to the time and place of its birth (see Adibekyan et al. 2011 and references therein). In order to check that this is indeed true also for our sample of relatively young transiting exoplanet host stars, we  \ncompare [X / Fe] 6 ratio versus [Fe / H] of our targets with abundances of nearby (distance < 200 pc) field stars of the Galactic thin disk selected from the Hypatia Catalogue (Hinkel et al. 2014). In particular, we considered targets with stellar parameters similar to those of our sample (i.e. 5000 < T e ff (K) < 6000, 4 . 4 < log g (dex) < 4 . 6). Figure 1 shows overall good agreement between the [X / Fe] versus [Fe / H] distribution of our targets and the nearby field stars with iron abundance close to zero. In the same figure, we overplot the abundances homogeneously derived within GAPS (except for Ce, taken from Delgado Mena et al. 2017) for older FGK transiting planet-hosting stars. The agreement is clear, with our targets compatible with solar values.  \nThe thin- and thick-disk stars can also be chemically distinguished by looking at their α -element content at a given iron abundance. Adibekyan et al. (2011) showed that for dwarf stars  \n∆  \nT  \nTable 3. Internal errors on abundance measurements due to uncertainties in stellar parameters.  \nTOI-5398  \n[Fe  \n[Fe  \n[C  \n[O  \ni  \ni  \n/  \nii  \n/  \ni  \nH]  \n/  \nH]  \n0.03  \nH]  \n0.05  \n/  \n[Na  \nH]  \n0.07  \ni  \n[Mg  \n[Al  \n[Si  \n[S  \ni  \ni  \ni  \n/  \n[Ca  \n[Ti  \n[Ti  \n[Cr  \n[Cr  \n[Ni  \n[Zn  \n[Y  \ni  \n/  \ni  \n/  \n/  \nH]  \n/  \nH]  \nH]  \nH]  \nH]  \n0.04  \ni  \n/  \n/  \nii  \ni  \nH]  \nH]  \n/  \n/  \nii  \ni  \nH]  \n0.01  \nH]  \n/  \n/  \ni  \nii  \n[Zr  \nH]  \n0.03  \nH]  \n/  \n/  \nii  \n[Ce  \n[Nd  \ni  \nH]  \nH]  \n0.00  \n/  \n/  \nH]  \n0.00  \nH]  \nii  \n/  \nH]  \n-  \n-  \ne  \nff =  \nff =  \nT  \ne  \n-  \n6000 K  \nlog  \n-  \n/  \n+  \n/  \n0.05  \n-  \n/  \n-  \n/  \n-  \n70 K  \n0.05  \n0.00  \n0.01  \n0.03  \n0.04  \n/  \n/  \n0.03  \n0.04  \n0.01  \n/  \n0.02  \n0.03  \n0.01  \n/  \n0.02  \n0.03  \n0.00  \n/  \n0.01  \n-  \n0.02  \n0.00  \n0.02  \n/  \n/  \n0.04  \n0.05  \n0.02  \n/  \n0.06  \n/  \n0.07  \n0.00  \n0.01  \n/  \n0.05  \n-  \n0.05  \n0.00  \n0.01  \n/  \n/  \n0.03  \n0.04  \n0.00  \n/  \n0.01  \n/  \n0.02  \n0.01  \n/  \n0.01  \n/  \n0.01  \n0.02  \n/  \n0.01  \n0.02  \n-  \n-  \n-  \n-  \n-  \ng  \n=  \n-  \n4  \n/  \n+  \n-  \n44 dex  \n.  \n0  \n10 dex  \n0.01  \n0.02  \n/  \n/  \n0.03  \n0.04  \n0.03  \n/  \n0.03  \n0.03  \n0.00  \n/  \n0.04  \n-  \n/  \n-  \n/  \n0.04  \n0.01  \n0.02  \n0.01  \n0.01  \n0.01  \n0.00  \n0.00  \n/  \n0.00  \n0.01  \n/  \n0.02  \n-  \n/  \n0.03  \n0.00  \n0.02  \n0.02  \n0.00  \n0.02  \n/  \n0.02  \n/  \n0.05  \n0.04  \n0.00  \n0.02  \n/  \n0.03  \n/  \n0.04  \n0.02  \n0.00  \n0.01  \n/  \n0.01  \n0.00  \n0.05  \n/  \n0.04  \n0.04  \n0.04  \n/  \n0.04  \n0.04  \n0.02  \n/  \n0.04  \n0.04  \n0.01  \n/  \n0.04  \n0.04  \n0.00  \n/  \n/  \n1  \n/  \n+  \n12 km  \n.  \n0  \n-  \n/  \n-  \n/  \n/  \ns  \n[Fe  \n10 km  \n0.02  \n...  \n0.02  \n...  \n0.00  \n0.00  \n-  \n-  \n/  \n-  \n/  \n0.01  \n0.01  \n0.00  \n0.01  \n0.00  \n0.00  \n0.00  \n/  \n0.00  \n0.00  \n/  \n0.00  \n0.00  \n-  \n-  \n/  \n-  \n/  \n-  \n/  \n-  \n/  \n-  \n/  \n-  \n/  \n-  \n/  \n-  \n/  \n-  \n/  \n-  \n/  \n0.02  \n0.00  \n0.01  \n0.00  \n0.01  \n0.02  \n0.00  \n0.02  \n0.01  \n0.00  \n0.05  \n0.04  \n0.02  \n0.01  \n0.01  \n-  \n-  \n-  \n-  \n-  \n=  \n0.09 dex  \n-  \n/  \n+  \n/  \n...  \n/  \n/  \n...  \n0.00  \n/  \n0.04  \n/  \n0.03  \n0.00  \n/  \n0.00  \n/  \n0.00  \n/  \n0.01  \n/  \n0.00  \n/  \n0.00  \n/  \n0.00  \n/  \n0.02  \n/  \n0.04  \n0.00  \n/  \n0.02  \n/  \n0.02  \n0.00  \n/  \n0.02  \n0.02  \n/  \n0.03  \n0.03  \n/  \n0.03  \n0.03  \n/  \n0.03  \n0.03  \n/  \n0.03  \n0.03  \nNotes. Measurements were done for TOI-1430 (for all elements but Fe i , Fe ii , Ti i and Ti ii which refer to TOI-5082) and TOI-5398 (for all elements but Fe i , Fe ii , Ti i and Ti ii which refer to TOI-1136), respectively the coolest and the warmest stars in our sample. The values in table refer to the di ff erences between the abundances obtained with ( -and + ) and without the uncertainties in stellar parameters.  \nthe two populations are clearly separated in terms of [ α / Fe] 7 from low metallicities up to super-solar metallicities. The [ α / Fe] versus [Fe / H] position obtained for our targets is depicted in Fig.2, together with the division by Adibekyan et al. (2012) separating thin- and thick- disk populations (yellow dashed line). The corresponding values taken from the Hypatia Catalogue for nearby dwarf stars with similar parameters as our targets are overplotted with small black dots (Hinkel et al. 2014). It is evident that our targets (filled blue dots) are all below the dashed line and around solar [ α / Fe] (and [Fe / H]) values, consistently with the position of Galactic thin stars.', '4.2. KinematicpropertiesinthecontextoftheGalacticdisk': 'After the chemical analysis performed in the previous sections, we now characterise the population of our targets in terms of their kinematics. We expect indeed that not only the chemistry, but also the dynamical history of stars, including stellar kinematics, may impact the distribution and the architecture of planets (see Magrini et al. 2022; Biazzo et al. 2022, and references therein). Using the mean radial velocities ( < V rad > ) from the HARPS-N spectra, parallaxes ( π ) and proper motions ( µα, µδ ) from Gaia DR3 (Gaia Collaboration et al. 2016, 2021), we computed space velocity components UVW with respect to the local standard of rest (LSR), correcting for the solar motion derived by Co¸skunoˇglu et al. (2011): ( U /circledot , V /circledot , W /circledot ) = (8.50, 13.38, 6.49) km / s. We then considered the general outline of Johnson & Soderblom (1987) in a right-hand coordinate system  \n.  \n0  \n10 dex  \n-  \n-  \n/  \n-  \n/  \n[Fe  \n/  \nH]  \nH]  \n=  \n/  \n/  \n/  \n-  \n-  \n-  \n-  \n-  \n-  \nlog  \ng  \n=  \n-  \n-  \n-  \n-  \n-  \n-  \n-  \n-  \n-  \n.  \nξ  \n=  \n-  \n.  \n∆  \n∆  \nξ  \n=  \n/  \ns  \n∆  \nFig. 1. [X / Fe] versus [Fe / H] for our sample (filled blue dots). Horizontal and vertical dashed lines mark the solar abundances. Filled grey dots represent the [X / Fe] values of old transiting planet host stars analysed by Biazzo et al. (2022), except for Ce for which we plot the values by Delgado Mena et al. (2017). Small black dots in the background represent the [X / Fe] distribution of nearby field stars (distance ≤ 200 pc) in the Galactic thin disk with similar parameters as those of our sample (Hinkel et al. 2014).  \n<!-- image -->  \n(i.e. with U pointing towards the Galactic centre, V towards the local direction of rotation in the plane of the Galaxy, and W towards the North Galactic Pole). The uncertainties in the space velocity components were obtained considering the prescriptions by Gagné et al. (2014) and taking into account the error contributions of V rad, π , µα , and µδ . Combining the errors on parallaxes, proper motions, and radial velocities, the resulting average uncertainties in the U , V , W velocities are around 0.15 km / s (for more detail see Table 4). In the left panel of Fig. 3, we show the location of our targets in the Boettlinger diagram in the ( U , V ) plane, where the boundary containing nearby ( < 150pc) associations younger than ∼ 1Gyr is shown with a dashed line (see Eggen 1996; Gagné et al. 2018). The locus containing young stars as defined by Francis & Anderson (2009) is also overplotted (dotted line). All our targets are close to or within the young and nearby Galactic thin disk boundary and all of them are within the locus defined by Francis & Anderson (2009). In the  \nright panel of Fig. 3 we plot the position of our targets in the Toomre diagram, which represents the quadratic addition of the radial and vertical kinetic energies ( v tot) versus the rotational component. In first approximation, stars with a total velocity lower than 50 km / s belong to the Galactic thin disk, while those with 70 < ∼ v tot < ∼ 180km / s are likely to be thick disk stars (Nissen 2004; Bensby et al. 2014). As expected, all our targets are part of the thin disk population, with v tot within or nearby the limit of ∼ 50km / s, while some of them, with v tot around 20 km / s, are very close to the Sun.  \nFinally, we also computed the thick-to-thin disk probability ratio, considering the prescriptions reported by Bensby et al. (2014) for the Gaussian distribution of random velocities of different stellar populations. In particular, to get the probabilities D and TD for the thin and thick disk membership, we considered the asymmetric drift, the velocity dispersion, and the fraction of each population given by Bensby et al. (2014). Our targets show  \nTable 4. Basic information and stellar kinematic properties of our targets.  \nNotes: columns list right ascensions, declinations, Gaia DR3 magnitudes, and parallaxes taken from Gaia Collaboration et al. (2016, 2021), spectral types obtained from the Pecaut & Mamajek (2013) table (version 2022) and our T e ff , UVW velocities, and thick-to-thin disk probability (TD / D) as derived in the present work. The last column lists the membership group from the literature.  \nReferences: a Desidera et al. (2023), b Dai et al. (2023), c Nardiello et al. in prep., d Mann et al. (2020), e Newton et al. (2023), f Hedges et al. (2021), g Carleo et al. (2024), h Mantovan et al. (2024a).  \nFig. 2. [ α / Fe] versus [Fe / H]. The overplotted grey points are targets from Biazzo et al. (2022). The dashed yellow line is the fiducial division into thin- and thick-field disk populations, as defined by Adibekyan et al. (2012). Small black dots in the background represent the [ α / Fe] distribution of nearby field stars representative of the nearby Galactic disk taken from the Hypatia Catalogue (Hinkel et al. 2014).  \n<!-- image -->  \nTD / Dratio of ∼ 0 . 01 -0 . 02, which is compatible with the membership to the Galactic thin disk (Table 4).', '4.3. Lithiumasachemicalclockandcomparisonwith gyrochronologicalandactivitydatingindicators': "Figure 4 shows the lithium abundance (left panel) as a function of the e ff ective temperature, where our targets are plotted as blue dots. In the same figure, we also plot the position of members of the Pleiades ( ∼ 100Myr, Sestito & Randich (2005); yellow crosses), the Hyades ( ∼ 600Myr, Cummings et al. (2017); red diamonds), and the M67 ( ∼ 4.5 Gyr, Pasquini et al. (2008); green triangles) clusters, together with the older transiting planethosting stars (grey squares) with detection of the lithium line and investigated by Biazzo et al. (2022). From this figure, two targets (namely, TOI-2048 and TOI-2076) are placed in the lower envelope of the Li distribution of the Pleiades, TOI-179 appears to be in between the lower envelope of the Pleiades and the Hyades cluster, other four targets (TOI-1136, TOI-1726, TOI5082 and TOI-5398) are very close to the Hyades position, TOI-  \n1430 seems to follow the Hyades sequence for low abundance values, while TOI-4515 is located between the Hyades and the older M67.  \nWe estimated the age of our targets from the lithium content considering a set of empirical isochrones in the EW Li versus T e ff plane calibrated using all the data for open clusters with age of 1-6000Myr observed as part of the Gaia -ESO Survey (GES; Randich et al. 2022, Gilmore et al. 2022) and implemented in the eagles (Empirical AGes from lithium Equivalent widthS; Je ff ries et al. 2023) code. This code models a set of stars by fitting the empirical isochrones and an intrinsic dispersion around those isochrones, for data in the range of 3000-6500K and -0.3 dex < [Fe / H] < 0.2 dex. A combined likelihood for the stars is multiplied by a prior probability to form a posterior probability distribution of age. Since it does not use lithium abundance, it is independent of the temperature scale or the adoption of any set of stellar atmosphere models or NLTE corrections.  \nTable 5 reports in the second column the most probable age and an asymmetric 68% confidence interval obtained through the eagles tool, while the right panel of Fig. 4 shows the lithium equivalent width as a function of the e ff ective temperature. In this figure, our data (blue dots) are plotted over a set of model isochrones in the 50 Myr-2 Gyr range, as derived by Je ff ries et al. (2023) through the fitting of the GES training data.  \nAll the targets in our sample have the log R ' HK index measured as part of this study as well as their rotation period known from the literature or estimated in this work from TESS data with the same approach used in Messina et al. (2022) (see Table 1). As an independent check of the age determination, we also estimated the ages of our targets using the calibrated agerotation and age-activity relations by Mamajek & Hillenbrand (2008). The derived gyrochronologic and chromospheric ages are listed in Cols. 3 and 4 in Table 5. Other than TOI-5543 (Li line not deteced), TOI-1430 and TOI-4515 (for which we could measure lower limit in age through the Li line), only two stars, namely TOI-5082 and TOI-179, show ages derived from Li slightly outside 1 σ with respect to both gyrochronological and activity ages. For the remaining five stars, the three age estimation methods give similar results. We performed a two-tailed Student's test for dependent samples and found that the null hypothesis, that is, that the Li ages come from the same population of gyrochronological and activity ages, cannot be rejected (p > 0.7).  \n<!-- image -->  \nFig. 3. Position of our planet-hosting stars in the Boettlinger and Toomre diagrams, where errors are within the symbol dimension. Open diamonds in both panels represent the position of young and nearby clusters analysed by Gagné et al. (2018). Left panel: The dashed line represents the boundary separating young disk and old disk stars in the ( U , V ) plane according to Eggen (1996), while the dotted line shows the velocity ellipsoid defined by Francis & Anderson (2009) for young stars. Right panel: The dashed lines indicate constant peculiar total velocities v tot = ( U 2 + V 2 + W 2 ) 1 / 2 = 20, 50, and 70 km / s.  \n<!-- image -->  \nWe remark here that Je ff ries et al. (2023) claim that for EW Li lower than 50 mÅ the eagles code can yield lower limit to ages. This is because the probability distribution becomes quite asymmetric as data points approach regions where EW Li becomes less sensitive to age because lithium starts to be depleted. This leads to a likelihood function poorly constrained for older ( > ∼ 1 Gyr) ages. Within our sample, we have four targets with EW Li < 50mÅ. This is the case of TOI-1430 and TOI4515 for which we could not constrain the age because of the very low EW Li value (0.6 mÅ) and the EW Li undefined mean value ( < 2 mÅ), respectively. However, the lower limits in ages derived for these two targets using the lithium line are consistent with the recent estimates by Nardiello et al. in prep. (see also ages by Orell-Miquel et al. 2023 and Carleo et al. 2024 of 700 ± 150 Myr and 1.2 ± 0.2 Gyr, respectively). Again, the position of these two targets in the log /epsilon1 (Li)T e ff diagram is consistent respectively with ages close to and older than the Hyades. Two other targets with EW Li lower than 50 mÅ are TOI-5082, for which no age was estimated in the literature, and TOI179, for which the peak of the posterior probability distribution is at ∼ 900Myr that is older than the recent determination performed by Desidera et al. (2023) and based on gyrochronology (400 ± 100Myr). In the latter case, we also remark that the position of the target in the log /epsilon1 (Li)T e ff diagram is indeed more consistent with the age values older than the Pleiades and closer to the Hyades clusters, in agreement with the gyrochronological and activity age estimates (Table 5). For the other targets with EW Li > 50mÅ, we find good agreement within 1 σ with the literature values. TOI-1726 is indeed member of the Ursa Major group ( ∼ 400Myr; Mann et al. 2020), TOI-2048 is member of the Group X ( ∼ 300Myr; Messina et al. 2022), TOI-2076 is a comoving system of TOI-1807 (300 ± 80Myr; Nardiello et al. 2022, Hedges et al. 2021), and, for TOI-1136, Dai et al. (2023) found an age of 700 ± 100 Myr from gyrochonology and activity indicators. In the case of TOI-5398, we  \nfind an age from lithium abundance which is younger than that one found by Mantovan et al. (2024a) through girochronology (650 ± 150Myr), but closer to their value of ∼ 370Myr derived through activity indicators. However, the higher value provided by the eagles code (see Table 5) and the position in the log /epsilon1 (Li)T e ff diagram near the Hyades cluster (Fig. 4, left panel) seems to be more consistent with the gyrochronological age estimated in Mantovan et al. (2024a).", '4.4. Mineralogicdiagnostics': 'Abundance ratios in planet-hosting stars are useful tools to obtain important information about the composition of their planets since they govern the distribution and formation of chemical species in the protoplanetary disks where they formed. Here in the following, we discuss how some geochemical ratios derived for planet-hosting stars may become useful indicators of the planetary mantle and core compositions, in particular for low-mass planets.', '4.4.1. The importance of the [Mg/Si] ratio': 'In Fig. 5 we show the [Mg / Si] ratio as a function of [Fe / H]. As demonstrated by Bond et al. (2010), Thiabaud et al. (2015) and Adibekyan et al. (2015), [Mg / Si] may play an important role in the internal structure and composition of planets. In particular, Adibekyan et al. (2015) found that this ratio can be used to constrain the chemical composition of low-mass planets and also that planets around stars belonging to the Galactic thick disk have higher values of this ratio when compared to their thin disk counterparts. This is because Mg, Si, and Fe elements have di ff erences in the production sites within the Galactic disk and this is reflected in the observed trends. Indeed, Adibekyan et al. (2015) found that low-mass planets are more prevalent around stars with high [Mg / Si], thus implying  \n<!-- image -->  \nFig. 4. Lithium abundance and equivalent width vs T e ff plots. Left: non-LTE Li abundance, log /epsilon1 (Li), versus e ff ective temperature, T e ff . The targets of this work (blue filled dots) are plotted together with the targets from Biazzo et al. (2022) in grey open squares and the members of the Pleiades ( ∼ 100 Myr; Sestito & Randich 2005), Hyades ( ∼ 600 Myr; Cummings et al. 2017), and M67 ( ∼ 4.5 Gyr; Pasquini et al. 2008) clusters (in yellow crosses, red diamonds, and green triangles, respectively). Right: lithium equivalent width ( EW Li) versus e ff ective temperature. Here, the empirical model isochrones by Je ff ries et al. (2023) at 50, 100, 150, 200, 300, 400, 500, 600, 700, 800, 900 Myr, 1, 1.5, and 2 Gyr are overplotted.  \n<!-- image -->  \nTable 5. Age estimates obtained from di ff erent methods.  \nNote: Col. 2 lists the most probable age from lithium equivalent width with the lower and upper bounds of the asymmetric 68 per cent confidence interval, as derived from the eagles code (Je ff ries et al. 2023). Columns 3 and 4 show the age estimates obtained from gyrochronology and chromospheric activity indicators, while Col. 5 shows age values from the literature.  \nReferences: a Desidera et al. (2023), b Dai et al. (2023), c Orell-Miquel et al. (2023), d Damasso et al. (2023), e Newton et al. (2023), f Nardiello et al. (2022), g Carleo et al. (2024), h Mantovan et al. (2024a).  \nthat the host stellar atmospheric abundances of refractory elements (such as Mg, Si, Fe) not only can they be considered a useful diagnostics of the composition of the initial protoplanetary disk, but also a proxy of the composition of very low-mass  \nplanets (Adibekyan et al. 2021). Although the relatively small sample analysed in this work is composed of only thin disk stars, we can conclude that our targets show values comparable with those by Adibekyan et al. (2015), with all low-mass planethosting stars (blue-filled squares) having [Mg / Si] around 0 dex and higher than those of the two targets hosting both high- and low-mass planets. TOI-4515, that hosts the only gas giant planet in our sample, shows a [Mg / Si] value around 0.02 dex, within the range between -0 . 05 and + 0 . 05dex found for higher-mass planet-hosting stars in Adibekyan et al. (2015) at similar [Fe / H].', '4.4.2. Key exoplanetary geochemical ratios': 'In Fig. 6 we show the distribution of our stars in the C / O 8 versus Mg / Si and Fe / Mgversus Mg / Si diagrams, compared with the sample of transiting planet-hosting stars studied by Biazzo et al. (2022). These elemental ratios are very useful mineralogical diagnostics because they govern the distribution and formation of chemical species in the protoplanetary disks, and, by extension, the composition of planets. On one hand, refractory elements (as traced, e.g. by the Mg / Si and Fe / Si ratios) with their high condensation temperatures ( T cond) condense close to the star. Their stellar ratios, therefore, remain roughly constant throughout the planet-forming disk, thus reflecting its composition and possibly the composition of the planetary atmosphere (see Brewer et al. 2016, and references therein). On the other hand, volatile elements (as traced, e.g. by C / O) have low T cond and their ratios are governed by ices and, therefore, their values in planetary composition can deviate significantly from those of the host stars (Thiabaud et al. 2015). If the C / Oratio is above a threshold value (typically of ∼ 0.8) there is almost no free oxygen available to form silicates and the geology will be dominated by carbonates, while below that threshold, planetesimal geology primarily consists of magnesium silicates. In the latter case, the exact compo-  \nFig. 5. [Mg / Si] versus [Fe / H] for our sample (blue points), where mean errors on [Mg / H], [Si / H], and [Fe / H] are similar to each other ( ∼ 0.090.10 dex; as a reference, typical errors in Mg, Si, and Fe abundances are represented in the bottom right corner of the plot). Overplotted in grey and orange are the results by Biazzo et al. (2022) and Adibekyan et al. (2015), respectively. As in Adibekyan et al. (2015), squares and circles represent targets hosting only planets with M p < 30 M E (low-mass planets) and M p > 30 M E (high-mass planets), respectively. Within our sample, we mark with diamonds those targets hosting both high- and lowmass planets and with crosses targets whose planets have not a mass estimate yet. Dashed lines represent the solar abundances.  \n<!-- image -->  \ntion of silicates in the mantle of a planet depends on the Mg distribution: for Mg / Si < ∼ 1, Mg forms orthopyroxene and Si is present as other silicate species such as feldspars or olivine; for 1 < ∼ Mg / Si < ∼ 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 (Bond et al. 2010, Thiabaud et al. 2015). Further, another diagnostic, such as the Fe / Mg ratio, can be considered as an indicator of the degree of core-mantle fractionation, thus being useful for determining the core size of low-mass planets. Values of Fe / Mg greater than ∼ 1 are interpreted as a propensity for the planets to have bigger (iron) cores (see, Wang et al. 2022, and references therein).  \nOur transiting planet host stars are all below the threshold of 0.8 in C / O and are mainly concentrated to values close to 0.4, with a mean value of < C / O >= 0.43 ± 0.12, thus consistent with silicate planets. Most of them are at the lower edge of the similar distribution done by Biazzo et al. (2022) for older stars. Only two targets show values slightly higher than the solar C / O ratio (yellow dot in the figure at C / O = 0.57 ± 0.04). This highlights how the remark by Turrini et al. (2022) on possible biases in the interpretation of planetary compositional models based on the use of solar values of C / O as reference stellar abundances stands also at younger evolutionary stages as those of our stars (see also Turrini et al. 2021; Biazzo et al. 2022; Jorge et al. 2022). Moreover, all our targets show Mg / Si between 1.0 and 1.3, with a mean value of 1.14 ± 0.07, only slightly lower than the solar reference of 1.17 ± 0.08 in Biazzo et al. (2022) and within the errors, while the Fe / Mgratio displayed in the lower panel of Fig. 6 shows a mean value of < Fe / Mg >= 0.83 ± 0.05, which is slightly higher, within the errors, than the solar value of 0.79 ± 0.07.  \nSummarising, the peak of the Mg / Si-C / Odistribution for our sample is consistent with mantles of Si which will take solid form as SiO -4 4 and SiO2 and Mg made of a mixture of pyroxene and olivine assemblages, analogously to what was obtained by, for example Suárez-Andrés et al. (2018) for solar-type planet-hosting stars and by Biazzo et al. (2022) for older transiting planet-hosting stars with similar atmospheric properties. At the same time, the peak of the Mg / Si-Fe / Mg distribution is close to the threshold where planetary cores start to be bigger (Wang et al. 2022).  \nFig. 6. C / O (upper panel) and Fe / Mg (lower panel) ratios versus Mg / Si. In grey, the results by Biazzo et al. (2022) are shown, while the yellow big dot in both panels represents the position of the Sun as derived by the same authors (see also Biazzo & GAPS Team 2023). Vertical lines represent Mg / Si = 1.0 (Thiabaud et al. 2015), while horizontal lines are plotted for C / O = 0.4, 0.8 (Suárez-Andrés et al. 2018), and for Fe / Mg = 0.9 (Wang et al. 2022).  \n<!-- image -->  \nAs already mentioned above, although the C / O ratio provides useful constraints on the potential formation locations of gaseous exoplanets, carbon and oxygen alone are not enough to determine where a given planet originated. Other ratios were therefore proposed in the literature, in particular those including refractory elements, which are present in solid phase throughout most of the disk and therefore they reflect the solid-to-gas accretion ratio of the material incorporated into the planetary envelope during the formation (see Turrini et al. 2021). Very recently, Chachan et al. (2023) proposed O / Si and C / Si as possible ratios useful to alleviate the degeneracies in the interpretation of atmospheric C / O ratios, making it easier to relate the planetary abundance ratios to the compositions of the solids in the disk. In Fig. 7 we show the C / Si versus O / Si ratios for our targets and those analysed by Biazzo et al. (2022). All targets are below the solar value in the C / Si ratio (with the exception of one target investigated by Biazzo et al. 2022; see that paper for de-  \ntails), consistently with the statements by Chachan et al. (2023) for planet-hosting stars with high content of refractory elements. At the same time, both samples of transiting planet-hosting stars show a wide distribution in O / Si, with a propensity for lowermass planet-hosting stars to possess higher values of O / Si. Considering the planet-hosting stars homogeneously analysed in this work and in Biazzo et al. (2022), the mean O / Si for both samples is indeed of 17 . 3 ± 3 . 0 and 12 . 7 ± 2 . 7 for stars hosting planets with M p < 30 M E and M p > 30 M E, respectively. By chance, all the five targets hosting only planets with M p < 30 M E appear to have super-solar O / Si ratio. This seems to be again consistent with the claim by Chachan et al. (2023) for which hot Jupiter host stars should show lower content of O / Si and host planets with silicate clouds and low water content. The possible relation between stellar O abundance and planetary mass will be discussed in Sect.4.5.  \nSimilar results were found by analysing the C / S ratio versus O / S, which were recently proposed by Crossfield (2023) as useful diagnostics of the formation of warm and hot giant planets. In the C / S vs O / S relation, shown in Fig.A.1, we used the same symbols as in Fig.7 to distinguish targets hosting low- and highmass planets. Similarly to the C / Si vs O / Si distribution, using sulphur S as refractory element we observe a tendency for lowermass planet-hosting stars to be present at higher values of O / S. The mean O / S values for stars hosting planets with M p < 30 M E and M p > 30 M E, are indeed 37.4 ± 3.2 and 30.7 ± 9.0 respectively.  \nWeare aware that these are presumably speculations because at least part of the trends could be due to selection e ff ects (e.g. stellar age and atmospheric parameters), but we think that these types of possible relationships that include multiple refractory and volatile elemental ratios should be studied widely and systematically in the next-coming and future big spectroscopic surveys.  \nFig. 7. C / Si versus O / Si for our stars (blue points) and those (in grey) analysed by Biazzo et al. (2022). Squares and circles represent targets hosting only planets with M p < 30 M E (low-mass planets) and M p > 30 M E (high-mass planets), respectively. Diamonds mark those targets hosting both higher- and lower-mass planets, while the cross symbol indicates that planet-hosting target without any planetary mass measurement. The typical uncertainty was estimated referring to the methodology used by Delgado Mena et al. (2010) and it is shown on the top left corner. Dashed lines represent the solar abundance ratios by Biazzo et al. (2022).  \n<!-- image -->', '4.4.3. Bulk composition estimates of low-mass planets': "Due to the above-mentioned limitations in the use of stellar volatile abundances to make assumptions on planetary composition, Wang et al. (2019) introduced the idea of using devolatilised host stellar abundances to constrain the bulk composition and interior modelling of hypothetical habitable-zone terrestrial exoplanets. This approach is limited to this kind of planets since it is based on the observations of bulk composition di ff erences and similarities of the Solar System's rocky bodies relative to the Sun. Among the major rock-forming elements analysed by the authors, only Ca and Al are not observed to be depleted in rocky bodies relative to the Sun, while for volatiles like O, S, and C the depletion is over 80%. Therefore, in rocky planets, elements like C and O are no longer valid indicators of the mantle oxidation state, which however is essential to understand the planetary interiors. As a consequence, the authors proposed the bulk (O -Mg -2Si) / Fe as a diagnostics of the oxygen fugacity, which is an indicator of the mantle oxidation state of a silicate terrestrial exoplanet, with a depletion of Mg, Fe, and Si of around 15% (see Wang et al. 2019 for details).  \nFigure 8 shows the distribution of our targets in the (O -Mg -2Si) / Fe versus Mg / Si and (O -Mg -2Si) / Fe versus Fe / Mgdiagrams, where we applied the devolatilisation model to only those stars hosting planets with mass lower than 10 M E (or with radius lower than 3 R E in case of TOI-2048 and TOI-5082 for which no planetary mass estimate was obtained by the time of this publication). In our sample, only four targets host planets satisfying the mass (or radius) requirements (namely TOI1136, TOI-1430, TOI-2048, and TOI-5082). Figure 8 considers the Mg / Si ratio as an indicator of the dominant mineral assemblages (i.e. olivine or pyroxene) in the mantle of a silicate planet. In the same figure, the Fe / Mgratio is used as a proxy of the coremantle fractionation degree, hence the core size (see details in Wang et al. 2022).  \nThe dashed lines divide our low-mass planets into di ff erent categories, in which di ff erent planetary interior properties may be expected, as indicated by Wang et al. (2022). Two of our four sample planets (orbiting TOI-1136 and TOI-5082) are below the unit line of (O -Mg -2Si) / Fe, implying a potential large iron core. These two planets with (O-Mg-2Si) / Fe lower than 1 show Fe / Mg values consistent within the errors with the Earth value (square symbol in Fig. 8), thus their mantle mineralogies would be more or less the same, as also confirmed by their similar values in the Mg / Si ratio. The iron in the low-mass planet around TOI-1430, with (O -Mg -2Si) / Fe > 1, may be more oxidised and locked in the mantle, thus developing comparably smaller core. TOI-2048 is missing in the plot because it has no oxygen abundance estimate (as shown in Table 2).", '4.5. Stellarabundancesversusplanetaryproperties': 'With the aim to look for possible relations between stellar abundances and planetary parameters, we analysed the [X / Fe] ratios for the volatile elements measured in the present work (namely, C, O, and S) as a function of the the total planetary mass of the planetary system ( M ptot = ∑ nplanets i M pi ), which we adopt as a proxy into the planet formation e ffi ciency of the original protoplanetary disks of the stars in the sample. As in Suárez-Andrés et al. (2017) and Biazzo et al. (2022), we do not find a significant trend with C and S, while a notable correlation seems to be present for the [O / Fe] ratio versus M ptot (shown in Fig. 9, top), with a decreasing trend towards greater masses. In particular, defining the significance p -value threshold α = 0 . 05  \nFig. 8. (O -Mg -2Si) / Fe versus Mg / Si (left panel) and (O -Mg -2Si) / Fe versus Fe / Mg (right panel) diagrams. The vertical and horizontal dashed lines (Wang et al. 2022) highlight the areas in which planets are expected to have di ff erent planetary interior properties, along with the positions of Mars (cross) and Earth (square) shown for reference (see Wang et al. 2022). Overplotted with diamonds are the terrestrial-type exoplanets analysed by Wang et al. (2022).  \n<!-- image -->  \n(where a statistically significant relation has p < α ), we find a Spearman statistical significance of p ∼ 0 . 01 for the oxygen decreasing trend with M ptot . For carbon and sulphur, the relation is not significant ( p ∼ 0 . 65 -0 . 93), with the [S / Fe] distribution showing higher dispersion.  \nThe observed trend between [O / Fe] and the total planetary mass has an intuitive physical interpretation when we look into the role of O in shaping the mass of heavy elements in stars and the amount of planet-forming material in protoplanetary disks. Specifically, O is the heavy element that contributes most mass to the stellar metallicity (about 45% in the case of the solar composition, Lodders 2010) and comparative studies between the solar composition and meteoritic abundances reveal that, mass-wise, it is a key component in the planetary building blocks already from the inner regions ( ∼ few au) of protoplanetary disks. In the Solar Nebula up to half the protosolar budget of O was incorporated into rocks and dust before the water snowline (Lodders 2010) while a large fraction of the remaining O condensed as water ice immediately after (Öberg & Bergin 2021; Turrini et al. 2021). Studies of the interstellar medium, circumstellar disks and polluted white dwarfs suggest that this key role of O is general to the Solar neighbourhood (Öberg & Bergin 2021; Trierweiler et al. 2023). Increases in the abundance of O should therefore translate into proportional increases in the mass of solid material that can support the growth of planetary bodies and enhance the planet formation e ffi ciency of protoplanetary disks.  \nTo test this we performed the following population synthesis study with our Monte Carlo implementation of the GroMiT code (Polychroni et al. 2024) previously used in Mantovan et al. (2024b). The code implements the prescriptions for the growth and migration of solid planets by pebble accretion from Johansen et al. (2019) and those for the growth and migration of gas-accreting giant planets by Tanaka et al. (2020), complementing them with the treatment of the solids-to-gas ratio in the pro-  \ntoplanetary disk from Turrini et al. (2023). We simulated the formation of synthetic planetary populations in two identical protoplanetary disks di ff ering only for their O abundances and compared the resulting planets. The first disk was characterised by solar composition and metallicity (1.4%; Asplund et al. 2009), while the second disk was enriched in O by 0.3 dex, resulting in a supersolar metallicity of 2.1%. Since the amount of O sequestered by rocks is controlled by the abundances of Fe, Mg and Si (Burrows & Sharp 1999; Fegley & Schaefer 2010) and these elements have the same abundances in both disks, the larger abundance of O in the O-enriched disk should a ff ect only the amount of solids beyond the water snowline.  \nTo correctly reproduce this e ff ect, we modified the prescription for the solid-to-gas ratio from (Turrini et al. 2023) between the two disks. Specifically, we assume both disks being characterised by the same solid-to-gas ratio in their inner and hotter regions dominated by rocks, with the O-enriched disk having a larger solid-to-gas ratio beyond the water snowline to represent its larger budget of water ice. As a result, in the solar disk the solid-to-gas ratio inside the water snowline is set to 0.7% (0.5x the disk metallicity, Turrini et al. 2023) and is 1.05% (0.75x the disk metallicity, Turrini et al. 2023) beyond the water snowline based on the case of the Solar System (see Turrini et al. 2021, for a detailed discussion). In the O-enriched disk, the solid-to-gas ratio inside the water snowline is kept to 0.7% and is increased to 1.58% beyond the water snowline.  \nFor each disk we simulated the growth and migration of 5 × 10 4 planetary seeds uniformly distributed in semimajor axis across the disk extension. The initial and boundary conditions of the two population synthesis simulations are summarised in Table 6. From the resulting planetary populations we extracted all planets with mass greater than about 5 M E and divided them over two classes based on those considered in Biazzo et al. (2022) and in this work: massive planets between 30 and 954 M E (roughly  \nTable 6. Initial parameters used to run the modified GroMiT code.Table 7. Results from the implemented GroMiT code for the O-enriched disk and the solar disk, and divided for the planet mass ranges (5-30 M E and 30-954 M E) and for their formation period ( < 1 Myr and > 1 Myr).  \n0.1-3 M J) and smaller planets with masses between 5 and 30 M E. The higher metallicity of the O-enriched disk (2.1% vs 1.4%) makes it 1.7 × more e ffi cient in forming planets above 5 M E than the solar disk, in line with previous studies on the role of disk metallicity on planet formation e ffi ciency (Savvidou & Bitsch 2023).  \nIn the O-enriched disk the formation of the smaller planets is boosted more (2.5 × with respect to the solar disk) than that of the more massive planets (1.4 × with respect to the solar disk). In the O-enriched disk the smaller planets represent 35% of the extracted population, while in the solar disk they represent only 24% (see Table 7). We further divided each class of planets into two sub-classes: early-born planets formed in the first 1 Myr and late-born planets formed beyond 1 Myr. In the O-enriched disk there is a global increase in the relative abundance of lateborn planets, with the abundance of late-born planets between 5 and 30 M E increasing twice as much as that of late-born planets above 30 M E. These results show that O-rich stars have a higher probability of forming planets between 5 and 30 M E and support the interpretation provided here and in Biazzo et al. (2022) to the observed trends between [O / Fe] and the total planetary mass. In particular, the large contribution of O to the mass fraction of heavy elements in stars and their disks appears e ff ective in promoting the formation of low-mass planets since the earlier stages of their circumstellar disk lifetime.  \nFig. 9. [O / Fe] as well as C / Oand Mg / Si elemental abundance ratios as a function of the total planetary mass of the planetary system. Filled and empty symbols represent the targets that within our sample host only one confirmed transiting planet or multi-planetary systems, respectively. The targets by Biazzo et al. (2022) are overplotted in grey circles. The dashed grey lines represent the linear fits on the data points.  \n<!-- image -->  \nIn the central and bottom plots of Fig. 9, we show the C / O and Mg / Si ratios calculated in the present work as a function of the total planetary mass, respectively. As done for the abundances of volatile elements, we searched for possible trends between C / O and Mg / Si ratios and the total planetary masses. We find hints of a relation describing a decreasing trend with a rank correlation coe ffi cient of ρ ∼ -0 . 27 and Spearman statistical significance p ∼ 0 . 12 for Mg / Si. Instead, for the C / O ratio we obtain p ∼ 0 . 007 and a rank correlation coe ffi cient of ρ ∼ 0 . 52 showing an increasing trend with M ptot . Some evidence of a downward trend of Mg / SiM ptot (higher values of Mg / Si for lower-mass planets) and a slightly increasing trend of C / OM ptot (higher values of C / O for higher-mass planets) were also reported by Mishenina et al. (2021), but no conclusion was drawn because of the large scatter of the elemental ratios as a function of planetary masses.', '5. Conclusions': "In this paper, we presented the analysis of a sample of ten GK TESS stars with T e ff > 5000K, v sin i < 10km / s, and ages ranging between 0.2 and 2 Gyr placed in the solar neighbourhood and hosting at least one confirmed transiting planet. This study was conducted in the framework of the GAPS sub-programme focussed on searching for exoplanets around young objects. We performed an accurate and homogeneous characterisation of relatively young planet-hosting stars from high-resolution spectra acquired with the HARPS-N spectrograph at the Telescopio Nazionale Galileo . We derived the kinematic properties, atmospheric parameters, and abundances of 18 elements. In the following, we illustrate our main results:  \n- -The chemical and kinematic analysis allowed us to confirm that our targets belong to the Galactic thin disk.\n- -We estimated stellar ages through the lithium line and the use of empirical isochrones, thus confirming the youth of our sample and the compatibility for most of the targets with the ages derived from gyrochronological and activity indicators, and those available in the literature.\n- -We find that all our targets show a C / Oratio smaller than 0.8, with most of them being below the solar value of ∼ 0.57 and in the lower edge of older stars that host transiting planets. This result, together with the mean Mg / Si ratio being around 1.14, suggests that our sample of planet-hosting stars is consistent with silicate mantles with Mg mostly distributed between olivines and pyroxenes. The Fe / Mg ratio indicates a propensity for bigger planetary cores.\n- -From C / Si and O / Si ratios, we observed a prevalence for transiting planet-hosting stars with a high O / Si ratio to host low-mass planets. Meanwhile, the C / Si ratio is sub-solar and consistent with a high content of refractory elements.\n- -We applied a devolatilisation model and used the bulk (OMg-2Si) / Fe as diagnostics of the oxygen fugacity for terrestrial planets. We find that two out of four very low-mass ( < 10 M E) planets show (O-Mg-2Si) / Fe values compatible with the possible presence of large iron cores and Fe / Mg and Mg / Si consistent with Earth's values.\n- -Concerning volatiles, we confirm the tendency previously found in old transiting planet hosts in which the oxygen element seems to be higher for stars hosting lower-mass companions. This possible trend seems to extend also to relatively young systems the recent suggestion that was reported in the literature according to which the formation of lowmass planets is favoured at higher stellar oxygen values.  \nThanks to JWST and the future launch of Ariel in 2029, it will be possible to probe exoplanet atmospheres of young transiting planets in detail, allowing for insight into the atmospheric composition and the possible presence of clouds and hazes. This will also allow us to get insight into the formation and migration mechanisms of young exoplanets that are expected to be better understood with the future findings of the dedicated ELT instruments (such as ANDES).  \nAcknowledgements. We are very grateful to the referee Haiyang Wang for providing insightful and stimulating comments that have helped us to improve the presentation of our results. This research has made use of the Exoplanet Follow-up Observation Programme (ExoFOP; DOI: 10.26134 / ExoFOP5) website, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Programme, and also of the INSPECT database, version 1.0 (www.inspect-stars.net). This work has made use of the SIMBAD database and the VizieR catalogue access tool operated at CDS (Strasbourg, France), and also of the Hypatia Catalog Database, an online compilation of stellar abundance  \ndata, which was supported by NASA's Nexus for Exoplanet System Science (NExSS) research coordination network and the Vanderbilt Initiative in DataIntensive Astrophysics (VIDA). This work has been financially supported by the PRIN-INAF 2019 'Planetary systems at young ages' (PLATEA) and by the grant INAF 2022 TRAME@JWST (TRacing the Accretion Metallicity rElationship with NIRSpec@JWST). DP acknowledges the support from the Istituto Nazionale di Oceanografia e Geofisica Sperimentale (OGS) and CINECA through the program 'HPC-TRES (High Performance Computing Training and Research for Earth Sciences)' award number 2022-05 as well as the support of the ASI-INAF agreement n 2021-5-HH.1-2022. GMa acknowledges support from CHEOPS ASI-INAF agreement n. 2019-29-HH.0.", 'References': "Adamow, M. M. 2017, in American Astronomical Society Meeting Abstracts,  \nVol. 230, American Astronomical Society Meeting Abstracts #230, 216.07 Adibekyan, V., Dorn, C., Sousa, S. G., et al. 2021, Science, 374, 330 Adibekyan, V. Z., Santos, N. C., Figueira, P., et al. 2015, A&A, 581, L2 Adibekyan, V. Z., Santos, N. C., Sousa, S. G., & Israelian, G. 2011, A&A, 535, L11  \n- Adibekyan, V. Z., Sousa, S. G., Santos, N. C., et al. 2012, A&A, 545, A32 Amarsi, A. M., Asplund, M., Collet, R., & Leenaarts, J. 2015, MNRAS, 454, L11\n- Amarsi, A. M., Nissen, P. E., & Skúladóttir, Á. 2019, A&A, 630, A104  \nAsplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481  \nBaratella, M., D'Orazi, V., Biazzo, K., et al. 2020a, A&A, 640, A123  \nBaratella, M., D'Orazi, V., Carraro, G., et al. 2020b, A&A, 634, A34  \n- Bensby, T., Feltzing, S., & Oey, M. S. 2014, A&A, 562, A71\n- Biazzo, K., D'Orazi, V., Desidera, S., et al. 2022, A&A, 664, A161  \nBiazzo, K. & GAPS Team. 2023, in Memorie della Societa Astronomica Italiana,  \nVol. 94, 189  \n- Bond, J. C., O'Brien, D. P., & Lauretta, D. S. 2010, The Astrophysical Journal, 715, 1050\n- Brewer, J. M., Fischer, D. A., Valenti, J. A., & Piskunov, N. 2016, The Astrophysical Journal Supplement Series, 225, 32\n- Burrows, A. & Sharp, C. M. 1999, ApJ, 512, 843\n- Ca ff au, E., Ludwig, H.-G., Ste ff en, M., et al. 2008, A&A, 488, 1031  \nCarleo, I., Malavolta, L., Desidera, S., et al. 2024, A&A, 682, A135  \nCarleo, I., Malavolta, L., Lanza, A. F., et al. 2020, A&A, 638, A5  \n- Castelli, F. & Kurucz, R. L. 2003, in Modelling of Stellar Atmospheres, ed. N. Piskunov, W. W. Weiss, & D. F. Gray, Vol. 210, A20\n- Chachan, Y., Knutson, H. A., Lothringer, J., & Blake, G. A. 2023, The Astrophysical Journal, 943, 112\n- Claret, A. 2019, Research Notes of the AAS, 3, 17  \nCo¸skunoˇglu, B., Ak, S., Bilir, S., et al. 2011, MNRAS, 412, 1237  \n- Cosentino, R., Lovis, C., Pepe, F., et al. 2012, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 8446, Ground-based and Airborne Instrumentation for Astronomy IV, ed. I. S. McLean, S. K. Ramsay, &H. Takami, 84461V  \nCovino, E., Esposito, M., Barbieri, M., et al. 2013, A&A, 554, A28  \nCrossfield, I. J. M. 2023, ApJ, 952, L18  \n- Cummings, J. D., Deliyannis, C. P., Maderak, R. M., & Steinhauer, A. 2017, AJ, 153, 128\n- Dai, F., Masuda, K., Beard, C., et al. 2023, AJ, 165, 33\n- Damasso, M., Locci, D., Benatti, S., et al. 2024, in press, arXiv e-prints, arXiv:2408.10629  \nDamasso, M., Locci, D., Benatti, S., et al. 2023, A&A, 672, A126  \n- Delgado Mena, E., Israelian, G., Hernández, J. I. G., et al. 2010, The Astrophysical Journal, 725, 2349  \nDelgado Mena, E., Tsantaki, M., Adibekyan, V. Z., et al. 2017, A&A, 606, A94  \nDesidera, S., Damasso, M., Gratton, R., et al. 2023, A&A, 675, A158  \nDoyle, A. P., Davies, G. R., Smalley, B., Chaplin, W. J., & Elsworth, Y. 2014,  \nMonthly Notices of the Royal Astronomical Society, 444, 3592  \nEggen, O. J. 1996, AJ, 111, 1615  \nFegley, B. & Schaefer, L. 2010, in Astrophysics and Space Science Proceedings,  \nVol. 16, Principles and Perspectives in Cosmochemistry, 347  \nFonte, S., Turrini, D., Pacetti, E., et al. 2023, MNRAS, 520, 4683  \nFrancis, C. & Anderson, E. 2009, New A, 14, 615  \n- Gagné, J., Lafrenière, D., Doyon, R., Malo, L., & Artigau, É. 2014, ApJ, 783, 121  \nGagné, J., Mamajek, E. E., Malo, L., et al. 2018, ApJ, 856, 23 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2021, A&A, 649, A1 Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1 Gardner, J. P., Mather, J. C., Abbott, R., et al. 2023, PASP, 135, 068001 Gilmore, G., Randich, S., Worley, C. C., et al. 2022, A&A, 666, A120 Hedges, C., Hughes, A., Zhou, G., et al. 2021, AJ, 162, 54 Hinkel, N. R., Timmes, F. X., Young, P. A., Pagano, M. D., & Turnbull, M. C.  \n2014, AJ, 148, 54  \nJe ff ries, R. D., Jackson, R. J., Wright, N. J., et al. 2023, Monthly Notices of the Royal Astronomical Society, 523, 802 Johansen, A., Ida, S., & Brasser, R. 2019, A&A, 622, A202 Johansson, S., Litzén, U., Lundberg, H., & Zhang, Z. 2003, ApJ, 584, L107 Johnson, D. R. H. & Soderblom, D. R. 1987, AJ, 93, 864 Jorge, D. M., Kamp, I. E. E., Waters, L. B. F. M., Woitke, P., & Spaargaren, R. J. 2022, A&A, 660, A85 Kolecki, J. R. & Wang, J. 2022, AJ, 164, 87 Kupka, F. G., Ryabchikova, T. A., Piskunov, N. E., Stempels, H. C., & Weiss, W. W. 2000, Baltic Astronomy, 9, 590 Lind, K., Asplund, M., & Barklem, P. S. 2009, A&A, 503, 541 Lind, K., Asplund, M., Barklem, P. S., & Belyaev, A. K. 2011, A&A, 528, A103 Lodders, K. 2010, in Astrophysics and Space Science Proceedings, Vol. 16, Principles and Perspectives in Cosmochemistry, 379 MacDougall, M. G., Petigura, E. A., Gilbert, G. J., et al. 2023, AJ, 166, 33 Madhusudhan, N., Amin, M. A., & Kennedy, G. M. 2014, ApJ, 794, L12 Magrini, L., Danielski, C., Bossini, D., et al. 2022, A&A, 663, A161 Magrini, L., Viscasillas Vázquez, C., Spina, L., et al. 2023, A&A, 669, A119 Maldonado, J., Micela, G., Baratella, M., et al. 2020, A&A, 644, A68 Mamajek, E. E. & Hillenbrand, L. A. 2008, ApJ, 687, 1264 Mann, A. W., Johnson, M. C., Vanderburg, A., et al. 2020, AJ, 160, 179 Mantovan, G., Malavolta, L., Desidera, S., et al. 2024a, A&A, 682, A129 Mantovan, G., Malavolta, L., Locci, D., et al. 2024b, A&A, 684, L17 Marfil, E., Tabernero, H. M., Montes, D., et al. 2021, A&A, 656, A162 Messina, S., Nardiello, D., Desidera, S., et al. 2022, A&A, 657, L3 Mishenina, T., Basak, N., Adibekyan, V., Soubiran, C., & Kovtyukh, V. 2021, MNRAS, 504, 4252 Nardiello, D., Malavolta, L., Desidera, S., et al. 2022, A&A, 664, A163 Newton, E. R., Rampalli, R., Kraus, A. L., et al. 2023, AJ, 165, 136 Nissen, P. E. 2004, in Origin and Evolution of the Elements, ed. A. McWilliam &M. Rauch, 154 Öberg, K. I. & Bergin, E. A. 2021, Phys. Rep., 893, 1 Öberg, K. I., Murray-Clay, R., & Bergin, E. A. 2011, ApJ, 743, L16 Orell-Miquel, J., Lampón, M., López-Puertas, M., et al. 2023, A&A, 677, A56 Pacetti, E., Turrini, D., Schisano, E., et al. 2022, ApJ, 937, 36 Pasquini, L., Biazzo, K., Bonifacio, P., Randich, S., & Bedin, L. R. 2008, A&A, 489, 677 Pecaut, M. J. & Mamajek, E. E. 2013, ApJS, 208, 9 Pepe, F., Mayor, M., Rupprecht, G., et al. 2002, The Messenger, 110, 9 Polychroni, D., Turrini, D., & Pirani, S. 2024, GroMiT: Planet Growth and Migration Track code Randich, S., Gilmore, G., Magrini, L., et al. 2022, A&A, 666, A121 Reddy, B. E., Lambert, D. L., Laws, C., Gonzalez, G., & Covey, K. 2002, MNRAS, 335, 1005 Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 9143, Space Telescopes and Instrumentation 2014: Optical, Infrared, and Millimeter Wave, ed. J. Oschmann, Jacobus M., M. Clampin, G. G. Fazio, & H. A. MacEwen, 914320 Savvidou, S. & Bitsch, B. 2023, A&A, 679, A42 Sestito, P. & Randich, S. 2005, A&A, 442, 615 Sneden, C. 1973, ApJ, 184, 839 Sousa, S. G., Adibekyan, V., Delgado-Mena, E., et al. 2021, A&A, 656, A53 Sousa, S. G., Santos, N. C., Adibekyan, V., Delgado-Mena, E., & Israelian, G. 2015, A&A, 577, A67 Spitoni, E., Recio-Blanco, A., de Laverny, P., et al. 2023, A&A, 670, A109 Strassmeier, K. G., Weber, M., Gruner, D., et al. 2023, A&A, 671, A7 Suárez-Andrés, L., Israelian, G., González Hernández, J. I., et al. 2017, A&A, 599, A96 Suárez-Andrés, L., Israelian, G., González Hernández, J. I., et al. 2018, A&A, 614, A84 Szentgyorgyi, A. H. & Furész, G. 2007, in Revista Mexicana de Astronomia y Astrofisica Conference Series, Vol. 28, Revista Mexicana de Astronomia y Astrofisica Conference Series, ed. S. Kurtz, 129-133 Tanaka, H., Murase, K., & Tanigawa, T. 2020, ApJ, 891, 143 Thiabaud, A., Marboeuf, U., Alibert, Y., Leya, I., & Mezger, K. 2015, A&A, 574, A138 Tinetti, G., Eccleston, P., Haswell, C., et al. 2021, arXiv e-prints, arXiv:2104.04824 Tokovinin, A., Fischer, D. A., Bonati, M., et al. 2013, PASP, 125, 1336 Trierweiler, I. L., Doyle, A. E., & Young, E. D. 2023, The Planetary Science Journal, 4, 136 Turrini, D., Codella, C., Danielski, C., et al. 2022, Experimental Astronomy, 53, 225 Turrini, D., Marzari, F., Polychroni, D., et al. 2023, A&A, 679, A55 Turrini, D., Schisano, E., Fonte, S., et al. 2021, ApJ, 909, 40 Wang, H. S., Liu, F., Ireland, T. R., et al. 2019, Monthly Notices of the Royal Astronomical Society, 482, 2222 Wang, H. S., Quanz, S. P., Yong, D., et al. 2022, Monthly Notices of the Royal Astronomical Society, 513, 5829 Yun, S., Lee, Y. S., Kim, Y. K., et al. 2024, arXiv e-prints, arXiv:2402.19321", 'Appendix A: Additional figure': 'Following the analysis in Sec.4.4.2, Fig.A.1 shows the relation between the diagnostic elemental ratios C / S and O / S, with the symbols used to distinguish targets hosting low- and high-mass planets as in Fig.7.  \nFig. A.1. C / S versus O / S for our stars (blue points) and those (in grey) analysed by Biazzo et al. (2022). Symbols as in Fig.7. Dashed lines represent the solar abundance ratios by Biazzo et al. (2022). The typical uncertainty was estimated referring to the methodology used by Delgado Mena et al. (2010) and it is represented on the top left corner of the plot.  \n<!-- image -->'}