# The GAPS programme at TNG. XL. A puffy and warm Neptune-sized planet and an outer Neptune-mass candidate orbiting the solar-type star TOI-1422

L. Naponiello, L. Mancini, M. Damasso, A. S. Bonomo, A. Sozzetti, D.

Nardiello, K. Biazzo, R. G. Stognone, J. Lillo-Box, A. F. Lanza, et al.

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

## ► To cite this version:

L. Naponiello, L. Mancini, M. Damasso, A. S. Bonomo, A. Sozzetti, et al.. The GAPS programme at TNG. XL. A puffy and warm Neptune-sized planet and an outer Neptune-mass candidate orbiting the solar-type star TOI-1422. Astronomy and Astrophysics - A&A, 2022, 667, 10.1051/0004-6361/202244079. insu-03860268

HAL Id: insu-03860268

https://insu.hal.science/insu-03860268

Submitted on 19 Nov 2022

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

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

## XL. A puffy and warm Neptune-sized planet and an outer Neptune-mass candidate orbiting the solar-type star TOI-1422

L. Naponiello<sup>1,2</sup>, L. Mancini<sup>1,3,4</sup>, M. Damasso<sup>3</sup>, A. S. Bonomo<sup>3</sup>, A. Sozzetti<sup>3</sup>, D. Nardiello<sup>6,7</sup>, K. Biazzo<sup>5</sup>, R. G. Stognone<sup>8</sup>, J. Lillo-Box<sup>9</sup>, A. F. Lanza<sup>14</sup>, E. Poretti<sup>18,19</sup>, J. J. Lissauer<sup>35,36</sup>, L. Zeng<sup>10,11</sup>, A. Bieryla<sup>32</sup>, G. Hébrard<sup>21,22</sup>, M. Basilicata<sup>1</sup>, S. Benatti<sup>15</sup>, A. Bignamini<sup>12</sup>, F. Borsa<sup>18</sup>, R. Claudi<sup>7</sup>, R. Cosentino<sup>19</sup>, E. Covino<sup>13</sup>, A. de Gurtubai<sup>19</sup>, X. Delfosse<sup>24</sup>, S. Desidera<sup>7</sup>, D. Dragomir<sup>33</sup>, J. D. Eastman<sup>32</sup>, Z. Essack<sup>30,31</sup>, A. F. M. Fiorenzano<sup>19</sup>, P. Giacobbe<sup>3</sup>, A. Harutyunyan<sup>19</sup>, N. Heidari<sup>25,26,23</sup>, C. Hellier<sup>39</sup>, J. M. Jenkins<sup>34</sup>, C. Knapic<sup>12</sup>, P.-C. König<sup>20,21</sup>, D. W. Latham<sup>32</sup>, A. Magazzù<sup>37</sup>, A. Maggio<sup>15</sup>, J. Maldonado<sup>15</sup>, G. Micela<sup>15</sup>, E. Molinari<sup>16</sup>, M. Molinaro<sup>12</sup>, E. H. Morgan<sup>28</sup>, C. Moutou<sup>27</sup>, V. Nascimbeni<sup>7</sup>, E. Pace<sup>2</sup>, I. Pagano<sup>14</sup>, M. Pedani<sup>19</sup>, G. Piotto<sup>17</sup>, M. Pinamonti<sup>3</sup>, E. V. Quintana<sup>29</sup>, M. Rainer<sup>18</sup>, G. R. Ricker<sup>31</sup>, S. Seager<sup>28,30,40</sup>, J. D. Twicken<sup>34,35</sup>, R. Vanderspek<sup>28</sup>, and J. N. Winn<sup>38</sup>

*(Affiliations can be found after the references)*

Received 21 May 2022 / Accepted 7 July 2022

## ABSTRACT

**Context.** Neptunes represent one of the main types of exoplanets and have chemical-physical characteristics halfway between rocky and gas giant planets. Therefore, their characterization is important for understanding and constraining both the formation mechanisms and the evolution patterns of planets.

**Aims.** We investigate the exoplanet candidate TOI-1422 b, which was discovered by the TESS space telescope around the high proper-motion G2 V star TOI-1422 (V = 10.6 mag), 155 pc away, with the primary goal of confirming its planetary nature and characterising its properties.

**Methods.** We monitored TOI-1422 with the HARPS-N spectrograph for 1.5 yr to precisely quantify its radial velocity (RV) variation. We analyse these RV measurements jointly with TESS photometry and check for blended companions through high-spatial resolution images using the AstraLux instrument.

**Results.** We estimate that the parent star has a radius of $R_* = 1.019_{-0.013}^{+0.014} R_\odot$, and a mass of $M_* = 0.981_{-0.065}^{+0.062} M_\odot$. Our analysis confirms the planetary nature of TOI-1422 b and also suggests the presence of a Neptune-mass planet on a more distant orbit, the candidate TOI-1422 c, which is not detected in TESS light curves. The inner planet, TOI-1422 b, orbits on a period of $P_b = 12.9972 \pm 0.0006$ days and has an equilibrium temperature of $T_{eq,b} = 867 \pm 17$ K. With a radius of $R_b = 3.96_{-0.11}^{+0.13} R_\odot$, a mass of $M_b = 9.0_{-2.0}^{+2.3} M_\odot$ and, consequently, a density of $\rho_b = 0.795_{-0.235}^{+0.290} \text{ g cm}^{-3}$, it can be considered a warm Neptune-sized planet. Compared to other exoplanets of a similar mass range, TOI-1422 b is among the most inflated, and we expect this planet to have an extensive gaseous envelope that surrounds a core with a mass fraction around 10% - 25% of the total mass of the planet. The outer non-transiting planet candidate, TOI-1422 c, has an orbital period of $P_c = 29.29_{-0.20}^{+0.21}$ days, a minimum mass, $M_c \sin i$, of $11.1_{-2.3}^{+2.6} M_\odot$, an equilibrium temperature of $T_{eq,c} = 661 \pm 13$ K and, therefore, if confirmed, could be considered as another warm Neptune.

**Key words.** techniques: photometric – planetary systems – techniques: spectroscopic – techniques: radial velocities – stars; individual: TOI-1422 – methods: data analysis

## 1. Introduction

Exoplanetary science has expanded quickly from the simple detection of new worlds to their in-depth characterization. This characterization is especially feasible for planets orbiting bright stars on a plane almost aligned to our line of sight, meaning that their radius and mass can be derived by transit photometry and radial velocity (RV) measurements, respectively. The population of known transiting planets has increased significantly in the last two decades, mainly thanks to dedicated ground-based surveys, which were then followed by surveys from space that turned out to be much more efficient, considering the total number of discoveries.

Thus far, the *Kepler* and the K2 space missions (Borucki et al. 2010; Howell et al. 2014) have had a very important impact

on the exoplanet field by discovering thousands of confirmed and candidate planets, many of which are not amenable to RV follow-up due to the faintness of their host stars. The Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2014), currently at the end of its first extended mission and with a second one already proposed, was designed to target nearby and bright stars over a large portion of the sky (around 85% sky coverage during the primary mission alone) because such stars are easier to follow up by means of RV, and result in refined measurements of their own exoplanet masses, atmospheres, sizes, and therefore, densities. The opportunity to use a large exoplanet sample such as that of *Kepler*, which is based on homogeneous data and has minimal pollution from false positives (< 10%, Fressin et al. 2013), has allowed us to distinguish between several distinct exoplanet regimes (Weiss & Marcy 2014; Buchhave et al. 2014;
Zeng et al. 2019): the terrestrial-like planets ($R_p < 1.7 R_⊕$), the gas dwarf planets with rocky cores and hydrogen–helium envelopes, the H$_2$O-dominated ices and fluid water worlds (both of the latter two classes have $1.7 R_⊕ < R_p < 3.9 R_⊕$) and the ice or gas giant planets ($R_p > 3.9 R_⊕$).

Planet occurrence around main-sequence stars has been investigated thanks to Doppler surveys (e.g. Cumming et al. 2008; Wright et al. 2012). In particular, the Keck Eta-Earth survey (Howard et al. 2010) and the CORALIE+HARPS survey (Mayor et al. 2011) first explored the domain of low-mass (3–30 $M_⊕$) close-in ($P_{\text{orb}} \sim 50$ days) planets. These planets turned out to be an order of magnitude more common than giant planets.

Other studies for determining the occurrence rates of planets, based on the *Kepler* sample, agree that for planets with less than a 1-yr orbital period, their mean number per star is higher within the radius range $1 R_⊕ < R_p < 4 R_⊕$ rather than the range $4 R_⊕ < R_p < 16 R_⊕$ (Howard et al. 2012; Fressin et al. 2013; Petigura et al. 2013). The subsequent and gradual refinement of parent-star properties (especially thanks to high-resolution stellar spectra) revealed a clear bimodality of the radius distribution of close-in ($P < 100$ days), small-sized ($R_p < 4.0 R_⊕$) planets orbiting bright, main-sequence solar-type stars (Petigura et al. 2017; Fulton & Petigura 2018; Van Eylen et al. 2018)<fnref n="1" />. These two quite distinct populations were identified as ‘super-Earths’ ($R_p < 1.5 R_⊕$) and ‘sub-Neptunes’ ($R_p = 2–3 R_⊕$), which are also represented in the intermediate region ($R_p = 1.5–2 R_⊕$) with fewer planets. However, it is better to stress that, since we do not know for sure what they are made of, the space of physical parameters ($R_p, M_p$), for which the previous terms apply, are not strictly defined.

The advantage of studying transiting planets is the possibility, in many cases, to measure both the planetary radius and mass, and therefore determine their density and bulk composition. Knowing the structural properties, one should be able to distinguish among the various scenarios of exoplanet formation and evolution. Unfortunately, theoretical models (e.g. Bitsch et al. 2019; Turbet et al. 2020) tell us that the mass-radius relationships for small planets present degeneracy due to the vastness of possible different compositions and amounts of rock, ice, and gas, especially in the transition between rocky super-Earths and Neptune-like planets (e.g. Miller-Ricci et al. 2009; Lozovsky et al. 2018). A detailed investigation of the mass-radius relation for small planets can be useful for throwing light on several open questions, such as the diversity of planet core masses and compositions, or where they form (in situ or beyond the snowline), and the existence of the radius gap. We refer the reader to the recent review by Biazzo et al. (2022) for an exhaustive discussion on this topic.

It is therefore clear how RV follow-up observations and planetary-mass measurements play an important role in understanding this process and why there is currently a tremendous effort in this field by many teams (e.g. KESPRINT: Gandolfi et al. 2018; HARPS-N consortium: Cloutier et al. 2020; NCORES: Armstrong et al. 2020; TESS-Keck Survey: Chontos et al. 2022; GAPS: Carleo et al. 2021) to confirm TESS small-planet candidates.

Probing the chemical composition of the atmosphere of a large number of sub-Neptune planets would also be helpful to unravel the skein. Various techniques (such as high-resolution spectroscopy, transmission, and emission spectroscopy) have

been implemented and applied successfully using the *Hubble* Space Telescope (HST) instruments or the high-resolution spectrographs mounted on large-class ground-based telescopes (e.g. CRIRES: Snellen et al. 2010; HARPS: Wyttenbach et al. 2015; LDSS3C: Diamond-Lowe et al. 2018; GIANO: Brogi et al. 2018; CARMENES: Casasayas-Barris et al. 2019; HARPS-N: Pino et al. 2020; ESPRESSO: Borsa et al. 2021). Unfortunately, these techniques for probing the planetary atmospheres are currently effectively applicable only to giant planets, as we know only a few sub-Neptune planets for which the transmission-spectrum signal can be detected with a sufficient signal-to-noise ratio (S/N) that allows us to discriminate between different atmospheric models. The featureless transmission spectra of GJ 436 b (Knutson et al. 2014) and GJ 1214 b (Kreidberg et al. 2014) are emblematic.

The situation should improve soon thanks to the *James Webb* Space Telescope (JWST; Barstow et al. 2015), which is about to go into operation, and with the next generation of space-based and large ground-based telescopes (Ariel: Tinetti et al. 2021; ELT: Ramsay et al. 2021; TMT: Skidmore et al. 2015). In the meantime, it is important that we continue to work to uncover new exoplanets, especially those of small size ($R_p < 5 R_⊕$) that orbit bright ($V < 11$ mag) main-sequence dwarf stars. This is currently possible thanks to the large number of planet candidates (more than 5000) that TESS is discovering at the present time. The recent detection of water vapour in the atmosphere of the super-Neptune TOI-674 b with the HST (Brande et al. 2022) is a successful example of this effort.

On the 6 November 2019, the TESS target star TIC 333473672 was officially named TOI-1422 (TESS Object of Interest; Guerrero et al. 2021), following the Data Validation Report Summary (DRS) produced by the TESS Science Processing Operations Center (SPOC; Jenkins et al. 2016) pipeline at the NASA Ames Research Center through the Transiting Planet Search (TPS; Jenkins 2002; Jenkins et al. 2010, 2020) and Data Validation (DV; Twicken et al. 2018, Li et al. 2019) modules. In particular, TOI-1422 01 was flagged as a potential planet with an orbital period of 13.0020 ± 0.0040 days, a transit depth of 1422 ± 94 ppm (parts per million), and a corresponding radius of $3.85 ± 0.90 R_⊕$, which is compatible with Neptune’s radius. The candidate passed all SPOC DV diagnostic tests, and, furthermore, all TIC (version 8) objects other than the target star were excluded as sources of the transit signal through the difference image centroid offsets (Twicken et al. 2018).

The long-term, multi-programme Global Architecture of Planetary Systems (GAPS; Covino et al. 2013; Poretti et al. 2016) exploits Doppler measurements taken with the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N; Cosentino et al. 2012) instrument at the Telescopio Nazionale Galileo (TNG) in La Palma (Spain). This high-resolution spectrograph (resolving power $R \approx 115\,000$) delivers the highest RV precision (~1 m s⁻<fnref n="1" />) currently achievable in the northern hemisphere. One of the aims of the GAPS programme is to confirm and obtain an accurate mass determination of planets having an intermediate-mass between super-Earths and super-Neptunes; for this reason, TOI-1422 was selected for RV follow-up observations, which started in June 2020.

In the present work, we report the results of our measurements and analyses that allowed us to confirm TOI-1422 as a new planetary system. The paper is organized as follows: Sect. 2 contains the details of the instruments, and the photometric and RV measurements; the results of our analyses are presented in Sect. 3

<footnote n="1">¹ For a possible explanation of Fulton’s gap, see Modirrousta-Galian et al. (2020).</footnote>
and discussed in Sect. 4; and we finally address the conclusions in Sect. 5.

## 2. Observations and data reduction

### 2.1. TESS photometry

Since late July 2018, TESS has observed more than 200 000 stars with its four wide-field optical Charged-Coupled Devices (CCD) cameras (24 × 96 degrees), each having a focal ratio of $f/1.4$ and a broad-band filter range between 600 and 1000 nm. The pre-selected target TIC 333473672 was observed in Sectors 16 and 17 between 11 September 2019 and 2 November 2019, and the first of a total of four transiting events were recorded on 19 September 2019. The two-minute cadence photometry of TOI-1422 from TESS spans a total of ≈ 50 days and to analyse it, we used the Presearch Data Conditioning Simple Aperture Photometry (PDC-SAP; Stumpe et al. 2012, 2014, Smith et al. 2012) light curve, which is provided by the TESS SPOC pipeline and retrieved through the Python package `lightkurve` (Lightkurve Collaboration 2018) from the Mikulski Archive for Space Telescopes (MAST). We jointly fitted the transit model and a Gaussian process (GP) using a simple (approximate) Matern kernel, which was implemented in the Python modelling tool `juliet<fnref n="2" />` (Espinoza et al. 2019) via `celerite` (Foreman-Mackey et al. 2017), of the form:

$$k(\tau_{i,j}) = \sigma_{GP}^2 M(\tau_{i,j}, \rho) + (\sigma_i^2 + \sigma_w^2) \delta_{i,j}, \quad (1)$$

