# The GAPS Programme at TNG. LI. Investigating the correlations between transiting system parameters and host chromospheric activity

R. Claudi, G. Bruno, L. Fossati, A. F. Lanza, A. Maggio, G. Micela, J.

Maldonado, S. Benatti, K. Biazzo, A. Bignamini, et al.

► To cite this version:

R. Claudi, G. Bruno, L. Fossati, A. F. Lanza, A. Maggio, et al. The GAPS Programme at TNG. LI. Investigating the correlations between transiting system parameters and host chromospheric activity. Astronomy & Astrophysics - A&A, 2024, 682, 10.1051/0004-6361/202347079. insu-04851594

HAL Id: insu-04851594

https://insu.hal.science/insu-04851594v1

Submitted on 20 Dec 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# The GAPS Programme at TNG

## LI. Investigating the correlations between transiting system parameters and host chromospheric activity

R. Claudi<fnref n="1" />,², G. Bruno³, L. Fossati⁴, A. F. Lanza³, A. Maggio⁵, G. Micela⁵, J. Maldonado⁵, S. Benatti⁵, K. Biazzo⁶, A. Bignamini⁷, L. Cabona<fnref n="1" />, I. Carleo⁸,⁹, C. Daniëlski¹⁰,¹¹, S. Desidera<fnref n="1" />, L. Malavolta<fnref n="1" />,¹², L. Mancini¹³,¹⁴,⁹, M. Montalto³, D. Nardiello<fnref n="1" />, M. Rainer¹⁵, G. Scandariato³, A. Sozzetti⁹, R. Cosentino¹⁶, E. Covino¹⁷, L. Di Fabrizio¹⁶, A. Ghedina¹⁶, V. Lorenzi¹⁶, E. Molinari¹⁸, M. Molinaro⁷, I. Pagano³, G. Piotto¹², and E. Poretti¹⁵

*(Affiliations can be found after the references)*

Received 2 June 2023 / Accepted 14 November 2023

## ABSTRACT

**Context.** Stellar activity is the most relevant types of astrophysical noise that affect the discovery and characterization of extrasolar planets. On the other hand, the amplitude of stellar activity could hint at an interaction between the star and a close-in giant planet. Progress has been made in recent years in understanding how to deal with stellar activity and search for observational evidence of star-planet interactions.

**Aims.** The aim of this work is to characterize the chromospheric activity of stars hosting short-period exoplanets by studying the correlations between the chromospheric emission (CE) in the Ca II H&K and the planetary parameters.

**Methods.** We measured CE in the Ca II H&K lines using more than 1900 high-resolution spectra of a sample composed of 76 targets, observed with the HARPS-N spectrograph between 2012 and 2020. We transformed the fluxes into bolometric- and photospheric-corrected chromospheric emission ratios, $R'_{HK}$. Furthermore, we completed the sample of hosts digging for data in previous works. Stellar parameters $T_{eff}$, B-V, and V were retrieved homogeneously from the Gaia DR3. Then, $M_\star$, $R_\star$, and ages were determined from isochrone fitting. We retrieved planetary data from the literature and catalogs. The search for correlations between the $log(R'_{HK})$ and planetary parameters have been performed through both Spearman's rank and its statistics as well as the more sophisticated Gaussian mixture model method.