where $\sigma_i$ is the error bar of the i-th data point, $\sigma_{GP}$ the amplitude of the GP in parts per million (ppm), $\sigma_w$ an added jitter term (in ppm), $\delta_{i,j}$ the Kronecker's delta, $k(\tau_{i,j})$ the element $i,j$ of the covariance matrix as a function of $\tau_{i,j} = |t_i - t_j|$, with $t_i$ and $t_j$ being the $i,j$ GP regressors (i.e. the observing times), while

$$M(\tau_{i,j}, \rho) = (1 + 1/\epsilon) e^{-[1-\epsilon]\sqrt{3\tau_{i,j}/\rho}} + (1 - 1/\epsilon) e^{-[1+\epsilon]\sqrt{3\tau_{i,j}/\rho}} \quad (2)$$

is the kernel with its characteristic time scale $\rho$. The parameter $\epsilon$ controls the quality of the approximation since, in the limit $\epsilon \to 0$, Eq. (2) becomes the Matern-3/2 function. In `juliet`, the possible polluting sources inside the TESS aperture<fnref n="3" /> (Fig. 1), which might result in a smaller transit depth compared to the real one, are taken into account with a dilution factor ($D$) that, in this case, has been neglected because the PDC-SAP photometry is already corrected for dilution from other objects contained within the aperture using the Create Optimal Apertures (COA) module (Bryson et al. 2010, 2020)<fnref n="4" />. In order to efficiently sample the whole plausible zone in the $(b, k)$ plane, where $b$ is the impact parameter and $k$ is the planet-to-star radius ratio, we used the $(r_1, r_2)$ parametrization described in Espinoza (2018). This is the same approach that we adopted for the modelling of the transits in the joint analysis with the RVs (see Sect. 3.3). Moreover, here we make use of the limb-darkening parametrizations of Kipping (2013) for two-parameter limb-darkening laws $(q_1, q_2 \to u_1, u_2)$.

The PDC-SAP light curve of TOI-1422 and its detrending are plotted in Fig. 2. We also analysed the SPOC SAP photometry (Twicken et al. 2010; Morris et al. 2020), which presents a small

<figure><img src="image_3.png" /><figcaption>Fig. 1. Target pixel file from the TESS observation of Sector 16, made with `tpfplotter` (Aller et al. 2020) and centred on TOI-1422, which is marked with a white cross. The SPOC pipeline aperture is shown by shaded red squares, and the Gaia satellite eDR3 catalogue (Brown et al. 2018; Prusti et al. 2016) is also overlaid with symbol sizes proportional to the magnitude difference with TOI-1422. The difference image centroid locates the source of the transits within $1.89 \pm 5$ arcsec of the target star's location, as reported by the TicOffset for the multi-sector DV report for this system.</figcaption></figure>

long-term variability that might be due to systematics, but no other feature or modulation can be discerned within the experimental uncertainties, aside from a possible single extra transit event, which is discussed at the end of Sect. 3.4, and a steep flux drop at the end of both the SAP and PDC-SAP light curves, which are probably due to high levels of background noise.

### 2.2. High-spatial resolution imaging – AstraLux

We observed TOI-1422 with the AstraLux high-spatial resolution camera (Hormuth et al. 2008), located at the 2.2 m telescope of the Calar Alto Observatory (CAHA, Almería, Spain) using the lucky-imaging technique. This technique obtains diffraction-limited images by acquiring thousands of short-exposure frames and selecting those with the highest Strehl ratio (Strehl 1902) to finally combine them into a co-added high-spatial resolution image. We observed this target on the night of 29 September 2021 under good weather conditions with a mean seeing of 1 arcsec, and obtained 50 000 frames with 20 ms exposure time in the Sloan Digital Sky Survey z filter (SDSSz), with a field of view windowed to $6 \times 6$ arcsec. The datacube was reduced by the instrument pipeline (Hormuth et al. 2008) and we selected the best quality 10% frames to produce the final high-resolution image. We obtained the sensitivity limits of the co-added image by using our own developed ASTRASENS package<fnref n="5" /> with the procedure described in Lillo-Box et al. (2012, 2014). The $5\sigma$ sensitivity curve is shown in Fig. 3. We could discard sources down to 0.2 arcsec with a magnitude contrast of $\Delta Z < 4$ mag, corresponding to a maximum contamination level of 2.5%. By using this high-spatial resolution image, we also estimated the probability of an undetected blended source. This probability (fully described in Lillo-Box et al. 2014) is called the blended source confidence (BSC). We used a python implementation of this approach (bsc, by J. Lillo-Box), which uses

<footnote n="2">² https://juliet.readthedocs.io</footnote>

<footnote n="3">³ `tpfplotter` is a python package developed by J. Lillo-Box and publicly available on www.github.com/jlillo/tpfplotter.</footnote>

<footnote n="4">⁴ Since the release of the light curve products from Year 2, the SPOC background estimation algorithm has been updated due to an over-correction bias, which was significant for dim and/or crowded targets. For this particular TOI, we estimated this over-correction to be negligible for the planetary radius estimation as it is significantly smaller than the transit depth uncertainty.</footnote>

<footnote n="5">⁵ https://github.com/jlillo/astrasens</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Light curve of TOI-1422 as collected by TESS in Sectors 16 and 17 with a 2-min cadence. **Top panel:** light curve from the PDC-SAP pipeline. The black line represents the best-fit model obtained through GP detrending, as detailed in Sect. 2.1. **Bottom panel:** residuals of the best-fit model in parts per million.</figcaption></figure>

the TRILEGAL<fnref n="6" /> galactic model (v1.6; Girardi et al. 2012) to retrieve a simulated source population of the region around the corresponding target<fnref n="7" />. This simulated population was used to compute the density of stars around the target position (radius $r = 1°$) and derive the probability of chance alignment at a given contrast magnitude and separation. When applied to the TOI-1422 location, we used a maximum contrast magnitude of $\Delta m_{b,max} = 6.97$ mag in the SDSSz passband, corresponding to the maximum contrast of a blended eclipsing binary that could mimic the observed transit depth of planet b (~ 1000 ppm). Thanks to our high-resolution image, we estimated the probability of an undetected blended source to be 0.28%. The probability of such an undetected source being an appropriate eclipsing binary was thus even lower and consequently, we could assume that the transit signal was not due to a blended eclipsing binary.

## 2.3. HARPS-N radial velocities

Between June 2020 and January 2022, we collected a total of 112 RV measurements of TOI-1422 with HARPS-N (Table A.1). The RVs were calculated using the TERRA pipeline (Anglada-Escudé & Butler 2012), version 1.8, through the YABI workflow interface (Hunter et al. 2012), which is maintained by the Italian center for Astronomical Archive (IA2). TERRA is an algorithm based on the template matching technique, and is preferred for the RVs retrieval in this paper over the standard Data Reduction Software (DRS) pipeline, which returned a slightly lower overall RV precision<fnref n="8" /> on this target. We used the RVs calculated using all the spectral orders, and from the full sample of RVs, we removed four points following Chauvenet's criterion. TERRA RVs have an average measurement error of 2.6 m s⁻¹, a root mean square error of 4.5 m s⁻¹, and a S/N ≈ 35, measured at a reference wavelength of 5500 Å. A long linear trend is evident in HARPS-N RV data, as we discuss in Sect. 3.2.

TOI-1422 was also observed with the SOPHIE instrument, a stabilized échelle spectrograph mounted at the 193-cm Telescope of Observatoire de Haute-Provence in France (Perruchot et al. 2008, Bouchy, F. et al. 2013). However, for signals of low

<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Blended source confidence (BSC) curve from the AstraLux SDSSz image (solid black line). The colour on each angular separation and contrast bin represent the probability of a source aligned at the location of the target, based on the TRILEGAL model. The horizontal dotted line shows the maximum contrast of a blended binary that is capable of imitating the planet's transit depth. The green region represents the regime that is not explored by the high-spatial resolution image. The BSC curve corresponds to the integration of $P_{\text{aligned}}$ over this region.</figcaption></figure>

semi-amplitudes such as those we discuss in this work, the RV measurements of SOPHIE, due to higher uncertainties compared to HARPS-N, do not increase the significance of the results presented in Sect. 5, and therefore have not been utilized.

# 3. System characterization

## 3.1. Parent star

From the co-added spectrum built from individual HARPS-N spectra extracted with the standard DRS pipeline, we derived the following atmospheric parameters of the planet's host star TOI-1422: effective temperature $T_{\text{eff}}$, surface gravity $\log g$, microturbulence velocity $\xi$, iron abundance [Fe/H],

<footnote n="6">⁶ http://stev.oapd.inaf.it/cgi-bin/trilegal</footnote>

<footnote n="7">⁷ This is done in python by using the astrobase implementation by Bhatti et al. (2020).</footnote>

<footnote n="8">⁸ For a comparison of the performances of TERRA vs. DRS see Perger et al. (2017).</footnote>
and rotational velocity $v \sin i_{\star}$. For $T_{\text{eff}}$, log g, $\xi$, and [Fe/H], we applied a method based on equivalent widths of iron lines taken from Biazzo et al. (2015) and the spectral analysis package MOOG (Sneden 1973; version 2017). The Castelli & Kurucz (2003) grid of model atmospheres was adopted. $T_{\text{eff}}$ and $\xi$ were derived by imposing that the abundance of Fe I was not dependent on the line excitation potentials and the reduced equivalent widths (i.e. EW/λ), respectively, while log g was obtained by imposing the Fe I/Fe II ionization equilibrium condition. The $v \sin i_{\star}$ was measured with the same MOOG code, by applying the spectral synthesis of three regions around 5400, 6200, and 6700 Å, and adopting the same grid of model atmosphere after fixing the macroturbulence velocity to the value of 3.4 km s⁻¹ from the relationship by Doyle et al. (2014). From these results, the star can be classified as a G2 V dwarf with a low projected rotation velocity $v \sin i_{\star}$ of 1.9 ± 0.8 km s⁻¹, implying a maximum rotation period of $27_{-8}^{+19}$ d at $1\sigma$. Analogously, using an empirical relation based on the Full Width at Half Maximum (FWHM)<fnref n="9" /> derived by the HARPS-N DRS, we find $v \sin i_{\star} \sim 2.2$ km s⁻¹.

The field of TOI-1422 was also observed in 2004, 2006, and 2007 during the WASP transit-search survey (Pollacco et al. 2006). A total of 20 000 photometric data points were obtained by observing the field every ~15 min on clear nights, over spans of ~120 days in each year. We searched the data for any rotational modulation using the methods from Maxted et al. (2011) and found no significant periodicity between 1 and 100 days, with a 95%-confidence upper limit on the amplitude of 2 mmag. The TESS light curve shows no modulation either (Sect. 2.1), confirming that the star is rather magnetically quiet over a period of ~100 days.

Moreover, the spectrum of TOI-1422 clearly shows a lithium line at $\lambda = 6707.8$ Å. We therefore estimated the lithium abundance log A(Li)NLTE by measuring the lithium EW and considering our $T_{\text{eff}}$, log g, $\xi$, and [Fe/H] previously derived together with the NLTE corrections by Lind et al. (2009). The value of the lithium abundance is listed in Table 1 and its position in a log A(Li)-$T_{\text{eff}}$ diagram is compatible with the M67 open cluster advanced age (~4.5 Gyr; see Pasquini et al. 2009) in agreement with the star's low activity level. The physical parameters of TOI-1422 are also displayed in Table 1 and were determined with the EXOFASTv2 Bayesian code (Eastman 2017; Eastman et al. 2019), by fitting the stellar spectral energy distribution (SED) and by employing the MESA Isochrones and Stellar Tracks (Dotter 2016) to more precisely constrain the stellar mass. In addition, in the table we report the stellar magnitudes used for the SED modelling, while the SED best fit is shown in Fig. 4. Gaussian priors were imposed on the Gaia eDR3 parallax (Gaia Collaboration 2021) as well as on the $T_{\text{eff}}$ and [Fe/H], as derived above from the analysis of the HARPS-N spectra. An upper limit was set on the V-band extinction, $A_V$, from reddening maps (Schlegel et al. 1998; Schlafly & Finkbeiner 2011).

## 3.2. RV and activity indexes periodogram analysis

We computed the generalized Lomb-Scargle (GLS) periodogram for the HARPS-N RVs and different stellar activity indexes<fnref n="10" />

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Spectral energy distribution computed for TOI-1422, where the black curve is the most likely atmospheric stellar model and the blue dots correspond to the model fluxes over each passband. The horizontal and vertical red error bars represent, respectively, the effective width of the passbands and the reported photometric measurement uncertainties (refer to the magnitudes in Table 1).</figcaption></figure>

using the Python package `astropy` v4.3.1 (Price-Whelan et al. 2018). The periodogram of the RVs shows the main peak around 29 days, and a significant peak at 13 days (TOI-1422 b transiting period), after correcting for a linear trend of ~4 m s⁻¹ yr⁻¹, observed in HARPS-N data. No index shows signs of the 29-day periodicity, but a linear trend is also present in the FWHM and log $R'_{\text{HK}}$ (see Fig. D.1), with the former correlating the most with the RVs, unveiling a moderate Spearman coefficient (Spearman 1904) of 0.41 (p-value 0.01%). Therefore, in order to explain the nature of the main peak in the RVs, we present the GLS periodogram of these coefficients posterior to the removal of their linear trends in Fig. 5 (see Fig. D.2 for a closer look at the RVs panel), but again no trace of the 29-day signal is found. We also performed a GP regression analysis, using a quasi-periodic model, of the log $R'_{\text{HK}}$ index corrected for the linear trend over the time series, and find no evidence of any particular periodic modulation in the posterior distribution of the periodic time-scale hyper-parameter. In short, there is no evidence pointing to a specific periodic rotation of the star TOI-1422, other than the tentative estimation from $v \sin i_{\star}$.

A query from the Gaia eDR3 archive returns astrometric excess noise and renormalized unit weight error (RUWE) values of 80 µas and 1.09, respectively, for TOI-1422. Thus, the star is astrometrically quiet. The analysis of Sect. 2.2 rules out the existence of obvious sub-arcsec stellar companions, and no co-moving objects are present in Gaia eDR3 data in a 600 arcsec radius. The linear trend seen in the RV data, along with a few activity indexes, can therefore be explained by long star magnetic activity, rather than by the presence of a companion<fnref n="11" />.

also analysed the chromospheric log $R'_{\text{HK}}$ index, and additional activity diagnostics derived from the spectroscopic lines He I, Na I, Ca I, Hα06 and Hα16 as defined in the code ACTIN (https://github.com/gomesdasilva/ACTIN v.1.3.9, Gomes da Silva et al. 2018) which has been used for the calculation. In particular, the two H-alpha indices have 1.6 and 0.6 Å band-pass width, respectively.

<footnote n="9">⁹ This relation was calibrated using a set of well-aligned transiting exoplanet systems, for which we could infer $v \sin i_{\star}$ as equal to their equatorial velocities. We estimate the equatorial velocities from the stellar radii and rotational period, and correlated these values directly to the FWHMs.</footnote>

<footnote n="10">¹⁰ The FWHM and the Bisector inverse span (BIS) are calculated using the cross correlation function (CCF) derived by the DRS pipeline. We</footnote>

<footnote n="11">¹¹ In this case, at a projected separation of 0.1 arcsec (~15 au at the distance of TOI-1422), the lower limit of the AstraLux imaging data, a maximum RV slope of the magnitude measured in this work would be produced by a companion of ~30 M<sub>Jup</sub> (i.e. either a very low-mass star or a massive sub-stellar companion).</footnote>
<figure><img src="image_7.png" /><figcaption>Fig. 5. GLS periodogram of HARPS-N RVs and of various activity indexes specified in the labels, after the removal of a linear trend (Sect. 3.2). The main peak of the RV GLS periodogram and that of the TOI-1422 b period are highlighted with a green and a red dashed line, respectively. They do not overlap with any of the peaks from the indexes, which in general do not suggest any clear stellar rotation period. The period corresponding to the highest peak in the RV GLS periodogram, and its False Alarm Probability (FAP), are written on the top of the first panel, while the horizontal dashed lines remark the 10% and 1% confidence levels (evaluated with the bootstrap method), respectively. The three peaks surrounding the RVs main frequency can all be explained as aliases of the 29-day signal due to the two highest frequencies of the window function (190 and 390.3 days, as shown in Fig. D.2 and Fig. D.3).</figcaption></figure>

### 3.3. RV and photometry joint analysis

A joint transit and RV analysis was carried out with juliet, which employs different Python tools: batman<fnref n="12" /> (Kreidberg 2015) for the modelling of transits, RadVel<fnref n="13" /> (Fulton et al. 2018) for the modelling of RVs, and stochastic processes, which are

treated as GPs with the packages george<fnref n="14" /> (Ambikasaran et al. 2015) and celerite<fnref n="15" />. The RV model that we used in juliet is the following:

$$M(t) = K(t) + \epsilon(t) + \bar{\mu} + A t + B, \quad (3)$$

<footnote n="12">¹² https://github.com/lkreidberg/batman</footnote>

<footnote n="13">¹³ https://radvel.readthedocs.io</footnote>

<footnote n="14">¹⁴ https://george.readthedocs.io</footnote>

<footnote n="15">¹⁵ https://celerite.readthedocs.io</footnote>
<table><caption>**Table 1.** TOI-1422 parameters.</caption><thead><tr><th>Parameter</th><th>Unit</th><th>Value</th><th>Source</th></tr></thead><tbody><tr><td colspan="4">Cross-identifications</td></tr><tr><td>TOI</td><td>...</td><td>TOI-1422</td><td>TOI catalogue</td></tr><tr><td>TIC ID</td><td>...</td><td>333473672</td><td>TIC</td></tr><tr><td>Tycho</td><td>...</td><td>3235-00524-1</td><td>Tycho</td></tr><tr><td>2MASS ID</td><td>...</td><td>J23365789+3938218</td><td>2MASS</td></tr><tr><td>Gaia ID</td><td>...</td><td>1920333449169516288</td><td>Gaia eDR3</td></tr><tr><td colspan="4">Astrometric properties</td></tr><tr><td>R.A.</td><td>J2016</td><td>354.240817</td><td>Gaia eDR3</td></tr><tr><td>Dec</td><td>J2016</td><td>+39.639275</td><td>Gaia eDR3</td></tr><tr><td>Parallax</td><td>mas</td><td>6.4418 ± 0.0138</td><td>Gaia eDR3</td></tr><tr><td>μ<sub>α</sub></td><td>mas yr<sup>-1</sup></td><td>-67.564 ± 0.015</td><td>Gaia eDR3</td></tr><tr><td>μ<sub>δ</sub></td><td>mas yr<sup>-1</sup></td><td>-31.180 ± 0.011</td><td>Gaia eDR3</td></tr><tr><td>Distance</td><td>pc</td><td>154.56<sup>+0.037</sup><sub>-0.027</sub></td><td>VizieR</td></tr><tr><td colspan="4">Photometric properties</td></tr><tr><td>B<sub>T</sub></td><td>mag</td><td>11.31 ± 0.07</td><td>Tycho</td></tr><tr><td>V<sub>T</sub></td><td>mag</td><td>10.62 ± 0.05</td><td>Tycho</td></tr><tr><td>J</td><td>mag</td><td>9.585 ± 0.022</td><td>2MASS</td></tr><tr><td>H</td><td>mag</td><td>9.275 ± 0.030</td><td>2MASS</td></tr><tr><td>K<sub>S</sub></td><td>mag</td><td>9.190 ± 0.022</td><td>2MASS</td></tr><tr><td>i'</td><td>mag</td><td>10.311 ± 0.075</td><td>APASS</td></tr><tr><td>W1</td><td>mag</td><td>9.161 ± 0.023</td><td>AllWISE</td></tr><tr><td>W2</td><td>mag</td><td>9.201 ± 0.020</td><td>AllWISE</td></tr><tr><td>W3</td><td>mag</td><td>9.161 ± 0.033</td><td>AllWISE</td></tr><tr><td>A<sub>V</sub></td><td>mag</td><td>&lt; 0.077</td><td>This work</td></tr><tr><td colspan="4">Stellar parameters</td></tr><tr><td>L<sub>*</sub></td><td>L<sub>☉</sub></td><td>1.116 ± 0.037</td><td>This work</td></tr><tr><td>M<sub>*</sub></td><td>M<sub>☉</sub></td><td>0.981<sup>+0.062</sup><sub>-0.065</sub></td><td>This work</td></tr><tr><td>R<sub>*</sub></td><td>R<sub>☉</sub></td><td>1.019<sup>+0.014</sup><sub>-0.013</sub></td><td>This work</td></tr><tr><td>T<sub>eff</sub></td><td>K</td><td>5840 ± 62</td><td>This work</td></tr><tr><td>log g<sub>*</sub></td><td>cgs</td><td>4.41 ± 0.11</td><td>This work</td></tr><tr><td>ξ</td><td>km s<sup>-1</sup></td><td>0.89 ± 0.07</td><td>This work</td></tr><tr><td>[Fe/H]</td><td>dex</td><td>-0.09 ± 0.07</td><td>This work</td></tr><tr><td>Spectral type <sup>(a)</sup></td><td></td><td>G2 V</td><td>This work</td></tr><tr><td>ρ<sub>*</sub></td><td>g cm<sup>-3</sup></td><td>1.3 ± 0.1</td><td>This work</td></tr><tr><td>ν sin i<sub>*</sub></td><td>km s<sup>-1</sup></td><td>1.7 ± 0.4</td><td>This work</td></tr><tr><td>log A(Li)<sub>NLTE</sub></td><td></td><td>1.97 ± 0.05</td><td>This work</td></tr><tr><td>log R'<sub>HK</sub></td><td>dex</td><td>-4.95 ± 0.03</td><td>This work</td></tr><tr><td>Age</td><td>Gyr</td><td>5.1<sup>+3.9</sup><sub>-3.1</sub></td><td>This work</td></tr></tbody></table>

**Notes.** <sup>(a)</sup>Spectral type defined according to the stellar spectral classification of Gray & Corbally (2009).

**References.** TESS Primary Mission TOI catalogue (Guerrero et al. 2021); TIC (Stassun et al. 2018; Stassun et al. 2019); Tycho (Høg et al. 2000); 2MASS (Skrutskie et al. 2006); Gaia eDR3 (Brown et al. 2021); AllWISE (Cutri et al. 2021); APASS (Henden et al. 2015); VizieR Online Data catalogue (Bailer-Jones et al. 2021).

where $\epsilon(t)$ is a noise model for the HARPS-N instrument, here assumed to be white-Gaussian noise, in other words $\epsilon(t) \approx N(0, \sigma(t)^2 + \sigma_w^2)$, with $\sigma(t)^2$ being the formal uncertainty of the RV point at time $t$, $\sigma_w^2$ being an added jitter term, and $N(\mu, \sigma^2)$ denoting a normal distribution with mean $\mu$ and variance $\sigma^2$. $K(t)$ is the Keplerian model of the RV star perturbations due to the orbiting planet, $\bar{\mu}$ is the systemic velocity linked to the instrument, and the coefficients $A, B$ (also referred to as RV slope and RV intercept) represent an additional linear trend used for modelling non-Keplerian signals with a period longer than the observation span. For a total number of data points $N$, we assumed the model likelihood to follow the likelihood of an

<figure><img src="image_8.png" /><figcaption>**Fig. 6.** GLS periodogram of the transiting one-planet model RV residuals. The main peak is highlighted in red and corresponds to a period of 29.2 days, with a FAP of 0.45% (evaluated with the bootstrap method), while the horizontal dashed lines show the 10% and 1% confidence levels.</figcaption></figure>

N-dimensional multi-variate Gaussian:

$$ \ln p(y|\theta) = -\frac{1}{2} [N \ln 2\pi + \ln |\Sigma| + r^T \Sigma^{-1} r], \quad (4) $$

where $y$ and $\theta$ are vectors containing, respectively, all the RV data points and instrumental parameters, while $r$ is the residual vector given by

$$ r(t_i) = y(t_i) - M(t). \quad (5) $$

The elements of the covariance matrix $\Sigma$ are:

$$ \Sigma(t_i, t_j) = k(\mu_i, \mu_j) + (\sigma_w^2 + \sigma_t^2)\delta_{i,t_j}, \quad (6) $$

with $k$ equal to any GP kernel model, or zero for a pure white-noise one. In order to estimate the Bayesian posteriors and evidence, $Z$, of different models, we used the dynamic nested sampling package *dynesty* (Speagle 2020), which adaptively allocates samples based on a posterior structure and, at the same time, estimates evidence and sampling from multi-modal distributions. In general, dynamic nested sampling algorithms sample a dynamic number of live points from the prior ‘volume’ and sequentially replace the point with the lowest likelihood with a new one, while updating the Bayesian evidence by the difference $\Delta Z$. Usually, the stopping criterion is a defined value of $\Delta Z$, below which the algorithm is said to have converged ($\Delta Z \approx 0.5$). However, here we used the default criterion described in Sect. 3.4 of Speagle (2020).

In order to reveal the transiting object suggested by the TESS light curve, we first ran the RV and photometry joint analysis with a simple one-planet model, using the parameters in the DVR produced by the SPOC pipeline as transit-related priors, both with a fixed null and uniformly-sampled eccentricity via the parametrization $S_1 = \sqrt{e} \sin \omega, S_2 = \sqrt{e} \cos \omega$, which is described in Eastman et al. (2013). All the priors are defined in Table B.1. In particular, we set Gaussian priors on both the limb-darkening coefficients (Claret 2017) and the star mean density $\rho_*$ (Sect. 3.1), which was implemented here instead of the scaled semi-major axis, $a/R_*$, because the latter can be recovered using Kepler’s third law using only the period of the respective planet, which is a direct result of any juliet run. In this way, from
<figure><img src="image_9.png" /><figcaption>**Fig. 7.** RV measurements of TOI-1422 versus time are shown on the top panel, while their residuals over the model fit are in the bottom panel. The circles with blue error bars are the RV data taken with HARPS-N. The large and small error bars indicate $\sigma_t$ and $\sigma_w$ (the added jitter term), respectively. In the top panel, the black line represents the two-planet model fit.</figcaption></figure>

the single value of $\rho_*$ we can evenly derive $a/R_*$ in the case of multiple planets.

The best one-planet RV model fit is found with $e = 0$ ($\Delta \ln Z_{e_b=0}^{e_b=0} = 0.7$), but the scatter of the residuals is higher than the average photon-noise uncertainties for this kind of star. In fact, the same peak of 29 days, which was found in the RV GLS periodogram, is also distinctly found in the residuals of the transiting one-planet model (see Fig. 6). Consequently, we proceeded to test two-planet models, whose priors are summed in Table B.2. Since they have comparable statistical significance ($\Delta \ln Z_{e_b \neq 0}^{e_b=0} = 0.4$), we use the results of the eccentric model for the rest of the paper. The two-planet eccentric model is plotted on top of the RVs in Fig. 7, along with its residuals. TOI-1422 b RV semi-amplitude and orbital period are found to be $K_b = 2.47_{-0.46}^{+0.50} \text{ m s}^{-1}$ and $P_b = 12.9972 \pm 0.0006$ days, respectively. The second planet, candidate TOI-1422 c, has an RV semi-amplitude of $K_c = 2.36_{-0.40}^{+0.42} \text{ m s}^{-1}$, orbital period of $P_c = 29.29_{-0.20}^{+0.21}$ days and $T_{0,c} = 2458776.6 \pm 4.6$ BJD (see the posteriors in Fig. C.1 and Table B.3). The eccentricities turn out to be $e_b = 0.04_{-0.03}^{+0.05}$ and $e_c = 0.14_{-0.10}^{+0.17}$, but it is important to note that when they are fixed to zero, the orbital parameters of TOI-1422 b and TOI-1422 c remain, within $1-\sigma$, compatible with those of the eccentric model.

<figure><img src="image_10.png" /><figcaption>**Fig. 8.** RV semi-amplitude *K* and orbital period *P*, along with their $1-\sigma$ error bars, for candidate planet *c*, as functions of the number of data points used for the two-planet eccentric model analysis with *juliet*.</figcaption></figure>

## 3.4. Results

TOI-1422 c's orbital period explains both the main peak found in the residuals of the one-planet model (Fig. 6) and in the RV GLS periodogram (Fig. 5); it is also in 9:4 orbital resonance with the first planet. The difference between the Bayesian evidence of the two-planet eccentric model and the one-planet model ($\Delta \ln Z_{1p}^{2p} = 5.1$) is barely above the *very strong* evidence threshold defined in Kass & Raftery (1995) ($\Delta \ln Z > 5$), so even if the existence of candidate planet *c* remains unproven, we believe the two-planet model is currently the better one to explain the 29-day signal observed in the RVs, due to the lack of evidence of star activity.

Furthermore, the two-planet analysis was replicated with different numbers of data points in order to understand how and if new measurements were impacting the significance of the

second planet detection. As shown in Fig. 8, both the RV semi-amplitude and the period seem to stabilize after ≈60 measurements, which matches the beginning of the second observation season, while the significance of the 29-day peak also grows (Fig. D.4). It is noteworthy to mention that the GLS periodogram of the residuals of the two-planet model does not show peaks below 50% FAP, and hence does not suggest the presence of additional detectable signals.

A phase-folded plot of both the transit and the RVs is shown in Fig. 9 for the eccentric two-planet model. The radius for TOI-1422 b was calculated with the transformations provided by Espinoza (2018) and, using the stellar radius from Sect. 3.1, its revised value turns out to be $R_b = 3.96_{-0.11}^{+0.13} R_⊕$. Using the stellar radius from Table 1, we derived the mass of both objects to be
<figure><img src="image_11.png" /><figcaption>**Fig. 9.** TESS light curve and RV curves phase-folded. Top panel: TOI-1422 b transit, compared to the best-fitting model. Bottom panels: HARPS-N RV data phase-folded to the period of planet b (middle) and candidate c (bottom), along with their residuals over the model. The red circles represent the average value of phased RV data points.</figcaption></figure>

$M_b = 9.0^{+2.3}_{-2.0} M_⊕$ and $M_c \sin i_c = 11.1^{+2.6}_{-2.3} M_⊕$. Their final parameters are reported in Table 2. An independent joint analysis of the HARPS-N RVs and TESS photometry, after the transits were normalized through a local linear fitting, was also performed with a DE-MCMC method (Eastman et al. 2013, 2019), following the same implementation as in Bonomo et al. (2014, 2015). The obtained results are consistent, within $1-\sigma$, with those reported in Table 2.

In order to evaluate possible transit time variations (TTVs) due to the influence of candidate TOI-1422 c over TOI-1422 b, we plot the four mid-transit times minus their expected values (based on the two-planet eccentric model) in Fig. 10, along with different TTV predictions made with the code described in Agol & Deck (2016). Unfortunately, there are not enough transits to draw any conclusions, as all the delays are compatible with zero within $1-\sigma$. Therefore, further precise monitoring of TOI-1422 b transits is encouraged in order to confirm the existence of TOI-1422 c and, overall, better characterize the planetary system.

### 3.5. Other transit events

In the search for TOI-1422 c transits, we found a possible single transit-like event around 2 458 756.35 BTJD days, as shown in Fig. 11, which cannot be related to either TOI-1422 b or TOI-1422 c. We fitted this potential transit using the light curve from the pipeline PATHOS (Nardiello et al. 2019) and retrieved a possible radius of $R_d = 2.82^{+0.38}_{-0.05} R_⊕$, which is compatible with the transit depth observed in the PDC-SAP and SAP light curves as well. The duration of the transit suggests an orbital period longer than that for TOI-1422 c, but this is very uncertain, while the lack of other transits in the TESS light curve suggests an