**Results.** We found that the distribution of $log(R'_{HK})$ for the transiting planet hosts is different from the distribution of field main-sequence and sub-giant stars. The $log(R'_{HK})$ of planetary hosts is correlated with planetary parameters proportional to the planetary radius to the power of $n (R_p^n)$, indicating a common origin for the correlations. The statistical analysis has also highlighted four clusters of host stars with different behavior in terms of their stellar activity with respect to the planetary surface gravity. Some of the host stars have a value of $log(R'_{HK})$ that is lower than the basal level of activity for main sequence stars. The planets of these systems are very close to filling their Roche lobe, suggesting that they evaporate through hydrodynamic escape under the strong irradiation of the host star, creating shrouds that absorb the core of the chromospheric resonance lines.

**Key words.** planets and satellites: fundamental parameters – planet-star interactions – stars: activity – planetary systems

## 1. Introduction

The majority of the 5000 exoplanets known thus far have been discovered thanks to transit and radial velocity (RV) methods<fnref n="1" /> (Schneider et al. 2011). The exoplanet population displays a tremendous diversity in terms of masses, radii, and system architectures. The broad diversity among the extrasolar planets is also reflected in the range of atmospheric compositions that has been emerging from several studies on transmission and emission spectroscopy data (for a review, see, e.g., Madhusudhan 2019). The atmospheric chemical compositions of these planets are expected to provide information on different scenarios of planetary formation and evolution. A severe obstacle in gathering this information is presented by the stellar noise, mainly due to the chromospheric activity (CA) of stars. The induced stellar variability due to the interaction between the turbulent stellar plasma and the existent magnetic field causes modifications to

the stellar spectral features that could be confused with planetary ones during a planetary transit or the observation of a phase curve (e.g., Ballerini et al. 2012; Robertson et al. 2015; Cegla et al. 2018; Barclay et al. 2021).

On the other hand, the stellar CA is not only a nuisance with the capability of hiding or mimicking the presence of a planetary companion and its atmospheric features; in fact, the enhancement of the CA of a stellar host is also considered to be an effect brought on by the presence of a hot-Jupiter orbiting very close to its host star via star-planet interactions (SPI) or it could be modulated or even affected by the presence of evaporating material shrouds lost by the close-in planet. The SPI can be induced by tidal and magnetic interactions (Cuntz et al. 2000). In the former case, planets can induce tidal bulges on the star with an interaction, depending on the semi-major axis of the planet, which could lead to an enhancement of the CA of the star ($\propto a_p^{-3}$). In the latter case, the SPI is driven by the planetary magnetic field that interacts with the stellar magnetic field

<footnote n="1">¹ http://exoplanet.eu/</footnote>
and the intensity of the interaction is $\propto a_p^{-2}$. Also, included in this case is a possible enhancement of the host CA. The first studies on the induced activity on planet hosts were performed by Bastian et al. (2000) in the radio wavelength range and by Saar & Cuntz (2001) in the optical, with poor results due to the insufficient sensitivity of their measurements. Later, Shkolnik et al. (2003, 2005, 2008), Kashyap et al. (2008), Scharf (2010), Gurdemir et al. (2012), Pillitteri et al. (2014), and Maggio et al. (2015) reported the detection of a significant increase of stellar activity in stars hosting hot Jupiters (hot gas-giant planets orbiting very close to their host stars). However, these studies were disputed by several authors, such as Poppenhaeger et al. (2010), Poppenhaeger & Schmitt (2011), Scandariato et al. (2013), and Miller et al. (2015), who attributed those detections to biases or selection effects. Furthermore, Canto Martins et al. (2011) searched for possible effects of SPI on the stellar chromosphere and did not find any significant correlations between the chromospheric activity (CA) indicator $log(R'_{HK})$ and planetary parameters.

Jupiter-like planets have a hydrogen-dominated atmosphere that evaporates when subjected to the EUV radiation and X-ray ($\lambda < 91.2$ nm) produced by their stellar host. In particular, the study of hot-Jupiters gave some evidence that mass loss can be a common phenomenon (Lammer et al. 2003; Lecavelier des Étangs et al. 2004; Lecavelier des Étangs 2007; Penz et al. 2008; Sanz-Forcada et al. 2010, 2011; Locci et al. 2019). A typical example is the extreme hot-Jupiter WASP-12 ($P_{orb} = 1.09$ d, Hebb et al. 2009), which is characterized by an exosphere surrounding the planet that fills its Roche lobe (Fossati et al. 2010; Haswell et al. 2012). The resulting mass loss seems to reinforce the circumstellar gas shroud that veils the host stars in some systems.

The most widely used CA indicator for F, G, and K stars is $log(R'_{HK})$, derived from the ratio between the Ca II H&K line-core chromospheric excess flux, evaluated by the Mount Wilson index, S, and the bolometric emission of the star (Noyes et al. 1984). Main sequence stars typically show $log(R'_{HK}) > -5.1$, a basal level exhibited by the quiet Sun. The determination of the $log(R'_{HK})$ index for the Sun has been discussed by, for instance, Egeland et al. (2017), while the recent work by Milbourne et al. (2019) shows that our star had a rotational modulation amplitude $\Delta log(R'_{HK}) \sim 0.015$, while the maximum variation in its chromospheric index was ~0.05 between July 2015 and October 2017, when the Sun was observed as a star by the HARPS-N solar telescope. Such amplitudes of variation appear to be broadly similar to those observed by Baliunas et al. (1995) in low-activity main sequence stars with a colour index comparable to the Sun over three decades. Therefore, they can be regarded as representative of the intrinsic variations in the $log(R'_{HK})$ index of quiet solar-like stars.

The anomalously low value of the WASP-12's Ca II H&K line-core flux ($log(R'_{HK}) \sim -5.5$) is interpreted as due to the presence of a circumstellar shroud of gas that is able to depress the stellar chromospheric flux to the observed level (Fossati et al. 2013). In subsequent works, other hosts of close-in planets were presented with values of $log(R'_{HK}) < -5.1$ and they all are likely viewed through shrouds of circumstellar gas coming from planetary atmosphere ablation (Staab et al. 2017; Haswell et al. 2020).

Following the study of Knutson et al. (2010) and considering a sample of 39 transiting planets, Hartman (2010) reported a correlation between the stellar CA ($log(R'_{HK})$) and the planetary surface gravity of transiting hot Jupiters (HJs) with a confidence level of 99.5%. Specifically, this author

found that the chromospheric emission index is lower for stars with planets with lower surface gravity. In the following years, Figueira et al. (2014) reviewed the $log(R'_{HK}) - log(g_p)$ correlation by collecting the available data on $log(R'_{HK})$ and stellar and planet properties and re-evaluated the correlation in the light of an extended dataset (108 transiting planets). A theoretical model was elaborated by Lanza (2014) to explain this correlation. In that model, the planetary material, ablated by the action of the UV radiation of the host star spreads towards the star along the magnetic field lines of the extended stellar corona in which the planet is lodged. The material is retained by the magnetic field lines forming prominence-like structures able to absorb at the core of chromospheric resonance lines (Mg II H&K, Ca II H&K). Strong absorption of the chromospheric resonance lines reduces the measured value of the $log(R'_{HK})$ indicator and the effect is stronger when the gravity of the planet is lower, explaining the correlation.

Exploiting a set of 41 transiting hosts, Fossati et al. (2015) improved the Lanza (2014) model by treating the sample of stars as a mixture of two subsets with intrinsic CA above or below the Vaughan-Preston gap (Vaughan & Preston 1980) observed for solar-like stars. Collier Cameron & Jardine (2018) highlighted an alternative explanation for the existence of the CA versus planetary surface gravity correlation, interpreting it as a direct consequence of the resulting distribution of orbital distances for transiting planets. Collier Cameron & Jardine (2018) found that the detection probabilities for transiting planets scale inversely with the increase of the planetary tidal migration rate, creating a statistical bias toward more massive planets being more likely to be observed around young and more active stars. Because the surface gravity of giant planets is mostly directly proportional to their mass, this hypothesis explains the correlation between CA and the surface gravity found by Hartman (2010) and confirmed by Figueira et al. (2014). However, the bimodality of the correlation proposed by Fossati et al. (2015) does not seem to be interpretable in the framework of this hypothesis.

In this paper, we evaluate $log(R'_{HK})$ for an extended sample of 76 transiting planet hosts observed over a time frame of about 9 yr by the Global Architecture of Planetary Systems (GAPS; Covino et al. 2013) collaboration. To this sample, we added original $S_{HK}$ data from other 60 transiting hosts drawn from the literature. Finally, we evaluated, in a uniform way, the $log(R'_{HK})$ value for both samples to review several correlations between $log(R'_{HK})$ and planetary parameters using this extended dataset.

The observations and data of this study are described in Sect. 2. In Sect. 3, we provide a description and characterization of the host stars. In Sect. 4, we explain how we obtained the chromospheric index for each star and the method used to evaluate the same for the data obtained by previous works. In Sect. 5, we discuss the time variation in the CA for the GAPS targets and we put the transiting hosts into the context of the other stars in Sect. 6. Finally, in Sects. 7 and 8, we present our analysis. In Sect. 9, we give our summary and conclusions.

## 2. Observations and data

The GAPS project is an Italian collaboration aimed at exploiting the HARPS-N spectrograph (Cosentino et al. 2012) at the Telescopio Nazionale Galileo (TNG) to search for new planets in known planetary systems and around other peculiar stars (of low metallicity and M stars) by means of high precision radial velocity time series. The project began in 2012 and completed its first phase in 2017, after 5 yr of observations. The second phase of GAPS started in 2017, triggered by
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** GAPS observations of transiting hosts considered in this work. Top: number of spectra per star of the GAPS sample of transiting hosts. Middle: number of observing nights per star. Bottom: S/N value per star.</figcaption></figure>

the implementation of the new TNG observing mode GIARPS (Claudi et al. 2016, 2017), which allows for the simultaneous use of HARPS-N and GIANO-B, covering both the optical (380–690 nm) and near-infrared (970–2945) wavelength ranges at the same time. In this second phase, GAPS has been devoted to studies of young planets (Carleo et al. 2018, 2020; Damasso et al. 2020; Benatti 2021; Nardiello et al. 2022) and planetary atmospheres (Borsa et al. 2019; Pino et al. 2020; Guilluy et al. 2020; Carleo et al. 2022), and this project is still ongoing. In both phases of GAPS, planetary systems with transiting planets were observed, in some cases collecting quite long time series of spectral data that have also proven useful in characterizing the CA level of such targets. In this study, we consider 76 transiting planet hosts observed in the past years by GAPS, to which we added 19 objects available from Knutson et al. (2010), and Staab et al. (2017), along with another 41 objects taken from Figueira et al. (2014). The complete sample contains 136 objects.

## 2.1. GAPS sample

The 76 stars in the original GAPS sample were regularly monitored with the spectrograph HARPS-N since 2012, collecting data for about five to 8 yr, in some cases. In the following, we discarded four objects (KELT-7, KELT-9, KELT-20, and WASP-33) for the reasons explained in Sect. 4. The journal of the observations for the remaining 72 objects is reported in Table A.1, including the date, $T_0$, in which the target has been observed for the first time, the total time span in days covered by GAPS data, and the number of available HARPS-N spectra. Figure 1 shows the distribution of spectra per star (top panel) and the distribution of the number of observing nights per star (central panel). For the transiting planets hosts, we have a total of 1910 spectra gathered in 867 observing nights, with a minimum of two spectra for WASP-3 and the highest number of spectra (123) for

HD189733. The median is 21 spectra per star. The high number of spectra obtained for some stars is due to the frequent observing cadence during the planetary transit night we observed to detect the Rossiter-McLaughlin effect on the radial velocity curves (e.g., Esposito et al. 2017). This is the case of WASP-11 with 26 spectra taken in one night, while the maximum number of nights has been devoted to XO-2N observed for 56 nights in about 3 yr.

The main aim of the GAPS observations is to obtain precise radial velocity (RV) measurements for new exoplanet detections. To achieve this goal, in our observations, we exploited the simultaneous reference technique using both optical fibers: the target fed the first fiber, while the second fiber was illuminated by a stable calibration source, either ThAr or Fabry Perot lamp<fnref n="2" />. If the stars are fainter than 11 mag, we used the second fiber instead for sky subtraction to avoid lamp pollution. The typical exposure time was 900 s to average out solar-like pulsations. The median value of the V magnitude of GAPS targets is 11.2 mag, so the median value of the reached signal-to-noise ratio (S/N) is 44 (see the bottom panel of Fig. 1). For the brighter tail of the target distribution, the single observations can reach S/N > 100 per pixel at 550 nm, to achieve a radial velocity precision of 1 m s⁻¹ or better, depending on the airmass of the observation and the weather conditions. The spectra were reduced through the HARPS-N data reduction software (DRS), including the radial velocity extraction based on the cross-correlation function (CCF) method (Pepe et al. 2002), where the observed spectrum is correlated with a binary mask matching the spectral type of each object. With a dedicated pipeline implemented at the INAF Trieste Observatory<fnref n="3" /> through the YABI workflow interface (Hunter et al. 2012), for the GAPS component of the sample, we also evaluated the CA index from the H&K lines of Ca II (S-index and log($R'_{HK}$)), through the procedure described in Sect. 4.

## 2.2. Objects from previous similar works

We gathered other transiting planet hosts for which S-indices and/or log($R'_{HK}$) values are available in the literature (see Table 1).

We did not take into consideration the set from Canto Martins et al. (2011) because the data are not homogeneous with the previous ones. Furthermore, those objects are not transiting their stars and the measured planetary masses are just minimum masses.

In summary, we selected from Knutson et al. (2010) 18 objects (with the measured S-index value) that are not included in the list of GAPS targets and another 40 objects (with no S-index) from Figueira et al. (2014) that are not included in either GAPS or Knutson's lists. Furthermore, we take into consideration the other two objects (WASP-72 and WASP-103) from Staab et al. (2017) for which the S-index values are given. Finally, the sample of 72 GAPS objects includes 43 stars already listed in previous works and 29 new entries. The complete target list considered in this study consists of 132 objects.

For homogeneity, we would like to evaluate the log($R'_{HK}$) from the original value of the S-index, using the stellar parameters that we obtained for all systems by the Gaia DR3 (see Sect. 3). Generally, the discovery papers reported in a few cases the complete time series of the $S_{HK}$. This was the case

<footnote n="2">² Since mid-2016, the Fabry Perot lamp has been operative, replacing the ThAr lamp due to a sharp decrease in the lamp flux.</footnote>

<footnote n="3">³ http://ia2.inaf.it</footnote>
<table><caption>**Table 1.** Compilations of $log(R'_{HK})$ measurements for planetary host stars.</caption><thead><tr><th>Compilation by</th><th>Obj.</th><th>log(R'_{HK})<sub>MAX</sub></th><th>Host</th><th>log(R'_{HK})<sub>MIN</sub></th><th>Host</th><th>Note</th></tr></thead><tbody><tr><td>Knutson et al. (2010)</td><td>50</td><td>-4.331</td><td>CoRoT-2</td><td>-5.500</td><td>WASP-12</td><td>a</td></tr><tr><td>Hartman (2010)</td><td>39</td><td>-4.331</td><td>CoRoT-2</td><td>-5.500</td><td>WASP-12</td><td>b</td></tr><tr><td>Figueira et al. (2014)</td><td>95</td><td>-4.100</td><td>WASP-59</td><td>-5.500</td><td>WASP-12</td><td>c</td></tr><tr><td>Fossati et al. (2015)</td><td>41</td><td>-4.100</td><td>WASP-59</td><td>-5.319</td><td>HAT-P-13</td><td>d</td></tr><tr><td>Canto Martins et al. (2011)</td><td>74</td><td>-4.42</td><td>HD128311</td><td>-5.22</td><td>HD4203</td><td>e</td></tr></tbody></table>

**Notes.** Columns indicate the compilation origin, number of objects, maximum and the minimum values, and the corresponding objects, from left to right. The last column gives explicative notes on each compilation. a: S-index available; b: subset of Knutson et al. (2010) compilation. No error is given; c: set of Hartman (2010) plus other multi planetary systems; d: same list of Figueira et al. (2014) but trimmed for systems with a ≤0.1 au and $M_p \ge 0.1 M_J$; e: no transiting systems.

of the following systems HAT-P-27, HAT-P-28, HAT-P-33, HAT-P-34, HAT-P-35, HAT-P-39, HAT-P-40, HAT-P-41, HAT-P-44, HAT-P-45, HAT-P-46, and HD97658 – for which we evaluate the mean value of $S_{HK}$. Instead, for the following systems, 55 Cnc, HAT-P-25, HAT-P-38, Kepler-19, Kepler-62, Kepler-68, WASP-22, and WASP-23, the bibliography reported only a single value of $S_{HK}$. For the remaining systems, when possible, we evaluated the $S_{HK}$ value from the available values of $log(R'_{HK})$ and the $B-V$ used by the paper authors. This was the case of: Kepler-10, Kepler-17, Kepler-78, WASP-5, WASP-15, WASP-58, WASP-70 A, and WASP-84. Not all the papers outline the value of $B-V$ used in evaluating the $log(R'_{HK})$, but most of them report the $T_{eff}$. We used the latter parameter to obtain the value of $B-V$ and then the value of $S_{HK}$ from Modern Mean Dwarf Stellar Color and Effective Temperature Sequence by Mamajek (2021, Pecaut & Mamajek 2013)$^4$. These systems are: Kepler-25, Kepler-48, Kepler-93, WASP-7, WASP-16, WASP-41, WASP-42, WASP-52, WASP-117, and Kepler-16. Due to the binarity of the latter system, the value of $S_{HK}$ obtained is pretty uncertain. Finally, for Kepler-20, we obtain the value of $B-V$ from Zacharias et al. (2012). All the hosts are listed in Table 2 with the bibliography used.

### 3. Stellar parameters

We used Gaia/DR3 (Gaia Collaboration 2016, 2018, 2023) to obtain the values of the G magnitude and Bailer-Jones et al. (2021) for the distance of all the targets in our sample$^5$. We found that only four stars of our sample (55 Cnc, HD 189733, HD 209458, and HD 97658) require the saturation correction described by Riello et al. (2021) for the G mag. Once applied, the correction for G is 0.002 mag for all four stars. The retrieved distances range from a minimum of 9.77 pc (GJ 436) to a maximum of 990.44 pc (Kepler-8) with a median of 221.12 pc (see the upper panel of Fig. 2). Thus, 85% of the sample is at a distance larger than 100 pc from the Sun.

We used the procedure exploited in developing the PLATO Input Catalog (PIC, Montalto et al. 2021) to retrieve the stellar parameters of the hosts homogeneously. We derived the $(B-V)_0$ colors from the intrinsic $G_{BP} - G_{RP}$ Gaia DR3 color by inverting the cubic equation reported in Table 2 of Evans et al. (2018). The reddening correction was obtained by interpolating the reddening map of Lallement et al. (2019, 2022), using the distances from Bailer-Jones et al. (2021). The effective temperatures of the stars were obtained from our custom-calibrated Gaia intrinsic color-effective temperature relationship used for the PIC.

The error in the color (median uncertainty of 3.6%) accounts for the reddening uncertainties as deduced from the reddening map, while the error on the temperature is derived from a Monte Carlo simulation, randomizing the input color within its uncertainties assuming normal error distribution and repeating the temperature calculation as described above. A normal random perturbation with a standard deviation equal to 200 K is added to the Monte Carlo temperature estimates to account for the expected difference between our temperature scale and literature results. All the details about this procedure are described in Montalto et al. (2021). In this way, we have a homogeneous determination of the $T_{eff}$ for most of the targets in the sample, with a trade-off of a median uncertainty of 3.5%. This is about four times the relative error for the determination of the $T_{eff}$ listed in the SWEET Catalog$^6$ (Santos et al. 2013; Sousa et al. 2018, 2021), but as these authors warn, these data are not homogeneous. In the catalog, only those systems with the SWFlag=1 can identify the stars with derived spectroscopic parameters in a homogeneous way by Sousa et al. (2021). Furthermore, we transformed the Gaia G magnitude to Johnson V magnitude using the relation by Evans et al. (2018). The middle and bottom panels of Fig. 2 show the distribution of the apparent V magnitude and of the $T_{eff}$ of the whole sample, respectively. The stars have apparent magnitude ranging $5.9 \le V \le 14.3$ mag, with a median of 11.6 mag, while the temperature range is $3461 \le T_{eff} \le 6781$ K. We evaluated the spectral type of each member of the sample by exploiting the table from Pecaut & Mamajek (2013). The sample comprises 49 F stars, 42 G stars, 39 K stars, and 2 M stars.

For a small subset of our sample (27 stars), we considered the metallicity homogeneously derived through the methodology described in Biazzo et al. (2022) and based on iron lines. For the other hosting stars, since no homogeneous metallicities were obtained in the literature, we considered the values listed in SWEET-Cat and derived through similar spectroscopic methods identified by SWFlag=1 (see, e.g., Sousa et al. 2021). Figure 3 shows the distribution of [Fe/H] spanning between -0.30 and 0.48 dex with a median value of 0.010 dex and a median error of 0.04 dex.

We retrieved the values of the mass, the radius, log g, and age of the stars through the PARAM web interface$^7$, which makes use of a Bayesian method, described in da Silva et al. (2006) applied on the theoretical isochrone set PARSEC (Bressan et al. 2012). The input parameters to the PARAM tool were the apparent magnitude V, the parallax, the $T_{eff}$, and [Fe/H].

<footnote>$^4$ https://www.pas.rochester.edu/~emamajek/</footnote>

<footnote>$^5$ Using SIMBAD (Wenger et al. 2000) and VIZIER (Ochsenbein et al. 2000).</footnote>

<footnote>$^6$ http://sweetcat.iastro.pt/</footnote>

<footnote>$^7$ http://stev.oapd.inaf.it/cgi-bin/param_1.3</footnote>
<table><caption>**Table 2.** For all systems, the table lists *Gaia* DR3 id., and (*B* − *V*), and *T*<sub>eff</sub> obtained by *Gaia* Collaboration (2018) as explained in the text.</caption><thead><tr><th>System</th><th><i>Gaia</i> DR3 Id.</th><th>(<i>B</i> − <i>V</i>)</th><th><i>T</i><sub>eff</sub><br>(K)</th><th><i>S</i><sub>MWO</sub></th><th>log(<i>R'</i><sub>HK</sub>)<sub>N</sub></th><th>log(<i>R'</i><sub>HK</sub>)<sub>M</sub></th><th><i>R'</i><sub>HK</sub> × 10<sup>5</sup><br>rms</th></tr></thead><tbody><tr><td>55Cnc</td><td>704967037090946688</td><td>0.841 ± 0.029</td><td>5214 ± 217</td><td>0.165<sup>(1)</sup></td><td>−5.0308</td><td>−4.8557</td><td>—</td></tr><tr><td>CoRoT-1</td><td>3105507886130792448</td><td>0.554 ± 0.021</td><td>6066 ± 205</td><td>0.124<sup>(2)</sup></td><td>−5.2761</td><td>−5.2975</td><td>—</td></tr><tr><td>CoRoT-2</td><td>4287820848378092672</td><td>0.821 ± 0.028</td><td>5265 ± 192</td><td>0.435<sup>(2)</sup></td><td>−4.4536</td><td>−4.2139</td><td>—</td></tr><tr><td>CoRoT-7</td><td>310726717757848576</td><td>0.824 ± 0.028</td><td>5256 ± 187</td><td>0.225<sup>(2)</sup></td><td>−4.8207</td><td>−4.6104</td><td>—</td></tr><tr><td>GJ3470<sup>(a)</sup></td><td>654687847820642560</td><td>1.917 ± 0.014</td><td>3558 ± 204</td><td>1.755 ± 0.321</td><td>−5.0625 ± 0.0843</td><td>−4.6696 ± 0.0852</td><td>0.158</td></tr><tr><td>GJ436<sup>(a)</sup></td><td>4017860992519744384</td><td>2.072 ± 0.014</td><td>3461 ± 208</td><td>0.735 ± 0.060</td><td>−5.3263 ± 0.0347</td><td>−5.3727 ± 0.0455</td><td>0.038</td></tr><tr><td>HAT-P-1</td><td>1928431764627661440</td><td>0.612 ± 0.023</td><td>5869 ± 199</td><td>0.151 ± 0.014</td><td>−5.0408 ± 0.1151</td><td>−4.8964 ± 0.1404</td><td>0.220</td></tr><tr><td>HAT-P-2<sup>(a,b)</sup></td><td>1380825667768742144</td><td>0.447 ± 0.019</td><td>6459 ± 202</td><td>0.181 ± 0.006</td><td>−4.7539 ± 0.0271</td><td>−4.6819 ± 0.0388</td><td>0.111</td></tr><tr><td>HAT-P-3</td><td>1510191594552968832</td><td>0.903 ± 0.031</td><td>5062 ± 185</td><td>0.220 ± 0.036</td><td>−4.8967 ± 0.1287</td><td>−4.6748 ± 0.1901</td><td>0.266</td></tr><tr><td>HAT-P-4</td><td>1291120362349158016</td><td>0.609 ± 0.023</td><td>5880 ± 203</td><td>0.147 ± 0.016</td><td>−5.0702 ± 0.1058</td><td>−4.8645 ± 0.1062</td><td>0.242</td></tr><tr><td>HAT-P-5</td><td>4606030169272920320</td><td>0.622 ± 0.023</td><td>5838 ± 207</td><td>0.148<sup>(2)</sup></td><td>−5.0683</td><td>−4.9209</td><td>—</td></tr><tr><td>HAT-P-6<sup>(a)</sup></td><td>1925321658551399040</td><td>0.378 ± 0.017</td><td>6570 ± 80<sup>(3)</sup></td><td>0.181 ± 0.032</td><td>−4.6872 ± 0.1006</td><td>−4.6464 ± 0.1525</td><td>0.648</td></tr><tr><td>HAT-P-7</td><td>2129256395211984000</td><td>0.469 ± 0.019</td><td>6376 ± 204</td><td>0.140 ± 0.004</td><td>−5.0318 ± 0.0390</td><td>−5.0585 ± 0.0744</td><td>0.078</td></tr><tr><td>HAT-P-8</td><td>1891507552826485632</td><td>0.521 ± 0.020</td><td>6181 ± 171</td><td>0.148 ± 0.004</td><td>−5.0072 ± 0.0276</td><td>−4.9462 ± 0.0395</td><td>0.062</td></tr><tr><td>HAT-P-9</td><td>898130030131443584</td><td>0.510 ± 0.020</td><td>6220 ± 204</td><td>0.141<sup>(2)</sup></td><td>−5.0567</td><td>−5.0381</td><td>—</td></tr><tr><td>HAT-P-11</td><td>2086512227851023872</td><td>1.075 ± 0.036</td><td>4681 ± 229</td><td>0.544 ± 0.010</td><td>−4.6808 ± 0.0081</td><td>−4.3017 ± 0.009</td><td>0.039</td></tr><tr><td>HAT-P-12</td><td>1499514786891168640</td><td>1.110 ± 0.037</td><td>4612 ± 213</td><td>0.362 ± 0.063</td><td>−4.9148 ± 0.0789</td><td>−4.5383 ± 0.0946</td><td>0.218</td></tr><tr><td>HAT-P-13</td><td>1014520826353577088</td><td>0.699 ± 0.025</td><td>5602 ± 193</td><td>0.141<sup>(2)</sup></td><td>−5.1400</td><td>−4.8945</td><td>—</td></tr><tr><td>HAT-P-14</td><td>1340497608486742400</td><td>0.452 ± 0.019</td><td>6439 ± 185</td><td>0.164 ± 0.003</td><td>−4.8479 ± 0.0184</td><td>−4.8152 ± 0.0295</td><td>0.058</td></tr><tr><td>HAT-P-15</td><td>179498266829041664</td><td>0.910 ± 0.031</td><td>5043 ± 204</td><td>0.182 ± 0.019</td><td>−5.0100 ± 0.0616</td><td>−4.7565 ± 0.0708</td><td>0.134</td></tr><tr><td>HAT-P-16</td><td>381592313648387200</td><td>0.548 ± 0.021</td><td>6085 ± 190</td><td>0.164 ± 0.060</td><td>−4.9157 ± 0.1950</td><td>−4.7965 ± 0.2226</td><td>1.036</td></tr><tr><td>HAT-P-17</td><td>1849786481032616960</td><td>0.824 ± 0.028</td><td>5257 ± 191</td><td>0.156 ± 0.042</td><td>−5.0629 ± 0.1392</td><td>−4.8963 ± 0.2363</td><td>0.391</td></tr><tr><td>HAT-P-18</td><td>1334573817793362560</td><td>1.045 ± 0.035</td><td>4744 ± 216</td><td>0.334 ± 0.079</td><td>−4.8553 ± 0.1295</td><td>−4.5282 ± 0.164</td><td>0.349</td></tr><tr><td>HAT-P-20</td><td>869913435026514688</td><td>1.203 ± 0.040</td><td>4436 ± 213</td><td>1.142 ± 0.107</td><td>−4.5609 ± 0.0389</td><td>−4.0656 ± 0.0409</td><td>0.259</td></tr><tr><td>HAT-P-21</td><td>770622651659107712</td><td>0.687 ± 0.025</td><td>5636 ± 200</td><td>0.306 ± 0.044</td><td>−4.5313 ± 0.0713</td><td>−4.3084 ± 0.0703</td><td>0.588</td></tr><tr><td>HAT-P-22</td><td>846946629987527168</td><td>0.812 ± 0.028</td><td>5290 ± 183</td><td>0.157 ± 0.016</td><td>−5.0590 ± 0.0702</td><td>−4.819 ± 0.0727</td><td>0.161</td></tr><tr><td>HAT-P-24</td><td>3167323052618369408</td><td>0.470 ± 0.019</td><td>6370 ± 197</td><td>0.146 ± 0.011</td><td>−4.9793 ± 0.1005</td><td>−4.9962 ± 0.2171</td><td>0.213</td></tr><tr><td>HAT-P-25</td><td>111322601672419712</td><td>0.854 ± 0.029</td><td>5182 ± 179</td><td>0.166<sup>(4)</sup></td><td>−5.0322</td><td>−4.8562</td><td>—</td></tr><tr><td>HAT-P-26</td><td>3668036348641580288</td><td>0.901 ± 0.031</td><td>5066 ± 246</td><td>0.172 ± 0.079</td><td>−5.0372 ± 0.2043</td><td>−4.8528 ± 0.2747</td><td>0.581</td></tr><tr><td>HAT-P-27<sup>(c)</sup></td><td>1159336403336463872</td><td>0.865 ± 0.030</td><td>5153 ± 202</td><td>0.231 ± 0.017<sup>(5)</sup></td><td>−4.8368</td><td>−4.6176</td><td>—</td></tr><tr><td>HAT-P-28</td><td>363702817083391232</td><td>0.821 ± 0.028</td><td>5265 ± 224</td><td>0.163 ± 0.013<sup>(6)</sup></td><td>−5.0318 ± 0.0690</td><td>−4.7908 ± 0.172</td><td>—</td></tr><tr><td>HAT-P-29</td><td>359058441314838528</td><td>0.595 ± 0.022</td><td>5925 ± 227</td><td>0.160 ± 0.020</td><td>−4.9705 ± 0.1778</td><td>−4.8274 ± 0.147</td><td>0.316</td></tr><tr><td>HAT-P-30</td><td>3096441729861716224</td><td>0.518 ± 0.020</td><td>6191 ± 195</td><td>0.142 ± 0.012</td><td>−5.0585 ± 0.2019</td><td>−5.0176 ± 0.0975</td><td>0.211</td></tr><tr><td>HAT-P-31</td><td>4583004815239650176</td><td>0.608 ± 0.022</td><td>5884 ± 231</td><td>0.130 ± 0.011</td><td>−5.2268 ± 0.1293</td><td>−5.0236 ± 0.1308</td><td>0.174</td></tr><tr><td>HAT-P-32</td><td>356348286886230272</td><td>0.469 ± 0.019</td><td>6376 ± 219</td><td>0.228 ± 0.010</td><td>−4.5879 ± 0.0327</td><td>−4.4472 ± 0.0397</td><td>0.194</td></tr><tr><td>HAT-P-33<sup>(a)</sup></td><td>893550942158776832</td><td>0.429 ± 0.018</td><td>6533 ± 226</td><td>0.172 ± 0.001<sup>(7)</sup></td><td>−4.7806 ± 0.0230</td><td>−4.7494 ± 0.084</td><td>—</td></tr><tr><td>HAT-P-34<sup>(a)</sup></td><td>1810218734058821632</td><td>0.423 ± 0.018</td><td>6555 ± 210</td><td>0.174 ± 0.002<sup>(8)</sup></td><td>−4.7656 ± 0.0290</td><td>−4.7349 ± 0.088</td><td>—</td></tr><tr><td>HAT-P-35</td><td>3094733088793018240</td><td>0.554 ± 0.021</td><td>6064 ± 182</td><td>0.129 ± 0.002<sup>(8)</sup></td><td>−5.2118 ± 0.0370</td><td>−5.1918 ± 0.087</td><td>—</td></tr><tr><td>HAT-P-36</td><td>1541532207133249920</td><td>0.706 ± 0.025</td><td>5582 ± 189</td><td>0.263 ± 0.028</td><td>−4.6378 ± 0.0671</td><td>−4.4364 ± 0.0699</td><td>0.364</td></tr><tr><td>HAT-P-38</td><td>324889227693298560</td><td>0.842 ± 0.029</td><td>5211 ± 141</td><td>0.144<sup>(9)</sup></td><td>−5.1290</td><td>−4.9821</td><td>—</td></tr><tr><td>HAT-P-39</td><td>3169391577586325120</td><td>0.470 ± 0.019</td><td>6371 ± 227</td><td>0.179 ± 0.002<sup>(10)</sup></td><td>−4.7795 ± 0.0250</td><td>−4.6761 ± 0.089</td><td>—</td></tr><tr><td>HAT-P-40</td><td>1962153854973972096</td><td>0.600 ± 0.022</td><td>5910 ± 188</td><td>0.145 ± 0.001<sup>(10)</sup></td><td>−5.0794 ± 0.0200</td><td>−4.8822 ± 0.076</td><td>—</td></tr><tr><td>HAT-P-41</td><td>4290415081653653632</td><td>0.428 ± 0.018</td><td>6534 ± 242</td><td>0.148 ± 0.002<sup>(10)</sup></td><td>−4.9237 ± 0.0360</td><td>−5.0068 ± 0.096</td><td>—</td></tr><tr><td>HAT-P-44</td><td>1506159170017461888</td><td>0.824 ± 0.028</td><td>5258 ± 212</td><td>0.125 ± 0.011<sup>(11)</sup></td><td>−5.2453 ± 0.0730</td><td>−5.1386 ± 0.213</td><td>—</td></tr><tr><td>HAT-P-45</td><td>4270547520706275328</td><td>0.483 ± 0.019</td><td>6322 ± 200</td><td>0.121 ± 0.008<sup>(11)</sup></td><td>−5.2540</td><td>−5.1684</td><td>—</td></tr><tr><td>HAT-P-46</td><td>4177560894633502720</td><td>0.651 ± 0.065</td><td>5826 ± 322</td><td>0.127 ± 0.004</td><td>−5.2756 ± 0.0490</td><td>−5.173 ± 0.182</td><td>—</td></tr><tr><td>HAT-P-47<sup>(a)</sup></td><td>131624808999798656</td><td>0.449 ± 0.019</td><td>6451 ± 215</td><td>0.177 ± 0.013</td><td>−4.7751 ± 0.1040</td><td>−4.71 ± 0.0564</td><td>0.247</td></tr><tr><td>HD 149026</td><td>1331356474971716992</td><td>0.566 ± 0.021</td><td>6024 ± 190</td><td>0.152<sup>(2)</sup></td><td>−5.0048</td><td>−4.8367</td><td>—</td></tr><tr><td>HD 17156</td><td>545560867790611072</td><td>0.584 ± 0.022</td><td>5961 ± 201</td><td>0.146 ± 0.022</td><td>−5.0654 ± 0.1362</td><td>−4.883 ± 0.1447</td><td>0.353</td></tr><tr><td>HD 189733</td><td>1827242816201846144</td><td>0.921 ± 0.031</td><td>5017 ± 186</td><td>0.521 ± 0.037</td><td>−4.4815 ± 0.0480</td><td>−4.1955 ± 0.0545</td><td>0.257</td></tr><tr><td>HD 209458</td><td>1779546757669063552</td><td>0.553 ± 0.021</td><td>6066 ± 180</td><td>0.156 ± 0.003</td><td>−4.9685 ± 0.0210</td><td>−4.8578 ± 0.0261</td><td>0.054</td></tr><tr><td>HD 80606</td><td>1019003226022657920</td><td>0.737 ± 0.026</td><td>5492 ± 203</td><td>0.149 ± 0.012</td><td>−5.0887 ± 0.1070</td><td>−4.9221 ± 0.1321</td><td>0.148</td></tr><tr><td>HD 97658</td><td>3997075206232885888</td><td>0.848 ± 0.029</td><td>5197 ± 231</td><td>0.175 ± 0.001<sup>(12)</sup></td><td>−4.9943 ± 0.0190</td><td>−4.8099 ± 0.122</td><td>—</td></tr><tr><td>KELT-3<sup>(a)</sup></td><td>806492023789218688</td><td>0.464 ± 0.019</td><td>6393 ± 198</td><td>0.189 ± 0.009</td><td>−4.7287 ± 0.0361</td><td>−4.6289 ± 0.0471</td><td>0.170</td></tr><tr><td>KELT-6</td><td>1464700950221781504</td><td>0.494 ± 0.020</td><td>6280 ± 215</td><td>0.151 ± 0.003</td><td>−4.9580 ± 0.0217</td><td>−4.9219 ± 0.0323</td><td>0.055</td></tr></tbody></table>

**Notes.** The Ca II H&K line strength is listed for both the original GAPS sample and the sample obtained by literature. We report in the log(<i>R'</i><sub>HK</sub>)<sub>N</sub> column the classical value (Noyes et al. 1984) of the activity as obtained in this work or in the literature while in the log(<i>R'</i><sub>HK</sub>)<sub>M</sub> the values corrected with the Mittag et al. (2013) procedure as described in Sect. 4.2. In the last column, the rms of the several measured time series are reported. (a) The R-indices obtained with the relation of Mittag et al. (2013) are derived in the color range of 0.44 < *B* − *V* < 1.6. Nevertheless, we give values for all the objects in our samples, and values for spectral types out of this color range should be taken with a pinch of salt. (b) Also known as HD147506. (c) Also known as WASP-40. (d) Also known as HAT-P-10.

**References.** (1) Wright et al. (2004); (2) Knutson et al. (2010); (3) Noyes et al. (2008); (4) Quinn et al. (2012); (5) Béky et al. (2011); (6) Buchhave et al. (2011); (7) Hartman et al. (2011); (8) Bakos et al. (2012); (9) Sato et al. (2012); (10) Hartman et al. (2012); (11) Hartman et al. (2014); (12) Howard et al. (2011); (13) Dumusque et al. (2014); (14) Winn et al. (2011); (15) Bonomo et al. (2012); (16) Malavolta et al. (2017); (17) Gautier et al. (2012); (18) Borucki et al. (2012); (19) Marcy et al. (2014); (20) Borucki et al. (2013); (21) Gilliland et al. (2013); (22) Pepe et al. (2013); (23) Triaud et al. (2010); (24) Albrecht et al. (2012); (25) Brown et al. (2012); (26) Anderson et al. (2011); (27) Triaud et al. (2011); (28) Maxted et al. (2011); (29) Lendl et al. (2012); (30) Hébrard et al. (2013); (31) Anderson et al. (2014); (32) Staab et al. (2017); (33) Lendl et al. (2014); (34) Tsantaki et al. (2014).
<table><caption>Table 2. continued.</caption><thead><tr><th>System</th><th>Gaia DR3 Id.</th><th>(B - V)</th><th>T<sub>eff</sub><br>(K)</th><th>S<sub>MWO</sub></th><th>log(R'<sub>HK</sub>)<sub>N</sub></th><th>log(R'<sub>HK</sub>)<sub>M</sub></th><th>R'<sub>HK</sub> × 10<sup>5</sup><br>rms</th></tr></thead><tbody><tr><td>Kepler-4</td><td>2132152916856093952</td><td>0.637 ± 0.023</td><td>5792 ± 204</td><td>0.168<sup>(2)</sup></td><td>-4.9437</td><td>-4.7219</td><td>-</td></tr><tr><td>Kepler-5</td><td>2079018300195390464</td><td>0.605 ± 0.022</td><td>5892 ± 240</td><td>0.148<sup>(2)</sup></td><td>-5.0584</td><td>-4.856</td><td>-</td></tr><tr><td>Kepler-6</td><td>2086636884980514304</td><td>0.715 ± 0.025</td><td>5556 ± 215</td><td>0.160<sup>(2)</sup></td><td>-5.0160</td><td>-4.7755</td><td>-</td></tr><tr><td>Kepler-7</td><td>210211781259036672</td><td>0.561 ± 0.021</td><td>6042 ± 202</td><td>0.155<sup>(2)</sup></td><td>-4.9801</td><td>-4.8148</td><td>-</td></tr><tr><td>Kepler-8</td><td>2116730994965905280</td><td>0.532 ± 0.021</td><td>6141 ± 225</td><td>0.153<sup>(2)</sup></td><td>-4.9764</td><td>-4.8897</td><td>-</td></tr><tr><td>Kepler-10</td><td>2132155017099178624</td><td>0.667 ± 0.024</td><td>5697 ± 189</td><td>0.164<sup>(13)</sup></td><td>-4.9782</td><td>-4.8045</td><td>-</td></tr><tr><td>Kepler-16</td><td>2133476355197071616</td><td>1.373 ± 0.015</td><td>4111 ± 186</td><td>0.913<sup>(14)</sup></td><td>-4.9293</td><td>-4.7535</td><td>-</td></tr><tr><td>Kepler-17</td><td>2086449761846310784</td><td>0.703 ± 0.025</td><td>5591 ± 201</td><td>0.420<sup>(15)</sup></td><td>-4.3611</td><td>-4.1517</td><td>-</td></tr><tr><td>Kepler-19</td><td>2051106987063242880</td><td>0.724 ± 0.026</td><td>5530 ± 216</td><td>0.164 ± 0.011<sup>(16)</sup></td><td>-4.9951 ± 0.0680</td><td>-4.8165 ± 0.169</td><td>-</td></tr><tr><td>Kepler-20</td><td>2102548708017562112</td><td>0.750 ± 0.026</td><td>5456 ± 215</td><td>0.183 ± 0.005<sup>(17)</sup></td><td>-4.9173 ± 0.0340</td><td>-4.7287 ± 0.131</td><td>-</td></tr><tr><td>Kepler-22</td><td>2127941757262806656</td><td>0.687 ± 0.025</td><td>5637 ± 188</td><td>0.149 ± 0.004<sup>(18)</sup></td><td>-5.0792 ± 0.0320</td><td>-4.913 ± 0.11</td><td>-</td></tr><tr><td>Kepler-25</td><td>2100451630105041152</td><td>0.502 ± 0.020</td><td>6251 ± 213</td><td>0.121<sup>(19)</sup></td><td>-5.2757</td><td>-5.1947</td><td>-</td></tr><tr><td>Kepler-48</td><td>2075112109039378688</td><td>0.894 ± 0.030</td><td>5082 ± 223</td><td>0.226<sup>(19)</sup></td><td>-4.8733</td><td>-4.6502</td><td>-</td></tr><tr><td>Kepler-62</td><td>2107001283431633408</td><td>0.953 ± 0.032</td><td>4945 ± 154</td><td>0.212 ± 0.002<sup>(20)</sup></td><td>-4.9622 ± 0.0340</td><td>-4.731 ± 0.104</td><td>-</td></tr><tr><td>Kepler-68</td><td>2129550445852902656</td><td>0.634 ± 0.023</td><td>5799 ± 196</td><td>0.139<sup>(21)</sup></td><td>-5.1500</td><td>-5.0069</td><td>-</td></tr><tr><td>Kepler-78</td><td>2078373642776670080</td><td>0.918 ± 0.031</td><td>5025 ± 178</td><td>0.467<sup>(22)</sup></td><td>-4.5297</td><td>-4.2495</td><td>-</td></tr><tr><td>Kepler-93</td><td>2052747119115620352</td><td>0.689 ± 0.025</td><td>5630 ± 230</td><td>0.165<sup>(19)</sup></td><td>-4.9801</td><td>-4.8031</td><td>-</td></tr><tr><td>Qatar-1</td><td>2244830490514284928</td><td>1.031 ± 0.035</td><td>4774 ± 211</td><td>0.597 ± 0.079</td><td>-4.5713 ± 0.0594</td><td>-4.2177 ± 0.065</td><td>0.368</td></tr><tr><td>Qatar-2</td><td>3620030644476623616</td><td>1.194 ± 0.040</td><td>4451 ± 186</td><td>1.122 ± 0.696</td><td>-4.5543 ± 0.2054</td><td>-4.0655 ± 0.2157</td><td>1.736</td></tr><tr><td>TrES-1</td><td>2098964849867337856</td><td>0.831 ± 0.029</td><td>5239 ± 235</td><td>0.229 ± 0.042</td><td>-4.8148 ± 0.1275</td><td>-4.6021 ± 0.1516</td><td>0.388</td></tr><tr><td>TrES-2</td><td>2131314401800665344</td><td>0.634 ± 0.023</td><td>5800 ± 198</td><td>0.156 ± 0.015</td><td>-5.0150 ± 0.1294</td><td>-4.8531 ± 0.1552</td><td>0.222</td></tr><tr><td>TrES-3</td><td>4609131509318715136</td><td>0.732 ± 0.026</td><td>5507 ± 214</td><td>0.305<sup>(2)</sup></td><td>-4.5663</td><td>-4.3571</td><td>-</td></tr><tr><td>TrES-4</td><td>4609062308806929152</td><td>0.507 ± 0.020</td><td>6232 ± 211</td><td>0.137 ± 0.008</td><td>-5.0888 ± 0.0757</td><td>-5.0365 ± 0.1156</td><td>0.134</td></tr><tr><td>WASP-1</td><td>2862548428079638912</td><td>0.570 ± 0.022</td><td>6011 ± 213</td><td>0.136 ± 0.071</td><td>-5.1429 ± 0.3168</td><td>-4.9851 ± 0.3438</td><td>1.186</td></tr><tr><td>WASP-2</td><td>1748596020745038208</td><td>0.892 ± 0.030</td><td>5087 ± 220</td><td>0.159<sup>(2)</sup></td><td>-5.0789</td><td>-4.9109</td><td>-</td></tr><tr><td>WASP-3</td><td>2095108312831835648</td><td>0.484 ± 0.019</td><td>6317 ± 218</td><td>0.147 ± 0.006</td><td>-4.9875 ± 0.0513</td><td>-4.9795 ± 0.0889</td><td>0.112</td></tr><tr><td>WASP-4</td><td>6535499658122055552</td><td>0.753 ± 0.026</td><td>5448 ± 210</td><td>0.194<sup>(2)</sup></td><td>-4.8750</td><td>-4.6819</td><td>-</td></tr><tr><td>WASP-5</td><td>6533307957785865856</td><td>0.679 ± 0.024</td><td>5661 ± 196</td><td>0.224<sup>(23)</sup></td><td>-4.7297</td><td>-4.5368</td><td>-</td></tr><tr><td>WASP-7</td><td>6681720724498802176</td><td>0.449 ± 0.019</td><td>6453 ± 196</td><td>0.447<sup>(24)</sup></td><td>-4.1655</td><td>-3.9722</td><td>-</td></tr><tr><td>WASP-10</td><td>1909762228985058944</td><td>1.083 ± 0.036</td><td>4666 ± 209</td><td>0.145 ± 0.074</td><td>-5.2953 ± 0.1494</td><td>-5.0891 ± 0.2144</td><td>0.285</td></tr><tr><td>WASP-11<sup>(d)</sup></td><td>123376685084303360</td><td>0.997 ± 0.033</td><td>4846 ± 208</td><td>0.295 ± 0.053</td><td>-4.8489 ± 0.0655</td><td>-4.56 ± 0.076</td><td>0.281</td></tr><tr><td>WASP-12</td><td>3435282862461427072</td><td>0.570 ± 0.022</td><td>6009 ± 218</td><td>0.112 ± 0.059</td><td>-5.5035 ± 0.4301</td><td>-5.3623 ± 0.5349</td><td>0.989</td></tr><tr><td>WASP-13</td><td>702295464353437952</td><td>0.591 ± 0.022</td><td>5940 ± 172</td><td>0.140 ± 0.008</td><td>-5.1172 ± 0.0905</td><td>-4.9306 ± 0.0966</td><td>0.127</td></tr><tr><td>WASP-14</td><td>1242084170974175232</td><td>0.447 ± 0.019</td><td>6458 ± 200</td><td>0.152 ± 0.004</td><td>-4.9130 ± 0.0293</td><td>-4.9334 ± 0.0555</td><td>0.077</td></tr><tr><td>WASP-15</td><td>6171009049549197824</td><td>0.491 ± 0.020</td><td>6289 ± 198</td><td>0.165 ± 0.015<sup>(23)</sup></td><td>-4.8704 ± 0.1040</td><td>-4.7925 ± 0.197</td><td>-</td></tr><tr><td>WASP-16</td><td>6283723285046532864</td><td>0.704 ± 0.025</td><td>5587 ± 178</td><td>0.145<sup>(25)</sup></td><td>-5.1102</td><td>-4.9466</td><td>-</td></tr><tr><td>WASP-17<sup>(a)</sup></td><td>6042793005779654656</td><td>0.345 ± 0.016</td><td>6561 ± 287</td><td>0.121<sup>(2)</sup></td><td>-5.0183</td><td>-4.8822</td><td>-</td></tr><tr><td>WASP-18</td><td>4955371367334610048</td><td>0.479 ± 0.019</td><td>6337 ± 238</td><td>0.116<sup>(2)</sup></td><td>-5.3271</td><td>-5.2571</td><td>-</td></tr><tr><td>WASP-19</td><td>5411736896952029568</td><td>0.775 ± 0.027</td><td>5386 ± 188</td><td>0.252<sup>(2)</sup></td><td>-4.7152</td><td>-4.5051</td><td>-</td></tr><tr><td>WASP-21</td><td>2831084391023184128</td><td>0.637 ± 0.023</td><td>5789 ± 214</td><td>0.160 ± 0.007</td><td>-4.9904 ± 0.0432</td><td>-4.8248 ± 0.048</td><td>0.099</td></tr><tr><td>WASP-22</td><td>5086537022856406272</td><td>0.581 ± 0.022</td><td>5972 ± 165</td><td>0.170 ± 0.010<sup>(26)</sup></td><td>-4.8985</td><td>-4.7533</td><td>-</td></tr><tr><td>WASP-23</td><td>5557345496687437696</td><td>0.895 ± 0.030</td><td>5080 ± 212</td><td>0.320 ± 0.040<sup>(27)</sup></td><td>-4.6895 ± 0.0960</td><td>-4.4376 ± 0.191</td><td>-</td></tr><tr><td>WASP-24</td><td>1153682508388170112</td><td>0.542 ± 0.021</td><td>6105 ± 179</td><td>0.159 ± 0.015</td><td>-4.9446 ± 0.1074</td><td>-4.8376 ± 0.143</td><td>0.260</td></tr><tr><td>WASP-26</td><td>2416782701664155008</td><td>0.604 ± 0.022</td><td>5897 ± 192</td><td>0.147 ± 0.009</td><td>-5.0687 ± 0.0777</td><td>-4.9361 ± 0.095</td><td>0.139</td></tr><tr><td>WASP-31</td><td>3545873945303316224</td><td>0.518 ± 0.020</td><td>6192 ± 238</td><td>0.150 ± 0.023</td><td>-4.9873 ± 0.1462</td><td>-4.9222 ± 0.2079</td><td>0.411</td></tr><tr><td>WASP-32</td><td>2546413408888429696</td><td>0.563 ± 0.021</td><td>6034 ± 164</td><td>0.175 ± 0.023</td><td>-4.8599 ± 0.1110</td><td>-4.7194 ± 0.1305</td><td>0.385</td></tr><tr><td>WASP-35</td><td>3211188618762023424</td><td>0.608 ± 0.022</td><td>5884 ± 180</td><td>0.157 ± 0.036</td><td>-4.9917 ± 0.1411</td><td>-4.8433 ± 0.1671</td><td>0.571</td></tr><tr><td>WASP-38</td><td>4453211899986180352</td><td>0.527 ± 0.020</td><td>6161 ± 200</td><td>0.152 ± 0.003</td><td>-4.9801 ± 0.0232</td><td>-4.901 ± 0.0325</td><td>0.052</td></tr><tr><td>WASP-39</td><td>3643098875168270592</td><td>0.781 ± 0.027</td><td>5371 ± 240</td><td>0.172 ± 0.035</td><td>-4.9778 ± 0.1726</td><td>-4.795 ± 0.2549</td><td>0.372</td></tr><tr><td>WASP-41</td><td>6160877943251930368</td><td>0.764 ± 0.027</td><td>5416 ± 191</td><td>0.261<sup>(28)</sup></td><td>-4.6853</td><td>-4.4754</td><td>-</td></tr><tr><td>WASP-42</td><td>6139698733664726784</td><td>0.955 ± 0.032</td><td>4940 ± 185</td><td>0.202<sup>(29)</sup></td><td>-4.9899</td><td>-4.7652</td><td>-</td></tr><tr><td>WASP-43</td><td>3767805209112436736</td><td>1.360 ± 0.015</td><td>4162 ± 204</td><td>1.618 ± 0.125</td><td>-4.6610 ± 0.0341</td><td>-4.0522 ± 0.0349</td><td>0.849</td></tr><tr><td>WASP-48</td><td>2141754578242371584</td><td>0.520 ± 0.020</td><td>6184 ± 267</td><td>0.169 ± 0.011</td><td>-4.8648 ± 0.0601</td><td>-4.724 ± 0.0713</td><td>0.192</td></tr><tr><td>WASP-49</td><td>2991284369063612928</td><td>0.728 ± 0.026</td><td>5517 ± 188</td><td>0.158 ± 0.012</td><td>-5.0330 ± 0.0633</td><td>-4.8587 ± 0.071</td><td>0.142</td></tr><tr><td>WASP-50</td><td>5160557726183065984</td><td>0.771 ± 0.027</td><td>5399 ± 180</td><td>0.161 ± 0.002</td><td>-5.0260 ± 0.0108</td><td>-4.8508 ± 0.0125</td><td>0.842</td></tr><tr><td>WASP-52</td><td>26660158787557546496</td><td>0.914 ± 0.031</td><td>5034 ± 191</td><td>0.585<sup>(30)</sup></td><td>-4.4180</td><td>-4.1318</td><td>-</td></tr><tr><td>WASP-54</td><td>3661983850663908608</td><td>0.541 ± 0.021</td><td>6111 ± 192</td><td>0.141 ± 0.010</td><td>-5.0825 ± 0.1189</td><td>-4.9593 ± 0.1665</td><td>0.173</td></tr><tr><td>WASP-57</td><td>6340377202215639936</td><td>0.668 ± 0.024</td><td>5695 ± 216</td><td>0.213 ± 0.046</td><td>-4.7592 ± 0.1804</td><td>-4.5697 ± 0.1939</td><td>0.639</td></tr><tr><td>WASP-58</td><td>2115245554756763392</td><td>0.593 ± 0.022</td><td>5932 ± 207</td><td>0.347<sup>(30)</sup></td><td>-4.3887</td><td>-4.1979</td><td>-</td></tr><tr><td>WASP-59</td><td>2841043011433969152</td><td>1.280 ± 0.014</td><td>4270 ± 193</td><td>0.707 ± 0.103</td><td>-4.8952 ± 0.0660</td><td>-4.3494 ± 0.0719</td><td>0.185</td></tr><tr><td>WASP-60</td><td>2868528637464028160</td><td>0.641 ± 0.023</td><td>5779 ± 209</td><td>0.147 ± 0.009</td><td>-5.0877 ± 0.0704</td><td>-4.8597 ± 0.0675</td><td>0.134</td></tr><tr><td>WASP-69</td><td>6910753016653587840</td><td>1.033 ± 0.035</td><td>4769 ± 228</td><td>0.671 ± 0.013</td><td>-4.5221 ± 0.0082</td><td>-4.1625 ± 0.0089</td><td>0.057</td></tr><tr><td>WASP-70A</td><td>6888469489370835840</td><td>0.650 ± 0.024</td><td>5750 ± 220</td><td>0.132<sup>(31)</sup></td><td>-5.2200</td><td>-5.0829</td><td>-</td></tr><tr><td>WASP-72</td><td>5065640460769428224</td><td>0.545 ± 0.021</td><td>6097 ± 164</td><td>0.124 ± 0.009<sup>(32)</sup></td><td>-5.2701 ± 0.1430</td><td>-5.1909 ± 0.205</td><td>-</td></tr><tr><td>WASP-76</td><td>2512326349403275520</td><td>0.532 ± 0.021</td><td>6141 ± 211</td><td>0.140 ± 0.009</td><td>-5.0802 ± 0.0696</td><td>-4.9708 ± 0.0855</td><td>0.158</td></tr><tr><td>WASP-80</td><td>4223507222112425344</td><td>1.505 ± 0.016</td><td>3971 ± 199</td><td>1.773 ± 0.087</td><td>-4.8242 ± 0.0210</td><td>-4.1707 ± 0.0213</td><td>0.073</td></tr><tr><td>WASP-84</td><td>3078836109158636928</td><td>0.843 ± 0.029</td><td>5210 ± 175</td><td>0.445<sup>(31)</sup></td><td>-4.4652</td><td>-4.2175</td><td>-</td></tr></tbody></table>
<table><caption>**Table 2.** continued.</caption><thead><tr><th>System</th><th>Gaia DR3 Id.</th><th>(B-V)</th><th>T<sub>eff</sub><br>(K)</th><th>S<sub>MWO</sub></th><th>log(R'<sub>HK</sub>)<sub>N</sub></th><th>log(R'<sub>HK</sub>)<sub>M</sub></th><th>R'<sub>HK</sub> × 10<sup>5</sup><br>rms</th></tr></thead><tbody><tr><td>WASP-103</td><td>4439085988769170432</td><td>0.601 ± 0.022</td><td>5908 ± 191</td><td>0.247 ± 0.030<sup>(32)</sup></td><td>-4.6106 ± 0.0990</td><td>-4.4066±</td><td>-</td></tr><tr><td>WASP-106</td><td>3788394461991295488</td><td>0.547 ± 0.021</td><td>6089 ± 214</td><td>0.182 ± 0.006</td><td>-4.8130 ± 0.0318</td><td>-4.6738 ± 0.0378</td><td>0.107</td></tr><tr><td>WASP-107</td><td>3578638842054261248</td><td>1.272 ± 0.016</td><td>4262 ± 194</td><td>1.191 ± 0.089</td><td>-4.6556 ± 0.0324</td><td>-4.0914 ± 0.0339</td><td>0.166</td></tr><tr><td>WASP-117</td><td>4746157737910069888</td><td>0.566 ± 0.021</td><td>6023 ± 182</td><td>0.160<sup>(33)</sup></td><td>-4.9497</td><td>-4.8242</td><td>-</td></tr><tr><td>WASP-127</td><td>3778075717162985600</td><td>0.638 ± 0.023</td><td>5788 ± 211</td><td>0.156 ± 0.004</td><td>-5.0157 ± 0.0280</td><td>-4.7914 ± 0.0271</td><td>0.064</td></tr><tr><td>XO-1</td><td>1316708918505350528</td><td>0.690 ± 0.025</td><td>5630 ± 236</td><td>0.182 ± 0.011</td><td>-4.8929 ± 0.0506</td><td>-4.7086 ± 0.0547</td><td>0.143</td></tr><tr><td>XO-2N</td><td>934346809278715776</td><td>0.827 ± 0.028</td><td>5248 ± 222</td><td>0.175 ± 0.031</td><td>-4.9833 ± 0.1272</td><td>-4.7398 ± 0.1346</td><td>0.291</td></tr><tr><td>XO-3<sup>(a)</sup></td><td>470650506779277952</td><td>0.325 ± 0.016</td><td>6781 ± 44<sup>(34)</sup></td><td>0.221 ± 0.036</td><td>-4.4877 ± 0.0802</td><td>-4.4124 ± 0.1164</td><td>0.789</td></tr><tr><td>XO-4</td><td>990291507088739072</td><td>0.488 ± 0.020</td><td>6300 ± 176</td><td>0.146 ± 0.055</td><td>-5.0001 ± 0.1775</td><td>-4.991 ± 0.2271</td><td>1.009</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Distributions of the stellar parameters of the GAPS sample. Upper panel: distribution of the distances of the star of the sample. Middle panel: distribution of the apparent V mag. Bottom panel: $T_{\text{eff}}$ distribution. The corresponding spectral types are indicated at the very top.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Distribution of the spectroscopic stellar parameter: metallicities of the hosts.</figcaption></figure>

In Fig. 4 (from top to bottom), the distribution of surface gravity, masses, and radii are reported. They range between $3.92 \le \log g \le 4.94$ with a median uncertainty of 1%, $0.27 \le M/M_{\odot} \le 1.51$, with a median relative error of 4%, and $0.26 \le R/R_{\odot} \le 2.12$ with an uncertainty of 5%, respectively.

The age of the sample is between $0.3 \le \text{age} \le 7.3$ Gyr, with 4.3 Gyr as the median age (see Fig. 5). The median

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Distribution of the stellar parameters: surface gravity (top panel), masses (middle panel), and radii (bottom panel).</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Distribution of the stellar ages. Top panel: all hosts. Bottom panel: age distribution, considering only those with $\Delta\text{age}/\text{age} \le 50\%$.</figcaption></figure>

relative uncertainty on age determination is 83%. Such a high uncertainty is generally given by two orders of problems tied both to the isochrones method. The age determination with the
<figure><img src="image_7.png" /><figcaption>Fig. 6. HIPPARCOS average MS by Wright (2004) and the transiting planet hosts. Light-green dashed lines indicate loci 0.45 and 2.0 mag above the MS.</figcaption></figure>

isochrones method introduces large uncertainties due to the poor knowledge of star luminosities caused by the errors in the values of the parallaxes; furthermore, a second effect is due to the long permanence of the star on the main sequence (MS). In our case, the parallax comes from the Gaia DR3 with a small uncertainty of 0.3%. Nevertheless, $T_{eff}$ is also an input parameter useful for evaluating the star's age. In this case (as already described), for most of the objects of our sample, we used the values of the $T_{eff}$ from the PIC, with a median error of 203 K. Such a high uncertainty on $T_{eff}$ has an impact on the choice of the right isochrones. Taking this into account, we selected all the ages with an uncertainty $\le 50\%$ (see bottom panel of Fig. 5). The peak of the distribution at 4.2 Gyr disappears, leaving a quite uniform distribution of ages.

To discern between evolved and unevolved stars in our sample, we evaluated the absolute visual magnitude for each host star taking into account their parallax. The parallax errors allow us to distinguish between M-S and slightly evolved objects. In Fig. 6, we compare the objects with the HIPPARCOS average M-S fitted by Wright (2004) in the (B-V) vs. $M_V$ plane. The two dashed curves indicate the portion of the plane above the M-S by 0.45 and 2.0 mag, where moderately evolved stars fall. Some objects (33) are in the sub-giant region. Objects with a height greater than 2.0 mag are considered evolved objects, but none of the stars fall in this region. It is worth noting that this is a selection effect that limits the number of evolved stars with respect to MS stellar hosts. In fact, the planet detection around a star and the determination of the planetary mass is made hard by the RV jitter induced by magnetic activity (e.g., Wright 2005) and stellar oscillations (e.g., Hatzes 1996; Beck et al. 2015). The latter ones are typically enhanced in amplitude and periods in evolved stars together with the increasing strength of surface convective motions. Both cause an increase in the RV jitter that may conceal the Keplerian signal of a planetary companion.

# 4. Chromospheric activity of the hosts

## 4.1. Extracting S-index and log $R'_{HK}$

To evaluate the S-index for each spectrum, we used the procedure described by Lovis et al. (2011). After the standard data reduction obtained by the DRS on-line at the telescope, this procedure performs background subtraction to correct for the diffuse light

due to the ThAr or Fabry Perot lamp when we observe in simultaneous reference mode. The correction for the background is essential for the spectrophotometric accuracy and precision in evaluating the S-index. The S-index is then evaluated in the same way and with the same spectral passbands as defined for Mt Wilson Observatory spectrophotometer (Duncan et al. 1991). In particular, the Mt Wilson index is defined as:

$$S = \alpha \frac{H + K}{R + V}, \qquad (1)$$

where H and K are the total flux in the CaII H&K line cores evaluated over triangular bands centered at 393.3664 nm (K) and 396.8470 nm (H), with an FWHM of 0.109 nm. Instead, V and R are the continuum bandpasses at 390.1070 nm and 400.1070 nm, respectively, with a width of 2.0 nm. The $\alpha$ value is a historical conversion factor discussed in depth by Hall et al. (2007), and fixed either at 2.40 or 2.30 (Duncan et al. 1991), or adjusted following the calibration by Baliunas et al. (1995). Lovis et al. (2011) found that this S-index is equal (or very close) to the ratio of the mean fluxes per wavelength interval ($\bar{X} = X/\Delta\lambda_X$, with X representing any bandpass, and $\Delta\lambda_X$ its width in wavelength). This reduces edge effects from finite sampling at the boundaries of the bands and it leads to a value of $\alpha \sim 1$ (Lovis et al. 2011). Together with the S value, its uncertainty is also computed considering the photon shot noise through error propagation. The value of the S-index obtained with HARPS (or HARPS-N) is indeed very close to the $S_{MWO}$ index, as can be seen by Eq. (3) in Lovis et al. (2011), with an observed difference between the two indices of 0.03.

The value of $R'_{HK}$ is the chromospheric excess flux of the CaII H&K lines, $F'_{HK}$, normalized by the bolometric flux of the star (Noyes et al. 1984):

$$R'_{HK} = \left( \frac{F'_{HK}}{\sigma T_{eff}^4} \right) = R_{HK} - R_{HK,phot}, \qquad (2)$$

where $R_{HK}$ is the stellar surface flux $F'_{HK}$ in the CaII H&K lines normalized by the bolometric stellar flux and $R_{HK,phot}$ is the photospheric flux contribution, $F'_{HK,phot}$, in the line center normalized by the stellar bolometric flux. The former is tied to the measured $S_{MWO}$ index via the following relation:

$$F_{HK} = \frac{F_{RV}}{\alpha} S_{MWO}, \qquad (3)$$

where $F_{RV}$ is the surface flux in both the continua, while $\alpha$ is the same factor of Eq. (1).

The procedure on YABI measures in an automatic way the value of the S-indices of the time series obtained by GAPS for each object. In the period between 2012 Sep. 28 (JD=2456199) and 2012 Oct. 25 (JD=2456228), the HARPS-N CCD had a failure in the red chip. We did not consider the spectra taken in that period. We exploited the S-index values obtained with the procedure on YABI and derived the $log(R'_{HK})_N$ values following Noyes et al. (1984). Some stars of the sample show high dispersion in the time series of $log(R'_{HK})_N$, due to the intrinsic activity cycle or flares. These variations can be occasionally observed during intense flares, but in most of the transient events, the ranges of the $log(R'_{HK})_N$ changes are comparable or smaller than those associated with the rotational modulation or the activity cycles. Nevertheless, the sparseness of the time series in a significant fraction of our targets prevents us from identifying the individual sources of variations as seen in the $log(R'_{HK})_N$ datasets
of the individual stars. To characterize the mean activity level of each of our targets, we decided to adopt the median of the observed log(R'_{HK}) timeseries of each star as the best indicator for that level. We computed also the arithmetic mean and compared it with the median to provide information on the deviation of the distribution of the log(R'_{HK}) measurements from a symmetric distribution that may suggest a preponderance of transient events. In general, these tend to increase the mean level of the log(R'_{HK}) with respect to its median. We found only three hosts (HAT-P-50, WASP-43, and WASP-50) with the average log(R'_{HK}) values differing more than 0.1 dex from the corresponding median. Moreover, the S/N of HAT-P-50 spectra is so low that we decided to exclude it from our analysis. The S-index and the corresponding log(R'_{HK})_N median values are reported for each object of the GAPS sample in Table 2. We discarded KELT-7, KELT-9, KELT-20, and WASP-33 from our sample because of their high temperatures and very blue colors that are very far from the bluer limit of the calibration range of log(R'_{HK}). Moreover, we evaluate the log(R'_{HK})_N also for all the hosts for which the S_{HK} values have been extracted by the literature (see Sect. 2.2, and Table 2).

## 4.2. The correction for the photospheric line flux

The log(R'_{HK})_N index obtained with the Lovis' procedure underestimates the value of the activity of the star, because the photospheric flux line contribution is not fully taken into consideration in the correction by Noyes et al. (1984). This choice results in an under-correction of the photospheric flux contribution by a factor of 2 (Mittag et al. 2013), thereby leading to an underestimation of the true flux excess in the Ca II H&K lines.

To convert the S-index into absolute CaII H&K chromospheric flux, Mittag et al. (2013) derived new color-dependent photospheric flux relation for the MS and evolved (subgiant and giant) stars in the color range between $0.44 \le (B-V) \le 1.6$. This will allow us to compare different compilations of the CA indicator from the literature.

Mittag et al. (2013) derived the following relation between the continuum flux and color index $B-V$ for the MS stars (with $0.44 \le B-V \le 1.60$) and subgiant stars (with $0.44 \le B-V \le 1.10$):

$$ \log \left( \frac{\mathcal{F}_{RV}}{\alpha} \right) = 8.25 - 1.67(B - V). \qquad (4) $$

With respect to the $R_{HK,phot}$ in Mittag et al. (2013), there are also the following relations with the color index for stars of different luminosity classes:

$$ \left\{ \begin{array}{ll} \text{MS stars} & \\ \log \mathcal{F}_{HK,phot} = 7.49 - 2.06(B-V) & 0.44 \le B-V < 1.28 \\ \log \mathcal{F}_{HK,phot} = 6.19 - 1.04(B-V) & 1.28 \le B-V < 1.60 \end{array} \right. \\ \text{Subgiant stars} \\ \log \mathcal{F}_{HK,phot} = 7.57 - 2.25(B-V) & 0.44 \le B-V < 1.10 \\ \text{Giant stars} \\ \log \mathcal{F}_{HK,phot} = 7.61 - 2.37(B-V) & 0.76 \le B-V < 1.18. \end{array} \right. \qquad (5) $$

This method allows for a simple and reliable conversion of the S-index and it will allow us to compare historical and new S-indices (Mittag et al. 2013).

Using the relations from Eqs. (2) to (5), we estimated the value of log(R'_{HK})_M for all the hosts (see Table 2). We use the suffix M to distinguish the value of log(R'_{HK}) obtained with the procedure described by Mittag et al. (2013) from that (log(R'_{HK})_N) obtained with the procedure described by Noyes et al. (1984).

We compared the values obtained with both methods to check if they are different in a significant way. We test the difference between the two mean log(R'_{HK}) values against the null hypothesis using the t-student statistic technique. It results that most of the means are significantly different with a confidence level of <0.001. The objects: HAT-P-6, HAT-P-24, HAT-P-30, WASP-3, WASP-12, WASP-31, and XO-4, have averages that cannot be considered different due to large errors that hamper the evaluation of the two log(R'_{HK}) values.

The obtained values of the CA allow us to find an empirical relation between log(R'_{HK})_N and log(R'_{HK})_M for the different evolutionary states of the hosts. The empirical relation will be useful to correct the classic index for the CA (log(R'_{HK})_N) for the photospheric line flux obtaining log(R'_{HK})_M when the S-indices are missing in the literature. We separated the sample into MS stars and slightly evolved (subgiant) stars as they come out from the comparison with the HIPPARCOS average MS (see Sect. 3 and Fig. 6) and for both the subsamples we search for a linear correlation between them. We take into consideration all the significant different values for all the objects in the MS with color ranging in the $0.44 \le B-V \le 1.60$ interval. For the subgiant stars instead, we considered those objects with color $0.44 \le B-V \le 1.10$. Figure 7 displays the comparison between the log(R'_{HK})_N and log(R'_{HK})_M for the MS objects (top panel) and sub-giant stars (bottom panel) and the correlation is described by the following:

$$ \left\{ \begin{array}{l} \text{Main sequence stars} \\ \log(R'_{HK})_M = (1.255 \pm 0.229) + (1.223 \pm 0.047) \times \log(R'_{HK})_N \\ \\ \text{Subgiant stars} \\ \log(R'_{HK})_M = (1.205 \pm 0.485) + (1.202 \pm 0.096) \times \log(R'_{HK})_N. \end{array} \right. \qquad (6) $$

In the following of this work, we will use the values of log(R'_{HK})_M, and consequently, we drop the subscript M. The distribution of the log(R'_{HK}) is shown in Fig. 8. Most of the transiting planet hosts are in the inactive and active region of the plot as they are defined by Wright (2004) and Henry et al. (1996), fewer are in the very active region, while some objects fall in the very inactive region with log(R'_{HK}) $\le -5.1$ (basal activity limit for MS stars). The upper and lower limits of our sample are respectively individuated by WASP-7 (log(R'_{HK}) = -3.9722) and GJ436 (log(R'_{HK}) = -5.3727). The median value of the sample is log(R'_{HK}) = -4.8099 held by HD 97658. This distribution follows the typical bias against very active stars of RV and (to a lesser degree) transit surveys.

# 5. R'_{HK} raw dispersion and chromospheric emission variability

For the targets observed by GAPS, we have long time series of the CA indicator, but their inhomogeneity does not allow a global analysis. Nevertheless, we used the recorded variation of CA to better characterize this sample. Information on the activity variability can be obtained via the standard deviation of the stellar activity indicator. The left panel of Fig. 9 shows the distribution
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Comparison between the values of $log(R'_{HK,Mittag})$ and $log(R'_{HK,Noyes})$ for the objects of the sample. The continuous blue line represents the best linear fit and the residuals of the fit are plotted beneath. The median error bars are indicated at the top of the plot. Top panel: comparison for the MS objects. Bottom panel: same as the top panel but for the sub-giant stars.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Distribution of the $log(R'_{HK})$ values of the considered sample of transiting objects. Plot also indicates the activity regions following Wright (2004) and Henry et al. (1996).</figcaption></figure>

of the standard deviation of $10^5 \times R'_{HK}$ (hereafter, $\sigma R_5$) as evaluated by the GAPS time series. The median of the distribution is 0.2 dex, and the distribution shows a tail of stars with high dispersion, possibly due to short-term variability (e.g., flares) or long-term magnetic cycles. The M2.5 dwarf GJ436 records the minimum value (0.038 dex) of the chromospheric emission (CE) variability. We observed this star on two distinct nights about 2 yr apart, so it could not be considered proof of the lowest dispersion level of our sample. Following other authors (e.g., Lovis et al. 2011; Gomes da Silva et al. 2021), we selected hosts observed on more than 20 nights on a time basis longer than 1000 days. The host with the lowest relative dispersion is HAT-P-7 (F6V) with an absolute $\sigma R_5 = 0.079$ dex and a relative dispersion for $S_{MW}$ of 2.0%. Comparing HAT-P-7 with the values for the quiet star $\tau$ Ceti: $S_{MW}/\sigma S_{MW} = 0.35\%$, obtained by Lovis et al. (2011) on a time bases of 7 yr of observations (157 nights), we can see that also HAT-P-7 could not be considered a long-term quiet star. So we have just a mere indication of the precision of our $R'_{HK}$ measurement.

Stellar activity compilations and catalogs measured for a high number of stars show that the higher envelope of activity dispersion decreases with decreasing activity level and evolutionary stage. The majority of stars have a variation in activity that ranges $0.01 \le \sigma R_5 \le 0.32$ dex (Gomes da Silva et al. 2021). In our sample, there are systems with a dispersion larger than 0.5 dex (see Table 3). Among them, Qatar-2 shows the highest amplitude of variation (1.736 dex) but because of the very short time span (less than 1 day), it could be due to a flare. All the others have longer laptime, instead.

In general, the activity variability is proportional to the level of stellar activity, and therefore, regarding age, young stars have greater activity variation than old ones. This behavior has been also highlighted by Scandariato et al. (2017) for a sample of low-activity early-type M dwarfs. The right panel of Fig. 9 reports the $\sigma R_5$ as a function of the $log(R'_{HK})$. The different colors represent the time span of the measurements. Considering only the red squares, corresponding to values with a time baseline longer than 100 days, it is possible to see that active stars have a greater spread of the CE variation. Nevertheless, some stars in the inactive zone ($-5.1 \le \log(R'_{HK}) \le -4.75$) exhibit a broad variation. They are WASP-1 (F9.5V), HAT-P-16 (F9V), XO-4 (F7V), HAT-P-26 (K2V), and WASP-35 (G1V). This situation seems to be corroborated by Gomes da Silva et al. (2021), who observed a high activity variation for stars of spectral types F and K with a low level of activity, as well as a larger spread of variation for the former spectral type than for the latter.

For the time being, we do not have the capacity to dig deeply into the activity variation behavior, and we leave it for future works; this is also because it is not the focal topic of this paper. In any case, the data indicate that the more active the stars, the more the activity varies. Furthermore, non-active stars could have high activity variation, in particular F stars.

**6. Activity of planetary hosts in a general context**

Figure 10 shows the distribution of the mean $log(R'_{HK})$ obtained for each host as a function of the $T_{eff}$ of the stars. In this plot, it is possible to note that F stars (blue circles) have a scattered distribution with a concentration around $log(R'_{HK}) = -5.0$. Differently, the G stars (green circles) range between $-5.1 \le \log(R'_{HK}) \le -4.1$, with a slightly high concentration around $\log(R'_{HK}) \sim -4.9$. Then, K stars (orange circles) show a scattered
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Standard deviation of $\sigma R_5$ for the GAPS transiting planet hosts: Left panel: distribution of $\sigma R_5$. Right panel: standard deviation of $R_5$ versus the log($R'_{HK}$) for the GAPS stars. The different colors indicate the time baselines of the measurements as reported in the legend. The position of the Sun is marked with a black star.</figcaption></figure>

<table><caption>**Table 3.** GAPS transiting systems with $\sigma R_5 > 0.5$.</caption><thead><tr><th>System</th><th>Time span (d)</th><th>σR<sub>5</sub></th></tr></thead><tbody><tr><td>Qatar-2</td><td>0.1713</td><td>1.736</td></tr><tr><td>WASP-1</td><td>1433.944</td><td>1.186</td></tr><tr><td>HAT-P-16</td><td>1849.973</td><td>1.036</td></tr><tr><td>XO-4</td><td>1164.830</td><td>1.009</td></tr><tr><td>WASP-12</td><td>1955.973</td><td>0.989</td></tr><tr><td>XO-3</td><td>1481.935</td><td>0.789</td></tr><tr><td>HAT-P-6</td><td>1859.903</td><td>0.648</td></tr><tr><td>WASP-57</td><td>426.063</td><td>0.639</td></tr><tr><td>HAT-P-21</td><td>1494.9679</td><td>0.588</td></tr><tr><td>HAT-P-26</td><td>0.1882</td><td>0.581</td></tr><tr><td>WASP-35</td><td>2661.838</td><td>0.571</td></tr></tbody></table>

distribution in the activity range of $-5.0 \le \log(R'_{HK}) \le -4.0$. There are only two M stars (red circles), too few to say anything about stars of this spectral type. We can just highlight their outlying position in the diagram. The outer atmosphere of these stars, due to their faintness in the visible wavelength range, still remains poorly understood. Despite this, there are several analyses of the activity of M stars (e.g., Maldonado et al. 2017, and reference therein). The FGK host stars' behavior seems to recall what Gomes da Silva et al. (2021) found by analyzing the distribution of the activity level of about 1650 stars for the different spectral types. In particular, for K stars, they found a distribution with three peaks: an inactive peak at -4.9 dex, an intermediate peak near -4.75 dex and an active peak near -4.5 dex.

We have too few K stars to recognize the peaks reported by Gomes da Silva et al. (2021), but it is possible to observe that also for our sample, the K stars show a large scatter.

## 6.1. Comparison with MS and evolved stars

To put the activity of host stars into context, we exploited the sample composed by Pace (2013) with a compilation of S-index measurements for about 2000 field stars from the literature. To remain homogeneous with the planetary host data, we search

<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Values of $\log(R'_{HK})$ as a function of the $T_{eff}$ of the star of sample. The color code is explained in the legend and the median error of both $T_{eff}$ and $\log(R'_{HK})$ are also indicated (in black).</figcaption></figure>

each field star in Gaia DR3 and eventually we evaluate the ($B-V$) color and the temperature for each star in the same way as described in Sect. 3. Since MS and evolved stars have very different activity distributions (Wright 2004; Mittag et al. 2013; Staab et al. 2017), to separate them in the field star sample, we used the absolute magnitude obtained from the XHIP catalog (Anderson & Francis 2012) and built an HR diagram (Fig. 11) using the average MS defined by Wright (2004) as in Sect. 3.

We computed the $\log(R'_{HK})$ values for both the MS and evolved stars identified in the Pace (2013) sample, following the procedure described in Sect. 4.2. The two left panels of Fig. 12 show the $\log(R'_{HK})$ measurements for the planetary hosts compared with MS field stars (top panels) or evolved stars (bottom panels). In the top left panel, the basal limit for MS stars is also indicated with a horizontal black line. The basal level of activity is the stellar analog of the quiet Sun reached by stars with no chromospherically active regions (e.g., Schrijver 1987; Schröder et al. 2012). The basal activity limit of $\log(R'_{HK}) > -5.1$ is applied only for dwarfs with solar metallicity. It is a little bit different ($-5.15$) for MS stars with super-solar metallicity
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** HR Diagram for the field stars in the Pace (2013) sample. The HIPPARCOS average MS by Wright (2004) and the loci 0.45 and 2 mag above it are indicated in orange and green, as in Fig. 6.</figcaption></figure>

(Henry et al. 1996; Wright 2004; Saar 2011). Furthermore, the basal level is also a slightly complicate function of the ($B-V$) of the star (Mittag et al. 2013, and references therein). Considering the polynomial model of the basal level described in Table 3 of Mittag et al. (2013), we superimposed it (the blue dashed curve) on the top-left panel of Fig. 12. The largest difference between the two descriptions of the basal level is visible in the red portion of the plot where the model has a quadratic dependence on the $B-V$ color. The presence of a relative maximum in the range $1.1 \le B-V \le 1.45$ has been justified by Mittag et al. (2013) with the absence of inactive stars in that color range. Instead, in the bluer range of the diagram, the difference between the two is inside the median measurement error. In any case, the polynomial description does not represent the basal limit in a better way than a constant value of $-5.1$ (e.g. Henry et al. 1996; Wright 2004; Staab et al. 2017; Haswell et al. 2020).

The right panels of Fig. 12 compare the histogram and cumulative distribution function of $\log(R'_{HK})$ for both the host and MS stars (top – right), and for host and slightly evolved stars (bottom – right). A comparison of the hosts with MS and evolved star populations shows that they actually come from different parent distributions in both cases (see Appendix B). This difference can be interpreted as due to observational selection biases that make planet detection more likely for less active hosts.

Among the MS stars in Pace (2013), there are very few stars (about 0.3%) below the basal activity level of $-5.1$, and about 1% in the complete sample. In our sample of transiting-planet hosts, ten stars (CoRoT-1, GJ436, HAT-P-35, HAT-P-44, HAT-P-45, HAT-P-46, Kepler-25, WASP-12, WASP-18, WASP-72), corresponding to 8% of our sample, show activity levels below the basal limit. Two of these, WASP-12, and WASP-72, are reported to be sub-giants having the $\Delta M_V > 0.45$ (0.74, and 1.26, respectively, see Sect. 3)<fnref n="8" />. Evolved stars can be inactive

stars well below the basal limit for MS stars (e.g., Wright 2004). The remaining host stars are all MS stars and in those cases, a mechanism that suppresses the star activity, making it lower than the basal limit should be at work. The activity suppression seems more common in our sample of planet-transiting hosts than in the field stars. This mechanism could also veil the CE of more active stars, but in this case, it will be very difficult to identify it. Hence, the activity suppression due to evaporated shrouds could be a more common phenomenon than it appears.

The presence of some planetary hosts below the MS activity basal limit has been interpreted as an indication that the presence of short-period massive planets mimics an activity suppression, but it is actually due to circumstellar absorption linked to planetary mass loss, which is sensitive to the surface gravity of the planet (Lanza 2014; Fossati et al. 2015). The same behavior of some MS stars has been considered, with success, which is indicative of undiscovered mass-losing planets (Haswell et al. 2020). A model that reproduces the observed correlation was proposed by Lanza (2014) and Fossati et al. (2015). On the other hand, since the transiting planet hosts are also distant stars (see the upper panel of Fig. 2 for the sample under study), a contribution to this suppression of the observed chromospheric index may be also given by the interstellar medium (Staab et al. 2017).

## 6.2. Considering whether the ISM absorption could bias the chromospheric activity level

The value of $\log(R'_{HK})$ of 8% of our sample of stars is lower than the activity basal limit for the MS. One of the causes that can contribute to justifying this low value of the $\log(R'_{HK})$ is the possible absorption of the HK flux by the interstellar medium (ISM).

Figure 13 shows the distribution of the transit hosts into the Galaxy. The red stars indicate the line of view of the hosts with an activity level lower than the basal one. In the plot, there is no evidence of a preferred direction. As we saw in Sect. 3, most of the stars we are considering have distances greater than 100 pc. For these systems, the ISM absorption has a significant influence on the measurement of the activity. As a matter of fact, beyond 100 pc the column density of the ISM Ca II is about $10^{12} \text{ cm}^{-2}$ (Welsh et al. 2010; Wyman & Redfield 2013; Fossati et al. 2017).

Following the work of Fossati et al. (2017), for each target we evaluate the column density of ISM Ca II. We utilized the ISM maps by Lallement et al. (2019, 2022) to evaluate $E(B-V)$ and successively the relation by Diplas & Savage (1994) to derive the column density of H I. The relation between the $N_{HII}$ and the ionic gas-phase abundance of ISM by Wakker & Mathis (2000) gave us the column density of Ca II. Finally, to obtain the correction of the $\log(R'_{HK})$ for the ISM absorption, we used the IDL procedure CORRECTION.PRO, described by Fossati et al. (2017, more details are given in Table A.2).

We note that with the only exception of HAT-P-44 and HAT-P-45, all the other eight MS systems with $\log(R'_{HK}) \le -5.1$ remain below the basal level after the application of the correction (see the systems highlighted in Table A.2). We can conclude that for these systems, the ISM absorption only partially explains the suppression of the CA index below the basal level, so that other mechanisms, such as absorption from evaporated shrouds, are at work.

<footnote n="8">⁸ We note in passing that the classification of WASP-12 as a subgiant, that is supported by the modeling of our HARPS-N spectra, may account for the decay rate of the orbit of its transiting planet due to the dissipation of tidal inertial gravity waves excited by the planet inside the star. According to the present theory of dynamic tide, those waves have a sufficient dissipation rate only if the core of WASP-12 is radiative, as generally predicted by stellar structure models in the case of a subgiant,</footnote>

<footnote>but not in the case of a MS star because the latter would have a convective core leading to a much smaller dissipation (e.g., Bailey & Goodman 2019; Barker 2020).</footnote>
<figure><img src="image_13.png" /><figcaption>Fig. 12. Comparison of the CA of host stars with the CA of MS and evolved stars. Top panel: distribution of the MS-field stars with $\Delta M_V \le 0.45$ above the average HIPPARCOS MS (cyano small filled circles) and the transiting planet hosts (left). In the plot, we also indicated two representations of the activity basal limit for MS stars. The black dashed line is the generally considered constant limit. The blue dashed line is the polynomial model of the basal level (Mittag et al. 2013). The comparison of the normalized histograms and the cumulative distribution functions of the two distributions in the left panel (right). Bottom panel: evolved field stars with $0.45 \le \Delta M_V \le 2.00$ above the average Hypparcos MS (small filled green circles) and the transiting planet hosts (of both samples) are on the left. A comparison of the normalized histograms and the cumulative distribution functions of the two distributions in the left panel is on the right.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>Fig. 13. Aitoff projection of the distribution of transit hosts in our sample in Galactic coordinates l and b of sight lines. The Galactic center is located at the center of the plot and the longitude increases counterclockwise. In the plot, the hosts with $\log(R'_{HK}) \le -5.1$ are indicated as red stars.</figcaption></figure>

## 7. Host activity and planetary parameters analysis

We searched for correlations between the $\log(R'_{HK})$ and planetary parameters, such as the mass ($M_P$), the radius ($R_P$), the density ($\rho_P$), and the orbital distance ($a_P$) of the planets. We checked also for correlation of the activity level with the planet's insolation ($S_*/S_0$; $S_0$ is the Earth's insolation at the top of the atmosphere), the equilibrium temperature ($T_{eq,P}$) of the planet assuming zero bond albedo. Considering also the star-planet interaction (SPI), we take into account the possibility that tidal interaction between the planet and its star can generate an enhancement of local turbulent velocity, resulting in an increment of the heating and a rise in stellar activity level (Cuntz et al. 2000). Hence, we searched for a correlation with the stellar gravitational perturbation by the planet (Cuntz et al. 2000):

$$ \frac{\Delta g_{\star}}{g_{\star}} = \frac{M_{P}}{M_{\star}} \frac{2R_{\star}^{3}}{a_{P}^{3}} $$

where $M_P$, $M_\star$, and $R_\star$ are the mass of the planet, the mass of the star, and the radius of the star, respectively; $g_\star$ is the surface gravity of the host star, and $a_\star$ is the orbital axis of the planet.

While the former set of variables consists of planetary quantities, the latter set contains quantities evaluated with both star and planet parameters. We used planetary data determined in a previous work (Bonomo et al. 2017), while for planets in our sample not considered in that work, we used the Exoplanet Encyclopedia<fnref n="9" /> (Schneider et al. 2011), searching the literature, to

<footnote n="9">⁹ Accessible from http://exoplanet.eu/</footnote>
<table><caption>**Table 4.** Statistical analysis of the correlation between the CA indicator <log(R'<sub>HK</sub>) > and different planetary parameters of the transiting planetary systems.</caption><thead><tr><th>Parameters</th><th>Condition</th><th>ν = N − 2</th><th>ρ<sup>(a)</sup></th><th>FAP</th></tr></thead><tbody><tr><td rowspan="2">M<sub>P</sub></td><td>None</td><td>140</td><td>0.029</td><td>73%</td></tr><tr><td>a<sub>p</sub> ≤ 0.1 au and M<sub>p</sub> ≥ 0.1 M<sub>J</sub></td><td>109</td><td>0.081</td><td>40%</td></tr><tr><td rowspan="2">R<sub>P</sub></td><td>None</td><td>140</td><td>−0.216</td><td>1.0%</td></tr><tr><td>a<sub>p</sub> ≤ 0.1 au and M<sub>p</sub> ≥ 0.1 M<sub>J</sub></td><td>109</td><td>−0.315</td><td>0.09%</td></tr><tr><td rowspan="2">ρ<sub>P</sub></td><td>None</td><td>140</td><td>0.186</td><td>3.0%</td></tr><tr><td>a<sub>p</sub> ≤ 0.1 au and M<sub>p</sub> ≥ 0.1 M<sub>J</sub></td><td>109</td><td>0.256</td><td>0.7%</td></tr><tr><td rowspan="2">a<sub>P</sub></td><td>None</td><td>137</td><td>−0.179</td><td>4.0%</td></tr><tr><td>a<sub>p</sub> ≤ 0.1 au and M<sub>p</sub> ≥ 0.1 M<sub>J</sub></td><td>109</td><td>−0.149</td><td>12%</td></tr><tr><td rowspan="2">log(g<sub>P</sub>)</td><td>None</td><td>140</td><td>0.204</td><td>1.5%</td></tr><tr><td>a<sub>p</sub> ≤ 0.1 au and M<sub>p</sub> ≥ 0.1 M<sub>J</sub></td><td>109</td><td>0.219</td><td>2.3%</td></tr><tr><td rowspan="2">S<sub>★</sub>/S<sub>0</sub></td><td>None</td><td>137</td><td>−0.213</td><td>1.0%</td></tr><tr><td>a<sub>p</sub> ≤ 0.1 au and M<sub>p</sub> ≥ 0.1 M<sub>J</sub></td><td>109</td><td>−0.322</td><td>0.06%</td></tr><tr><td rowspan="2">T<sub>eq,P</sub></td><td>None</td><td>137</td><td>−0.213</td><td>1.0%</td></tr><tr><td>a<sub>p</sub> ≤ 0.1 au and M<sub>p</sub> ≥ 0.1 M<sub>J</sub></td><td>109</td><td>−0.334</td><td>0.04%</td></tr><tr><td rowspan="2">Δg<sub>★</sub>/g<sub>★</sub></td><td>None</td><td>140</td><td>0.007</td><td>94%</td></tr><tr><td>a<sub>p</sub> ≤ 0.1 au and M<sub>p</sub> ≥ 0.1 M<sub>J</sub></td><td>108</td><td>−0.047</td><td>63%</td></tr></tbody></table>

**Notes.** In the second column, the condition of the sample is given; N is the total number of planets in the sample that fulfils the condition. In the last two columns we reported the Spearman's coefficient, ρ, and the value of the probability of the null hypothesis. (a) Rank correlation of two populations; the false alarm probability (FAP, in the last column) denotes the two-sided significance of its deviation from 0 by random chance, i.e., small values indicate a significant correlation.

obtain the mass, and radius and all the other useful parameters (a<sub>p</sub>, log(g<sub>P</sub>), etc.). We are aware that the derivation of the planetary parameters obtained in this way has a degree of inhomogeneity, which can produce biased results, and we stress that the investigation that follows should be seen with this caveat in mind.

The analysis of all the possible correlations has been performed using the Spearman's rank coefficient (Press et al. 1992), as was already done on smaller samples of transiting planets by Hartman (2010) and Figueira et al. (2014). For each correlation coefficient, we obtained the false alarm probability (FAP) or the probability of having a larger or equal correlation coefficient under the hypothesis that the data pairs are uncorrelated (our null hypothesis). The tests were performed considering the whole sample of host-planet pairs with no restrictions and a second sub-set responding to the condition: a<sub>p</sub> ≤ 0.1, and M<sub>p</sub> ≥ 0.1 M<sub>J</sub>. In Table 4, we summarize the results.

Among the planetary parameters, the correlation between log(R'<sub>HK</sub>) and M<sub>P</sub> has a value of Spearman's correlation coefficient that is not different from zero with a FAP of 73% considering the whole sample. This seems to contradict the statistical explanation of Collier Cameron & Jardine (2018) of the correlation reported by Hartman (2010) and Figueira et al. (2014) between the CE and the surface gravity of close-orbiting planets based on the statistical bias towards younger systems produced by the faster inward migration of more massive planets because of their stronger tidal interaction with their host stars. The correlation with the orbital semi-major axis a<sub>p</sub> has a FAP = 4%. This is the only SPI correlation that has been considered so far by further dedicated studies (e.g. Figueira et al. 2014), to the best knowledge of the authors. Instead, the correlation with the R<sub>P</sub> (FAP = 1.0% and FAP = 0.09%) have a confidence level greater than 99%, while log g<sub>P</sub> (FAP = 1.5% whole sample and 2.3% the

limited one), and the density of the planet ρ<sub>P</sub> (FAP = 3% and FAP = 0.7%), have all a confidence in the range between 97% and 99%.

Concerning the parameter containing stellar and planetary quantities, the correlation with insolation S<sub>★</sub>/S<sub>0</sub>, and the planetary T<sub>eq,P</sub> show also a confidence greater than 99%. This is not the same result, however, for the quantity Δg<sub>★</sub>/g<sub>★</sub> (FAP = 94% and 63%). In Fig. 14, we reported the value of log(R'<sub>HK</sub>) of the object as a function of the ratio, η, between the semi-major axis, a<sub>p</sub>, of the planet and the value of the Roche limit of the system (Faber et al. 2005; Ford & Rasio 2006):

$$a_{\text{Roche}} = 2.16 R_{\text{P}} \left( \frac{M_{\star}}{M_{\text{P}}} \right)^{1/3}, \qquad (7)$$

where a<sub>Roche</sub> is the critical separation at which the planet fills its Roche's lobe and it begins to lose mass. This limit is a function of the planetary radius R<sub>p</sub>, and of the stellar (M<sub>★</sub>) and planetary (M<sub>P</sub>) masses.

The analysis of the correlation between η and log(R'<sub>HK</sub>) does not give a significant correlation (Fig. 14). The Spearman's rank coefficient results ρ = 0.02 with a FAP ~ 83%, but it is worth noting that all the hosts with log(R'<sub>HK</sub>) lower than the basal limit have planets with a mean density (1.7×10<sup>3</sup> kg m<sup>-3</sup>) that is lower than the average density (2.2×10<sup>3</sup> kg m<sup>-3</sup>) of the planets with the ratio η in the interval [1, 10]. Actually, in the same range of η values, there is a high spread in the activity level of the hosts. The spread could be due to the dispersion of the host age and/or to the spread of the planetary companions' densities. In fact, in both cases, they are the stellar and planetary parameters with higher standard deviations by their mean values. There is not a clear indication of which parameter can account for it. Another reason that accounts for the spread in this region of the plot could
<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Activity level as a function of the ratio $\eta$ between the semi-major axis and the Roche lobe radius. The hosts with $\log(R'_{HK}) < -5.1$ are clearly identified.</figcaption></figure>

be the presence of matter evaporated by the planet that, depending on its quantity, could veil (with different levels of efficiency) the CE of the stars.

In the region of the plot for $\eta > 10$, the distribution of the activity levels has a reduced spread, and it is possible to define a border-line ($\log(R'_{HK}) = 0.27 \log(\eta) - 5.41$) that separates the space in the diagram between a region that is populated by the hosts of transiting planets and one that is empty. Because we are considering only transiting planets, the population of the $\eta > 10$ zone in the $\eta-\log(R'_{HK})$ plane could be the result of a selection effect. Planets observed with the RV method (e.g., the sample discussed by Canto Martins et al. 2011) could help to populate this zone and unveil whether there is a selection effect present. For those planets, the value of the radius is unknown and the mass obtained is the minimum mass, making it impossible to evaluate the $\eta$ value for those systems using measured quantities.

To conclude this section, we can say that there is no evidence for a strong single correlation between the activity level of the planet host star and one of the planetary and stellar parameters. Instead, a slight correlation seems to arise when we considered planetary parameters depending on the radius of the planet (surface gravity, density, insolation). All the parameters that have a correlation with $\log(R'_{HK})$ with a confidence value greater than 98% are $\propto R_p^n$ suggesting a common origin. The same is true for the insolation and the $T_{eq,P}$ that show a negative correlation with the activity. The analysis of these correlations confirms that close-in transiting planets swell under the irradiation from their host stars corroborating the hypothesis already discussed in preceding works (e.g., Lanza 2014; Fossati et al. 2013, 2015) that the planetary atmosphere evaporates under the action of the stellar irradiation and the matter lost by the planet forms shrouds that absorb at the core of the chromospheric resonance lines (e.g., Mg II h&k, Ca II H&K), making the observed chromospheric emission correspondingly lower. Due to the role of the planetary radius in these correlations, the mechanism at work could be the hydrodynamic escape to which some of these planets (those with the $\log(R'_{HK})$ lower than the basal value) should be subject. In this direction, we get also an indication that hosts with $\log(R'_{HK})$ lower than the basal limit for MS stars are very close to filling their Roche lobe even if other unknown phenomena could have their effects.

On the other hand, it is worth noting that the analyzed correlations with the possible star-planet interaction ($a_p$ and $\Delta g_\star/g_\star$)

are not significant ones. This seems to indicate that the higher activity level of some of the considered transit hosts in this study is not due to the SPI, but it appears natural.

## 8. Machine learning classification

In order to test the Lanza (2014) model of planetary material absorbing at the core of chromospheric resonance lines, Fossati et al. (2015) applied a cluster-weighted model (Ingrassia et al. 2012) to distinguish different sub-populations in their data. This model reproduces the bimodal distribution observed in the CA level of solar-like stars in the field, the so-called *Vaughan-Preston gap* (e.g., Wright 2004). Specifically, the existence of two linear regressions can be considered to be the natural outcome of the underlying bimodal distribution of the chromospheric index. These authors explored the likelihood of one to three components in their data set and found two components to be the most likely representation.

Following these results, we also classified our data based on the $\log(R'_{HK})$ and the inverse of the planetary surface gravity, $g_p^{-1}$. In this analysis, we only considered systems with main sequence stars and excluded those with an eccentric orbit and planetary masses in the Neptunian range. For the latter, in particular, there would not be enough planetary material available to create the absorbing torus predicted by Lanza (2014). Moreover, as the effect of planetary material on stellar chromospheric lines is only expected for close-in planets, we excluded systems with $a_p > 0.1$ AU.

Our data set is larger than that of Fossati et al. (2015)'s by almost a factor of 3, and we visually identified four possible components: the same high- and low-activity components they highlighted, a third component with average activity and high inverse of the surface gravity $g_p^{-1}$, namely, low planetary-surface gravity and a high-activity component with low correlation with $g_p^{-1}$. We hypothesize the third population might represent a "saturated" regime for stellar irradiation-driven planet evaporation, where planetary mass loss $\dot{M}_p$ becomes proportional to the square root of the stellar UV flux (Murray-Clay et al. 2009), rather than being proportional to the UV flux as assumed in the energy-limited evaporation model by Lanza (2014). Therefore, for this subset, a lower evaporation rate for a given $g_p^{-1}$ than assumed by Lanza (2014) would hold and this would imply a lower absorption by the evaporated matter than predicted by his model. The fourth population might instead represent a set of stars with an intrinsically high level of activity and a strong enough stellar wind to blow away all evaporated matter, making the measured stellar chromospheric emission virtually unaffected by planet evaporation.

In light of the already available information, as well as our aforementioned hypotheses, we implemented a supervised machine-learning classification scheme to determine a statistically-motivated separation between the subsets. The details of its training and implementation can be found in Appendix C. This resulted in the grouping presented in Fig. 15, where each data point is represented by a marker and colour corresponding to its predicted cluster membership, and where white-filled data points are those with an assignment probability <70%: a similar threshold was used by Fossati et al. (2015) to identify points that might likely belong to a different cluster. The high- and low-activity, low $g_p^{-1}$ populations (recalling those of Fossati et al. 2015) were then used to carry out linear fits in the $\log(R'_{HK}) - g_p^{-1}$ space, following the Lanza (2014) model. This required including the data point uncertainties both in the classification and fits: to do this, we first conservatively expanded
<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Clustering in the $g_p^{-1} - \log(R'_{HK})$ space provided by our decision tree implementation. Each cluster is represented by a different color and marker. Linear fits of $\log(R'_{HK})$ as a function of $g_p^{-1}$, as discussed in the text, are also shown. The portion of parameter space and the linear relationships found by Fossati et al. (2015) are shown with the gray region and black lines, respectively.</figcaption></figure>

the error bars of the high-activity systems ($\log(R'_{HK}) > -4.75)$ to 0.25 dex, based on the values shown in Fig. 3 of Wright (2004). This allowed us to take into account realistic variations due to the activity cycle. We then bootstrapped our data set 1000 times, on the basis of each data point's Gaussian uncertainty in $\log(R'_{HK})$. At every iteration, the previously-trained classification was carried out on the new data set, and linear fits were performed on the synthetic data set. For every realization, we only considered points with a >70% probability of being assigned to a given cluster, to then fit the populations with the linear model $\log R'_{HK} = \log R'_{HK}^{(0)} + \gamma g_p^{-1}$. We obtained:

- $\gamma = -336.16 \pm 181.66$, $\log R'_{HK}^{(0)} = -4.19 \pm 0.10$ for the high-activity cluster;

- $\gamma = -54.26 \pm 44.56$, $\log R'_{HK}^{(0)} = -4.87 \pm 0.04$ for the low-activity cluster.

Both populations are in $3\sigma$ agreement with Fossati et al. (2015): the parameter space investigated by these authors and the trends they identified are shown in Fig. 15. The $\gamma$ coefficients of both populations are, however, less significant than those found by Fossati et al., while the one on the high-activity population is more significant than the one found by Daniëlski et al. (2022). Such differences might relate to different factors, such as the kind of classification algorithm, the details of the conversion between the S and the $\log(R'_{HK})$ activity indexes, and, most importantly, the size of the sample. All are aspects that play a key role in the classification process and identified correlations. Another relevant aspect is the probability of assignment to the high- or low-activity cluster or even to the high-activity and higher $g_p^{-1}$ group, of single data points. For these reasons, we highlight the need to perform further investigations which use a larger homogeneous sample, to shed light on how planetary systems populate this parameter space. This might reduce the reliance on training information and allow for the use of unsupervised algorithms to increase the robustness of the clustering and, therefore, the estimates of the correlation coefficients.

## 9. Summary and conclusions

In this work, we have presented the characterization of a part of the hosts of transiting planets of the GAPS sample, mainly

analyzing their stellar chromospheric activity (CA). In the process, we added to the observed data others found in the literature obtained for other hosts of transiting planets, enlarging the data set on which the cited previous works were based.

For each of the 132 stars in the sample, we obtained the stellar parameters ($T_{eff}$, $B-V$, $V$ in a homogeneous way, exploiting the Gaia DR3 data release (Gaia Collaboration 2016, 2018, 2023). Subsequently, the complete set was compared with the activity levels of the field stars. To do this, we used the Mittag et al. (2013) procedure for calculating $\log(R'_{HK})$ values. We report an activity index between -5.373 and -3.972, with a median value of -4.810. Somewhat more than half of the sample (58%) are inactive stars and this was expected, as most of our sample was selected from radial velocity follow-up campaigns, which privilege inactive stellar hosts. This study allowed us to determine a relationship between the $\log(R'_{HK})$ values calculated with the classical procedure (Noyes et al. 1984) and those obtainable with Mittag's procedure, in the event that the S-index value is not available for the given object. This relationship has been determined both for MS and slightly evolved objects.

The analysis of the time series of $\log(R'_{HK})$ for the GAPS objects observed for more than 20 nights showed that some of the targets have a dispersion value $\sigma R_5 > 0.5$ dex. In one case, Qatar-2, the high dispersion value (1.736 dex) and the low time-lapse of the time series seem to indicate a flare event. In general, as expected, the activity dispersion pattern is proportional to the activity level of the star, but in the cases of the F stars WASP-1, HAT-P-16, XO-4, and in that of the K star HAT-P-26, although these stars have a $\log(R'_{HK})$ value included in the zone of inactive stars ($-5.1 \le \log(R'_{HK}) \le -4.75$), they show a large variation in activity. This behavior confirms the high activity variation for K and F stars with low activity levels observed by Gomes da Silva et al. (2021).

We contrasted also the transiting planet hosts (both GAPS and those taken from the literature) with the field objects of both MS and slightly evolved by Pace (2013) after that we obtained also for field stars the stellar parameters by Gaia/DR3 in order to compare homogeneous samples. Kolmogorov-Smirnov tests show that the population of transiting objects is different from both populations (stars of MS and sub-giants) of the field objects. In addition, the population of transiting objects shows
a higher percentage of objects (8% of the sample corresponding to 11 stars) with an activity level lower than the basal activity level for main sequence stars than in the case of field objects (1%). While three of these objects are subgiants, the other eight are stars of MS, for which there is a process of suppression of chromospheric line core fluxes that we have ascertained not to be due to the interstellar medium.

We analyzed the presence or the absence of possible correlations between the level of CA of the host stars and the stellar and planetary parameters of the systems considered. In particular, we tested the correlations of the activity level with the following planetary parameters: $M_P$, $R_P$, log $g_P$, $\rho_P$ and $a_P$, the insolation at the top of the planet's atmosphere, the $T_{eq,p}$, the equilibrium temperature of the planet and the stellar equilibrium tide. The analysis has been performed for the whole sample and for a limited sub-set to $M_P \ge 0.1 M_J$ and $a_P \le 0.1$ au. We did not find significant correlations with $a_P$, the only planetary parameter linked to the existence of star-planet interaction. Furthermore, the lack of a correlation with $M_P$ seems to contradict the Collier Cameron & Jardine (2018)'s explanation of the correlation firstly reported by Hartman (2010). Conversely, all correlations with FAP $\le 1\%$ are with planetary parameters proportional to $R_P^n$, indicating a common origin. We conclude that some planets evaporate by means of hydrodynamic escape under the strong irradiation of the host star and the evaporated matter forms a system of shrouds that absorbs the core of chromospheric resonance lines, as described by Lanza (2014). This is confirmed by the fact that all the systems with activity lower than the basal level host planets that are very close to filling their Roche lobe.

An analysis made with a machine learning technique supported the presence of four different clusters of hosts with different behavior concerning the activity of stars as a function of the planetary surface gravity. Two of these clusters were already individuated by Fossati et al. (2015). The third cluster, with average activity and low planetary surface gravity, has been interpreted as due to a "saturated" regime for stellar irradiation-driven evaporation. The fourth cluster likely represents stars with a high level of activity and a strong stellar wind able to remove the shrouds of evaporated matter. So far, these interpretations are qualitative ones and they are the subject of ongoing and future works.

**Acknowledgements.** The authors are grateful to an anonymous referee for a careful reading of the manuscript and a very useful report that helped them to remarkably improve their work and its presentation. The GAPS project in Italy acknowledges the support by INAF through the Progetti Premiali funding scheme of the Italian Ministry of Education, University, and Research. Based on observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the island of La Palma by the Fundacion Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias This research has made use of the VizieR catalog access tool, CDS, Strasbourg, France (DOI: 10.26093/cds/vizier) and of SIMBAD database, operated at CDS, Strasbourg, France. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular, the institutions participating in the Gaia Multilateral Agreement.

## References

* Albrecht, S., Winn, J. N., Butler, R. P., et al. 2012, ApJ, 744, 189

* Alspach, D., & Sorenson, H. 1972, IEEE Trans. Automatic Control, 17, 439

* Anderson, E., & Francis, C. 2012, Astron. Lett., 38, 331

* Anderson, D. R., Collier Cameron, A., Gillon, M., et al. 2011, A&A, 534, A16

* Anderson, D. R., Collier Cameron, A., Delrez, L., et al. 2014, MNRAS, 445, 1114

* Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, AJ, 161, 147

* Bailey, A., & Goodman, J. 2019, MNRAS, 482, 1872

* Bakos, G. Á., Hartman, J. D., Torres, G., et al. 2012, AJ, 144, 19

* Baliunas, S. L., Donahue, R. A., Soon, W. H., et al. 1995, ApJ, 438, 269

* Ballerini, P., Micela, G., Lanza, A. F., & Pagano, I. 2012, A&A, 539, A140

* Barclay, T., Kostov, V. B., Colón, K. D., et al. 2021, AJ, 162, 300

* Barker, A. J. 2020, MNRAS, 498, 2270

* Bastian, T. S., Dulk, G. A., & Leblanc, Y. 2000, ApJ, 545, 1058

* Beck, P. G., Kambe, E., Hillen, M., et al. 2015, A&A, 573, A138

* Béky, B., Bakos, G. Á., Hartman, J., et al. 2011, ApJ, 734, 109

* Benatti, S. 2021, in Posters from the TESS Science Conference II (TSC2), 21

* Biazzo, K., D'Orazi, V., Desidera, S., et al. 2022, A&A, 664, A161

* Bonomo, A. S., Hébrard, G., Santerne, A., et al. 2012, A&A, 538, A96

* Bonomo, A. S., Desidera, S., Benatti, S., et al. 2017, A&A, 602, A107

* Borsa, F., Rainer, M., Bonomo, A. S., et al. 2019, A&A, 631, A34

* Borucki, W. J., Koch, D. G., Batalha, N., et al. 2012, ApJ, 745, 120

* Borucki, W. J., Agol, E., Fressin, F., et al. 2013, Science, 340, 587

* Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127

* Brown, D. J. A., Cameron, A. C., Anderson, D. R., et al. 2012, MNRAS, 423, 1503

* Buchhave, L. A., Bakos, G. Á., Hartman, J. D., et al. 2011, ApJ, 733, 116

* Canto Martins, B. L., Das Chagas, M. L., Alves, S., et al. 2011, A&A, 530, A73

* Carleo, I., Benatti, S., Lanza, A. F., et al. 2018, A&A, 613, A50

* Carleo, I., Malavolta, L., Lanza, A. F., et al. 2020, A&A, 638, A5

* Carleo, I., Giacobbe, P., Guilluy, G., et al. 2022, AJ, 164, 101

* Cegla, H. M., Watson, C. A., Shelyag, S., et al. 2018, ApJ, 866, 55

* Claudi, R., Benatti, S., Carleo, I., et al. 2016, SPIE Conf. Ser., 9908, 99081A

* Claudi, R., Benatti, S., Carleo, I., et al. 2017, Eur. Phys. J. Plus, 132, 364

* Collier Cameron, A., & Jardine, M. 2018, MNRAS, 476, 2542

* Cosentino, R., Lovis, C., Pepe, F., et al. 2012, SPIE Conf. Ser., 8446, 84461V

* Covino, E., Esposito, M., Barbieri, M., et al. 2013, Science, 345, A54, A28

* Cuntz, M., Saar, S. H., & Musielak, Z. E. 2000, ApJ, 533, L151

* da Silva, L., Girardi, L., Pasquini, L., et al. 2006, A&A, 458, 609

* Damasso, M., Lanza, A. F., Benatti, S., et al. 2020, A&A, 642, A133

* Daniëlski, C., Brucalassi, A., Benatti, S., et al. 2022, Exp. Astron., 53, 473

* Diplas, A., & Savage, B. D. 1994, ApJ, 427, 274

* Dumusque, X., Bonomo, A. S., Haywood, R. D., et al. 2014, ApJ, 789, 154

* Duncan, D. K., Vaughan, A. H., Wilson, O. C., et al. 1991, ApJS, 76, 383

* Egeland, R., Soon, W., Baliunas, S., et al. 2017, ApJ, 835, 25

* Esposito, M., Covino, E., Desidera, S., et al. 2017, A&A, 601, A53

* Evans, D. W., Riello, M., De Angeli, F., et al. 2018, A&A, 616, A4

* Faber, J. A., Rasio, F. A., & Willems, B. 2005, Icarus, 175, 248

* Figueira, P., Oshagh, M., Adibekyan, V. Z., & Santos, N. C. 2014, A&A, 572, A51

* Ford, E. B., & Rasio, F. A. 2006, ApJ, 638, L45

* Fossati, L., Haswell, C. A., Froning, C. S., et al. 2010, ApJ, 714, L222

* Fossati, L., Ayres, T. R., Haswell, C. A., et al. 2013, ApJ, 766, L20

* Fossati, L., Ingrassia, S., & Lanza, A. F. 2015, ApJ, 812, L35

* Fossati, L., Marcelja, S. E., Staab, D., et al. 2017, A&A, 601, A104

* Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1

* Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1

* Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1

* Gautier, Thomas N., I., Charbonneau, D., Rowe, J. F., et al. 2012, ApJ, 749, 15

* Geron, A. 2017, Hands-on Machine Learning with Scikit-Learn and TensorFlow (Sebastopol, CA: O'Reilly Media, Inc.)

* Gilliland, R. L., Marcy, G. W., Rowe, J. F., et al. 2013, ApJ, 766, 40

* Gomes da Silva, J., Santos, N. C., Adibekyan, V., et al. 2021, A&A, 646, A77

* Gray, R. O., Corbally, C. J., Garrison, R. F., et al. 2006, AJ, 132, 161

* Guilluy, G., Andretta, V., Borsa, F., et al. 2020, A&A, 639, A49

* Gurdemir, L., Redfield, S., & Cuntz, M. 2012, PASA, 29, 141

* Hall, J. C., Lockwood, G. W., & Skiff, B. A. 2007, AJ, 133, 862

* Hartman, J. D. 2010, ApJ, 717, L138

* Hartman, J. D., Bakos, G. Á., Torres, G., et al. 2011, ApJ, 742, 59

* Hartman, J. D., Bakos, G. Á., Béky, B., et al. 2012, AJ, 144, 139

* Hartman, J. D., Bakos, G. Á., Torres, G., et al. 2014, AJ, 147, 128

* Haswell, C. A., Fossati, L., Ayres, T., et al. 2012, ApJ, 760, 79

* Haswell, C. A., Staab, D., Barnes, J. R., et al. 2020, Nat. Astron., 4, 408

* Hatzes, A. P. 1996, PASP, 108, 839

* Hebb, L., Collier-Cameron, A., Loeillet, B., et al. 2009, ApJ, 693, 1920

* Hébrard, G., Collier Cameron, A., Brown, D. J. A., et al. 2013, A&A, 549, A134

* Henry, T. J., Soderblom, D. R., Donahue, R. A., & Baliunas, S. L. 1996, AJ, 111, 439

* Howard, A. W., Johnson, J. A., Marcy, G. W., et al. 2011, ApJ, 730, 10
Hunter, A. A., Macgregor, A. B., Szabo, T. O., Wellington, C. A., & Bellgard, M. I. 2012, Source Code Biol. Med., 7, 1

Ingrassia, S., Minotti, S. C., & Vittadini, G. 2012, J. Classif., 29, 363

Kashyap, V. L., Drake, J. J., & Saar, S. H. 2008, ApJ, 687, 1339

Knutson, H. A., Howard, A. W., & Isaacson, H. 2010, ApJ, 720, 1569

Lallement, R., Babusiaux, C., Vergely, J. L., et al. 2019, A&A, 625, A135

Lallement, R., Vergely, J. L., Babusiaux, C., & Cox, N. L. J. 2022, A&A, 661, A147

Lammer, H., Selsis, F., Ribas, I., et al. 2003, ApJ, 598, L121

Lanza, A. F. 2014, A&A, 572, L6

Lecavelier des Étangs, A. 2007, A&A, 461, 1185

Lecavelier des Étangs, A., Vidal-Madjar, A., McConnell, J. C., & Hébrard, G. 2004, A&A, 418, L1

Lendl, M., Anderson, D. R., Collier-Cameron, A., et al. 2012, A&A, 544, A72

Lendl, M., Triaud, A. H. M. J., Anderson, D. R., et al. 2014, A&A, 568, A81

Locci, D., Cecchi-Pestellini, C., & Micela, G. 2019, A&A, 624, A101

Lovis, C., Dumusque, X., Santos, N. C., et al. 2011, arXiv e-prints [arXiv:1107.5325]

Madhusudhan, N. 2019, ARA&A, 57, 617

Maggio, A., Pillitteri, I., Scandariato, G., et al. 2015, ApJ, 811, L2

Malavolta, L., Borsato, L., Granata, V., et al. 2017, AJ, 153, 224

Maldonado, J., Scandariato, G., Stelzer, B., et al. 2017, A&A, 598, A27

Marcy, G. W., Isaacson, H., Howard, A. W., et al. 2014, ApJS, 210, 20

Maxted, P. F. L., Anderson, D. R., Collier Cameron, A., et al. 2011, PASP, 123, 547

Milbourne, T. W., Haywood, R. D., Phillips, D. F., et al. 2019, ApJ, 874, 107

Miller, B. P., Gallo, E., Wright, J. T., & Pearson, E. G. 2015, ApJ, 799, 163

Mittag, M., Schmitt, J. H. M. M., & Schröder, K. P. 2013, A&A, 549, A117

Montalto, M., Piotto, G., Marrese, P. M., et al. 2021, A&A, 653, A98

Murray-Clay, R. A., Chiang, E. I., & Murray, N. 2009, ApJ, 693, 23

Nardiello, D., Malavolta, L., Desidera, S., et al. 2022, A&A, 664, A163

Noyes, R. W., Bakos, G. Á., Torres, G., et al. 2008, ApJ, 673, L79

Noyes, R. W., Hartmann, L. W., Baliunas, S. L., Duncan, D. K., & Vaughan, A. H. 1984, ApJ, 279, 763

Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23

Pace, G. 2013, A&A, 551, L8

Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9

Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, J. Mach. Learn. Res., 12, 2825

Penz, T., Micela, G., & Lammer, H. 2008, A&A, 477, 309

Pepe, F., Mayor, M., Galland, F., et al. 2002, A&A, 388, 632

Pepe, F., Cameron, A. C., Latham, D. W., et al. 2013, Nature, 503, 377

Pillitteri, I., Wolk, S. J., Lopez-Santiago, J., et al. 2014, ApJ, 785, 145

Pino, L., Désert, J.-M., Brogi, M., et al. 2020, ApJ, 894, L27

Poppenhaeger, K., Robrade, J., & Schmitt, J. H. M. M. 2010, A&A, 515, A98

Poppenhaeger, K., & Schmitt, J. H. M. M. 2011, ApJ, 735, 59

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in C. The art of scientific computing (Cambridge: University Press)

Quinn, S. N., Bakos, G. Á., Hartman, J., et al. 2012, ApJ, 745, 80

Riello, M., De Angeli, F., Evans, D. W., et al. 2021, A&A, 649, A3

Robertson, P., Roy, A., & Mahadevan, S. 2015, ApJ, 805, L22

Saar, S. H. 2011, ASP Conf. Ser., 448, 1239

Saar, S. H., & Cuntz, M. 2001, MNRAS, 325, 55

Santos, N. C., Sousa, S. G., Mortier, A., et al. 2013, A&A, 556, A150

Sanz-Forcada, J., Ribas, I., Micela, G., et al. 2010, A&A, 511, L8

Sanz-Forcada, J., Micela, G., Ribas, I., et al. 2011, A&A, 532, A6

Sato, B., Hartman, J. D., Bakos, G. Á., et al. 2012, PASJ, 64, 97

Scandariato, G., Maggio, A., Lanza, A. F., et al. 2013, A&A, 552, A7

Scandariato, G., Maldonado, J., Affer, L., et al. 2017, A&A, 598, A28

Scharf, C. A. 2010, ApJ, 722, 1547

Schneider, J., Dedieu, C., Le Sidaner, P., Savalle, R., & Zolotukhin, I. 2011, A&A, 532, A79

Schrijver, C. J. 1987, A&A, 172, 111

Schröder, K. P., Mittag, M., Pérez Martínez, M. I., Cuntz, M., & Schmitt, J. H. M. M. 2012, A&A, 540, A130

Shkolnik, E., Walker, G. A. H., & Bohlender, D. A. 2003, ApJ, 597, 1092

Shkolnik, E., Walker, G. A. H., Rucinski, S. M., Bohlender, D. A., & Davidge, T. J. 2005, AJ, 130, 799

Shkolnik, E., Bohlender, D. A., Walker, G. A. H., & Collier Cameron, A. 2008, ApJ, 676, 628

Sorenson, H., & Alspach, D. 1971, Automatica, 7, 465

Sousa, S. G., Adibekyan, V., Delgado-Mena, E., et al. 2018, A&A, 620, A58

Sousa, S. G., Adibekyan, V., Delgado-Mena, E., et al. 2021, A&A, 656, A53

Staab, D., Haswell, C. A., Smith, G. D., et al. 2017, MNRAS, 466, 738

Triaud, A. H. M. J., Collier Cameron, A., Queloz, D., et al. 2010, A&A, 524, A25

Triaud, A. H. M. J., Queloz, D., Hellier, C., et al. 2011, A&A, 531, A24

Tsantaki, M., Sousa, S. G., Santos, N. C., et al. 2014, A&A, 570, A80

Vaughan, A. H., & Preston, G. W. 1980, PASP, 92, 385

Wakker, B. P., & Mathis, J. S. 2000, ApJ, 544, L107

Welsh, B. Y., Lallement, R., Vergely, J. L., & Raimond, S. 2010, A&A, 510, A54

Wenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9

Winn, J. N., Albrecht, S., Johnson, J. A., et al. 2011, ApJ, 741, L1

Wright, J. T. 2004, AJ, 128, 1273

Wright, J. T. 2005, PASP, 117, 657

Wright, J. T., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2004, ApJS, 152, 261

Wyman, K., & Redfield, S. 2013, ApJ, 773, 96

Zacharias, N., Finch, C. T., Girard, T. M., et al. 2012, VizieR Online Data Catalog: I/322A

¹ INAF – Osservatorio Astronomico di Padova, vicolo Osservatorio, 5, 35122 Padova, Italy
e-mail: riccardo.claudi@oapd.inaf.it

² Dipartimento di Matematica e Fisica, Università Roma Tre, Via della Vasca Navale 84, 00146 Roma, Italy

³ INAF-Osservatorio Astrofisico di Catania, via S. Sofia, 78, 95123 Catania, Italy

⁴ Space Research Institute, Austrian Academy of Sciences, Schmiedl-Strasse 6, 8042 Graz, Austria

⁵ INAF – Osservatorio Astronomico di Palermo, Piazza del Parla-mento 1, 90134 Palermo, Italy

⁶ INAF – Osservatorio Astronomico di Roma, Via Frascati, 33, 00078 Monte Porzio (RM), Italy

⁷ INAF – Osservatorio Astronomico di Trieste, via Tiepolo 11, 34143 Trieste, Italy

⁸ Instituto de Astrofísica de Canarias (IAC), 38205 La Laguna, Tenerife, Spain

⁹ INAF – Osservatorio Astrofisico di Torino, via Osservatorio 20, 10025 Pino Torinese, Italy

¹⁰ Instituto de Astrofísica de Andalucía, CSIC, Glorieta de la Astronomía, 18080 Granada, Spain

¹¹ Sorbonne Universités, UPMC Université Paris 6 et CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, France

¹² Dipartimento di Fisica e Astronomia Galileo Galilei, Università di Padova, Vicolo dell'Osservatorio 2, 35122 Padova (Pd), Italy

¹³ Dipartimento di Fisica, Università di Roma 'Tor Vergata', Via della Ricerca Scientifica 1, 00133 Rome, Italy

¹⁴ Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany

¹⁵ INAF – Osservatorio Astronomico di Brera, via E. Bianchi 46, 23807 Merate, Italy

¹⁶ Fundación Galileo Galilei – INAF, Rambla José Ana Fernandez Pérez 7, 38712 Breña Baja, TF, Spain

¹⁷ INAF – Osservatorio Astronomico di Capodimonte, Salita Moiariello 16, 80131, Napoli, Italy

¹⁸ INAF – Osservatorio Astronomico di Cagliari, Via della Scienza 5, 09047, Selargius (CA), Italy
Appendix A: Additional tables

<table><caption>Table A.1. Journal of GAPS Observations.</caption><thead><tr><th>System</th><th>T<sub>0</sub></th><th>Time span (d)</th><th>Spectra n.</th><th>System</th><th>T<sub>0</sub></th><th>Time span (d)</th><th>Spectra n.</th></tr></thead><tbody><tr><td>GJ436</td><td>2018-04-01</td><td>698.1226</td><td>47</td><td>TRES-1</td><td>2013-04-23</td><td>1534.9043</td><td>20</td></tr><tr><td>GJ3470</td><td>2018-01-13</td><td>744.2591</td><td>82</td><td>TRES-2</td><td>2013-04-24</td><td>1533.9317</td><td>15</td></tr><tr><td>HAT-P-1</td><td>2012-08-08</td><td>2587.8794</td><td>55</td><td>TRES-4</td><td>2013-03-11</td><td>1190.853</td><td>19</td></tr><tr><td>HAT-P-2</td><td>2013-03-12</td><td>2526.9941</td><td>29</td><td>WASP-1</td><td>2012-08-21</td><td>1433.9444</td><td>13</td></tr><tr><td>HAT-P-3</td><td>2013-01-30</td><td>791.6893</td><td>38</td><td>WASP-3</td><td>2012-08-07</td><td>25.949</td><td>2</td></tr><tr><td>HAT-P-4</td><td>2013-02-01</td><td>2565.9444</td><td>30</td><td>WASP-10</td><td>2012-11-12</td><td>0.1557</td><td>13</td></tr><tr><td>HAT-P-6</td><td>2012-08-27</td><td>1859.9033</td><td>18</td><td>WASP-11</td><td>2012-10-02</td><td>2301.7128</td><td>50</td></tr><tr><td>HAT-P-7</td><td>2013-04-24</td><td>2310.8561</td><td>21</td><td>WASP-12</td><td>2012-09-07</td><td>1955.9731</td><td>71</td></tr><tr><td>HAT-P-8</td><td>2012-09-04</td><td>1423.1446</td><td>14</td><td>WASP-13</td><td>2012-12-05</td><td>2577.9926</td><td>41</td></tr><tr><td>HAT-P-11</td><td>2019-07-07</td><td>0.1814</td><td>26</td><td>WASP-14</td><td>2013-01-30</td><td>1229.7274</td><td>15</td></tr><tr><td>HAT-P-12</td><td>2015-03-14</td><td>41.9912</td><td>37</td><td>WASP-21</td><td>2013-07-11</td><td>1540.8895</td><td>10</td></tr><tr><td>HAT-P-14</td><td>2013-03-25</td><td>1654.6531</td><td>14</td><td>WASP-24</td><td>2013-01-30</td><td>1281.6053</td><td>14</td></tr><tr><td>HAT-P-15</td><td>2014-10-10</td><td>403.0163</td><td>34</td><td>WASP-26</td><td>2012-09-05</td><td>1492.8679</td><td>12</td></tr><tr><td>HAT-P-16</td><td>2012-08-22</td><td>1849.9729</td><td>21</td><td>WASP-31</td><td>2013-01-13</td><td>1166.7565</td><td>14</td></tr><tr><td>HAT-P-17</td><td>2013-10-13</td><td>2159.2268</td><td>53</td><td>WASP-32</td><td>2012-09-06</td><td>1492.019</td><td>12</td></tr><tr><td>HAT-P-18</td><td>2013-06-11</td><td>689.2777</td><td>28</td><td>WASP-35</td><td>2012-09-06</td><td>2661.8378</td><td>40</td></tr><tr><td>HAT-P-20</td><td>2012-11-21</td><td>1136.8736</td><td>42</td><td>WASP-38</td><td>2013-03-24</td><td>1041.062</td><td>18</td></tr><tr><td>HAT-P-21</td><td>2013-01-03</td><td>1494.9679</td><td>59</td><td>WASP-39</td><td>2015-05-04</td><td>1099.1804</td><td>61</td></tr><tr><td>HAT-P-22</td><td>2013-01-03</td><td>2628.7719</td><td>68</td><td>WASP-43</td><td>2012-12-27</td><td>1156.8908</td><td>54</td></tr><tr><td>HAT-P-24</td><td>2012-11-21</td><td>1260.6462</td><td>17</td><td>WASP-48</td><td>2013-04-24</td><td>1533.9315</td><td>12</td></tr><tr><td>HAT-P-26</td><td>2015-03-27</td><td>0.1882</td><td>27</td><td>WASP-49</td><td>2013-12-30</td><td>1024.2611</td><td>9</td></tr><tr><td>HAT-P-29</td><td>2012-08-27</td><td>2684.7157</td><td>56</td><td>WASP-50</td><td>2012-09-30</td><td>1190.6492</td><td>43</td></tr><tr><td>HAT-P-30</td><td>2013-01-03</td><td>1386.0594</td><td>15</td><td>WASP-54</td><td>2013-01-13</td><td>1217.7109</td><td>16</td></tr><tr><td>HAT-P-31</td><td>2013-04-28</td><td>1142.8728</td><td>15</td><td>WASP-57</td><td>2015-04-02</td><td>426.0631</td><td>20</td></tr><tr><td>HAT-P-32</td><td>2017-11-17</td><td>731.139</td><td>32</td><td>WASP-59</td><td>2016-09-12</td><td>0.1925</td><td>19</td></tr><tr><td>HAT-P-36</td><td>2013-02-22</td><td>0.0866</td><td>9</td><td>WASP-60</td><td>2013-10-20</td><td>0.2334</td><td>21</td></tr><tr><td>HAT-P-47</td><td>2019-11-21</td><td>0.2101</td><td>31</td><td>WASP-69</td><td>2019-07-24</td><td>0.1742</td><td>25</td></tr><tr><td>HAT-P-50</td><td>2016-12-15</td><td>0.2822</td><td>40</td><td>WASP-76</td><td>2017-10-26</td><td>762.1462</td><td>110</td></tr><tr><td>HD17156</td><td>2012-08-22</td><td>1521.8131</td><td>25</td><td>WASP-80</td><td>2019-08-09</td><td>43.0984</td><td>32</td></tr><tr><td>HD80606</td><td>2020-02-02</td><td>0.4319</td><td>9</td><td>WASP-106</td><td>2016-02-01</td><td>0.0239</td><td>3</td></tr><tr><td>HD189733</td><td>2017-05-31</td><td>505.9913</td><td>123</td><td>WASP-107</td><td>2019-02-08</td><td>86.0099</td><td>38</td></tr><tr><td>HD209458</td><td>2018-07-08</td><td>423.1412</td><td>105</td><td>WASP-127</td><td>2019-02-20</td><td>0.2575</td><td>25</td></tr><tr><td>KELT-3</td><td>2018-01-19</td><td>0.1836</td><td>17</td><td>XO-1</td><td>2013-03-25</td><td>1243.8207</td><td>13</td></tr><tr><td>KELT-6</td><td>2014-02-10</td><td>854.8237</td><td>73</td><td>XO-2N</td><td>2012-11-21</td><td>2263.9112</td><td>73</td></tr><tr><td>QATAR-1</td><td>2012-09-03</td><td>1857.105</td><td>27</td><td>XO-3</td><td>2012-09-28</td><td>1481.9351</td><td>21</td></tr><tr><td>QATAR-2</td><td>2014-04-27</td><td>0.1713</td><td>14</td><td>XO-4</td><td>2012-10-26</td><td>1164.83</td><td>17</td></tr></tbody></table>
<table><caption>**Table A.2.** Correction ($\Delta \log R'_{HK}$) and corrected value of $\log R'_{HK}$ for the ISM absorption for the systems of the sample. The highlighted systems have still $\log(R'_{HK}) \le -5.1$.</caption><thead><tr><th>System</th><th>Δ log R'HK</th><th>log R'HK</th><th>System</th><th>Δ log R'HK</th><th>log R'HK</th><th>System</th><th>Δ log R'HK</th><th>log R'HK</th></tr></thead><tbody><tr><td>55Cnc</td><td>0.036</td><td>-4.8197</td><td>HAT-P-45</td><td>0.07</td><td>-5.0994</td><td>WASP-14</td><td>0.042</td><td>-4.8914</td></tr><tr><td><b>CoRoT-1</b></td><td>0.044</td><td>-5.2535</td><td><b>HAT-P-46</b></td><td>0.06</td><td>-5.1160</td><td>WASP-15</td><td>0.052</td><td>-4.7405</td></tr><tr><td>CoRoT-2</td><td>0.075</td><td>-4.1389</td><td>HAT-P-47</td><td>0.06</td><td>-4.6470</td><td>WASP-16</td><td>0.052</td><td>-4.8946</td></tr><tr><td>CoRoT-7</td><td>0.048</td><td>-4.5624</td><td>HD149026</td><td>0.04</td><td>-4.7937</td><td>WASP-17</td><td>0.063</td><td>-4.8192</td></tr><tr><td>GJ3470</td><td>0.040</td><td>-4.6296</td><td>HD17156</td><td>0.04</td><td>-4.8430</td><td><b>WASP-18</b></td><td>0.034</td><td>-5.2231</td></tr><tr><td><b>GJ436</b></td><td>0.028</td><td>-5.3447</td><td>HD189733</td><td>0.05</td><td>-4.1505</td><td>WASP-19</td><td>0.057</td><td>-4.4481</td></tr><tr><td>HAT-P-1</td><td>0.046</td><td>-4.8504</td><td>HD209458</td><td>0.04</td><td>-4.8178</td><td>WASP-21</td><td>0.052</td><td>-4.7728</td></tr><tr><td>HAT-P-2</td><td>0.045</td><td>-4.6369</td><td>HD80606</td><td>0.04</td><td>-4.8841</td><td>WASP-22</td><td>0.045</td><td>-4.7083</td></tr><tr><td>HAT-P-3</td><td>0.048</td><td>-4.6268</td><td>HD97658</td><td>0.04</td><td>-4.7709</td><td>WASP-23</td><td>0.051</td><td>-4.3866</td></tr><tr><td>HAT-P-4</td><td>0.043</td><td>-4.8215</td><td>KELT-3</td><td>0.05</td><td>-4.5809</td><td>WASP-24</td><td>0.053</td><td>-4.7846</td></tr><tr><td>HAT-P-5</td><td>0.048</td><td>-4.8729</td><td>KELT-6</td><td>0.05</td><td>-4.8749</td><td>WASP-26</td><td>0.039</td><td>-4.8971</td></tr><tr><td>HAT-P-6</td><td>0.063</td><td>-4.5834</td><td>Kepler-4</td><td>0.05</td><td>-4.6719</td><td>WASP-31</td><td>0.043</td><td>-4.8792</td></tr><tr><td>HAT-P-7</td><td>0.043</td><td>-5.0155</td><td>Kepler-5</td><td>0.06</td><td>-4.7970</td><td>WASP-32</td><td>0.050</td><td>-4.6694</td></tr><tr><td>HAT-P-8</td><td>0.048</td><td>-4.8982</td><td>Kepler-6</td><td>0.06</td><td>-4.7205</td><td>WASP-35</td><td>0.045</td><td>-4.7983</td></tr><tr><td>HAT-P-9</td><td>0.042</td><td>-4.9961</td><td>Kepler-7</td><td>0.06</td><td>-4.7558</td><td>WASP-38</td><td>0.047</td><td>-4.8540</td></tr><tr><td>HAT-P-11</td><td>0.053</td><td>-4.2487</td><td>Kepler-8</td><td>0.06</td><td>-4.8347</td><td>WASP-39</td><td>0.055</td><td>-4.7400</td></tr><tr><td>HAT-P-12</td><td>0.053</td><td>-4.4853</td><td>Kepler-10</td><td>0.05</td><td>-4.7575</td><td>WASP-41</td><td>0.057</td><td>-4.4184</td></tr><tr><td>HAT-P-13</td><td>0.042</td><td>-4.8525</td><td>Kepler-16</td><td>0.06</td><td>-4.6955</td><td>WASP-42</td><td>0.061</td><td>-4.7042</td></tr><tr><td>HAT-P-14</td><td>0.048</td><td>-4.7672</td><td>Kepler-17</td><td>0.07</td><td>-4.0817</td><td>WASP-43</td><td>0.050</td><td>-4.0022</td></tr><tr><td>HAT-P-15</td><td>0.077</td><td>-4.6795</td><td>Kepler-19</td><td>0.05</td><td>-4.7635</td><td>WASP-48</td><td>0.051</td><td>-4.6730</td></tr><tr><td>HAT-P-16</td><td>0.049</td><td>-4.7475</td><td>Kepler-20</td><td>0.05</td><td>-4.6767</td><td>WASP-49</td><td>0.046</td><td>-4.8127</td></tr><tr><td>HAT-P-17</td><td>0.050</td><td>-4.8463</td><td>Kepler-22</td><td>0.05</td><td>-4.8650</td><td>WASP-50</td><td>0.047</td><td>-4.8038</td></tr><tr><td>HAT-P-18</td><td>0.059</td><td>-4.4692</td><td><b>Kepler-25</b></td><td>0.04</td><td>-5.1557</td><td>WASP-52</td><td>0.059</td><td>-4.0728</td></tr><tr><td>HAT-P-20</td><td>0.050</td><td>-4.0156</td><td>Kepler-48</td><td>0.06</td><td>-4.5892</td><td>WASP-54</td><td>0.042</td><td>-4.9173</td></tr><tr><td>HAT-P-21</td><td>0.051</td><td>-4.2574</td><td>Kepler-62</td><td>0.06</td><td>-4.6710</td><td>WASP-57</td><td>0.061</td><td>-4.5087</td></tr><tr><td>HAT-P-22</td><td>0.042</td><td>-4.7770</td><td>Kepler-68</td><td>0.04</td><td>-4.9659</td><td>WASP-58</td><td>0.057</td><td>-4.1409</td></tr><tr><td>HAT-P-24</td><td>0.046</td><td>-4.9502</td><td>Kepler-78</td><td>0.06</td><td>-4.1905</td><td>WASP-59</td><td>0.059</td><td>-4.2904</td></tr><tr><td>HAT-P-25</td><td>0.071</td><td>-4.7852</td><td>Kepler-93</td><td>0.05</td><td>-4.7541</td><td>WASP-60</td><td>0.055</td><td>-4.8047</td></tr><tr><td>HAT-P-26</td><td>0.051</td><td>-4.8018</td><td>QATAR-1</td><td>0.06</td><td>-4.1587</td><td>WASP-69</td><td>0.053</td><td>-4.1095</td></tr><tr><td>HAT-P-27</td><td>0.056</td><td>-4.5616</td><td>QATAR-2</td><td>0.06</td><td>-4.0035</td><td>WASP-70A</td><td>0.043</td><td>-5.0399</td></tr><tr><td>HAT-P-28</td><td>0.057</td><td>-4.7338</td><td>TRES-1</td><td>0.06</td><td>-4.5441</td><td><b>WASP-72</b></td><td>0.036</td><td>-5.1549</td></tr><tr><td>HAT-P-29</td><td>0.061</td><td>-4.7664</td><td>TRES-2</td><td>0.05</td><td>-4.8071</td><td>WASP-76</td><td>0.042</td><td>-4.9288</td></tr><tr><td>HAT-P-30</td><td>0.035</td><td>-4.9826</td><td>Tres-3</td><td>0.06</td><td>-4.3001</td><td>WASP-80</td><td>0.051</td><td>-4.1197</td></tr><tr><td>HAT-P-31</td><td>0.046</td><td>-4.9776</td><td>TRES-4</td><td>0.04</td><td>-4.9925</td><td>WASP-84</td><td>0.048</td><td>-4.1695</td></tr><tr><td>HAT-P-32</td><td>0.065</td><td>-4.3822</td><td>WASP-1</td><td>0.05</td><td>-4.9381</td><td>WASP-103</td><td>0.066</td><td>-4.3406</td></tr><tr><td>HAT-P-33</td><td>0.046</td><td>-4.7034</td><td>WASP-2</td><td>0.06</td><td>-4.8529</td><td>WASP-106</td><td>0.046</td><td>-4.6278</td></tr><tr><td>HAT-P-34</td><td>0.053</td><td>-4.6819</td><td>WASP-3</td><td>0.06</td><td>-4.9245</td><td>WASP-107</td><td>0.050</td><td>-4.0414</td></tr><tr><td><b>HAT-P-35</b></td><td>0.033</td><td>-5.1588</td><td>WASP-4</td><td>0.05</td><td>-4.6329</td><td>WASP-117</td><td>0.040</td><td>-4.7842</td></tr><tr><td>HAT-P-36</td><td>0.050</td><td>-4.3864</td><td>WASP-5</td><td>0.05</td><td>-4.4888</td><td>WASP-127</td><td>0.041</td><td>-4.7504</td></tr><tr><td>HAT-P-38</td><td>0.060</td><td>-4.9221</td><td>WASP-7</td><td>0.06</td><td>-3.9172</td><td>XO-1</td><td>0.051</td><td>-4.6576</td></tr><tr><td>HAT-P-39</td><td>0.054</td><td>-4.6221</td><td>WASP-10</td><td>0.06</td><td>-5.0311</td><td>XO-2N</td><td>0.050</td><td>-4.6898</td></tr><tr><td>HAT-P-40</td><td>0.057</td><td>-4.8252</td><td>WASP-11</td><td>0.07</td><td>-4.4940</td><td>XO-3</td><td>0.074</td><td>-4.3384</td></tr><tr><td>HAT-P-41</td><td>0.055</td><td>-4.9518</td><td><b>WASP-12</b></td><td>0.04</td><td>-5.3203</td><td>XO-4</td><td>0.053</td><td>-4.9380</td></tr><tr><td>HAT-P-44</td><td>0.043</td><td>-5.0956</td><td>WASP-13</td><td>0.04</td><td>-4.8916</td><td></td><td></td><td></td></tr></tbody></table>

## Appendix B: Comparison between CA index distributions of field and transit host stars

The comparison of the four statistical moments of each distribution (see Table B.1) reveals that all are leptokurtic (kurt(X) > 0), and asymmetric distributions towards the higher values of $\log(R'_{HK})$ (skew(X)>0). In this contest, the host stars are more dispersed with respect to field stars, including both M-S and evolved field stars ($\sigma_{host} > \sigma_{all} > \sigma_{M-S} > \sigma_{sub-giant}$). On the contrary, the host stars distribution is more symmetric with respect to the others (kurt$_{host} < kurt_{M-S} < kurt_{all} < kurt_{sub-giant}$).

<table><caption>**Table B.1.** Comparison of the values of the distribution moments of the different populations of objects.</caption><thead><tr><th></th><th>Hosts</th><th>Field M-S stars</th><th>Field Sub-Giants</th><th>Field stars (all)</th></tr></thead><tbody><tr><td>Mean (μ)</td><td>-4.7429</td><td>-4.5862</td><td>-4.7307</td><td>-4.6231</td></tr><tr><td>Variance (σ)</td><td>0.088</td><td>0.063</td><td>0.055</td><td>0.065</td></tr><tr><td>Skewness (skew(X))</td><td>0.015</td><td>0.605</td><td>3.739</td><td>0.790</td></tr><tr><td>Kurtosis (kurt(X))</td><td>0.626</td><td>0.542</td><td>1.577</td><td>0.724</td></tr></tbody></table>
A Kolmogorov-Smirnov test yields very low probabilities for the null hypothesis: the probability that the CA distributions for M-S field stars and planet hosts are drawn from the same parent population (Figure 12 top-right panel) is Prob = 1.29 × 10⁻¹³ (D=0.35, N_e = 120), while for evolved stars (Figure 12 bottom-right panel) we found Prob = 4.4 × 10⁻² (D=0.1, N_e = 103). Considering instead the whole distributions for field stars and planet hosts, we obtained Prob = 4.14 × 10⁻⁹ (D=0.3, N_e = 122) that they are drawn from the same parent population. It is worth noting that the field stars' bimodality indicated by the different average values of log(R'_{HK}) for M-S (−4.5862) and sub-giant (−4.7307), is expected (Wright 2004; Gray et al. 2006). In conclusion, the M-S and evolved stars' CA distribution comes from a different parent distribution of the transit host stars.

## Appendix C: Details of machine learning analysis

<figure><img src="image_17.png" /><figcaption>**Fig. C.1.** Training clustering provided to the algorithm in the log(R'_{HK}) − q_p⁻¹ space. Black points were not assigned to any cluster and were used as test data.</figcaption></figure>

In machine learning investigations, data points are represented by vectors containing different “features”, and classification methods divide *n*-dimensional samples into clusters according to different possible metrics. Another fundamental distinction is between unsupervised and supervised methods: the former do not rely on prior, or “training”, information, contrarily to the latter. Following Fossati et al. (2015), we first attempted a classification without prior information and assumed our data set could be represented by a finite number of Gaussian distributions. As shown by Daniëlski et al. (2022), this can be done with a Gaussian Mixture Model (GMM, Sorenson & Alspach 1971; Alspach & Sorenson 1972).

Before running the classification, we prepared our data set by normalizing it, in order to assign equal weight to each feature. This is usually done by removing the data mean and by scaling it to unit variance:

$$Z = \frac{X - \mu}{\sigma}, \qquad (C.1)$$

where Z is the transformed data set, X is the raw data set, and μ and σ are its mean and standard deviation across every feature or parameter.

Then we adopted the GMM implementation in SCIKIT-LEARN (Pedregosa et al. 2011) and tested it by imposing an increasing number of sub-clusters. Even if the GMM was able to reproduce Fossati et al. (2015)'s result, we noticed its tendency to prefer three mixtures on our larger data set, where the entire sample of Fossati et al. (2015) is included in a single cluster. This can probably be explained by the different way clusters are formed by a GMM and by the cluster-weighted model adopted in their study. While in the former case, the data points are assumed to populate the log(R'_{HK}) − q_p⁻¹ parameter space according to a number of multivariate Gaussian distributions, each with a center and a covariance matrix, the latter associates them to linear relationships following Lanza (2014)'s model.

We then preferred to use a supervised learning algorithm. After testing several algorithms, we chose a Decision Tree approach, which learns information from data using simple decision rules (e.g. Geron 2017). Once again, we relied on SCIKIT-LEARN's implementation of this algorithm and chose a maximum depth of 4, that is, the one that allowed us to reproduce a close clustering to the one of Fossati et al. (2015)'s on their data.

We trained the algorithm after visual inspection on how four clusters might be identified in the data set (see Sect. 8 in the main text). We assigned 94% of the data point to different clusters, as shown by the colours used in Figure C.1, leaving out a few points at the edge of the clusters. A cross-validation test (e.g. Geron 2017) yielded an accuracy of (91 ± 6)% for our combination of algorithms and training data. Then, we ran the Decision Tree on the whole data set and obtained the classification shown in Figure 15 in the main text. We remark that we could not use a classical “score test” to measure the accuracy of the prediction, as we worked in a framework of multi-label, and not binary, classification (see SCIKIT-LEARN's documentation of the Decision Tree algorithm).