<table><caption>**Table 2.** Best-fit median values, with upper and lower 68% credibility bands as errors, of the fitted and derived parameters for TOI-1422 b and TOI-1422 c, as extracted from the posterior distribution of the two-planet eccentric model (Table B.3 and Fig. C.1).</caption><thead><tr><th></th><th>TOI-1422 b</th><th>TOI-1422 c</th></tr></thead><tbody><tr><td colspan="3"><strong>Transit and orbital parameters</strong></td></tr><tr><td>K (m s<sup>-1</sup>)</td><td>2.47<sup>+0.50</sup><sub>-0.46</sub></td><td>2.36<sup>+0.42</sup><sub>-0.40</sub></td></tr><tr><td>P<sub>orb</sub> (days)</td><td>12.9972 ± 0.0006</td><td>29.29<sup>+0.21</sup><sub>-0.20</sub></td></tr><tr><td>T<sub>0</sub> (BJD)</td><td>2 458 745.9205<sup>+0.0012</sup><sub>-0.0011</sub></td><td>2 458 776.6<sup>+4.6</sup><sub>-4.5</sub></td></tr><tr><td>T<sub>14</sub> (hours)</td><td>4.52 ± 0.16</td><td>—</td></tr><tr><td>R<sub>p</sub>/R<sub>★</sub></td><td>0.0356<sup>+0.0007</sup><sub>-0.0005</sub></td><td>—</td></tr><tr><td>b</td><td>0.19<sup>+0.11</sup><sub>-0.10</sub></td><td>—</td></tr><tr><td>i (deg)</td><td>89.52<sup>+0.26</sup><sub>-0.28</sub></td><td>—</td></tr><tr><td>a/R<sub>★</sub></td><td>22.72<sup>+0.31</sup><sub>-0.40</sub></td><td>39.05<sup>+0.50</sup><sub>-0.73</sub></td></tr><tr><td>q<sub>1</sub></td><td>0.28<sup>+0.11</sup><sub>-0.08</sub></td><td>—</td></tr><tr><td>q<sub>2</sub></td><td>0.30<sup>+0.05</sup><sub>-0.05</sub></td><td>—</td></tr><tr><td>√e sin ω</td><td>0.018<sup>+0.108</sup><sub>-0.095</sub></td><td>0.120<sup>+0.221</sup><sub>-0.233</sub></td></tr><tr><td>√e cos ω</td><td>-0.149<sup>+0.153</sup><sub>-0.128</sub></td><td>-0.070<sup>+0.349</sup><sub>-0.304</sub></td></tr><tr><td colspan="3"><strong>Derived parameters</strong></td></tr><tr><td>M<sub>p</sub> (M<sub>⊕</sub>)</td><td>9.0<sup>+2.3</sup><sub>-2.0</sub></td><td>—</td></tr><tr><td>M<sub>p</sub> sin i (M<sub>⊕</sub>)</td><td>—</td><td>11.1<sup>+2.6</sup><sub>-2.3</sub></td></tr><tr><td>R<sub>p</sub> (R<sub>⊕</sub>)</td><td>3.96<sup>+0.13</sup><sub>-0.11</sub></td><td>—</td></tr><tr><td>ρ<sub>p</sub> (g cm<sup>-3</sup>)</td><td>0.795<sup>+0.290</sup><sub>-0.235</sub></td><td>—</td></tr><tr><td>log g<sub>p</sub> (cgs)</td><td>2.75<sup>+0.08</sup><sub>-0.14</sub></td><td>—</td></tr><tr><td>a (AU)</td><td>0.108 ± 0.003</td><td>0.185 ± 0.006</td></tr><tr><td>T<sub>eq</sub><sup>(3)</sup> (K)</td><td>867 ± 17</td><td>661 ± 13</td></tr><tr><td>u<sub>1</sub></td><td>0.32<sup>+0.12</sup><sub>-0.10</sub></td><td>—</td></tr><tr><td>u<sub>2</sub></td><td>0.21<sup>+0.10</sup><sub>-0.08</sub></td><td>—</td></tr><tr><td>e</td><td>0.04<sup>+0.05</sup><sub>-0.03</sub></td><td>0.14<sup>+0.17</sup><sub>-0.10</sub></td></tr><tr><td>ω (deg)</td><td>153<sup>+20</sup><sub>-56</sub></td><td>99<sup>+63</sup><sub>-64</sub></td></tr></tbody></table>

Notes. (‡) This is the equilibrium temperature for a zero Bond albedo and uniform heat redistribution to the night side.

<figure><img src="image_12.png" /><figcaption>**Fig. 10.** Residuals for the mid-transit timings of TOI-1422 b versus a linear ephemeris, with $1-\sigma$ error bars, are plotted in black. The green circles, red diamonds, and blue stars represent TTV predictions in the cases of null, average, or maximum eccentricities, respectively, with the error bars showing the uncertainty due to $T_{0,c}$ (see Table 2). The points have been slightly shifted on the x-axis to allow for more visibility.</figcaption></figure>

orbital period between 17 and 22, or longer than, 35 days, thus incompatible with that of TOI-1422 c. PATHOS is a PSF-based approach to TESS data that minimizes the dilution effects in crowded environments, and here it is utilized to extract high-precision photometry of TOI-1422 to independently confirm the
<figure><img src="image_13.png" /><figcaption>**Fig. 11.** PDC-SAP and PATHOS light curves. **Top panel:** TOI-1422 b transits highlighted in red in the PDC-SAP light curve, and the expected TOI-1422 c transits, with their uncertainties, highlighted in blue. A single planetary-transit event is also marked with a vertical line green, and is discussed at the end of Sect. 3.4. **Bottom panel:** Single transit-like event as seen in the PATHOS light curve and the corresponding fit.</figcaption></figure>

presence of this transit even after the application of a differ-
ent neighbour-subtraction technique. Neither the single transit
nor TOI-1422 b transits show correlation with the X,Y pixels
and the sky background signal (Fig. D.5), and the single transit
depth also does not change with different photometric aper-
tures (Fig. D.6). Nevertheless, the three-planet model for the
joint transit-RV analysis is not statistically significant and the
lack of other transits makes the suggestion of another candidate
impossible to justify.

However, no transit compatible with the expected T<sub>0,c</sub> and P<sub>c</sub>
evaluated with the RV and photometry joint analysis, was found
in the SPOC (both SAP and PDC-SAP) light curves, even though
a small part of the supposed transiting window was missed by
TESS. When we take into account both the time-span of the
TESS light curve and TOI-1422 c expected (non-grazing) tran-
sit duration, the probability that such transits would have been
missed can be estimated to be around 1% and 7%, with 1σ
and 3σ uncertainty, respectively, on T<sub>0,c</sub>. Other than misaligned
orbits, another possible explanation for the lack of TOI-1422 c
transits is that despite its mass, which is greater than that of
planet b, its size could be much smaller (similar to the high-
density sub-Neptune, BD+20594b of Espinoza et al. 2016), as
any object with a radius approximately below 2.8R<sub>⊕</sub> might be
disguised in the light curve noise (as proven by the, so far
undetected and uncertain, single transit-like event). Ultimately,
it remains unknown if candidate planet c is transiting or not,
so further high-precision long photometric follow-up observa-
tions will be important to clear up this possibility, along with the
nature of the single transit event. The new TESS observations
of this target, during Sector 57, are definitely welcome as they
might shed some light on both matters.

4. Discussion

4.1. Orbital resonance

As we have seen, candidate c is within 1-σ, in 9:4 orbital
resonance with planet b. This is likely coincidental since the
resonance is fifth-order, and thus very weak, unless one of the
planets is quite eccentric<fnref n="16" /> or the mutual inclination is high.
The exact 9:4 (or 2.25) resonance is within uncertainty, perhaps
only because the uncertainty of the orbital period of TOI-1422 c
is large compared to the tight period uncertainties of transit-
ing planets. As a matter of fact, period ratios a little above two
have been found within many exoplanetary systems (Winn &
Fabrycky 2015), but it is also possible that the 9:4 resonance is
actually the result of a resonant chain of three planets in first-
order 3:2 resonances among each other, with the middle one yet
to be seen. If that is the case, since the period ratios of Kepler
planets near first-order resonances are usually slightly wide of
resonance, the likely orbital period for this unknown exoplanet
would be slightly more than 19.5 days, and thus compatible with
the observed single transit discussed in Sect. 3.4. Given this
orbital period and assuming that an RV semi-amplitude roughly
up to 2 m s⁻¹ might be hidden in the residuals of the two-planet
model, this middle object should not have a mass higher than
≈8 M<sub>⊕</sub>, or a density higher than ≈2 g cm⁻³.

<footnote n="16">¹⁶ We note that, even with the *e*’s suggested by the eccentric fits, which are unusually high compared to most multi-transiting planetary systems according to Xie *et al.* (2016), the 9:4 would not be as strong as a first-order resonance.</footnote>
<figure><img src="image_14.png" /><figcaption>Fig. 12. Planetary masses and radii of the known transiting exoplanets (values taken from the Transiting Extrasolar Planet catalogue, TEP-Cat, which is available at http://www.astro.keele.ac.uk/jkt/tepcat/ catalogue; Southworth 2010, 2011) with equilibrium temperature $T_{eq}$ between 600 and 1000 K and host star radius between 0.6 and 1.5 $R_{\odot}$. Different lines correspond to different mass fractions of relatively cold hydrogen envelopes. The ice giants of the Solar System are displayed in filled black circles. TOI-1422 b is on the low-density envelope of planets with precise mass and/or radius estimations ($\sigma M_{\text{p}}/M_{\text{p}} \le 30\%$; $\sigma R_{\text{p}}/R_{\text{p}} \le 10\%$), one of the reasons that make it potentially valuable for transit spectroscopy.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Fig. 13. Transmission spectroscopy observations (TSM) values with the JWST over the equilibrium temperature for planets with a measured mass in the radius range $2.75 < R_{\oplus} < 4.0$, including TOI-1422 b (green star). Filled black dots and empty squares identify the sample of planets around stars with $T_{\text{eff}} > 5400$ K and $T_{\text{eff}} < 5400$ K, respectively.</figcaption></figure>

## 4.2. Mass-radius diagram and internal structure of planet b

TOI-1422 b is one of the puffier planets with a density of ~0.8 g cm⁻³, which is close to that of Saturn and, therefore, lower than most exoplanets in this mass range. It lies towards the upper-left corner of the mass-radius diagram (Fig. 12), making it very similar to Kepler-36 c (Vissapragada et al. 2020) and especially to Kepler-11 e (Lissauer et al. 2013), which even shares the same kind of host star but is on a longer orbit. On one hand, it has a similar radius compared to Neptune and Uranus in our solar system, but on the other hand, its mass is only about 50–60% that of our ice giants. Thus, an extensive gaseous envelope, surrounding a massive core, is expected to be found in TOI-1422 b. More precisely, the mass fraction of this envelope is expected to be around 10–25% of the total mass of the planet (using the equations of state from Becker et al. 2014), suggesting that the atmosphere has not been blown away by the stellar wind. The nature of this extensive envelope as well as its core requires further investigation. For this purpose, we assess the expected S/N of the JWST/NIRISS measurements<fnref n="17" /> of TOI-1422 b transits compared to planets of similar sizes, by evaluating the transmission spectroscopy metric (TSM) defined in Kempton et al. (2018):

$$ \text{TSM} = (\text{Scale factor}) \times \frac{R_{\oplus}^3 T_{\text{eq}}}{M_{\oplus} R_{\star}^2} \times 10^{-0.2 J}, \quad (7) $$

where the scale factor is a dimensionless normalization constant, equal to 1.28 for planets with $2.75 < R_{\oplus} < 4.0$, and J is the apparent magnitude of the host star in the J band (a filter that is near the middle of the NIRISS bandpass). As a result (Fig. 13), TOI-1422 b ranks fourth<fnref n="18" /> among Neptunes ($2.75 < R_{\oplus} < 4.0$) orbiting G-F dwarfs ($T_{\text{eff}} > 5400$K), but being the one with the lowest density, it is definitely an interesting candidate for atmospheric characterization by the JWST.

## 5. Conclusions

In this paper, we have confirmed the planetary nature of the TESS transiting planet TOI-1422 b, which turns out to be a low-density and warm Neptune-sized planet orbiting an astrometrically, and overall magnetically, quiet G2 V star. Therefore, TOI-1422 b is the latest addition to the low-populated range of exoplanets with the size of Neptune, but with Saturn-like density. In order to well constrain the mass of TOI-1422 b, a long RV monitoring with more than a hundred observations was necessary with the HARPS-N instrument at the TNG in La Palma, which resulted in fully characterized orbital and physical parameters of this new planetary system. On top of that, our RV measurements also suggest the presence in the system of a possibly non-transiting, heavier candidate planet, TOI-1422 c, in a weak 9:4 orbital resonance with its inner brother, which will require further study to validate.

*Acknowledgements.* We acknowledge the use of public TESS data from pipelines at the TESS Science Office and at the TESS Science Processing Operations Center. Resources supporting this work were provided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center for the production of the SPOC data products. The research has made use of the SIMBAD database, operated at CDS, Strasbourg, France, NASA’s Astrophysics Data System and the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. This research was also partly supported by the Department of Energy under awards DE-NA0003904 and DE-FOA0002633*

<footnote n="17">¹⁷ From a 10-h observing programme assuming a cloud-free, solar-metallicity, H₂-dominated atmosphere.</footnote>

<footnote n="18">¹⁸ Following TOI-561 c (Lacedelli et al. 2020, 2022), HD 136352 c (Kane et al. 2020; Delrez et al. 2021) and HD 63935 b (Scarsdale et al. 2021)</footnote>
(to S.B.J., principal investigator, and collaborator Li Zeng) with Harvard University and by the Sandia Z Fundamental Science Program. This research represents the authors' views and not those of the Department of Energy. The work is 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 work has also 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. L.M. acknowledges support from the "Fondi di Ricerca Scientifica d'Ateneo 2021" of the University of Rome "Tor Vergata". J.L-B. acknowledges financial support received from "la Caixa" Foundation (ID 100010434) and from the European Union's Horizon 2020 research and innovation programme under the Marie Slodowska-Curie grant agreement No 847648, with fellowship code LCF/BQ/PI20/11760023. This research has also been partly funded by the Spanish State Research Agency (AEI) Projects No. PID2019-107061GB-C61 and No. MDM-2017-0737 Unidad de Excelencia "María de Maeztu"-Centro de Astrobiología (INTA-CSIC). We acknowledge financial contribution from the agreement ASI-NAF no. 2018-16-HH.0. D.D. acknowledges support from the TESS Guest Investigator Program grants 80NSSC21K0108 and 80NSSC22K0185. L.N. acknowledges the support of the ARIEL ASI-NAF agreement 2021-5-HH.0.

## References

Agol, E., & Deck, K. 2016, ApJ, 818, 177

Aller, A., Lillo-Box, J., Jones, D., Miranda, L. F., & Barceló Forteza, S. 2020, A&A, 635, A128

Ambikasaran, S., Foreman-Mackey, D., Greengard, L., Hogg, D. W., & O'Neil, M. 2015, IEEE Trans. Pattern Anal. Mach. Intell., 38, 252

Anglada-Escudé, G., & Butler, R. P. 2012, ApJS, 200, 15

Armstrong, D. J., Lopez, T. A., Adibekyan, V., et al. 2020, Nature, 583, 39

Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, VizieR Online Data Catalog: I/352

Barstow, J. K., Aigrain, S., Irwin, P. G. J., Kendrew, S., & Fletcher, L. N. 2015, MNRAS, 448, 2546

Becker, A., Lorenzen, W., Fortney, J. J., et al. 2014, ApJS, 215, 21

Bhatti, W., Bouma, L., Joshua, John., & Price-Whelan, A. 2020, https://doi.org/10.5281/zenodo.3723832

Biazzo, K., Gratton, R., Desidera, S., et al. 2015, A&A, 583, A135

Biazzo, K., Bozza, V., Mancini, L., & Sozzetti, A. 2022, in Astrophysics and Space Science Library, Demographics of Exoplanetary Systems, Lecture Notes of the 3rd Advanced School on Exoplanetary Science, eds. K. Biazzo, V. Bozza, L. Mancini, & A. Sozzetti, 466, 143

Bitsch, B., Raymond, S. N., & Izidoro, A. 2019, A&A, 624, A109

Bonomo, A. S., Sozzetti, A., Lovis, C., et al. 2014, A&A, 572, A2

Bonomo, A. S., Sozzetti, A., Santerne, A., et al. 2015, A&A, 575, A85

Borsa, F., Allart, R., Casasayas-Barris, N., et al. 2021, A&A, 645, A24

Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977

Bouchy, F., Díaz, R. F., Hébrard, G., et al. 2013, A&A, 549, A49

Brande, J., Crossfield, I. J. M., Kreidberg, L., et al. 2022, AJ, accepted [arXiv:2201.04197]

Brogi, M., Giacobbe, P., Guilluy, G., et al. 2018, A&A, 615, A16

Brown, A. G. A., Vallenari, A., Prusti, T., et al. 2018, A&A, 616, A1

Brown, A. G. A., Vallenari, A., Prusti, T., et al. 2021, A&A, 650, C3

Bryson, S. T., Jenkins, J. M., Klaus, T. C., et al. 2020, Kepler Data Processing Handbook: Target and Aperture Definitions: Selecting Pixels for Kepler Downlink, Kepler Science Document KSCI-19081-003

Bryson, S. T., Jenkins, J. M., Klaus, T. C., et al. 2010, SPIE Conf. Ser., 7740, 77401D

Buchhave, L. A., Bizzarro, M., Latham, D. W., et al. 2014, Nature, 590, 593

Carleo, I., Desidera, S., Nardiello, D., et al. 2021, A&A, 645, A71

Casasayas-Barris, N., Pallé, E., Yan, F., et al. 2019, A&A, 628, A9

Castelli, F., & Kurucz, R. L. 2003, in Modelling of Stellar Atmospheres, eds. N. Piskunov, W. W. Weiss, & D. F. Gray, 210, A20

Chontos, A., Murphy, J. M. A., MacDougall, M. G., et al. 2022, AJ, 163, 297

Claret, A. 2017, A&A, 600, A30

Cloutier, R., Eastman, J. D., Rodriguez, J. E., et al. 2020, AJ, 160, 3

Cosentino, R., Lovis, C., Pepe, F., et al. 2012, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Harps-N: the new planet hunter at TNG, 8446, 84461V

Covino, E., Esposito, M., Barbieri, M., et al. 2013, A&A, 554, A28

Cumming, A., Butler, R. P., Marcy, G. W., et al. 2008, PASP, 120, 531

Cutri, R. M., Wright, E. L., Conrow, T., et al. 2021, VizieR Online Data Catalog: II/328

Delrez, L., Ehrenreich, D., Alibert, Y., et al. 2021, Nat. Astron., 5, 775

Diamond-Lowe, H., Berta-Thompson, Z., Charbonneau, D., & Kempton, E. M. R. 2018, AJ, 156, 42

Dotter, A. 2016, ApJS, 222, 8

Doyle, A. P., Davies, G. R., Smalley, B., Chaplin, W. J., & Elsworth, Y. 2014, MNRAS, 444, 3592

Eastman, J. 2017, EXOFASTv2: Generalized publication-quality exoplanet modeling code, Astrophysics Source Code Library [record ascl:1710.003]

Eastman, J., Gaudí, B. S., & Agol, E. 2013, PASP, 125, 83

Eastman, J. D., Rodriguez, J. E., Agol, E., et al. 2019, PASP, submitted [arXiv:1907.09480]

Espinoza, N. 2018, RNAAS, 2, 209

Espinoza, N., Brahm, R., Jordán, A., et al. 2016, ApJ, 830, 43

Espinoza, N., Kossakowski, D., & Brahm, R. 2019, MNRAS, 490, 2262

Foreman-Mackey, D., Agol, E., Ambikasaran, S., & Angus, R. 2017, AJ, 154, 220

Fressin, F., Torres, G., Charbonneau, D., et al. 2013, ApJ, 766, 81

Fulton, B. J., & Petigura, E. A. 2018, AJ, 156, 264

Fulton, B. J., Petigura, E. A., Blunt, S., & Sinukoff, E. 2018, PASP, 130, 044504

Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1

Gandolfi, D., Barragán, O., Livingston, J. H., et al. 2018, A&A, 619, A10

Girardi, L., Barbieri, M., Groenewegen, M. A. T., et al. 2012, Astrophys. Space Sc. Proc., 26, 165

Gomes da Silva, J., Figueira, P., Santos, N., & Faria, J. 2018, J. Open Source Softw., 3, 667

Gray, R. O., & Corbally, C. 2009, Stellar Spectral Classification (Princeton University Press)

Guerrero, N. M., Seager, S., Huang, C. X., et al. 2021, ApJS, 254, 39

Henden, A. A., Levine, S., Terrell, D., & Welch, D. L. 2015, in American Astronomical Society Meeting Abstracts, 225, 336.16

Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27

Hormuth, F., Brandner, W., Hippler, S., & Henning, T. 2008, J. Phys. Conf. Ser., 131, 012051

Howard, A. W., Marcy, G. W., Johnson, J. A., et al. 2010, Science, 330, 653

Howard, A. W., Marcy, G. W., Bryson, S. T., et al. 2012, ApJS, 201, 15

Howell, S. B., Sobeck, C., Haas, M., et al. 2014, PASP, 126, 398

Hunter, A., Macgregor, A. B., Szabo, T., Wellington, C., & Bellgard, M. I. 2012, Source Code Biol. Med., 7, 1

Jenkins, J. M. 2002, ApJ, 575, 493

Jenkins, J. M., Chandrasekaran, H., McCauliff, S. D., et al. 2010, SPIE Conf. Ser., 7740, 77400D

Jenkins, J. M., Twicken, J. D., McCauliff, S., et al. 2016, SPIE Conf. Ser., 9913, 99133E

Jenkins, J. M., Tenenbaum, P., Seader, S., et al. 2020, Kepler Data Processing Handbook: Transiting Planet Search, ed. J. M. Jenkins. Kepler Science Document KSCI-19081-003, 9

Kane, S. R., Yalçınkaya, S., Osborn, H. P., et al. 2020, AJ, 160, 129

Kass, R. E., & Raftery, A. E. 1995, J. Am. Stat. Assoc., 90, 773

Kempton, E. M.-R., Bean, J. L., Louie, D. R., et al. 2018, PASP, 130, 114401

Kipping, D. M. 2013, MNRAS, 435, 2152

Knutson, H. A., Benneke, B., Deming, D., & Homeier, D. 2014, Nature, 505, 66

Kreidberg, L. 2015, PASP, 127, 1161

Kreidberg, L., Bean, J. L., Désert, J.-M., et al. 2014, Nature, 505, 69

Lacedelli, G., Malavolta, L., Borsato, L., et al. 2020, MNRAS, 501, 4148

Lacedelli, G., Wilson, T. G., Malavolta, L., et al. 2022, MNRAS, 511, 4551

Lind, K., Asplund, M., & Barklem, P. S. 2009, A&A, 503, 541L

Lissauer, J. J., Jontof-Hutter, D., Rowe, J. F., et al. 2019, PASP, 131, 1

Lightkurve Collaboration (Cardoso, J. V. d. M., et al.) 2018, Lightkurve: Kepler and TESS time series analysis in Python, Astrophysics Source Code Library [record ascl:1812.013]

Lillo-Box, J., Barrado, D., & Bouy, H. 2012, A&A, 546, A10

Lillo-Box, J., Barrado, D., & Bouy, H. 2014, A&A, 566, A103

Lissauer, J. J., Jontof-Hutter, D., Rowe, J. F., et al. 2013, ApJ, 770, 131

Lozovsky, M., Helled, R., Dorn, C., & Venturini, J. 2018, ApJ, 866, 49

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

Mayor, M., Marmier, M., Lovis, C., et al. 2011, ArXiv e-prints, [arXiv:1109.2497]

Miller-Ricci, E., Seager, S., & Sasselov, D. 2009, ApJ, 690, 1056

Modirrousta-Galian, D., Locci, D., & Micela, G. 2020, ApJ, 891, 158

Morris, R. L., Twicken, J. D., Smith, J. C., et al. 2020, Kepler Data Processing Handbook: Photometric Analysis, ed. J. M. Jenkins, Kepler Science Document KSCI-19081-003, id. 6

Nardiello, D., Borsato, L., Piotto, G., et al. 2019, MNRAS, 490, 3806

Pasquini, L., Biazzo, K., Bonifacio, P., Randich, S., & Bedin, L. R. 2009, A&A, 489, 677
Perger, M., García-Piquer, A., Ribas, I., et al. 2017, A&A, 598, A26

Perruchot, S., Kohler, D., Bouchy, F., et al. 2008, in Ground-based and Air-borne Instrumentation for Astronomy II, eds. I. S. McLean, & M. M. Casali, International Society for Optics and Photonics (SPIE), 7014, 235

Petigura, E. A., Howard, A. W., Marcy, G. W., et al. 2017, AJ, 154, 107

Petigura, E. A., Marcy, G. W., & Howard, A. W. 2013, ApJ, 770, 69

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

Pollacco, D. L., Skillen, I., Collier Cameron, A., et al. 2006, PASP, 118, 1407

Poretti, E., Boccato, C., Claudi, R., et al. 2016, Mem. Soc. Astron. Italiana, 87, 141

Price-Whelan, A. M., Sipőcz, B. M., Günther, H. M., et al. 2018, AJ, 156, 123

Prusti, T., de Bruijne, J. H. J., Brown, A. G. A., et al. 2016, A&A, 595, A1

Ramsay, S., Cirasuolo, M., Amico, P., et al. 2021, The Messenger, 182, 3

Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, Proceedings of the SPIE, 9143, 914320

Scarsdale, N., Murphy, J. M. A., Batalha, N. M., et al. 2021, AJ, 162, 215

Schlaflay, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525

Skidmore, W., TMT International Science Development Teams, & Science Advisory Committee, T. 2015, Res. Astron. Astrophys., 15, 1945

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Smith, J. C., Stumpe, M. C., Cleve, J. E. V., et al. 2012, PASP, 124, 1000

Sneden, C. 1973, AJ, 184, 839

Snellen, I. A. G., de Kok, R. J., de Mooij, E. J. W., & Albrecht, S. 2010, Nature, 465, 1049

Southworth, J. 2010, MNRAS, 408, 1689

Southworth, J. 2011, MNRAS, 417, 2166

Speagle, J. S. 2020, MNRAS, 493, 3132

Spearman, C. 1904, Am. J. Psychol., 15, 72

Stassun, K. G., Oelkers, R. J., Pepper, J., et al. 2018, AJ, 156, 102

Stassun, K. G., Oelkers, R. J., Paegert, M., et al. 2019, AJ, 158, 138

Strehl, K. 1902, Astron. Nachr., 158, 89

Stumpe, M. C., Smith, J. C., Van Cleve, J. E., et al. 2012, PASP, 124, 985

Stumpe, M. C., Smith, J. C., Catanzarite, J. H., et al. 2014, PASP, 126, 100

Tinetti, G., Eccleston, P., Haswell, C., et al. 2021, ArXiv e-prints, [arXiv:2104.04824]

Trifonov, T. 2019, The Exo-Striker: Transit and radial velocity interactive fitting tool for orbital analysis and N-body simulations, Astrophysics Source Code Library [record ascl:1906.004]

Turbet, M., Bolmont, E., Ehrenreich, D., et al. 2020, A&A, 638, A41

Twicken, J. D., Catanzarite, J. H., Clarke, B. D., et al. 2018, PASP, 130, 064502

Twicken, J. D., Clarke, B. D., Bryson, S. T., et al. 2010, in Software and Cyberinfrastructure for Astronomy, eds. N. M. Radziwill, & A. Bridger, International Society for Optics and Photonics (SPIE), 7740, 749

Van Eylen, V., Agentoft, C., Lundkvist, M. S., et al. 2018, MNRAS, 479, 4786

Vissapragada, S., Jontof-Hutter, D., Shporer, A., et al. 2020, AJ, 159, 108

Weiss, L. M., & Marcy, G. W. 2014, ApJ, 783, L6

Winn, J. N., & Fabrycky, D. C. 2015, ARA&A, 53, 409

Wright, J. T., Marcy, G. W., Howard, A. W., et al. 2012, ApJ, 753, 160

Wytenbach, A., Ehrenreich, D., Lovis, C., Udry, S., & Pepe, F. 2015, A&A, 577, A62

Xie, J.-W., Dong, S., Zhu, Z., et al. 2016, Proc. Natl. Acad. Sci., 113, 11431

Zeng, L., Jacobsen, S. B., Sasselov, D. D., et al. 2019, Proc. Natl. Acad. Sci. U.S.A., 116, 9723

<footnote n="1">¹ Department of Physics, University of Rome “Tor Vergata”, Via della Ricerca Scientifica 1, 00133 Rome, Italy e-mail: luca.naponiello@unifi.it</footnote>

<footnote n="2">² Department of Physics and Astronomy, University of Florence, Largo Enrico Fermi 5, 50125 Firenze, Italy</footnote>

<footnote n="3">³ INAF – Turin Astrophysical Observatory, via Osservatorio 20, 10025 Pino Torinese, Italy</footnote>

<footnote n="4">⁴ Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany</footnote>

<footnote n="5">⁵ INAF – Osservatorio Astronomico di Roma, via Frascati 33, 00040 Monte Porzio Catone (RM), Italy</footnote>

<footnote n="6">⁶ Aix-Marseille Université, CNRS, CNES, LAM, Marseille, France</footnote>

<footnote n="7">⁷ INAF – Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy</footnote>

<footnote n="8">⁸ Dipartimento di Fisica, Università degli Studi di Torino, via Pietro Giuria 1, 10125 Torino, Italy</footnote>

<footnote n="9">⁹ Centro de Astrobiología, Depto. de Astrofísica, ESAC campus, 28692, Villanueva de la Cañada (Madrid), Spain</footnote>

<footnote n="10">¹⁰ Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA</footnote>

<footnote n="11">¹¹ Department of Earth and Planetary Sciences, Harvard University, 20 Oxford Street, Cambridge, MA 02138, USA</footnote>

<footnote n="12">¹² INAF – Osservatorio Astronomico di Trieste, via Tiepolo 11, 34143 Trieste, Italy</footnote>

<footnote n="13">¹³ INAF – Osservatorio Astronomico di Capodimonte, Salita Moiariello 16, 80131 Naples, Italy</footnote>

<footnote n="14">¹⁴ INAF – Osservatorio Astrofísico di Catania, via S. Sofia 78, 95123 Catania, Italy</footnote>

<footnote n="15">¹⁵ INAF – Osservatorio Astronomico di Palermo, Piazza del Parla-mento 1, 90134 Palermo, Italy</footnote>

<footnote n="16">¹⁶ INAF – Osservatorio Astronomico di Cagliari, via della Scienza 5, 09047 Selargius (CA), Italy</footnote>

<footnote n="17">¹⁷ Dipartimento di Fisica e Astronomia “Galileo Galilei” – Università di Padova, Vicolo dell’Osservatorio 2, 35122 Padova, Italy</footnote>

<footnote n="18">¹⁸ INAF – Osservatorio Astronomico di Brera, Via E. Bianchi 46, 23807 Merate (LC), Italy</footnote>

<footnote n="19">¹⁹ Fundación Galileo Galilei – INAF, Rambla José Ana Fernandez Pérez 7, 38712 Breña Baja, TF, Spain</footnote>

<footnote n="20">²⁰ European Southern Observatory, Karl-Schwarzschild-Straße 2, 85748 Garching, Germany</footnote>

<footnote n="21">²¹ Institut d'Astrophysique de Paris, CNRS, UMR 7095 & Sorbonne Université, UPMC Paris 6, 98bis Bd Arago, 75014 Paris, France</footnote>

<footnote n="22">²² Observatoire de Haute-Provence, CNRS, Université d'Aix-Marseille, 04870 Saint-Michel-l'Observatoire, France</footnote>

<footnote n="23">²³ Aix-Marseille Université, CNRS, CNES, LAM, Marseille, France</footnote>

<footnote n="24">²⁴ Université Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France</footnote>

<footnote n="25">²⁵ Department of Physics, Shahid Beheshti University, Tehran, Iran</footnote>

<footnote n="26">²⁶ Laboratoire J.-L. Lagrange, OCA, Université de Nice-Sophia Antipolis, CNRS, Campus Valrose, 06108 Nice Cedex 2, France</footnote>

<footnote n="27">²⁷ IRAP, Université de Toulouse, CNRS, UPS, CNES, 31400 Toulouse, France</footnote>

<footnote n="28">²⁸ Department of Physics and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA</footnote>

<footnote n="29">²⁹ NASA Goddard Space Flight Center, 8800 Greenbelt Road, Greenbelt, MD 20771, USA</footnote>

<footnote n="30">³⁰ Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA</footnote>

<footnote n="31">³¹ Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA</footnote>

<footnote n="32">³² Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA</footnote>

<footnote n="33">³³ Department of Physics and Astronomy, University of New Mexico, 210 Yale Blvd NE, Albuquerque, NM 87106, USA</footnote>

<footnote n="34">³⁴ NASA Ames Research Center, Moffett Field, CA 94035, USA</footnote>

<footnote n="35">³⁵ SETI Institute, Mountain View, CA 94043, USA</footnote>

<footnote n="36">³⁶ Space Science & Astrobiology Division MS 245-3, NASA Ames Research Center, Moffett Field, CA 94035, USA</footnote>

<footnote n="37">³⁷ INAF-Telescopio Nazionale Galileo, Apartado 565, 38700, Santa Cruz de La Palma, Spain</footnote>

<footnote n="38">³⁸ Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544, USA</footnote>

<footnote n="39">³⁹ Astrophysics Group, Keele University, Staffordshire ST5 5BG, UK</footnote>

<footnote n="40">⁴⁰ Department of Aeronautics and Astronautics, MIT, 77 Massachusetts Avenue, Cambridge, MA 02139, USA</footnote>
## Appendix A: HARPS-N RV datapoints

<table><caption>Table A.1. HARPS-N RV data points and activity indexes (used in Fig. 5) obtained with the TERRA reduction pipeline between June 8, 2020 and January 21, 2022. The four lines in bold highlight the RV data points that have been removed because they do not fit Chauvenet's criterion.</caption><thead><tr><th>BJD<sub>UTC</sub><br/>-2457000 [days]</th><th>RV<br/>[m s<sup>-1</sup>]</th><th>±lσ<sub>RV</sub></th><th>FWHM</th><th>BIS<sup>(†)</sup></th><th>Exp.<sup>(*)</sup><br/>[sec]</th><th>S/N</th><th>I<sub>Ca ii</sub></th><th>I<sub>Hα06</sub></th><th>I<sub>He i</sub></th><th>I<sub>Na i</sub></th><th>I<sub>Ca i</sub></th><th>I<sub>Hα16</sub></th></tr></thead><tbody><tr><td>2008.69349316</td><td>-7.95</td><td>1.95</td><td>7174.86</td><td>-7.56</td><td>900</td><td>42.8</td><td>0.549907</td><td>0.084427</td><td>0.097028</td><td>0.180378</td><td>0.527652</td><td>0.683496</td></tr><tr><td>2009.71160774</td><td>-3.78</td><td>2.25</td><td>7181.28</td><td>-6.09</td><td>900</td><td>34.3</td><td>0.555074</td><td>0.078803</td><td>0.097526</td><td>0.179596</td><td>0.522912</td><td>0.63963</td></tr><tr><td>2026.66187880</td><td>-4.51</td><td>2.31</td><td>7164.91</td><td>0.894</td><td>900</td><td>44.2</td><td>0.554022</td><td>0.082134</td><td>0.100171</td><td>0.187731</td><td>0.505988</td><td>0.687073</td></tr><tr><td>2027.72062839</td><td>4.586</td><td>2.2</td><td>7153.54</td><td>-4.7</td><td>900</td><td>36.4</td><td>0.558749</td><td>0.083288</td><td>0.107005</td><td>0.187897</td><td>0.530667</td><td>0.687811</td></tr><tr><td>2028.69301741</td><td>3.085</td><td>2.06</td><td>7178.27</td><td>-6.91</td><td>900</td><td>39.3</td><td>0.587418</td><td>0.078941</td><td>0.098429</td><td>0.181754</td><td>0.53219</td><td>0.657997</td></tr><tr><td>2037.70925340</td><td>-1.04</td><td>2.78</td><td>7190.11</td><td>1.317</td><td>900</td><td>27.1</td><td>0.588273</td><td>0.087717</td><td>0.100143</td><td>0.181427</td><td>0.533928</td><td>0.670313</td></tr><tr><td>2038.71004788</td><td>-6.86</td><td>2.35</td><td>7174.7</td><td>-12.8</td><td>900</td><td>33.8</td><td>0.589738</td><td>0.08641</td><td>0.102117</td><td>0.18419</td><td>0.531931</td><td>0.671326</td></tr><tr><td>2039.71645615</td><td>-2.51</td><td>2.17</td><td>7176.79</td><td>-5.84</td><td>900</td><td>40.5</td><td>0.590585</td><td>0.078021</td><td>0.099401</td><td>0.182083</td><td>0.530063</td><td>0.658114</td></tr><tr><td>2040.70404346</td><td>-0.98</td><td>2.65</td><td>7182.67</td><td>-18.1</td><td>900</td><td>31.6</td><td>0.593109</td><td>0.085039</td><td>0.099033</td><td>0.182093</td><td>0.528559</td><td>0.650582</td></tr><tr><td>2050.69789104</td><td>-1.68</td><td>1.94</td><td>7172.47</td><td>-1.01</td><td>900</td><td>40.2</td><td>0.594075</td><td>0.087503</td><td>0.099635</td><td>0.183098</td><td>0.534144</td><td>0.666092</td></tr><tr><td>2051.69861636</td><td>-4.63</td><td>1.99</td><td>7177.72</td><td>-11.5</td><td>900</td><td>38.6</td><td>0.591485</td><td>0.076251</td><td>0.099171</td><td>0.182165</td><td>0.525447</td><td>0.654697</td></tr><tr><td>2054.72044480</td><td>-1.22</td><td>3.15</td><td>7162.98</td><td>-4.59</td><td>900</td><td>30.1</td><td>0.564703</td><td>0.094168</td><td>0.100699</td><td>0.182831</td><td>0.523618</td><td>0.663779</td></tr><tr><td>2056.61577511</td><td>-3.21</td><td>6.73</td><td>7165.8</td><td>-6.58</td><td>900</td><td>12.7</td><td>0.584758</td><td>0.089195</td><td>0.098718</td><td>0.179674</td><td>0.533928</td><td>0.670313</td></tr><tr><td>2069.66437348</td><td>1.928</td><td>1.46</td><td>7167.86</td><td>-1.31</td><td>900</td><td>52.5</td><td>0.562065</td><td>0.080459</td><td>0.103277</td><td>0.183022</td><td>0.530759</td><td>0.690607</td></tr><tr><td>2070.70250746</td><td>-3.5</td><td>1.57</td><td>7172.32</td><td>-3.47</td><td>900</td><td>47.8</td><td>0.561804</td><td>0.078625</td><td>0.09964</td><td>0.18108</td><td>0.529827</td><td>0.668034</td></tr><tr><td>2071.7080293</td><td>-0.88</td><td>2.12</td><td>7195.44</td><td>-4.27</td><td>900</td><td>39.4</td><td>0.589073</td><td>0.083154</td><td>0.099452</td><td>0.180534</td><td>0.525844</td><td>0.65915</td></tr><tr><td>2072.70309705</td><td>-2.2</td><td>1.86</td><td>7185.56</td><td>-14.3</td><td>900</td><td>42.3</td><td>0.557986</td><td>0.082279</td><td>0.098576</td><td>0.180502</td><td>0.525855</td><td>0.654666</td></tr><tr><td>2075.63970105</td><td>-3.37</td><td>1.9</td><td>7169.71</td><td>0.109</td><td>900</td><td>45.6</td><td>0.594915</td><td>0.083487</td><td>0.100564</td><td>0.181682</td><td>0.527055</td><td>0.685465</td></tr><tr><td>2076.71092255</td><td>-6.68</td><td>2.31</td><td>7174.73</td><td>-9.8</td><td>900</td><td>33.1</td><td>0.554061</td><td>0.07931</td><td>0.103664</td><td>0.184892</td><td>0.530864</td><td>0.665231</td></tr><tr><td>2078.68129258</td><td>-6.1</td><td>1.94</td><td>7175.62</td><td>-16.9</td><td>900</td><td>44.2</td><td>0.584749</td><td>0.082137</td><td>0.100369</td><td>0.183009</td><td>0.534477</td><td>0.67444</td></tr><tr><td>2079.69077171</td><td>-4.53</td><td>1.95</td><td>7176.88</td><td>-8.31</td><td>900</td><td>45.3</td><td>0.560252</td><td>0.079819</td><td>0.100953</td><td>0.183346</td><td>0.530948</td><td>0.673179</td></tr><tr><td>2091.68437135</td><td>8.017</td><td>1.77</td><td>7185.69</td><td>-8.84</td><td>900</td><td>41.1</td><td>0.549246</td><td>0.082682</td><td>0.099392</td><td>0.181514</td><td>0.526944</td><td>0.682402</td></tr><tr><td>2092.62561704</td><td>-1.48</td><td>2.43</td><td>7172.71</td><td>-7.71</td><td>900</td><td>34.1</td><td>0.581979</td><td>0.086567</td><td>0.100783</td><td>0.183421</td><td>0.530712</td><td>0.666342</td></tr><tr><td>2093.65087003</td><td>7.487</td><td>7.04</td><td>7164.43</td><td>-16.6</td><td>900</td><td>13.4</td><td>0.564012</td><td>0.107488</td><td>0.101191</td><td>0.188208</td><td>0.50388</td><td>0.653611</td></tr><tr><td>2094.60798244</td><td>-1.1</td><td>1.7</td><td>7168.72</td><td>-6.22</td><td>900</td><td>42.4</td><td>0.5533</td><td>0.084541</td><td>0.100049</td><td>0.184004</td><td>0.527844</td><td>0.70821</td></tr><tr><td>2095.61169299</td><td>-0.2</td><td>2.3</td><td>7152.3</td><td>3.952</td><td>900</td><td>36.4</td><td>0.585249</td><td>0.077758</td><td>0.105133</td><td>0.187703</td><td>0.525175</td><td>0.697218</td></tr><tr><td>2096.62962752</td><td>-0.29</td><td>1.86</td><td>7166.3</td><td>-11.5</td><td>900</td><td>46.4</td><td>0.541437</td><td>0.080363</td><td>0.09865</td><td>0.182500</td><td>0.530522</td><td>0.677679</td></tr><tr><td>2097.66412402</td><td>-2.95</td><td>2.64</td><td>7180.27</td><td>-4.25</td><td>900</td><td>33.7</td><td>0.584758</td><td>0.086406</td><td>0.102655</td><td>0.182256</td><td>0.52911</td><td>0.64144</td></tr><tr><td>2099.63975353</td><td>-8.99</td><td>2.53</td><td>7178.63</td><td>-1</td><td>900</td><td>32.3</td><td>0.511441</td><td>0.083238</td><td>0.09935</td><td>0.183486</td><td>0.52374</td><td>0.664088</td></tr><tr><td>2106.65854394</td><td>-2.12</td><td>2.29</td><td>7175.27</td><td>-17.8</td><td>900</td><td>34.3</td><td>0.542303</td><td>0.089511</td><td>0.101622</td><td>0.184597</td><td>0.529016</td><td>0.670419</td></tr><tr><td>2110.71566068</td><td>-6.18</td><td>3.41</td><td>7186.25</td><td>5.57</td><td>900</td><td>25.3</td><td>0.512902</td><td>0.082166</td><td>0.101099</td><td>0.18416</td><td>0.530842</td><td>0.67291</td></tr><tr><td>2111.55269319</td><td>-6.16</td><td>1.92</td><td>7162.99</td><td>-2.12</td><td>900</td><td>41</td><td>0.525014</td><td>0.080415</td><td>0.100104</td><td>0.185605</td><td>0.529595</td><td>0.68975</td></tr><tr><td>2112.58796961</td><td>-7.39</td><td>1.88</td><td>7169.58</td><td>-11.6</td><td>900</td><td>40</td><td>0.534984</td><td>0.083924</td><td>0.099656</td><td>0.182888</td><td>0.53337</td><td>0.680635</td></tr><tr><td>2119.69014147</td><td>11.98</td><td>2.85</td><td>7190.32</td><td>-9.13</td><td>900</td><td>25.8</td><td>0.545065</td><td>0.087481</td><td>0.098395</td><td>0.182445</td><td>0.534123</td><td>0.626241</td></tr><tr><td>2120.67679558</td><td>-2.05</td><td>2.87</td><td>7208.68</td><td>-6.74</td><td>900</td><td>27.8</td><td>0.584749</td><td>0.086578</td><td>0.102032</td><td>0.186226</td><td>0.523501</td><td>0.627348</td></tr><tr><td>2125.54442015</td><td>-5.53</td><td>3.13</td><td>7169.97</td><td>-3.52</td><td>900</td><td>28.1</td><td>0.555249</td><td>0.085601</td><td>0.101368</td><td>0.186582</td><td>0.523009</td><td>0.665638</td></tr><tr><td>2126.56835804</td><td>-6.72</td><td>2.67</td><td>7174.3</td><td>-16.4</td><td>1200</td><td>32.1</td><td>0.567202</td><td>0.083452</td><td>0.099046</td><td>0.183561</td><td>0.527745</td><td>0.651399</td></tr><tr><td>2127.63415799</td><td>0.931</td><td>2.18</td><td>7180.78</td><td>-4.92</td><td>900</td><td>35.3</td><td>0.539365</td><td>0.085257</td><td>0.098857</td><td>0.18364</td><td>0.524582</td><td>0.655533</td></tr><tr><td>2130.59628332</td><td>-1.27</td><td>3.36</td><td>7187.06</td><td>-7.56</td><td>900</td><td>23.7</td><td>0.576932</td><td>0.083983</td><td>0.098713</td><td>0.183139</td><td>0.531644</td><td>0.651165</td></tr><tr><td>2134.62011272</td><td>0.717</td><td>3.58</td><td>7200.6</td><td>-10.9</td><td>900</td><td>22.7</td><td>0.580358</td><td>0.091434</td><td>0.100845</td><td>0.184286</td><td>0.525551</td><td>0.661751</td></tr><tr><td>2137.55086797</td><td>-6.08</td><td>1.75</td><td>7184.04</td><td>-1.32</td><td>900</td><td>39.9</td><td>0.587977</td><td>0.087623</td><td>0.097656</td><td>0.18263</td><td>0.528299</td><td>0.65772</td></tr><tr><td>2153.56077137</td><td>-4.39</td><td>7.54</td><td>7144.15</td><td>-24</td><td>900</td><td>14.5</td><td>0.561764</td><td>0.077853</td><td>0.102889</td><td>0.189554</td><td>0.527182</td><td>0.670189</td></tr><tr><td>2156.53992496</td><td>-0.87</td><td>2.45</td><td>7194.74</td><td>-1.33</td><td>900</td><td>34.9</td><td>0.596947</td><td>0.089431</td><td>0.103586</td><td>0.187334</td><td>0.531839</td><td>0.639791</td></tr><tr><td><strong>2157.58172574</strong></td><td><strong>15.06</strong></td><td><strong>6.49</strong></td><td><strong>7189.21</strong></td><td><strong>8.093</strong></td><td><strong>900</strong></td><td><strong>15.9</strong></td><td><strong>0.554652</strong></td><td><strong>0.093102</strong></td><td><strong>0.098102</strong></td><td><strong>0.18062</strong></td><td><strong>0.525648</strong></td><td><strong>0.651384</strong></td></tr><tr><td>2169.31333487</td><td>0.713</td><td>5.22</td><td>7202.4</td><td>-9.03</td><td>900</td><td>19.2</td><td>0.569926</td><td>0.086208</td><td>0.102926</td><td>0.186679</td><td>0.532988</td><td>0.706365</td></tr><tr><td>2170.34375539</td><td>-3.22</td><td>2.75</td><td>7181.57</td><td>-5.96</td><td>900</td><td>26.3</td><td>0.573886</td><td>0.092071</td><td>0.106226</td><td>0.189905</td><td>0.524569</td><td>0.677237</td></tr><tr><td>2171.31709456</td><td>5.22</td><td>2.1</td><td>7199.72</td><td>-2.35</td><td>900</td><td>40</td><td>0.567555</td><td>0.080705</td><td>0.097921</td><td>0.182549</td><td>0.526484</td><td>0.643617</td></tr><tr><td>2172.31291994</td><td>-1.16</td><td>1.78</td><td>7183.03</td><td>-13.1</td><td>900</td><td>51.7</td><td>0.602331</td><td>0.082135</td><td>0.100467</td><td>0.184348</td><td>0.521872</td><td>0.655749</td></tr><tr><td>2189.39389821</td><td>-5.12</td><td>2.18</td><td>7187.35</td><td>0.742</td><td>900</td><td>34.9</td><td>0.553612</td><td>0.087636</td><td>0.106345</td><td>0.186232</td><td>0.524683</td><td>0.636639</td></tr><tr><td>2190.35888875</td><td>-5.56</td><td>2.1</td><td>7181.28</td><td>-1.38</td><td>900</td><td>37.8</td><td>0.553642</td><td>0.08599</td><td>0.10436</td><td>0.185399</td><td>0.52397</td><td>0.649955</td></tr><tr><td>2192.33613311</td><td>-1.58</td><td>1.8</td><td>7190.22</td><td>-4.64</td><td>900</td><td>41.2</td><td>0.594328</td><td>0.086471</td><td>0.099245</td><td>0.182607</td><td>0.521139</td><td>0.661821</td></tr><tr><td>2212.37354384</td><td>7.092</td><td>3.59</td><td>7203.08</td><td>-9.36</td><td>900</td><td>23.9</td><td>0.593524</td><td>0.09262</td><td>0.102766</td><td>0.183774</td><td>0.521448</td><td>0.630614</td></tr><tr><td><strong>2213.42063681</strong></td><td><strong>16.78</strong></td><td><strong>4.42</strong></td><td><strong>7211.7</strong></td><td><strong>3.168</strong></td><td><strong>900</strong></td><td><strong>19</strong></td><td><strong>0.596428</strong></td><td><strong>0.099179</strong></td><td><strong>0.101745</strong></td><td><strong>0.181668</strong></td><td><strong>0.520206</strong></td><td><strong>0.654342</strong></td></tr></tbody></table>

Notes. (*) Duration of each individual exposure. (†) Bisector spans; error bars are twice those of RVs.
<table><caption>Table A.1. continued.</caption><thead><tr><th>BJD<sub>UTC</sub><br/>-2457000 [days]</th><th>RV<br/>[m s<sup>-1</sup>]</th><th>±1σRV</th><th>FWHM</th><th>BIS<sup>(†)</sup></th><th>Exp.<sup>(*)</sup><br/>[sec]</th><th>S/N</th><th>I<sub>Ca</sub><br/>u</th><th>I<sub>Ha06</sub></th><th>I<sub>He</sub><br/>i</th><th>I<sub>Na</sub><br/>i</th><th>I<sub>Ca</sub><br/>u</th><th>I<sub>Ha06</sub></th></tr></thead><tbody><tr><td>2216.40383735</td><td>-8.94</td><td>2.19</td><td>7203.82</td><td>-10.3</td><td>900</td><td>33.9</td><td>0.594565</td><td>0.089359</td><td>0.100295</td><td>0.182115</td><td>0.519823</td><td>0.629985</td></tr><tr><td>2235.36719701</td><td>6.639</td><td>2.74</td><td>7177.77</td><td>0.334</td><td>900</td><td>25.5</td><td>0.595733</td><td>0.091527</td><td>0.101648</td><td>0.182908</td><td>0.520753</td><td>0.626604</td></tr><tr><td>2236.31768047</td><td>5.56</td><td>2.2</td><td>7188.86</td><td>-8.92</td><td>900</td><td>38.7</td><td>0.567033</td><td>0.087864</td><td>0.100285</td><td>0.182334</td><td>0.518575</td><td>0.610068</td></tr><tr><td>2237.32364626</td><td>4.933</td><td>3.21</td><td>7181.9</td><td>-17.9</td><td>1200</td><td>28.2</td><td>0.559051</td><td>0.086238</td><td>0.104913</td><td>0.186654</td><td>0.525845</td><td>0.664309</td></tr><tr><td><strong>2239.31264031</strong></td><td><strong>21.6</strong></td><td><strong>2.85</strong></td><td><strong>7219.49</strong></td><td><strong>-1.592</strong></td><td><strong>900</strong></td><td><strong>31.5</strong></td><td><strong>0.591982</strong></td><td><strong>0.088633</strong></td><td><strong>0.099598</strong></td><td><strong>0.180553</strong></td><td><strong>0.524177</strong></td><td><strong>0.630278</strong></td></tr><tr><td>2240.31312271</td><td>7.368</td><td>2.31</td><td>7208.34</td><td>-2.51</td><td>900</td><td>37.6</td><td>0.594946</td><td>0.081703</td><td>0.10278</td><td>0.183971</td><td>0.522467</td><td>0.662409</td></tr><tr><td>2244.31706817</td><td>6.052</td><td>2.84</td><td>7195.35</td><td>0.352</td><td>900</td><td>32.5</td><td>0.594812</td><td>0.084492</td><td>0.100398</td><td>0.18777</td><td>0.524697</td><td>0.655945</td></tr><tr><td>2245.31980936</td><td>0.237</td><td>2.24</td><td>7209.07</td><td>-5.47</td><td>900</td><td>34.9</td><td>0.589965</td><td>0.087524</td><td>0.103149</td><td>0.183372</td><td>0.519622</td><td>0.663637</td></tr><tr><td><strong>2246.32103509</strong></td><td><strong>27.1</strong></td><td><strong>7.78</strong></td><td><strong>7248.46</strong></td><td><strong>34.51</strong></td><td><strong>1200</strong></td><td><strong>12.7</strong></td><td><strong>0.594927</strong></td><td><strong>0.090211</strong></td><td><strong>0.100742</strong></td><td><strong>0.184404</strong></td><td><strong>0.521228</strong></td><td><strong>0.74003</strong></td></tr><tr><td>2412.68868670</td><td>-3.7</td><td>2.54</td><td>7163.71</td><td>-3.37</td><td>900</td><td>32.8</td><td>0.550281</td><td>0.090724</td><td>0.100293</td><td>0.183955</td><td>0.523314</td><td>0.6709</td></tr><tr><td>2413.66378329</td><td>-1.5</td><td>2.16</td><td>7179.72</td><td>-6.74</td><td>900</td><td>34.4</td><td>0.558588</td><td>0.088191</td><td>0.099293</td><td>0.1792</td><td>0.526702</td><td>0.652927</td></tr><tr><td>2414.31760521</td><td>2.937</td><td>2.15</td><td>7166.9</td><td>4.331</td><td>900</td><td>37.1</td><td>0.555402</td><td>0.084723</td><td>0.098655</td><td>0.181032</td><td>0.530929</td><td>0.662741</td></tr><tr><td>2416.63114629</td><td>2.289</td><td>1.95</td><td>7191.8</td><td>-11.7</td><td>900</td><td>48</td><td>0.584258</td><td>0.083699</td><td>0.099042</td><td>0.180173</td><td>0.526418</td><td>0.68478</td></tr><tr><td>2417.65137234</td><td>3.912</td><td>1.73</td><td>7191.95</td><td>-8.14</td><td>900</td><td>40.8</td><td>0.583166</td><td>0.087257</td><td>0.098779</td><td>0.18156</td><td>0.52857</td><td>0.664774</td></tr><tr><td>2418.65486069</td><td>4.286</td><td>3.33</td><td>7195.12</td><td>-13.2</td><td>900</td><td>21.6</td><td>0.587824</td><td>0.10659</td><td>0.100434</td><td>0.180994</td><td>0.498526</td><td>0.662161</td></tr><tr><td>2427.7318121</td><td>1.618</td><td>2.29</td><td>7191.43</td><td>-0.36</td><td>900</td><td>34.2</td><td>0.595541</td><td>0.089463</td><td>0.100339</td><td>0.183087</td><td>0.530979</td><td>0.649579</td></tr><tr><td>2428.66018094</td><td>0.03</td><td>3.07</td><td>7187.2</td><td>-11.9</td><td>900</td><td>28.1</td><td>0.587265</td><td>0.089537</td><td>0.098082</td><td>0.181868</td><td>0.526651</td><td>0.668995</td></tr><tr><td>2430.66988978</td><td>-0.96</td><td>2.26</td><td>7186.32</td><td>-3.5</td><td>900</td><td>36</td><td>0.558113</td><td>0.085708</td><td>0.102817</td><td>0.184679</td><td>0.529491</td><td>0.668573</td></tr><tr><td>2431.71360583</td><td>-5.64</td><td>2.64</td><td>7190.54</td><td>-22.7</td><td>900</td><td>32.1</td><td>0.550955</td><td>0.07889</td><td>0.100986</td><td>0.181203</td><td>0.527072</td><td>0.693437</td></tr><tr><td>2443.66087057</td><td>8.667</td><td>3.5</td><td>7188.34</td><td>-4.35</td><td>900</td><td>25.2</td><td>0.590039</td><td>0.092999</td><td>0.109664</td><td>0.179042</td><td>0.526005</td><td>0.696157</td></tr><tr><td>2444.57982389</td><td>9.299</td><td>3.32</td><td>7180.37</td><td>-2.39</td><td>900</td><td>26.8</td><td>0.589859</td><td>0.091057</td><td>0.098188</td><td>0.179557</td><td>0.522807</td><td>0.666047</td></tr><tr><td>2445.59632825</td><td>10.81</td><td>1.86</td><td>7178.02</td><td>-8.87</td><td>900</td><td>36.9</td><td>0.557645</td><td>0.085087</td><td>0.10075</td><td>0.181234</td><td>0.52886</td><td>0.678235</td></tr><tr><td>2446.60873425</td><td>5.874</td><td>1.87</td><td>7182.03</td><td>2.719</td><td>900</td><td>45</td><td>0.552972</td><td>0.085681</td><td>0.101935</td><td>0.183638</td><td>0.525628</td><td>0.682197</td></tr><tr><td>2447.61736597</td><td>1.855</td><td>2.1</td><td>7181.21</td><td>-2</td><td>900</td><td>39.9</td><td>0.585974</td><td>0.086557</td><td>0.101268</td><td>0.182021</td><td>0.529423</td><td>0.6903</td></tr><tr><td>2448.59166621</td><td>2.524</td><td>1.79</td><td>7183.26</td><td>-6.53</td><td>900</td><td>46.5</td><td>0.550409</td><td>0.081722</td><td>0.099201</td><td>0.181717</td><td>0.528627</td><td>0.676518</td></tr><tr><td>2449.61666267</td><td>-2.34</td><td>1.85</td><td>7181.85</td><td>-9.02</td><td>900</td><td>46.7</td><td>0.55726</td><td>0.081361</td><td>0.100242</td><td>0.183494</td><td>0.527072</td><td>0.693437</td></tr><tr><td>2450.61247483</td><td>-5.57</td><td>2.17</td><td>7179.71</td><td>-10.5</td><td>900</td><td>38.1</td><td>0.561505</td><td>0.087524</td><td>0.101564</td><td>0.182982</td><td>0.529629</td><td>0.696181</td></tr><tr><td>2454.66732062</td><td>-2.11</td><td>5</td><td>7157.53</td><td>-13.6</td><td>900</td><td>16.4</td><td>0.534792</td><td>0.105417</td><td>0.106316</td><td>0.188685</td><td>0.531533</td><td>0.672621</td></tr><tr><td>2455.68136136</td><td>1.242</td><td>2.34</td><td>7201.87</td><td>-15</td><td>900</td><td>37.6</td><td>0.559782</td><td>0.085698</td><td>0.097922</td><td>0.182038</td><td>0.527513</td><td>0.642527</td></tr><tr><td>2456.70358427</td><td>1.194</td><td>1.71</td><td>7193.3</td><td>-2.59</td><td>900</td><td>39.1</td><td>0.539261</td><td>0.083335</td><td>0.100678</td><td>0.18192</td><td>0.526882</td><td>0.654302</td></tr><tr><td>2457.65444875</td><td>0.232</td><td>1.53</td><td>7186.93</td><td>-4.04</td><td>900</td><td>45.9</td><td>0.544456</td><td>0.082158</td><td>0.103781</td><td>0.185289</td><td>0.52979</td><td>0.678245</td></tr><tr><td>2458.62416725</td><td>1.797</td><td>2.1</td><td>7191.62</td><td>-2.73</td><td>900</td><td>42.5</td><td>0.571849</td><td>0.083688</td><td>0.10082</td><td>0.181891</td><td>0.525159</td><td>0.670374</td></tr><tr><td>2459.58868774</td><td>0</td><td>2.42</td><td>7166.31</td><td>-9.47</td><td>900</td><td>32</td><td>0.534143</td><td>0.084143</td><td>0.101018</td><td>0.183983</td><td>0.531303</td><td>0.671843</td></tr><tr><td>2460.68008744</td><td>2.859</td><td>2</td><td>7199.43</td><td>-10</td><td>900</td><td>39</td><td>0.561505</td><td>0.085291</td><td>0.09952</td><td>0.182178</td><td>0.530929</td><td>0.641539</td></tr><tr><td>2461.61247483</td><td>-3.43</td><td>2.01</td><td>7177.71</td><td>-14.6</td><td>900</td><td>39.8</td><td>0.531283</td><td>0.081882</td><td>0.101482</td><td>0.18371</td><td>0.523212</td><td>0.626268</td></tr><tr><td>2462.56999290</td><td>-1.49</td><td>2.41</td><td>7193.46</td><td>-10.2</td><td>900</td><td>34</td><td>0.53251</td><td>0.091757</td><td>0.1007</td><td>0.182059</td><td>0.53047</td><td>0.678952</td></tr><tr><td>2463.59087617</td><td>5.604</td><td>5.45</td><td>7198.52</td><td>-24.4</td><td>900</td><td>17.8</td><td>0.536789</td><td>0.095036</td><td>0.099161</td><td>0.183787</td><td>0.528055</td><td>0.686225</td></tr><tr><td>2465.55285533</td><td>1.953</td><td>1.95</td><td>7189.61</td><td>-10.6</td><td>900</td><td>37.5</td><td>0.525385</td><td>0.081987</td><td>0.102391</td><td>0.184078</td><td>0.528833</td><td>0.659363</td></tr><tr><td>2472.64437339</td><td>9.088</td><td>3.09</td><td>7206.96</td><td>-15.6</td><td>900</td><td>28.8</td><td>0.582423</td><td>0.088182</td><td>0.102267</td><td>0.18471</td><td>0.526636</td><td>0.671635</td></tr><tr><td>2473.58211801</td><td>-2.46</td><td>3.71</td><td>7182.04</td><td>-22.5</td><td>900</td><td>24</td><td>0.548713</td><td>0.090743</td><td>0.100974</td><td>0.184532</td><td>0.522132</td><td>0.670676</td></tr><tr><td>2475.57083577</td><td>3.505</td><td>2.35</td><td>7184.54</td><td>-16.3</td><td>1800</td><td>32.2</td><td>0.580897</td><td>0.092011</td><td>0.101234</td><td>0.182594</td><td>0.528772</td><td>0.656172</td></tr><tr><td>2476.56282582</td><td>7.827</td><td>1.91</td><td>7184.92</td><td>-6.25</td><td>900</td><td>44</td><td>0.565952</td><td>0.083328</td><td>0.098685</td><td>0.182411</td><td>0.526736</td><td>0.673459</td></tr><tr><td>2477.54247641</td><td>5.777</td><td>1.95</td><td>7190.45</td><td>-13.2</td><td>900</td><td>40.1</td><td>0.584375</td><td>0.084401</td><td>0.105974</td><td>0.188686</td><td>0.528026</td><td>0.684554</td></tr><tr><td>2478.53938306</td><td>3.887</td><td>2.42</td><td>7199.21</td><td>1.209</td><td>900</td><td>33.4</td><td>0.542519</td><td>0.091376</td><td>0.099644</td><td>0.180595</td><td>0.522094</td><td>0.662569</td></tr><tr><td>2479.53060547</td><td>0.796</td><td>2.4</td><td>7200.54</td><td>-6.9</td><td>900</td><td>35.5</td><td>0.578876</td><td>0.086324</td><td>0.099936</td><td>0.181075</td><td>0.527574</td><td>0.654195</td></tr><tr><td>2481.49382164</td><td>3.548</td><td>2.82</td><td>7186.71</td><td>-6.58</td><td>900</td><td>31.2</td><td>0.53899</td><td>0.091801</td><td>0.101227</td><td>0.183502</td><td>0.527366</td><td>0.676158</td></tr><tr><td>2513.44982261</td><td>1.106</td><td>2.07</td><td>7176.01</td><td>-3.44</td><td>900</td><td>41.6</td><td>0.561363</td><td>0.088328</td><td>0.100405</td><td>0.183709</td><td>0.523996</td><td>0.678603</td></tr><tr><td>2513.50218232</td><td>-4.57</td><td>2.69</td><td>7172.76</td><td>-4.86</td><td>900</td><td>33.2</td><td>0.566953</td><td>0.08974</td><td>0.10596</td><td>0.188059</td><td>0.524511</td><td>0.666513</td></tr><tr><td>2515.41917706</td><td>-2.34</td><td>1.93</td><td>7184.92</td><td>-5.97</td><td>900</td><td>39</td><td>0.576153</td><td>0.089897</td><td>0.101247</td><td>0.183188</td><td>0.526739</td><td>0.676359</td></tr><tr><td>2515.44453511</td><td>-3.02</td><td>1.76</td><td>7178.72</td><td>-6.73</td><td>900</td><td>49.2</td><td>0.567076</td><td>0.088259</td><td>0.098628</td><td>0.181831</td><td>0.520853</td><td>0.679953</td></tr><tr><td>2516.53077561</td><td>-4.93</td><td>2.92</td><td>7190.8</td><td>1.972</td><td>900</td><td>31.9</td><td>0.602149</td><td>0.092247</td><td>0.099896</td><td>0.182444</td><td>0.527717</td><td>0.666103</td></tr><tr><td>2563.36649887</td><td>0.845</td><td>1.64</td><td>7187.23</td><td>-3.3</td><td>900</td><td>54.5</td><td>0.596593</td><td>0.088658</td><td>0.098414</td><td>0.180042</td><td>0.522949</td><td>0.675468</td></tr><tr><td>2566.34542939</td><td>-5.47</td><td>2.15</td><td>7179.09</td><td>-6.14</td><td>900</td><td>45.4</td><td>0.601739</td><td>0.091592</td><td>0.099575</td><td>0.18055</td><td>0.525212</td><td>0.678527</td></tr><tr><td>2575.40497038</td><td>3.93</td><td>2.12</td><td>7181.99</td><td>-8.36</td><td>900</td><td>40.1</td><td>0.566264</td><td>0.087701</td><td>0.098051</td><td>0.182169</td><td>0.523518</td><td>0.671318</td></tr><tr><td>2579.38524573</td><td>5.059</td><td>2.18</td><td>7191.98</td><td>-5.62</td><td>900</td><td>38.3</td><td>0.562049</td><td>0.089366</td><td>0.104764</td><td>0.183719</td><td>0.525329</td><td>0.674917</td></tr><tr><td>2580.41266409</td><td>10.09</td><td>2.97</td><td>7187.22</td><td>4.876</td><td>900</td><td>19.9</td><td>0.572011</td><td>0.107006</td><td>0.107628</td><td>0.185394</td><td>0.52275</td><td>0.687077</td></tr><tr><td>2584.33980519</td><td>0.821</td><td>2.95</td><td>7183.34</td><td>-0.433</td><td>900</td><td>29.2</td><td>0.56581</td><td>0.092026</td><td>0.104575</td><td>0.185035</td><td>0.514433</td><td>0.687819</td></tr><tr><td>2588.34483394</td><td>1.991</td><td>5.05</td><td>7181.58</td><td>9.001</td><td>900</td><td>19.2</td><td>0.560173</td><td>0.094073</td><td>0.110324</td><td>0.186892</td><td>0.524911</td><td>0.680451</td></tr><tr><td>2601.31675212</td><td>5.918</td><td>2.03</td><td>7196.44</td><td>-6.76</td><td>900</td><td>45.2</td><td>0.565656</td><td>0.087291</td><td>0.10094</td><td>0.17769</td><td>0.523886</td><td>0.667702</td></tr></tbody></table>

Notes. (*) Duration of each individual exposure. (†) Bisector spans; error bars are twice those of RVs.
## Appendix B: Priors and posteriors

<table><caption>**Table B.1.** Prior volume for the parameters of the one-planet model fit of Sect. 2.3 processed with *juliet*. $\mathcal{U}(a, b)$ indicates a uniform distribution between $a$ and $b$; $\mathcal{L}(a, b)$ a log-normal distribution, $\mathcal{N}(a, b)$ a normal distribution, and $\mathcal{T}(a, b)$ a truncated normal distribution (where lower possible value equals zero) with mean $a$ and standard deviation $b$.</caption><thead><tr><th>Parameter</th><th>Prior distribution</th></tr></thead><tbody><tr><td colspan="2"><strong>Keplerian Parameters:</strong></td></tr><tr><td>ρ<sub>★</sub> [kg/m<sup>3</sup>]</td><td>N(1300, 100)</td></tr><tr><td>T<sub>0,b</sub> [BJD]</td><td>N(2458745.921, 0.003)</td></tr><tr><td>P<sub>b</sub> [days]</td><td>N(12.998, 0.002)</td></tr><tr><td>e<sub>b</sub><sup>*</sup></td><td>0</td></tr><tr><td>ω<sub>b</sub><sup>*</sup></td><td>90</td></tr><tr><td colspan="2"><strong>Transit Parameters:</strong></td></tr><tr><td>R<sub>p</sub>/R<sub>★</sub></td><td>U(0.0, 1.0)</td></tr><tr><td>D</td><td>1.0</td></tr><tr><td>q<sub>1</sub></td><td>N(0.31, 0.30)</td></tr><tr><td>q<sub>2</sub></td><td>N(0.25, 0.10)</td></tr><tr><td colspan="2"><strong>Light curve GP Hyperparameters:</strong></td></tr><tr><td>σ<sub>TESS</sub> [ppt]</td><td>L(10<sup>-3</sup>, 10)</td></tr><tr><td>ρ<sub>TESS</sub> [days]</td><td>L(10<sup>-1</sup>, 10)</td></tr><tr><td colspan="2"><strong>RV parameters:</strong></td></tr><tr><td>K<sub>b</sub> [m/s]</td><td>U(0.0, 10.0)</td></tr><tr><td>σ<sub>HARPS-N</sub> [m s<sup>-1</sup>]</td><td>U(0, 10)</td></tr><tr><td>A [m s<sup>-1</sup> days<sup>-1</sup>]</td><td>U(-1, 1)</td></tr><tr><td>B [m s<sup>-1</sup>]</td><td>U(-20, 20)</td></tr></tbody></table>

**Notes.** (*) In the case of non-null eccentricity, the priors were set as follows: $(\sqrt{e} \sin \omega, \sqrt{e} \cos \omega)$ in $\mathcal{U}(-1.0, 1.0)$.

<table><caption>**Table B.2.** Prior volume for the parameters of the two-planet model fit of Sect. 2.3 processed with *juliet*.</caption><thead><tr><th>Parameter</th><th>Prior distribution</th></tr></thead><tbody><tr><td colspan="2"><strong>Keplerian Parameters:</strong></td></tr><tr><td>ρ<sub>★</sub> [kg/m<sup>3</sup>]</td><td>N(1300, 100)</td></tr><tr><td>T<sub>0,b</sub> [BJD]</td><td>N(2458745.921, 0.003)</td></tr><tr><td>P<sub>b</sub> [days]</td><td>N(12.998, 0.002)</td></tr><tr><td>T<sub>0,c</sub> [BJD]</td><td>N(2458740, 2458790)</td></tr><tr><td>P<sub>c</sub> [days]</td><td>U(1, 100)</td></tr><tr><td>(e<sub>b</sub>, e<sub>c</sub>)<sup>*</sup></td><td>0</td></tr><tr><td>(ω<sub>b</sub>, ω<sub>c</sub>)<sup>*</sup></td><td>90</td></tr><tr><td colspan="2"><strong>Transit Parameters:</strong></td></tr><tr><td>R<sub>p</sub>/R<sub>★</sub></td><td>U(0.0, 1.0)</td></tr><tr><td>D</td><td>1.0</td></tr><tr><td>q<sub>1</sub></td><td>N(0.31, 0.30)</td></tr><tr><td>q<sub>2</sub></td><td>N(0.25, 0.10)</td></tr><tr><td colspan="2"><strong>Light curve GP Hyperparameters:</strong></td></tr><tr><td>σ<sub>TESS</sub> [ppt]</td><td>L(10<sup>-3</sup>, 10)</td></tr><tr><td>ρ<sub>TESS</sub> [days]</td><td>L(10<sup>-1</sup>, 10)</td></tr><tr><td colspan="2"><strong>RV parameters:</strong></td></tr><tr><td>K<sub>b</sub> [m/s]</td><td>U(0.0, 10.0)</td></tr><tr><td>K<sub>c</sub> [m/s]</td><td>U(0, 10)</td></tr><tr><td>σ<sub>HARPS-N</sub> [m s<sup>-1</sup>]</td><td>U(0, 10)</td></tr><tr><td>A [m s<sup>-1</sup> days<sup>-1</sup>]</td><td>U(-1, 1)</td></tr><tr><td>B [m s<sup>-1</sup>]</td><td>U(-20, 20)</td></tr></tbody></table>

**Notes.** (*) In the case of non-null eccentricity, the priors were set as follows: $(\sqrt{e} \sin \omega, \sqrt{e} \cos \omega)$ in $\mathcal{U}(-1.0, 1.0)$.
<table><caption>**Table B.3.** Posterior's result for the parameters of the two-planet eccentric model fit of Sect. 2.3 processed with *juliet*.</caption><thead><tr><th>Parameter</th><th>Value (±1σ)</th></tr></thead><tbody><tr><td colspan="2"><strong>Keplerian Parameters:</strong></td></tr><tr><td>ρ<sub>★</sub> [kg/m<sup>3</sup>]</td><td>1312<sup>+55</sup><sub>-68</sub></td></tr><tr><td>a<sub>b</sub>/R<sub>★</sub></td><td>22.72<sup>+0.31</sup><sub>-0.40</sub></td></tr><tr><td>a<sub>c</sub>/R<sub>★</sub></td><td>39.05<sup>+0.50</sup><sub>-0.73</sub></td></tr><tr><td>T<sub>0,b</sub> [BJD]</td><td>2458745.9205<sup>+0.0012</sup><sub>-0.0011</sub></td></tr><tr><td>P<sub>b</sub> [days]</td><td>12.9972 ± 0.0006</td></tr><tr><td>T<sub>0,c</sub> [BJD]</td><td>2458776.6<sup>+4.5</sup><sub>-4.6</sub></td></tr><tr><td>P<sub>c</sub> [days]</td><td>29.29<sup>+0.21</sup><sub>-0.20</sub></td></tr><tr><td colspan="2"><strong>Transit Parameters:</strong></td></tr><tr><td>R<sub>Pb</sub>/R<sub>★</sub></td><td>0.0356<sup>+0.0007</sup><sub>-0.0005</sub></td></tr><tr><td>q<sub>1</sub></td><td>0.28<sup>+0.11</sup><sub>-0.08</sub></td></tr><tr><td>q<sub>2</sub></td><td>0.30<sup>+0.05</sup><sub>-0.05</sub></td></tr><tr><td>b<sub>b</sub></td><td>0.19<sup>+0.11</sup><sub>-0.10</sub></td></tr><tr><td>i<sub>b</sub> [deg]</td><td>89.52<sup>+0.26</sup><sub>-0.28</sub></td></tr><tr><td colspan="2"><strong>Light curve GP Hyperparameters:</strong></td></tr><tr><td>σ<sub>TESS</sub> [ppt]</td><td>0.19<sup>+0.03</sup><sub>-0.02</sub></td></tr><tr><td>ρ<sub>TESS</sub> [days]</td><td>0.76<sup>+0.19</sup><sub>-0.15</sub></td></tr><tr><td colspan="2"><strong>RV parameters:</strong></td></tr><tr><td>K<sub>b</sub> [m/s]</td><td>2.47<sup>+0.50</sup><sub>-0.46</sub></td></tr><tr><td>K<sub>c</sub> [m/s]</td><td>2.36<sup>+0.42</sup><sub>-0.40</sub></td></tr><tr><td>σ<sub>HARPS-N</sub> [m s<sup>-1</sup>]</td><td>2.93<sup>+0.35</sup><sub>-0.32</sub></td></tr><tr><td>A [m s<sup>-1</sup> days<sup>-1</sup>]</td><td>0.0110 ± 0.0015</td></tr><tr><td>B [m s<sup>-1</sup>]</td><td>-9.1 ± 1.3</td></tr></tbody></table>
## Appendix C: Corner plots

<figure><img src="image_16.png" /><figcaption>**Fig. C.1.** Corner plot for the posterior distribution of the joint transit and RV analysis of Sect. 3.3 in the case of two planets, elaborated with *juliet*.</figcaption></figure>
## Appendix D: Additional plots

<figure><img src="image_17.png" /><figcaption>**Fig. D.1.** FWHM and log $R'_{HK}$ are plotted over time respectively in the upper and lower panel, along with their linear trends (orange line) and average value (dashed grey line).</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Fig. D.2.** Close-up look of the RV GLS periodogram, executed with the publicly available tool Exo-Striker (Trifonov 2019; https://github.com/3fon3fonov/exostriker) after the removal of a linear trend. The two vertical blue lines, around the 29-day signal (indicated by a vertical yellow line), show the main peak aliases due to the two highest frequencies of the window function, in the upper and bottom panels. The three horizontal dotted lines represent the 10%, 1%, and 0.1% FAP levels.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Fig. D.3.** Window function of the HARPS-N RV measurements, as evaluated with Exo-Striker. The two highest peaks, excluding the 1-day peak and frequencies close to zero, are indicated by the respective labels.</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>**Fig. D.4.** Unnormalized GLS power for a different number of HARPS-N observations. The power of the 29-day signal increases with more observations. The vertical dashed red and green lines indicate TOI-1422 b and TOI-1422 c orbital periods, respectively, while the horizontal dashed lines signal the 10% and 1% confidence levels, respectively (evaluated with the bootstrap method).</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>**Fig. D.5.** TOI-1422 b transits, as seen with PATHOS, folded on the first row of the left column and the single transit event on the right one, with X/Y and the sky background in the following rows, showing no correlation with the transits.</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Fig. D.6.** Single transit depth from PATHOS in different apertures, with the three rows showing the transit depth at an aperture radius of 2, 3, and 4 pixels, respectively.</figcaption></figure>