# TOI-663: A newly discovered multi-planet system with three transiting mini-Neptunes orbiting an early M star

M. Cointepas, F. Bouchy, J. M. Almenara, X. Bonfils, N. Astudillo-Defru, H.

Knierim, M. Stalport, L. Mignon, N. Grieves, J. Bean, et al.

► To cite this version:

M. Cointepas, F. Bouchy, J. M. Almenara, X. Bonfils, N. Astudillo-Defru, et al., TOI-663: A newly discovered multi-planet system with three transiting mini-Neptunes orbiting an early M star. Astronomy and Astrophysics - A&A, 2024, 685, pp.A19. 10.1051/0004-6361/202346899. hal-04564852

HAL Id: hal-04564852

https://hal.science/hal-04564852

Submitted on 30 Apr 2024

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

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

<figure><img src="image_1.png" /></figure>
# TOI-663: A newly discovered multi-planet system with three transiting mini-Neptunes orbiting an early M star

M. Cointepas, F. Bouchy, J. M. Almenara, X. Bonfils, N. Astudillo-Defru, H.

Knierim, M. Stalport, L. Mignon, N. Grieves, J. Bean, et al.

► To cite this version:

M. Cointepas, F. Bouchy, J. M. Almenara, X. Bonfils, N. Astudillo-Defru, et al., TOI-663: A newly discovered multi-planet system with three transiting mini-Neptunes orbiting an early M star. Astronomy & Astrophysics - A&A, 2024, 685, pp.A19. 10.1051/0004-6361/202346899. hal-04564852

HAL Id: hal-04564852

https://hal.science/hal-04564852v1

Submitted on 30 Apr 2024

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# TOI-663: A newly discovered multi-planet system with three transiting mini-Neptunes orbiting an early M star

M. Cointepas¹,², F. Bouchy², J. M. Almenara¹©, X. Bonfils¹©, N. Astudillo-Defru³©,
H. Knierim⁷©, M. Stalport²², L. Mignon², N. Grieves², J. Bean³²©, M. Brady³²©, J. Burt¹⁷©,
B. L. Canto Martins⁹©, K. A. Collins⁴©, K. I. Collins⁵, X. Delfosse¹, J. R. de Medeiros⁹©,
B.-O. Demory³⁴, C. Dorn⁸, T. Forveille¹©, A. Fukui²⁴,¹⁰©, T. Gan⁶©, Y. Gómez Maqueo Chew³⁵©, S. Halverson¹⁷©,
R. Helled⁷©, I. Helm⁵©, T. Hirano²⁸,²⁹©, K. Horne²¹©, S. B. Howell¹⁶©, K. Isogai²⁵,²⁶©, D. Kasper³²©,
K. Kawauchi²⁷,²⁶©, J. H. Livingston²⁸,²⁹,³⁰©, B. Massey¹⁹©, R. A. Matson¹⁸©, F. Murgas¹⁰,¹¹©, N. Narita²⁴,²⁸,¹⁰©,
E. Palle¹⁰,¹¹©, H. M. Relles⁴©, L. Sabin³⁶, N. Schanche³⁸, R. P. Schwarz⁴©, A. Seifahrt³²©, A. Shporer¹²©,
G. Stefansson¹⁵,∗¹⁶©, J. Sturmer³³, M. Tamura³¹,²⁸,²⁹©, T.-G. Tan²³©, J. D. Twicken³⁷,¹⁶©, N. Watanabe²⁶©,
R. D. Wells³⁴, F. P. Wilkin²⁰©, G. R. Ricker¹²©, S. Seager¹²,¹³,¹⁴©, J. N. Winn¹⁵©, and J. M. Jenkins¹⁶©

*(Affiliations can be found after the references)*

Received 14 May 2023 / Accepted 2 November 2023

## ABSTRACT

We present the detection of three exoplanets orbiting the early M dwarf TOI-663 (TIC 54962195; V = 13.7 mag, J = 10.4 mag, $R_*$ = 0.512 ± 0.015 $R_⊕$, $M_*$ = 0.514 ± 0.012 $M_⊕$, d = 64 pc). TOI-663 b, c, and d, with respective radii of 2.27 ± 0.10 $R_⊕$, 2.26 ± 0.10 $R_⊕$, and 1.92 ± 0.13 $R_⊕$ and masses of 4.45 ± 0.65 $M_⊕$, 3.65 ± 0.97 $M_⊕$, and <5.2 $M_⊕$ at 99%, are located just above the radius valley that separates rocky and volatile-rich exoplanets. The planet candidates are identified in two TESS sectors and are validated with ground-based photometric follow-up, precise radial-velocity measurements, and high-resolution imaging. We used the software package *juliet* to jointly model the photometric and radial-velocity datasets, with Gaussian processes applied to correct for systematics. The three planets discovered in the TOI-663 system are low-mass mini-Neptunes with radii significantly larger than those of rocky analogs, implying that volatiles, such as water, must predominate. In addition to this internal structure analysis, we also performed a dynamical analysis that confirmed the stability of the system. The three exoplanets in the TOI-663 system, similarly to other sub-Neptunes orbiting M dwarfs, have been found to have lower densities than planets of similar sizes orbiting stars of different spectral types.

**Key words.** techniques: photometric – techniques: radial velocities – planets and satellites: detection – stars: low-mass

## 1. Introduction

Small exoplanets discovered around other stars could be similar or very different to those of our Solar System in terms of composition, formation history, and atmospheric properties. They are also difficult to study due to the small transit and radial-velocity (RV) signals they produce. But the improved sensitivity of recent surveys, such as that conducted with the Transiting Exoplanet Survey Satellite (TESS), have enabled the detection of smaller exoplanets. TESS was launched in 2018 to conduct an all-sky survey and discover transiting planets around the nearest and brightest stars (Ricker et al. 2015).

When looking for small and cool exoplanets, low-mass stars have distinct advantages and are particularly interesting. The transit depth and RV semi-amplitude caused by small exoplanets orbiting M dwarfs are much greater than those caused by similar planets orbiting larger stars of earlier spectral types, making such planets easier to detect. More importantly, such systems are excellent candidates for atmospheric characterization via transmission or thermal emission spectroscopy (e.g., Kempton et al. 2018; Batalha et al. 2018). To determine which planets are predominantly rocky or volatile-rich and

whether they could have a significant atmosphere, precise radius measurements from transit surveys must be combined with mass measurements from dynamical observations, such as precise RV measurements. Exoplanets with precise mass, radius, and bulk density measurements orbiting M dwarfs are thus important objects for better understanding the so-called radius valley separating the population of super-Earths and mini-Neptunes (e.g., Fulton et al. 2017; Mayo et al. 2018; Cloutier & Menou 2020; Hardegree-Ullman et al. 2020).

The radius valley around M dwarfs is located at slightly smaller radii than around Sun-like stars (Fulton et al. 2017; Cloutier & Menou 2020) and marks the transition between rocky super-Earths and mini-Neptunes that host H/He envelopes (Owen & Wu 2017). A variety of physical mechanisms have been proposed to explain the emergence of the radius valley, including photo-evaporation, core-powered mass loss (Jin & Mordasini 2018; Gupta & Schlichting 2020; Wyatt et al. 2020), or a direct outcome of the planet formation process wherein super-Earths form at later times in either a gas-poor or completely gas-depleted environment (Lee & Connors 2021; Lopez & Rice 2018). Each of these mechanisms predicts a unique dependence of the radius valley on the orbital period such that measuring the radius valley's slope provides direct insight into the origin

<footnote>* NASA Sagan Fellow.</footnote>
# TOI-663: A newly discovered multi-planet system with three transiting mini-Neptunes orbiting an early M star

M. Cointepas<sup>1,2</sup>, F. Bouchy<sup>2</sup>, J. M. Almenara<sup>1</sup>, X. Bonfils<sup>1</sup>, N. Astudillo-Defru<sup>3</sup>,
H. Knierim<sup>7</sup>, M. Stalport<sup>22</sup>, L. Mignon<sup>2</sup>, N. Grieves<sup>2</sup>, J. Bean<sup>32</sup>, M. Brady<sup>32</sup>, J. Burt<sup>17</sup>,
B. L. Canto Martins<sup>9</sup>, K. A. Collins<sup>4</sup>, K. I. Collins<sup>5</sup>, X. Delfosse<sup>1</sup>, J. R. de Medeiros<sup>9</sup>,
B.-O. Demory<sup>34</sup>, C. Dorn<sup>8</sup>, T. Forveille<sup>1</sup>, A. Fukui<sup>24,10</sup>, T. Gan<sup>6</sup>, Y. Gómez Maqueo Chew<sup>35</sup>, S. Halverson<sup>17</sup>,
R. Helled<sup>7</sup>, I. Helm<sup>5</sup>, T. Hirano<sup>28,29</sup>, K. Horne<sup>21</sup>, S. B. Howell<sup>16</sup>, K. Isogai<sup>25,26</sup>, D. Kasper<sup>32</sup>,
K. Kawauchi<sup>27,26</sup>, J. H. Livingston<sup>28,29,30</sup>, B. Massey<sup>19</sup>, R. A. Matson<sup>18</sup>, F. Murgas<sup>10,11</sup>, N. Narita<sup>24,28,10</sup>,
E. Palle<sup>10,11</sup>, H. M. Relles<sup>4</sup>, L. Sabin<sup>36</sup>, N. Schanche<sup>38</sup>, R. P. Schwarz<sup>4</sup>, A. Seifahrt<sup>32</sup>, A. Shporer<sup>12</sup>,
G. Stefansson<sup>15,*</sup>, J. Sturmer<sup>33</sup>, M. Tamura<sup>31,28,29</sup>, T.-G. Tan<sup>23</sup>, J. D. Twicken<sup>37,16</sup>, N. Watanabe<sup>26</sup>,
R. D. Wells<sup>34</sup>, F. P. Wilkin<sup>20</sup>, G. R. Ricker<sup>12</sup>, S. Seager<sup>12,13,14</sup>, J. N. Winn<sup>15</sup>, and J. M. Jenkins<sup>16</sup>

*(Affiliations can be found after the references)*

Received 14 May 2023 / Accepted 2 November 2023

## ABSTRACT

We present the detection of three exoplanets orbiting the early M dwarf TOI-663 (TIC 54962195; $V = 13.7$ mag, $J = 10.4$ mag, $R_\star = 0.512 \pm 0.015 R_⊕$, $M_\star = 0.514 \pm 0.012 M_⊕$, $d = 64$ pc). TOI-663 b, c, and d, with respective radii of $2.27 \pm 0.10 R_⊕$, $2.26 \pm 0.10 R_⊕$, and $1.92 \pm 0.13 R_⊕$ and masses of $4.45 \pm 0.65 M_⊕$, $3.65 \pm 0.97 M_⊕$, and $<5.2 M_⊕$ at 99%, are located just above the radius valley that separates rocky and volatile-rich exoplanets. The planet candidates are identified in two TESS sectors and are validated with ground-based photometric follow-up, precise radial-velocity measurements, and high-resolution imaging. We used the software package *juliet* to jointly model the photometric and radial-velocity datasets, with Gaussian processes applied to correct for systematics. The three planets discovered in the TOI-663 system are low-mass mini-Neptunes with radii significantly larger than those of rocky analogs, implying that volatiles, such as water, must predominate. In addition to this internal structure analysis, we also performed a dynamical analysis that confirmed the stability of the system. The three exoplanets in the TOI-663 system, similarly to other sub-Neptunes orbiting M dwarfs, have been found to have lower densities than planets of similar sizes orbiting stars of different spectral types.

**Key words.** techniques: photometric – techniques: radial velocities – planets and satellites: detection – stars: low-mass

## 1. Introduction

Small exoplanets discovered around other stars could be similar or very different to those of our Solar System in terms of composition, formation history, and atmospheric properties. They are also difficult to study due to the small transit and radial-velocity (RV) signals they produce. But the improved sensitivity of recent surveys, such as that conducted with the Transiting Exoplanet Survey Satellite (TESS), have enabled the detection of smaller exoplanets. TESS was launched in 2018 to conduct an all-sky survey and discover transiting planets around the nearest and brightest stars (Ricker et al. 2015).

When looking for small and cool exoplanets, low-mass stars have distinct advantages and are particularly interesting. The transit depth and RV semi-amplitude caused by small exoplanets orbiting M dwarfs are much greater than those caused by similar planets orbiting larger stars of earlier spectral types, making such planets easier to detect. More importantly, such systems are excellent candidates for atmospheric characterization via transmission or thermal emission spectroscopy (e.g., Kempton et al. 2018; Batalha et al. 2018). To determine which planets are predominantly rocky or volatile-rich and

whether they could have a significant atmosphere, precise radius measurements from transit surveys must be combined with mass measurements from dynamical observations, such as precise RV measurements. Exoplanets with precise mass, radius, and bulk density measurements orbiting M dwarfs are thus important objects for better understanding the so-called radius valley separating the population of super-Earths and mini-Neptunes (e.g., Fulton et al. 2017; Mayo et al. 2018; Cloutier & Menou 2020; Hardegree-Ullman et al. 2020).

The radius valley around M dwarfs is located at slightly smaller radii than around Sun-like stars (Fulton et al. 2017; Cloutier & Menou 2020) and marks the transition between rocky super-Earths and mini-Neptunes that host H/He envelopes (Owen & Wu 2017). A variety of physical mechanisms have been proposed to explain the emergence of the radius valley, including photo-evaporation, core-powered mass loss (Jin & Mordasini 2018; Gupta & Schlichting 2020; Wyatt et al. 2020), or a direct outcome of the planet formation process wherein super-Earths form at later times in either a gas-poor or completely gas-depleted environment (Lee & Connors 2021; Lopez & Rice 2018). Each of these mechanisms predicts a unique dependence of the radius valley on the orbital period such that measuring the radius valley's slope provides direct insight into the origin

<footnote>* NASA Sagan Fellow.</footnote>
of the radius valley. Cloutier & Menou (2020) discovered that the radius valley separating rocky and gaseous planets exhibits a negative slope with insolation around low-mass stars, as opposed to the positive slope found around solar-type stars. This supports models of direct terrestrial planet formation in a gas-poor environment. The different slopes of models of gas-depleted formation and thermally driven mass loss naturally carve out a subset of the period-radius parameter space, within which the two models make opposing predictions about the bulk composition of planets. Characterizing the bulk composition of transiting planets allows us to determine whether the M dwarf radius valley originated directly from planet formation or as a result of post-formation atmospheric escape.

Multiple transiting planet systems are particularly interesting: they provide a unique opportunity for direct comparative planetology (Millholland et al. 2017) because the exoplanets formed within the same protoplanetary disk and evolved around the same host star. These systems also allow studies of planet-planet interactions (Barros et al. 2015) and planetary formation and migration (Delisle 2017). The discovery and accurate characterization of a system with multiple transiting planets thus represents a significant addition to the current knowledge about planetary systems formation. Furthermore, transiting multi-planet systems may contain additional long-period planets not yet detected. Deep searches to find additional exoplanets have already been conducted for a few stars, and both photometry and RV follow-up of transiting systems have successfully detected additional transiting exoplanets (e.g., the Trappist-1 system; Gillon et al. 2017; Luger et al. 2017 and the LHS 1140 system; Lillo-Box et al. 2020).

Here we present the discovery of three new exoplanets orbiting around the same M1-type star. This discovery was made using TESS, multiple ground-based photometry facilities, including ExTrA and the Las Cumbres Observatory (LCO) network, and precise RV measurements from the High Accuracy Radial velocity Planet Searcher (HARPS), MAROON-X, the InfraRed Doppler (IRD), and the Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO). Section 2 presents an analysis of the stellar properties of TOI-663. Section 3 describes the observations and data used in this study, including photometry from space and from the ground, RV measurements, and high-resolution imaging. Section 4 presents the global analysis of the available data, which is used to constrain the planetary properties. Section 5 presents the discussion and conclusion.

## 2. Stellar parameters

TOI-663 is an M-type star located at a distance of 64.23 ± 0.18 pc (Gaia Collaboration 2018, Lindegren et al. 2018, and Bailer-Jones et al. 2018). The astrometry, photometry, and stellar parameters are reported in Table 1.

TOI-663 mass and radius are derived from the stellar parallax and $K_s$-band magnitude, used to compute the absolute $K_s$-band magnitude $M_{K_s}$, and the empirically derived M dwarf mass-luminosity and radius-luminosity relations from Mann et al. (2019) and Mann et al. (2015), respectively. The analysis by Mann includes stars with close-to-solar metallicities. We used SpecMatch-Emp<fnref n="1" /> (Yee et al. 2017) and HARPS high-resolution spectra of TOI-663 to estimate the stellar metallicity, for which we found a value around 0.07 ± 0.12 dex (but the five best matches are spread over ~0.5 dex so the metallicity remains

<table><caption>**Table 1.** Stellar parameters for TOI-663 (TIC 54962195, UCAC4 408-053407, 2MASS J10401596-0830385, WISE J104015.85-083039.4, and APASS 3855837).</caption><thead><tr><th>Parameter</th><th>Value</th><th>Refs</th></tr></thead><tbody><tr><td><strong>Astrometry</strong></td><td></td><td></td></tr><tr><td>Right ascension (J2015.5), α</td><td>10<sup>h</sup>40<sup>m</sup>15.8<sup>s</sup></td><td>1,2</td></tr><tr><td>Declination (J2015.5), δ</td><td>-08°30′39.9″</td><td>1,2</td></tr><tr><td>Parallax, (mas)</td><td>15.639 ± 0.023</td><td>1,2</td></tr><tr><td>Distance, d (pc)</td><td>64.23 ± 0.18</td><td>1,2</td></tr><tr><td>Proper motion RA (mas yr<sup>-1</sup>)</td><td>-141.994 ± 0.030</td><td>1,2</td></tr><tr><td>Proper motion D (mas yr<sup>-1</sup>)</td><td>-83.131 ± 0.029</td><td>1,2</td></tr><tr><td><strong>Photometry</strong></td><td></td><td></td></tr><tr><td>V (mag)</td><td>13.667 ± 0.023</td><td>3</td></tr><tr><td>TESS magnitude (mag)</td><td>11.7602 ± 0.0074</td><td>3</td></tr><tr><td>J (mag)</td><td>10.448 ± 0.026</td><td>4</td></tr><tr><td>H (mag)</td><td>9.828 ± 0.022</td><td>4</td></tr><tr><td>K<sub>s</sub> (mag)</td><td>9.607 ± 0.021</td><td>4</td></tr><tr><td><strong>Stellar parameters</strong></td><td></td><td></td></tr><tr><td>Spectral type</td><td>M1V</td><td>5</td></tr><tr><td>M<sub>Ks</sub> (mag)</td><td>5.571 ± 0.022</td><td>8</td></tr><tr><td>Effective temperature, T<sub>eff</sub> (K)</td><td>3681 ± 70</td><td>7,8</td></tr><tr><td>Surface gravity, log g (cgs)</td><td>4.7266 ± 0.0085</td><td>3</td></tr><tr><td>Metallicity, (Fe/H) (dex)</td><td>0.07 ± 0.12</td><td>7,8</td></tr><tr><td>Stellar radius, R<sub>⋆</sub> (R<sub>☉</sub>)</td><td>0.512 ± 0.015</td><td>1,6,8</td></tr><tr><td>Stellar mass, M<sub>⋆</sub> (M<sub>☉</sub>)</td><td>0.514 ± 0.012</td><td>1,6,8</td></tr><tr><td>log R'<sub>HK</sub></td><td>-4.81 ± 0.15</td><td>8</td></tr></tbody></table>

**References.** (1) Gaia Collaboration (2018), (2) Lindegren et al. (2018), (3) Stassun et al. (2019), (4) Cutri et al. (2003), (5) Pecaut & Mamajek (2013), (6) Mann et al. (2015, 2019), (7) Yee et al. (2017), (8) this work.

poorly constrained), confirmed the stellar radius and obtained a more precise effective temperature than the one computed from the Mann et al. (2015, 2019) empirical relations.

We determine the stellar radius of TOI-663 independently by modeling the spectral energy distribution (SED) with stellar atmosphere and evolution models. To construct the SED, we used the magnitudes from Gaia Data Release 2 (DR2; Evans et al. 2018; Maíz Apellániz & Weiler 2018), the 2-Micron All-Sky Survey (2MASS; Skrutskie et al. 2006; Cutri et al. 2003), and Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010; Cutri & et al. 2013). We modeled these magnitudes using the procedure described in Díaz et al. (2014). We used informative priors for the effective temperature ($T_{eff}$), and metallicity ([Fe/H]) from the analysis of the HARPS co-added spectra with SpecMatch-Emp (Yee et al. 2017), and for the distance from Gaia DR2 (Gaia Collaboration 2018; Bailer-Jones et al. 2018). Non-informative priors were used for the rest of parameters. We used the PHOENIX/BT-Settl stellar atmosphere models (Allard et al. 2012). The data with the maximum a posteriori stellar atmosphere model is shown in Fig. 1. We obtained a stellar radius of 0.504 ± 0.010 $R_☉$, consistent with the estimation from Mann et al. (2015) presented in Table 1. We decided to use the stellar radius estimated from the empirical relation for the rest of the paper.

We tried to determine the stellar rotation period of TOI-663. We studied activity indices from the HARPS spectra to see if there was some periodicity in the $S_{index}$ of the sodium lines or in the $S_{index}$ of the $H\alpha$ line. After removing any long-term variation, we did not find a significant peak in the periodograms. The flux around the calcium lines is very low in the HARPS

<footnote n="1">¹ https://github.com/samuelyeewl/specmatch-emp</footnote>
of the radius valley. Cloutier & Menou (2020) discovered that the radius valley separating rocky and gaseous planets exhibits a negative slope with insolation around low-mass stars, as opposed to the positive slope found around solar-type stars. This supports models of direct terrestrial planet formation in a gas-poor environment. The different slopes of models of gas-depleted formation and thermally driven mass loss naturally carve out a subset of the period-radius parameter space, within which the two models make opposing predictions about the bulk composition of planets. Characterizing the bulk composition of transiting planets allows us to determine whether the M dwarf radius valley originated directly from planet formation or as a result of post-formation atmospheric escape.

Multiple transiting planet systems are particularly interesting: they provide a unique opportunity for direct comparative planetology (Millholland et al. 2017) because the exoplanets formed within the same protoplanetary disk and evolved around the same host star. These systems also allow studies of planet-planet interactions (Barros et al. 2015) and planetary formation and migration (Delisle 2017). The discovery and accurate characterization of a system with multiple transiting planets thus represents a significant addition to the current knowledge about planetary systems formation. Furthermore, transiting multi-planet systems may contain additional long-period planets not yet detected. Deep searches to find additional exoplanets have already been conducted for a few stars, and both photometry and RV follow-up of transiting systems have successfully detected additional transiting exoplanets (e.g., the Trappist-1 system; Gillon et al. 2017; Luger et al. 2017 and the LHS 1140 system; Lillo-Box et al. 2020).

Here we present the discovery of three new exoplanets orbiting around the same M1-type star. This discovery was made using TESS, multiple ground-based photometry facilities, including ExTrA and the Las Cumbres Observatory (LCO) network, and precise RV measurements from the High Accuracy Radial velocity Planet Searcher (HARPS), MAROON-X, the InfraRed Doppler (IRD), and the Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO). Section 2 presents an analysis of the stellar properties of TOI-663. Section 3 describes the observations and data used in this study, including photometry from space and from the ground, RV measurements, and high-resolution imaging. Section 4 presents the global analysis of the available data, which is used to constrain the planetary properties. Section 5 presents the discussion and conclusion.

## 2. Stellar parameters

TOI-663 is an M-type star located at a distance of 64.23 ± 0.18 pc (Gaia Collaboration 2018, Lindegren et al. 2018, and Bailer-Jones et al. 2018). The astrometry, photometry, and stellar parameters are reported in Table 1.

TOI-663 mass and radius are derived from the stellar parallax and $K_s$-band magnitude, used to compute the absolute $K_s$-band magnitude $M_{K_s}$, and the empirically derived M dwarf mass-luminosity and radius-luminosity relations from Mann et al. (2019) and Mann et al. (2015), respectively. The analysis by Mann includes stars with close-to-solar metallicities. We used SpecMatch-Emp<fnref n="1" /> (Yee et al. 2017) and HARPS high-resolution spectra of TOI-663 to estimate the stellar metallicity, for which we found a value around 0.07 ± 0.12 dex (but the five best matches are spread over ~0.5 dex so the metallicity remains

<table><caption>**Table 1.** Stellar parameters for TOI-663 (TIC 54962195, UCAC4 408-053407, 2MASS J10401596-0830385, WISE J104015.85-083039.4, and APASS 3855837).</caption><thead><tr><th>Parameter</th><th>Value</th><th>Refs</th></tr></thead><tbody><tr><td><strong>Astrometry</strong></td><td></td><td></td></tr><tr><td>Right ascension (J2015.5), α</td><td>10<sup>h</sup>40<sup>m</sup>15.8<sup>s</sup></td><td>1,2</td></tr><tr><td>Declination (J2015.5), δ</td><td>-08°30′39.9″</td><td>1,2</td></tr><tr><td>Parallax, (mas)</td><td>15.639 ± 0.023</td><td>1,2</td></tr><tr><td>Distance, d (pc)</td><td>64.23 ± 0.18</td><td>1,2</td></tr><tr><td>Proper motion RA (mas yr<sup>-1</sup>)</td><td>-141.994 ± 0.030</td><td>1,2</td></tr><tr><td>Proper motion D (mas yr<sup>-1</sup>)</td><td>-83.131 ± 0.029</td><td>1,2</td></tr><tr><td><strong>Photometry</strong></td><td></td><td></td></tr><tr><td>V (mag)</td><td>13.667 ± 0.023</td><td>3</td></tr><tr><td>TESS magnitude (mag)</td><td>11.7602 ± 0.0074</td><td>3</td></tr><tr><td>J (mag)</td><td>10.448 ± 0.026</td><td>4</td></tr><tr><td>H (mag)</td><td>9.828 ± 0.022</td><td>4</td></tr><tr><td>K<sub>s</sub> (mag)</td><td>9.607 ± 0.021</td><td>4</td></tr><tr><td><strong>Stellar parameters</strong></td><td></td><td></td></tr><tr><td>Spectral type</td><td>M1V</td><td>5</td></tr><tr><td>M<sub>Ks</sub> (mag)</td><td>5.571 ± 0.022</td><td>8</td></tr><tr><td>Effective temperature, T<sub>eff</sub> (K)</td><td>3681 ± 70</td><td>7,8</td></tr><tr><td>Surface gravity, log g (cgs)</td><td>4.7266 ± 0.0085</td><td>3</td></tr><tr><td>Metallicity, (Fe/H) (dex)</td><td>0.07 ± 0.12</td><td>7,8</td></tr><tr><td>Stellar radius, R<sub>⋆</sub> (R<sub>☉</sub>)</td><td>0.512 ± 0.015</td><td>1,6,8</td></tr><tr><td>Stellar mass, M<sub>⋆</sub> (M<sub>☉</sub>)</td><td>0.514 ± 0.012</td><td>1,6,8</td></tr><tr><td>log R'<sub>HK</sub></td><td>-4.81 ± 0.15</td><td>8</td></tr></tbody></table>

**References.** (1) Gaia Collaboration (2018), (2) Lindegren et al. (2018), (3) Stassun et al. (2019), (4) Cutri et al. (2003), (5) Pecaut & Mamajek (2013), (6) Mann et al. (2015, 2019), (7) Yee et al. (2017), (8) this work.

poorly constrained), confirmed the stellar radius and obtained a more precise effective temperature than the one computed from the Mann et al. (2015, 2019) empirical relations.

We determine the stellar radius of TOI-663 independently by modeling the spectral energy distribution (SED) with stellar atmosphere and evolution models. To construct the SED, we used the magnitudes from Gaia Data Release 2 (DR2; Evans et al. 2018; Maíz Apellániz & Weiler 2018), the 2-Micron All-Sky Survey (2MASS; Skrutskie et al. 2006; Cutri et al. 2003), and Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010; Cutri & et al. 2013). We modeled these magnitudes using the procedure described in Díaz et al. (2014). We used informative priors for the effective temperature ($T_{eff}$), and metallicity ([Fe/H]) from the analysis of the HARPS co-added spectra with SpecMatch-Emp (Yee et al. 2017), and for the distance from Gaia DR2 (Gaia Collaboration 2018; Bailer-Jones et al. 2018). Non-informative priors were used for the rest of parameters. We used the PHOENIX/BT-Settl stellar atmosphere models (Allard et al. 2012). The data with the maximum a posteriori stellar atmosphere model is shown in Fig. 1. We obtained a stellar radius of 0.504 ± 0.010 $R_{⊙}$, consistent with the estimation from Mann et al. (2015) presented in Table 1. We decided to use the stellar radius estimated from the empirical relation for the rest of the paper.

We tried to determine the stellar rotation period of TOI-663. We studied activity indices from the HARPS spectra to see if there was some periodicity in the $S_{index}$ of the sodium lines or in the $S_{index}$ of the $H\alpha$ line. After removing any long-term variation, we did not find a significant peak in the periodograms. The flux around the calcium lines is very low in the HARPS

<footnote n="1">¹ https://github.com/samuelyeewl/specmatch-emp</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** SED of TOI-663. The solid line is the maximum a posteriori PHOENIX/BT-Settl interpolated synthetic spectrum, red circles are the absolute photometric observations, and gray open circles are the result of integrating the synthetic spectrum in the observed bandpasses.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Periodogram of the ASAS and ZTF photometry and HARPS RV measurements.</figcaption></figure>

spectra but using the ESPRESSO spectra, we were able to measure the value of $\log R'_{HK} = -4.783 \pm 0.063$, derived from the calcium doublet of the ESPRESSO combined spectra. This indicates that the star has likely a rotation period estimated to $P_{rot} = 30 \pm 3.1$ days (Astudillo-Defru et al. 2017a). According to this rotation period and the stellar radius, we expect the $v\sin(i)$ to be less than a km s$^{-1}$, which is consistent with the values of full width at half maximum (FWHM) of the cross-correlation function (CCF) from the ESPRESSO spectra. No sign of periodicity were found in the All Sky Automated Survey (ASAS; Pojmanski 2002) and the Zwicky Transient Facility (ZTF; Bellm et al. 2019) photometry as well as in the HARPS dataset (see Fig. 2). Unfortunately, this difficulty to find the rotation period of the star prevented us from determining the age of the system using gyrochronology.

# 3. Observations

## 3.1. TESS photometry

TOI-663 (TIC 54962195) was observed in two TESS sectors (9 and 35), one during the primary mission (February–March 2019) and one during the first extended mission (February–March 2021) with a two-minute cadence. The TESS Science Processing Operations Center (SPOC; Jenkins et al. 2016) at

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** TESS target pixel file image of TOI-663 in Sector 9 (created with tfpplotter; Aller et al. 2020). The electron counts are color-coded. The red bordered pixels are used in the SPOC SAP. The size of the red circles indicates the Gaia DR2 magnitudes of all nearby stars.</figcaption></figure>

NASA Ames processed the TESS photometric data and the resulting Presearch Data Conditioning Simple Aperture Photometry (PDCSAP; Smith et al. 2012; Stumpe et al. 2012, 2014) light curve of TOI-663 was corrected for dilution in the TESS aperture by known contaminating sources. It is critical to ensure that no visually close-by targets that could affect the depth of the transit are present in the 21'' TESS pixel and to check for contaminating eclipsing binaries. Figure 3 shows a plot of the target pixel file and the aperture mask that is used for the simple aperture photometry (SAP). We can see that one star overlaps the TESS aperture but its faintness results in minimal dilution of the TESS light curve. The SPOC conducted a transit search of the sector 9 light curve on 2019 April 25 and the combined light curve for sectors 9 and 35 on 2021 May 27 with an adaptive, noise-compensating matched filter (Jenkins 2002; Jenkins et al. 2010, 2020). The sector 9 search produced threshold crossing events with orbital periods of 2.60d and 4.70d; the combined search produced a third threshold crossing event with period of 7.10d. An initial limb-darkened transit model was fitted (Li et al. 2019) and a suite of diagnostic tests were conducted to help make or break the planetary nature of each signal (Twicken et al. 2018). The transit signatures passed all diagnostic tests presented in the SPOC Data Validation reports, and all TIC objects other than the target star were excluded as the transit source in each case. The TESS Science Office reviewed the vetting information and issued alerts for TOI-663.01 and TOI-663.02 on 2019 May 07, and TOI-663.03 on 2021 July 08 (Guerrero et al. 2021).

Canto Martins et al. (2020) conducted a search for rotation period in TESS TOIs using fast Fourier transform, Lomb-Scargle, and wavelet techniques, accompanied by a rigorous visual inspection, and found that TOI-663 does not show any sign of chromospheric activity throughout the observed time ranges.

## 3.2. Ground-based photometry

The TESS pixel scale is ~21" pixel$^{-1}$ and photometric apertures typically extend out to roughly 1 arc-minute, generally causing multiple stars to blend in the TESS aperture. To definitionally exclude the presence of another star causing the signal
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** SED of TOI-663. The solid line is the maximum a posteriori PHOENIX/BT-Settl interpolated synthetic spectrum, red circles are the absolute photometric observations, and gray open circles are the result of integrating the synthetic spectrum in the observed bandpasses.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Periodogram of the ASAS and ZTF photometry and HARPS RV measurements.</figcaption></figure>

spectra but using the ESPRESSO spectra, we were able to measure the value of $\log R'_{HK} = -4.783 \pm 0.063$, derived from the calcium doublet of the ESPRESSO combined spectra. This indicates that the star has likely a rotation period estimated to $P_{rot} = 30 \pm 3.1$ days (Astudillo-Defru et al. 2017a). According to this rotation period and the stellar radius, we expect the $v\sin(i)$ to be less than a km s$^{-1}$, which is consistent with the values of full width at half maximum (FWHM) of the cross-correlation function (CCF) from the ESPRESSO spectra. No sign of periodicity were found in the All Sky Automated Survey (ASAS; Pojmanski 2002) and the Zwicky Transient Facility (ZTF; Bellm et al. 2019) photometry as well as in the HARPS dataset (see Fig. 2). Unfortunately, this difficulty to find the rotation period of the star prevented us from determining the age of the system using gyrochronology.

# 3. Observations

## 3.1. TESS photometry

TOI-663 (TIC 54962195) was observed in two TESS sectors (9 and 35), one during the primary mission (February–March 2019) and one during the first extended mission (February–March 2021) with a two-minute cadence. The TESS Science Processing Operations Center (SPOC; Jenkins et al. 2016) at

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** TESS target pixel file image of TOI-663 in Sector 9 (created with tpfplotter; Aller et al. 2020). The electron counts are color-coded. The red bordered pixels are used in the SPOC SAP. The size of the red circles indicates the Gaia DR2 magnitudes of all nearby stars.</figcaption></figure>

NASA Ames processed the TESS photometric data and the resulting Presearch Data Conditioning Simple Aperture Photometry (PDCSAP; Smith et al. 2012; Stumpe et al. 2012, 2014) light curve of TOI-663 was corrected for dilution in the TESS aperture by known contaminating sources. It is critical to ensure that no visually close-by targets that could affect the depth of the transit are present in the 21'' TESS pixel and to check for contaminating eclipsing binaries. Figure 3 shows a plot of the target pixel file and the aperture mask that is used for the simple aperture photometry (SAP). We can see that one star overlaps the TESS aperture but its faintness results in minimal dilution of the TESS light curve. The SPOC conducted a transit search of the sector 9 light curve on 2019 April 25 and the combined light curve for sectors 9 and 35 on 2021 May 27 with an adaptive, noise-compensating matched filter (Jenkins 2002; Jenkins et al. 2010, 2020). The sector 9 search produced threshold crossing events with orbital periods of 2.60d and 4.70d; the combined search produced a third threshold crossing event with period of 7.10d. An initial limb-darkened transit model was fitted (Li et al. 2019) and a suite of diagnostic tests were conducted to help make or break the planetary nature of each signal (Twicken et al. 2018). The transit signatures passed all diagnostic tests presented in the SPOC Data Validation reports, and all TIC objects other than the target star were excluded as the transit source in each case. The TESS Science Office reviewed the vetting information and issued alerts for TOI-663.01 and TOI-663.02 on 2019 May 07, and TOI-663.03 on 2021 July 08 (Guerrero et al. 2021).

Canto Martins et al. (2020) conducted a search for rotation period in TESS TOIs using fast Fourier transform, Lomb-Scargle, and wavelet techniques, accompanied by a rigorous visual inspection, and found that TOI-663 does not show any sign of chromospheric activity throughout the observed time ranges.

## 3.2. Ground-based photometry

The TESS pixel scale is ~21" pixel$^{-1}$ and photometric apertures typically extend out to roughly 1 arc-minute, generally causing multiple stars to blend in the TESS aperture. To definitionally exclude the presence of another star causing the signal
<table><caption>**Table 2.** Summary of ground-based follow-up observations.</caption><thead><tr><th>Telescope/observatory</th><th>UTC date</th><th>Filter</th><th>Coverage</th><th>TOI</th><th>Use</th></tr></thead><tbody><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2020-01-09</td><td>0.85–1.55 µm</td><td>No ingress</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (3 telescopes)/La Silla</td><td>2020-03-22</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-02-15</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-02-28</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-03-13</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-05-04</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (3 telescopes)/La Silla</td><td>2022-03-04</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2022-03-30</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2022-04-20</td><td>0.85–1.55 µm</td><td>No ingress</td><td>.01</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2022-05-03</td><td>0.85–1.55 µm</td><td>No ingress</td><td>.01</td><td>TTV</td></tr><tr><td>GMU</td><td>2021-02-07</td><td>R</td><td>Partial</td><td>.01</td><td>None</td></tr><tr><td>LCO-McD</td><td>2021-02-08</td><td>zs</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>LCO-CTIO</td><td>2020-02-05</td><td>ip</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>LCO-CTIO</td><td>2020-02-18</td><td>ip</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>LCO-SAAO</td><td>2021-01-28</td><td>zs</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>LCO-SAAO</td><td>2021-02-23</td><td>ip</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>LCO-SSO</td><td>2020-04-15</td><td>zs</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>LCO-SSO</td><td>2020-04-28</td><td>zs</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>LCO-SSO</td><td>2020-05-11</td><td>zs</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>LCO-SSO</td><td>2021-01-18</td><td>zs</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>MuSCAT2/Teide Observatory</td><td>2020-01-17</td><td>g, r, i, zs</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>PEST</td><td>2020-04-28</td><td>Rc</td><td>Full</td><td>.01</td><td>None</td></tr><tr><td>SAINT-EX/San Pedro Mártir Observatory</td><td>2020-02-05</td><td>zs</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-02-15</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-03-06</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-03-20</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-04-22</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>ExTrA (1 telescopes)/La Silla</td><td>2021-06-08</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2022-04-09</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>LCO-CTIO</td><td>2020-03-01</td><td>ip</td><td>Full</td><td>.02</td><td>Global + TTV</td></tr><tr><td>LCO-CTIO</td><td>2020-03-15</td><td>gp</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>LCO-SAAO</td><td>2020-06-07</td><td>zs</td><td>Full</td><td>.02</td><td>Global + TTV</td></tr><tr><td>LCO-SAAO</td><td>2021-02-06</td><td>zs</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>MuSCAT/Okayama Observatory</td><td>2020-03-24</td><td>g, r, zs</td><td>Full</td><td>.02</td><td>Global + TTV</td></tr><tr><td>OAAlbanya</td><td>2020-03-19</td><td>Ic</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>PEST</td><td>2020-04-26</td><td>Rc</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>ExTrA (1 telescopes)/La Silla</td><td>2022-04-13</td><td>0.85–1.55 µm</td><td>Full</td><td>.03</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2022-04-20</td><td>0.85–1.55 µm</td><td>Full</td><td>.03</td><td>Global + TTV</td></tr></tbody></table>

in the TESS data, and to monitor for potential transit timing variations relative to a linear ephemeris, we conducted ground-based photometric follow-up observations of the field around TOI-663 as part of the TESS Follow-up Observing Program<fnref n="2" /> Sub Group 1 (TFOP; Collins 2019). All the follow-up photometry is summarized in Table 2.

### 3.2.1. ExTrA photometry

ExTrA (Bonfils et al. 2015) is a near-infrared (0.85 to 1.55 µm) multi-object spectrograph fed by three 60-cm telescopes located at La Silla Observatory in Chile. Seven full transits (and three partial transits) of TOI-663.01, six full transits of TOI-663.02, and two full transits of TOI-663.03 were observed using the

ExTrA telescopes (see Table 2). We used 8" aperture fibers and the low-resolution mode (R ~ 20) of the spectrograph with an exposure time of 60 seconds. Five fiber positioners are used at the focal plane of each telescope to select light from the target and four comparison stars chosen with 2MASS J-magnitudes (Skrutskie et al. 2006) and Gaia effective temperatures (Gaia Collaboration 2018) similar to the target. The resulting ExTrA data were analyzed using a custom data reduction software to produce a synthetic photometry in a 0.85–1.55 micron bandpass, described in more detail in Cointepas et al. (2021).

### 3.2.2. LCOGT

We observed nine and four full predicted transit windows of TOI-663.01 and TOI-663.02, respectively, using the

<footnote n="2">² https://tess.mit.edu/followup</footnote>
<table><caption>**Table 2.** Summary of ground-based follow-up observations.</caption><thead><tr><th>Telescope/observatory</th><th>UTC date</th><th>Filter</th><th>Coverage</th><th>TOI</th><th>Use</th></tr></thead><tbody><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2020-01-09</td><td>0.85–1.55 µm</td><td>No ingress</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (3 telescopes)/La Silla</td><td>2020-03-22</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-02-15</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-02-28</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-03-13</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-05-04</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (3 telescopes)/La Silla</td><td>2022-03-04</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2022-03-30</td><td>0.85–1.55 µm</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2022-04-20</td><td>0.85–1.55 µm</td><td>No ingress</td><td>.01</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2022-05-03</td><td>0.85–1.55 µm</td><td>No ingress</td><td>.01</td><td>TTV</td></tr><tr><td>GMU</td><td>2021-02-07</td><td>R</td><td>Partial</td><td>.01</td><td>None</td></tr><tr><td>LCO-McD</td><td>2021-02-08</td><td>zs</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>LCO-CTIO</td><td>2020-02-05</td><td>ip</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>LCO-CTIO</td><td>2020-02-18</td><td>ip</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>LCO-SAAO</td><td>2021-01-28</td><td>zs</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>LCO-SAAO</td><td>2021-02-23</td><td>ip</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>LCO-SSO</td><td>2020-04-15</td><td>zs</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>LCO-SSO</td><td>2020-04-28</td><td>zs</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>LCO-SSO</td><td>2020-05-11</td><td>zs</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>LCO-SSO</td><td>2021-01-18</td><td>zs</td><td>Full</td><td>.01</td><td>Global + TTV</td></tr><tr><td>MuSCAT2/Teide Observatory</td><td>2020-01-17</td><td>g, r, i, zs</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>PEST</td><td>2020-04-28</td><td>Rc</td><td>Full</td><td>.01</td><td>None</td></tr><tr><td>SAINT-EX/San Pedro Mártir Observatory</td><td>2020-02-05</td><td>zs</td><td>Full</td><td>.01</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-02-15</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-03-06</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-03-20</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2021-04-22</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>ExTrA (1 telescopes)/La Silla</td><td>2021-06-08</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2022-04-09</td><td>0.85–1.55 µm</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>LCO-CTIO</td><td>2020-03-01</td><td>ip</td><td>Full</td><td>.02</td><td>Global + TTV</td></tr><tr><td>LCO-CTIO</td><td>2020-03-15</td><td>gp</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>LCO-SAAO</td><td>2020-06-07</td><td>zs</td><td>Full</td><td>.02</td><td>Global + TTV</td></tr><tr><td>LCO-SAAO</td><td>2021-02-06</td><td>zs</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>MuSCAT/Okayama Observatory</td><td>2020-03-24</td><td>g, r, zs</td><td>Full</td><td>.02</td><td>Global + TTV</td></tr><tr><td>OAAlbanya</td><td>2020-03-19</td><td>Ic</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>PEST</td><td>2020-04-26</td><td>Rc</td><td>Full</td><td>.02</td><td>TTV</td></tr><tr><td>ExTrA (1 telescopes)/La Silla</td><td>2022-04-13</td><td>0.85–1.55 µm</td><td>Full</td><td>.03</td><td>Global + TTV</td></tr><tr><td>ExTrA (2 telescopes)/La Silla</td><td>2022-04-20</td><td>0.85–1.55 µm</td><td>Full</td><td>.03</td><td>Global + TTV</td></tr></tbody></table>

in the TESS data, and to monitor for potential transit timing variations relative to a linear ephemeris, we conducted ground-based photometric follow-up observations of the field around TOI-663 as part of the TESS Follow-up Observing Program<fnref n="2" /> Sub Group 1 (TFOP; Collins 2019). All the follow-up photometry is summarized in Table 2.

### 3.2.1. ExTrA photometry

ExTrA (Bonfils et al. 2015) is a near-infrared (0.85 to 1.55 µm) multi-object spectrograph fed by three 60-cm telescopes located at La Silla Observatory in Chile. Seven full transits (and three partial transits) of TOI-663.01, six full transits of TOI-663.02, and two full transits of TOI-663.03 were observed using the

ExTrA telescopes (see Table 2). We used 8" aperture fibers and the low-resolution mode (R ~ 20) of the spectrograph with an exposure time of 60 seconds. Five fiber positioners are used at the focal plane of each telescope to select light from the target and four comparison stars chosen with 2MASS J-magnitudes (Skrutskie et al. 2006) and Gaia effective temperatures (Gaia Collaboration 2018) similar to the target. The resulting ExTrA data were analyzed using a custom data reduction software to produce a synthetic photometry in a 0.85–1.55 micron bandpass, described in more detail in Cointepas et al. (2021).

### 3.2.2. LCOGT

We observed nine and four full predicted transit windows of TOI-663.01 and TOI-663.02, respectively, using the

<footnote n="2">² https://tess.mit.edu/followup</footnote>
Las Cumbres Observatory Global Telescope (LCOGT; Brown et al. 2013) 1.0 m network nodes. We used the TESS Transit Finder, which is a customized version of the Tapir software package (Jensen 2013), to schedule our transit observations. The 1 m telescopes are equipped with 4096 × 4096 SINISTRO cameras having an image scale of 0.389 per pixel, resulting in a 26' × 26' field of view. The images were calibrated by the standard LCOGT BANZAI pipeline (McCully et al. 2018). Differential photometric data were extracted using AstroImageJ (Collins et al. 2017). We detected transit-like signals in all 13 light curves using photometric apertures that exclude flux from all known neighbors of TOI-663, confirming the event on-target relative to all known Gaia DR3 and TICv8 stars.

### 3.2.3. MuSCAT

We observed one transit of TOI-663.02 on 2020 March 24 UT with the multiband imager MuSCAT (Narita et al. 2015) on the National Astronomical Observatory of Japan (NAOJ) 188 cm telescope at the Okayama Astro-Complex in Okayama, Japan. MuSCAT has three optical channels each equipped with a 1k × 1k CCD camera with a pixel scale of 0.361 pixel⁻¹, enabling simultaneous *g*-, *r*-, and *z*<sub>s</sub>-band imaging. The observation was performed with exposure times of 30, 15, and 30 s for the *g*, *r*, and *z*<sub>s</sub> bands, respectively. We used optical diffusers to obtain slightly defocused images. We performed image calibration and aperture photometry to extract relative light curves using a custom pipeline described in Fukui et al. (2011).

### 3.2.4. MuSCAT2

We observed one transit of TOI-663.01 on 2020 January 17 UT with the multi-band imager MuSCAT2 (Narita et al. 2019) on the 1.52m TCS telescope at the Teide Observatory in the Canary Islands, Spain. MuSCAT2 has four optical channels each equipped with a 1k × 1k CCD camera with a pixel scale of 0.44 pixel⁻¹, enabling simultaneous *g*-, *r*-, *i*-, and *z*<sub>s</sub>-band imaging. The observation was performed with exposure times of 60, 30, 30, and 15 s for the *g*, *r*, *i*, and *z*<sub>s</sub> bands, respectively. We performed image calibration and aperture photometry to extract relative light curves using a dedicated MuSCAT2 pipeline described in Parviainen et al. (2020).

### 3.2.5. GMU

TOI-663.01 was observed on the night of 2021 February 07 in the *R* band with 90 second exposures with the George Mason University (GMU) 0.8m Ritchey-Chretien telescope using the automated software control from Reefe et al. (2022). Intermittent clouds were present throughout the night, and a detection of the transit was not obtained. Data were reduced with a custom Python code and plate-solved with astrometry.net. Aperture photometry, light curve extraction, and detrending were performed with AstroImageJ (Collins et al. 2017). The results are available at ExoFOP at the NASA Exoplanet Archive.

### 3.2.6. OAAlbanya

TOI-663.02 was observed on the night of 2020 March 19 UT in the *Ic* band using an exposure of 120 s and with a good focus. The observation was performed from the 0.4 m Schmidt-Cassegrain telescope from Observatori Astronòmic Albanyà, at an altitude of 230 m in the northeast of Spain. The night was mostly transparent during the window of observation, with some

minor clouds at the out of transit phase of the observation. Data reduction was performed with aperture differential photometry using the AstroimageJ suite. Results showed a tentative signal coherent with an on-time transit, with high RMS related to the predicted depth, but consistent with an event neverthe-less. The dataset was detrended using the airmass values and the movement of the *Y*-axis of the image sequence to reduce systematics.

### 3.2.7. PEST

We observed TOI-663.01 and TOI-663.02 on 2020 April 28 UT and 2020 April 26 UT, respectively, with the Perth Exoplanet Survey Telescope (PEST) in the *Rc* band. PEST is located near Perth, Australia. At the time, the 0.3 m telescope was equipped with a 1530 × 1020 SBIG ST-8XME camera with an image scale of 1.2 pixel⁻¹ resulting in a 31' × 21' field of view. A custom pipeline based on C-Munipack<fnref n="3" /> was used to calibrate the images and extract the differential photometry.

### 3.2.8. SAINT-EX

Time series observations of the transit of TOI-663.01 were acquired with the SAINT-EX Observatory (Search And characterIsatioN of Transiting EXoplanets) on the night of 2020 February 05 UT. SAINT-EX is a 1-m F/8 Ritchey-Chretien telescope that resides within the Observatorio Astronómico Nacional located at the Sierra de San Pedro Mártir in Baja California, Mexico (31.04342 N, 115.45476 W) at an altitude of 2780 m (Demory et al. 2020). The observations were acquired with a deep-depleted Andor Ikon CCD in the *z'* filter and consisted of 396 data points, which cover a full transit of TOI-663.01. Each observation had an exposure time of 2.0 s spanning 5 h of continuous observations. The data were reduced and the aperture photometry was done, using the custom pipeline PRINCE done on site, as detailed in Demory et al. (2020). For the final light curve, the aperture of 14 pixels was chosen. We then analyzed the light curve with AstroImageJ (Collins et al. 2017) to correct for airmass and total counts in the aperture thus reducing systematics.

## 3.3. Radial velocities

### 3.3.1. HARPS

We obtained 40 spectra of TOI-663 with the High Accuracy Radial velocity Planet Searcher (HARPS; Mayor et al. 2003) echelle spectrograph at the ESO 3.6m telescope at La Silla Observatory in Chile. The HARPS optical spectrograph has a resolution of *R* = 115 000 and is stabilized in pressure and temperature, in order to achieve a sub-m s⁻¹ precision. The observations span 463 days and were taken between UT 2020 March 02 and UT 2021 June 07 (ProgID 105.20AK.002, 106.215E.004, 1102.C-0339). The second fiber was on sky during the observations.

These data were reduced with the Data Reduction Software<fnref n="4" />. Radial velocities were computed through maximum likelihood between individual spectra and the median stellar spectrum, where spectral zones containing tellurics were flagged. Details of the procedure are given in Astudillo-Defru et al. (2017b). The HARPS RVs of TOI-663 have a median uncertainty per data point of 4.5 m s⁻¹.

<footnote n="3">³ http://c-munipack.sourceforge.net</footnote>

<footnote n="4">⁴ http://www.eso.org/sci/facilities/lasilla/instruments/harps/doc/DRS.pdf</footnote>
Las Cumbres Observatory Global Telescope (LCOGT; Brown et al. 2013) 1.0 m network nodes. We used the TESS Transit Finder, which is a customized version of the Tapir software package (Jensen 2013), to schedule our transit observations. The 1 m telescopes are equipped with 4096 × 4096 SINISTRO cameras having an image scale of 0.389 per pixel, resulting in a 26' × 26' field of view. The images were calibrated by the standard LCOGT BANZAI pipeline (McCully et al. 2018). Differential photometric data were extracted using AstroImageJ (Collins et al. 2017). We detected transit-like signals in all 13 light curves using photometric apertures that exclude flux from all known neighbors of TOI-663, confirming the event on-target relative to all known Gaia DR3 and TICv8 stars.

### 3.2.3. MuSCAT

We observed one transit of TOI-663.02 on 2020 March 24 UT with the multiband imager MuSCAT (Narita et al. 2015) on the National Astronomical Observatory of Japan (NAOJ) 188 cm telescope at the Okayama Astro-Complex in Okayama, Japan. MuSCAT has three optical channels each equipped with a 1k × 1k CCD camera with a pixel scale of 0.361 pixel⁻¹, enabling simultaneous *g*-, *r*-, and *z*-*s*-band imaging. The observation was performed with exposure times of 30, 15, and 30 s for the *g*, *r*, and *z*-*s* bands, respectively. We used optical diffusers to obtain slightly defocused images. We performed image calibration and aperture photometry to extract relative light curves using a custom pipeline described in Fukui et al. (2011).

### 3.2.4. MuSCAT2

We observed one transit of TOI-663.01 on 2020 January 17 UT with the multi-band imager MuSCAT2 (Narita et al. 2019) on the 1.52m TCS telescope at the Teide Observatory in the Canary Islands, Spain. MuSCAT2 has four optical channels each equipped with a 1k × 1k CCD camera with a pixel scale of 0.44 pixel⁻¹, enabling simultaneous *g*-, *r*-, *i*-, and *z*-*s*-band imaging. The observation was performed with exposure times of 60, 30, 30, and 15 s for the *g*, *r*, *i*, and *z*-*s* bands, respectively. We performed image calibration and aperture photometry to extract relative light curves using a dedicated MuSCAT2 pipeline described in Parviainen et al. (2020).

### 3.2.5. GMU

TOI-663.01 was observed on the night of 2021 February 07 in the *R* band with 90 second exposures with the George Mason University (GMU) 0.8m Ritchey-Chretien telescope using the automated software control from Reefe et al. (2022). Intermittent clouds were present throughout the night, and a detection of the transit was not obtained. Data were reduced with a custom Python code and plate-solved with astrometry.net. Aperture photometry, light curve extraction, and detrending were performed with AstroImageJ (Collins et al. 2017). The results are available at ExoFOP at the NASA Exoplanet Archive.

### 3.2.6. OAAlbanya

TOI-663.02 was observed on the night of 2020 March 19 UT in the *Ic* band using an exposure of 120 s and with a good focus. The observation was performed from the 0.4 m Schmidt-Cassegrain telescope from Observatori Astronòmic Albanyà, at an altitude of 230 m in the northeast of Spain. The night was mostly transparent during the window of observation, with some

minor clouds at the out of transit phase of the observation. Data reduction was performed with aperture differential photometry using the AstroimageJ suite. Results showed a tentative signal coherent with an on-time transit, with high RMS related to the predicted depth, but consistent with an event neverthe-less. The dataset was detrended using the airmass values and the movement of the *Y*-axis of the image sequence to reduce systematics.

### 3.2.7. PEST

We observed TOI-663.01 and TOI-663.02 on 2020 April 28 UT and 2020 April 26 UT, respectively, with the Perth Exoplanet Survey Telescope (PEST) in the *Rc* band. PEST is located near Perth, Australia. At the time, the 0.3 m telescope was equipped with a 1530 × 1020 SBIG ST-8XME camera with an image scale of 1.2 pixel⁻¹ resulting in a 31' × 21' field of view. A custom pipeline based on C-Munipack<fnref n="3" /> was used to calibrate the images and extract the differential photometry.

### 3.2.8. SAINT-EX

Time series observations of the transit of TOI-663.01 were acquired with the SAINT-EX Observatory (Search And characterIsatioN of Transiting EXoplanets) on the night of 2020 February 05 UT. SAINT-EX is a 1-m F/8 Ritchey-Chretien telescope that resides within the Observatorio Astronómico Nacional located at the Sierra de San Pedro Mártir in Baja California, Mexico (31.04342 N, 115.45476 W) at an altitude of 2780 m (Demory et al. 2020). The observations were acquired with a deep-depleted Andor Ikon CCD in the *z'* filter and consisted of 396 data points, which cover a full transit of TOI-663.01. Each observation had an exposure time of 2.0 s spanning 5 h of continuous observations. The data were reduced and the aperture photometry was done, using the custom pipeline PRINCE done on site, as detailed in Demory et al. (2020). For the final light curve, the aperture of 14 pixels was chosen. We then analyzed the light curve with AstroImageJ (Collins et al. 2017) to correct for airmass and total counts in the aperture thus reducing systematics.

## 3.3. Radial velocities

### 3.3.1. HARPS

We obtained 40 spectra of TOI-663 with the High Accuracy Radial velocity Planet Searcher (HARPS; Mayor et al. 2003) echelle spectrograph at the ESO 3.6m telescope at La Silla Observatory in Chile. The HARPS optical spectrograph has a resolution of *R* = 115 000 and is stabilized in pressure and temperature, in order to achieve a sub-m s⁻¹ precision. The observations span 463 days and were taken between UT 2020 March 02 and UT 2021 June 07 (ProgID 105.20AK.002, 106.215E.004, 1102.C-0339). The second fiber was on sky during the observations.

These data were reduced with the Data Reduction Software<fnref n="4" />. Radial velocities were computed through maximum likelihood between individual spectra and the median stellar spectrum, where spectral zones containing tellurics were flagged. Details of the procedure are given in Astudillo-Defru et al. (2017b). The HARPS RVs of TOI-663 have a median uncertainty per data point of 4.5 m s⁻¹.

<footnote n="3">³ http://c-munipack.sourceforge.net</footnote>

<footnote n="4">⁴ http://www.eso.org/sci/facilities/lasilla/instruments/harps/doc/DRS.pdf</footnote>
### 3.3.2. ESPRESSO

We obtained ten spectra of TOI-663 with ESPRESSO (Pepe et al. 2021) at the 8.2 m ESO Very Large Telescope (VLT) array, at the Paranal Observatory in Chile. The observations were obtained from 2022 April 05 to 2022 April 18 for an ESPRESSO program designed to test if the density of sub-Neptunes orbiting M-dwarfs is different than sub-Neptunes orbiting more massive stars (Program ID: 109.2391.001; PI: Grieves). We observed TOI-663 in HR2x1 mode (1 UT, R ~ 140 000) over a spectral range from ~380 to ~780 nm. ESPRESSO is contained in a temperature and pressure controlled vacuum vessel to minimize the night spectral drift.

ESPRESSO observations can be carried out with simultaneous calibration using a Fabry-Pérot etalon on a second calibration fiber. However, given the relatively faint magnitude of TOI-663 (V = 13.67) we did not use the Fabry-Pérot etalon for simultaneous calibration, and the calibration fiber was pointed at the sky for background calibration that is used as a sky subtraction calibration in the data processing. The ESPRESSO RVs were processed using the Data & Analysis Center for Exoplanets (DACE), a facility dedicated to extra-solar planets data visualization, exchange, and analysis hosted by the University of Geneva<fnref n="5" />. The RV and stellar activity indicators provided in DACE are computed by cross-correlating the calibrated spectra with stellar templates for the specific spectral type, using the latest version (3.0.0) of the publicly available ESPRESSO Data Reduction Software<fnref n="6" />. The ten ESPRESSO RVs of TOI-663 have a median uncertainty of 1.58 m s⁻¹.

### 3.3.3. MAROON-X

We obtained 14 spectra of TOI-663 with MAROON-X (Seifahrt et al. 2018) at Gemini North in two campaigns from UT 2021 February 22 to 2021 March 04 and UT 2021 April 19 to 2021 April 30 under Prog-ID GN-2021A-Q-230. MAROON-X is a fairly new Extreme Precision Radial Velocity spectrograph operating in the red-optical wavelength regime, delivering R ≈ 85 000 from 500–920 nm in two arms, dubbed blue (500–670 nm) and red (650–920 nm). Spectra were taken with a fixed exposure time of 20 min and showed an average peak signal-to-noise ratio (S/N) of 100 in the red arm. The data were reduced by the instrument team using a custom python3 data reduction pipeline to produce optimum extracted and wavelength calibrated one-dimensional spectra. The RV analysis was performed using *serval* (Zechmeister et al. 2018), a template matching RV retrieval code in a custom python3 implementation. On average, the RV uncertainty per data point was 2.3 m s⁻¹ for the blue and 1.7 m s⁻¹ for the red arm. MAROON-X uses a stabilized Fabry-Perot etalon for wavelength and drift calibration (Seifahrt et al. 2022); it can deliver 30 cm s⁻¹ on-sky RV precision over short timescales (Trifonov et al. 2021) but suffers from inter-run RV offsets, which in the case of our observations had a prior of 1.5 ± 0.5 m s⁻¹ for the blue and 1.6 ± 0.5 m s⁻¹ for the red channel.

### 3.3.4. IRD

We observed TOI-663 with the IRD instrument (Tamura et al. 2012; Kotani et al. 2018) on the Subaru 8.2 m telescope. IRD simultaneously covers the Y, J, and H bands with the spectral resolving power of $\lambda/\Delta\lambda \approx 70 000$. A total of six near-infrared

spectra were secured between UT 2021 February 02 and June 08, each with an integration time of 1200–1500 s. One frame was interrupted by bad weather, which we discarded in the subsequent analysis. We reduced the raw frames with the standard reduction procedure, and measured the relative RVs using the custom RV-analysis pipeline for IRD (Hirano et al. 2020). Except for the frame affected by bad weather, the extracted one-dimensional spectra had the typical S/Ns of 51–73 per pixel at 1000 nm, and derived RVs had the statistical errors of 3.5–4.8 m s⁻¹.

### 3.4. High-resolution imaging

As part of our standard process for validating transiting exoplanets to assess the possible contamination of bound or unbound companions on the derived planetary radii (Ciardi et al. 2015), we observed the TOI 663 with near-infrared adaptive optics (AO) imaging at Keck and with optical speckle imaging at Gemini. Gaia DR3 is also used to provide additional constraints on the presence of undetected stellar companions as well as wide companions.

#### 3.4.1. Near-infrared adaptive optics imaging

The Keck Observatory observations were made with the NIRC2 instrument on Keck-II behind the natural guide star AO system (Wizinowich et al. 2000) on 2019 June 10 UT in the standard three-point dither pattern that is used with NIRC2 to avoid the left-lower quadrant of the detector, which is typically noisier than the other three quadrants. The dither pattern step size was 3" and was repeated twice, with each dither offset from the previous dither by 0.5". NIRC2 was used in the narrow-angle mode with a full field of view of ~10" and a pixel scale of approximately 0.0099442" per pixel. The Keck observations were made in the narrowband Br-γ filter ($\lambda_o = 2.1686$; $\Delta\lambda = 0.0326$ µm) and in the Jcont filter ($\lambda_o = 1.2132$; $\Delta\lambda = 0.0198$ µm), both with an integration time of 10 seconds for a total of 90 seconds on target per filter.

The science frames were flat-fielded and sky-subtracted. The flat fields were generated from a median average of dark subtracted flats taken on-sky. The flats were normalized such that the median value of the flats is unity. The sky frames were generated from the median average of the dithered science frames; each science image was then sky-subtracted and flat-fielded. The reduced science frames were combined into a single combined image using an intra-pixel interpolation that conserves flux, shifts the individual dithered frames by the appropriate fractional pixels, and median-coadds the frames. The final resolutions of the combined dithers were determined from the FWHM of the point spread functions: 0.052" at 2 µm and 0.064" at 1.2 µm.

The sensitivities of the final combined AO image were determined by injecting simulated sources azimuthally around the primary target every 20° at separations of integer multiples of the central source's FWHM (Furlan et al. 2017). The brightness of each injected source was scaled until standard aperture photometry detected it with $5\sigma$ significance. The resulting brightness of the injected sources relative to TOI-663 set the contrast limits at that injection location. The final $5\sigma$ limit at each separation was determined from the average of all of the determined limits at that separation and the uncertainty on the limit was set by the rms dispersion of the azimuthal slices at a given radial distance. The Keck data have a sensitivity close-in of $\Delta.2$ mag at 0.06" (Br-γ) and $\delta.4$ mag at 0.07" (Jcont), and deeper sensitivity at wider separations ($\delta.7$ mag at $\ge 0.5$"); the final sensitivity curves

<footnote n="5">⁵ https://dace.unige.ch/dashboard/index.html</footnote>

<footnote n="6">⁶ https://www.eso.org/sci/software/pipelines/espresso/espresso-pipe-recipes.html</footnote>
### 3.3.2. ESPRESSO

We obtained ten spectra of TOI-663 with ESPRESSO (Pepe et al. 2021) at the 8.2 m ESO Very Large Telescope (VLT) array, at the Paranal Observatory in Chile. The observations were obtained from 2022 April 05 to 2022 April 18 for an ESPRESSO program designed to test if the density of sub-Neptunes orbiting M-dwarfs is different than sub-Neptunes orbiting more massive stars (Program ID: 109.2391.001; PI: Grieves). We observed TOI-663 in HR2x1 mode (1 UT, R ~ 140 000) over a spectral range from ~380 to ~780 nm. ESPRESSO is contained in a temperature and pressure controlled vacuum vessel to minimize the night spectral drift.

ESPRESSO observations can be carried out with simultaneous calibration using a Fabry-Pérot etalon on a second calibration fiber. However, given the relatively faint magnitude of TOI-663 (V = 13.67) we did not use the Fabry-Pérot etalon for simultaneous calibration, and the calibration fiber was pointed at the sky for background calibration that is used as a sky subtraction calibration in the data processing. The ESPRESSO RVs were processed using the Data & Analysis Center for Exoplanets (DACE), a facility dedicated to extra-solar planets data visualization, exchange, and analysis hosted by the University of Geneva<fnref n="5" />. The RV and stellar activity indicators provided in DACE are computed by cross-correlating the calibrated spectra with stellar templates for the specific spectral type, using the latest version (3.0.0) of the publicly available ESPRESSO Data Reduction Software<fnref n="6" />. The ten ESPRESSO RVs of TOI-663 have a median uncertainty of 1.58 m s⁻¹.

### 3.3.3. MAROON-X

We obtained 14 spectra of TOI-663 with MAROON-X (Seifahrt et al. 2018) at Gemini North in two campaigns from UT 2021 February 22 to 2021 March 04 and UT 2021 April 19 to 2021 April 30 under Prog-ID GN-2021A-Q-230. MAROON-X is a fairly new Extreme Precision Radial Velocity spectrograph operating in the red-optical wavelength regime, delivering R ≈ 85 000 from 500–920 nm in two arms, dubbed blue (500–670 nm) and red (650–920 nm). Spectra were taken with a fixed exposure time of 20 min and showed an average peak signal-to-noise ratio (S/N) of 100 in the red arm. The data were reduced by the instrument team using a custom python3 data reduction pipeline to produce optimum extracted and wavelength calibrated one-dimensional spectra. The RV analysis was performed using *serval* (Zechmeister et al. 2018), a template matching RV retrieval code in a custom python3 implementation. On average, the RV uncertainty per data point was 2.3 m s⁻¹ for the blue and 1.7 m s⁻¹ for the red arm. MAROON-X uses a stabilized Fabry-Perot etalon for wavelength and drift calibration (Seifahrt et al. 2022); it can deliver 30 cm s⁻¹ on-sky RV precision over short timescales (Trifonov et al. 2021) but suffers from inter-run RV offsets, which in the case of our observations had a prior of 1.5 ± 0.5 m s⁻¹ for the blue and 1.6 ± 0.5 m s⁻¹ for the red channel.

### 3.3.4. IRD

We observed TOI-663 with the IRD instrument (Tamura et al. 2012; Kotani et al. 2018) on the Subaru 8.2 m telescope. IRD simultaneously covers the Y, J, and H bands with the spectral resolving power of $\lambda/\Delta\lambda \approx 70 000$. A total of six near-infrared

spectra were secured between UT 2021 February 02 and June 08, each with an integration time of 1200–1500 s. One frame was interrupted by bad weather, which we discarded in the subsequent analysis. We reduced the raw frames with the standard reduction procedure, and measured the relative RVs using the custom RV-analysis pipeline for IRD (Hirano et al. 2020). Except for the frame affected by bad weather, the extracted one-dimensional spectra had the typical S/Ns of 51–73 per pixel at 1000 nm, and derived RVs had the statistical errors of 3.5–4.8 m s⁻¹.

## 3.4. High-resolution imaging

As part of our standard process for validating transiting exoplanets to assess the possible contamination of bound or unbound companions on the derived planetary radii (Ciardi et al. 2015), we observed the TOI 663 with near-infrared adaptive optics (AO) imaging at Keck and with optical speckle imaging at Gemini. Gaia DR3 is also used to provide additional constraints on the presence of undetected stellar companions as well as wide companions.

### 3.4.1. Near-infrared adaptive optics imaging

The Keck Observatory observations were made with the NIRC2 instrument on Keck-II behind the natural guide star AO system (Wizinowich et al. 2000) on 2019 June 10 UT in the standard three-point dither pattern that is used with NIRC2 to avoid the left-lower quadrant of the detector, which is typically noisier than the other three quadrants. The dither pattern step size was 3" and was repeated twice, with each dither offset from the previous dither by 0.5". NIRC2 was used in the narrow-angle mode with a full field of view of ~10" and a pixel scale of approximately 0.0099442" per pixel. The Keck observations were made in the narrowband Br-γ filter ($\lambda_o = 2.1686$; $\Delta\lambda = 0.0326$ µm) and in the Jcont filter ($\lambda_o = 1.2132$; $\Delta\lambda = 0.0198$ µm), both with an integration time of 10 seconds for a total of 90 seconds on target per filter.

The science frames were flat-fielded and sky-subtracted. The flat fields were generated from a median average of dark subtracted flats taken on-sky. The flats were normalized such that the median value of the flats is unity. The sky frames were generated from the median average of the dithered science frames; each science image was then sky-subtracted and flat-fielded. The reduced science frames were combined into a single combined image using an intra-pixel interpolation that conserves flux, shifts the individual dithered frames by the appropriate fractional pixels, and median-coadds the frames. The final resolutions of the combined dithers were determined from the FWHM of the point spread functions: 0.052" at 2 µm and 0.064" at 1.2 µm.

The sensitivities of the final combined AO image were determined by injecting simulated sources azimuthally around the primary target every 20° at separations of integer multiples of the central source's FWHM (Furlan et al. 2017). The brightness of each injected source was scaled until standard aperture photometry detected it with $5\sigma$ significance. The resulting brightness of the injected sources relative to TOI-663 set the contrast limits at that injection location. The final $5\sigma$ limit at each separation was determined from the average of all of the determined limits at that separation and the uncertainty on the limit was set by the rms dispersion of the azimuthal slices at a given radial distance. The Keck data have a sensitivity close-in of $\Delta.2$ mag at 0.06" (Br-γ) and $\delta.4$ mag at 0.07" (Jcont), and deeper sensitivity at wider separations ($\delta.7$ mag at $\ge 0.5$"); the final sensitivity curves

<footnote n="5">⁵ https://dace.unige.ch/dashboard/index.html</footnote>

<footnote n="6">⁶ https://www.eso.org/sci/software/pipelines/espresso/espresso-pipe-recipes.html</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Keck near-infrared AO imaging and sensitivity curve for TOI-663 taken in the Br-γ filter. The image reaches a contrast of ~7.7 magnitudes fainter than the host star within 0.′5. Inset: image of the central portion of the data.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** 5σ speckle imaging contrast curves of Gemini South Zorro in both filters as a function of the angular separation out to 1.2 arcsec, the end of speckle coherence. The inset shows the reconstructed 832 nm image with a 1 arcsec scale bar. The star, TOI-663, was found to have no close companions to within the angular and contrast levels achieved.</figcaption></figure>

for the Keck Br-γ imaging is shown in (Fig. 4). No close-in (≤1′′) stellar companions were detected by Keck in agreement with the optical speckle imaging.

### 3.4.2. Optical speckle imaging

TOI-663 was observed on 2019 June 10 and 2020 January 12 UT using the 'Alopeke and Zorro speckle instruments on Gemini South and North, respectively. Zorro and 'Alopeke (Scott et al. 2021) both provide simultaneous speckle imaging in two bands (562 and 832 nm) with output data products including a reconstructed image and robust contrast limits on companion detections (Howell et al. 2011). Both datasets gave consistent results and Fig. 5 shows the Gemini South Zorro resultant 5-sigma contrast curves and the reconstructed 832-nm speckle image. We find that TOI-663 is a single star with no companion brighter than about 5–9 magnitudes, respectively, from 0.1 to 1.2 arcsec. At the distance of TOI-663 (64 pc), these angular limits correspond to spatial limits of 6.4–77 AU.

### 3.4.3. Gaia assessment

In addition to the high-resolution imaging, we used Gaia to identify any wide stellar companions that may be bound members

of the system. Typically, these stars are already in the TESS Input Catalog and their flux dilution to the transit has already been accounted for in the transit fits and associated derived parameters. There are no additional widely separated companions identified by Gaia that have the same distance and proper motion as TOI-663 (see also Mugrauer & Michel 2020, 2021).

Additionally, the Gaia DR3 astrometry provides additional information on the possibility of inner companions that may have gone undetected by either Gaia or the high-resolution imaging. The Gaia renormalized unit weight error (RUWE) is a metric, similar to a reduced chi-square, where values that are ≤1.4 indicate that the Gaia astrometric solution is consistent with the star being single, whereas RUWE values ≥1.4 may indicate an astrometric excess noise, possibly caused by the presence of an unseen companion (e.g., Ziegler et al. 2020). TOI-663 has a Gaia DR3 RUWE value of 1.6 – while this value is slightly higher than the nominal value of 1.4, the 1.4 RUWE boundary is more of a guideline as opposed to a hard cutoff. Additionally, TOI-663 does not appear in the Gaia DR3 binary table; as such, the Gaia astrometric fit is still consistent with a single star model. However, this RUWE value could be linked to the long-term RV trend seen in the HARPS dataset, which could be associated with a bound companion.

## 4. Analysis and modeling

We used the software package juliet (Espinoza et al. 2019) to model the photometric and RV data. The algorithm is built on many publicly available tools for the modeling of transits (batman; Kreidberg 2015), RVs (radvel; Fulton et al. 2018), and Gaussian processes (GPs; george; Ambikasaran et al. 2015; celerite, Foreman-Mackey et al. 2017). In order to compare different models, juliet efficiently computes the Bayesian evidence (ln Z) using dynasty (Speagle 2020), a python package to estimate Bayesian posteriors and evidence using nested sampling methods. Instead of starting with an initial parameter vector centered on a likelihood maximum discovered through optimization techniques, nested sampling algorithms sample directly from the given priors. Throughout our various analyses, we ensured that we had enough live points given the number of free parameters to avoid missing peaks in the parameter space. Here we decided to conduct two independent analyses of our data. We started the analysis using only TESS photometry, and the resulting planet parameters were used as priors in a subsequent RV analysis. As the result of these two analyses, we were able to constrain the priors of some parameters and use them for a joint analysis of all our data.

### 4.1. TESS photometry

First, using juliet, we modeled the TESS PDCSAP light curve where our planet candidates was initially detected. The transit model fits the stellar density $ρ_★$ along with the planetary and jitter parameters. Except for the stellar density, which we calculated in Sect. 2, we used ExoFOP priors for the orbital parameters. We adopted a few parametrization modifications when dealing with the transit photometry. We assigned a quadratic limb-darkening law for TESS, as shown to be appropriate for space-based missions (Espinoza & Jordán 2015), which was then parameterized with the uniform sampling scheme ($q_1, q_2$) introduced by Kipping (2013). Additionally, rather than fitting directly for the planet-to-star radius ratio ($p = R_p/R_★$) and the impact parameter of the orbit ($b = a/R_★ \cos i$), juliet used the parameterization introduced in Espinoza (2018) and fit for the parameters $r_1$ and
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Keck near-infrared AO imaging and sensitivity curve for TOI-663 taken in the Br-γ filter. The image reaches a contrast of ~7.7 magnitudes fainter than the host star within 0.′5. Inset: image of the central portion of the data.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** 5σ speckle imaging contrast curves of Gemini South Zorro in both filters as a function of the angular separation out to 1.2 arcsec, the end of speckle coherence. The inset shows the reconstructed 832 nm image with a 1 arcsec scale bar. The star, TOI-663, was found to have no close companions to within the angular and contrast levels achieved.</figcaption></figure>

for the Keck Br-γ imaging is shown in (Fig. 4). No close-in (≤1′′) stellar companions were detected by Keck in agreement with the optical speckle imaging.

### 3.4.2. Optical speckle imaging

TOI-663 was observed on 2019 June 10 and 2020 January 12 UT using the 'Alopeke and Zorro speckle instruments on Gemini South and North, respectively. Zorro and 'Alopeke (Scott et al. 2021) both provide simultaneous speckle imaging in two bands (562 and 832 nm) with output data products including a reconstructed image and robust contrast limits on companion detections (Howell et al. 2011). Both datasets gave consistent results and Fig. 5 shows the Gemini South Zorro resultant 5-sigma contrast curves and the reconstructed 832-nm speckle image. We find that TOI-663 is a single star with no companion brighter than about 5–9 magnitudes, respectively, from 0.1 to 1.2 arcsec. At the distance of TOI-663 (64 pc), these angular limits correspond to spatial limits of 6.4–77 AU.

### 3.4.3. Gaia assessment

In addition to the high-resolution imaging, we used Gaia to identify any wide stellar companions that may be bound members

of the system. Typically, these stars are already in the TESS Input Catalog and their flux dilution to the transit has already been accounted for in the transit fits and associated derived parameters. There are no additional widely separated companions identified by Gaia that have the same distance and proper motion as TOI-663 (see also Mugrauer & Michel 2020, 2021).

Additionally, the Gaia DR3 astrometry provides additional information on the possibility of inner companions that may have gone undetected by either Gaia or the high-resolution imaging. The Gaia renormalized unit weight error (RUWE) is a metric, similar to a reduced chi-square, where values that are ≤1.4 indicate that the Gaia astrometric solution is consistent with the star being single, whereas RUWE values ≥1.4 may indicate an astrometric excess noise, possibly caused by the presence of an unseen companion (e.g., Ziegler et al. 2020). TOI-663 has a Gaia DR3 RUWE value of 1.6 – while this value is slightly higher than the nominal value of 1.4, the 1.4 RUWE boundary is more of a guideline as opposed to a hard cutoff. Additionally, TOI-663 does not appear in the Gaia DR3 binary table; as such, the Gaia astrometric fit is still consistent with a single star model. However, this RUWE value could be linked to the long-term RV trend seen in the HARPS dataset, which could be associated with a bound companion.

## 4. Analysis and modeling

We used the software package juliet (Espinoza et al. 2019) to model the photometric and RV data. The algorithm is built on many publicly available tools for the modeling of transits (batman; Kreidberg 2015), RVs (radvel; Fulton et al. 2018), and Gaussian processes (GPs; george; Ambikasaran et al. 2015; celerite, Foreman-Mackey et al. 2017). In order to compare different models, juliet efficiently computes the Bayesian evidence (ln Z) using dynasty (Speagle 2020), a python package to estimate Bayesian posteriors and evidence using nested sampling methods. Instead of starting with an initial parameter vector centered on a likelihood maximum discovered through optimization techniques, nested sampling algorithms sample directly from the given priors. Throughout our various analyses, we ensured that we had enough live points given the number of free parameters to avoid missing peaks in the parameter space. Here we decided to conduct two independent analyses of our data. We started the analysis using only TESS photometry, and the resulting planet parameters were used as priors in a subsequent RV analysis. As the result of these two analyses, we were able to constrain the priors of some parameters and use them for a joint analysis of all our data.

### 4.1. TESS photometry

First, using juliet, we modeled the TESS PDCSAP light curve where our planet candidates was initially detected. The transit model fits the stellar density $ρ_★$ along with the planetary and jitter parameters. Except for the stellar density, which we calculated in Sect. 2, we used ExoFOP priors for the orbital parameters. We adopted a few parametrization modifications when dealing with the transit photometry. We assigned a quadratic limb-darkening law for TESS, as shown to be appropriate for space-based missions (Espinoza & Jordán 2015), which was then parameterized with the uniform sampling scheme ($q_1, q_2$) introduced by Kipping (2013). Additionally, rather than fitting directly for the planet-to-star radius ratio ($p = R_p/R_★$) and the impact parameter of the orbit ($b = a/R_★ \cos i$), juliet used the parameterization introduced in Espinoza (2018) and fit for the parameters $r_1$ and
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Modeling of the TESS data. Two left panels: TESS photometry time series from Sectors 9 and 35 along with the maximum a posteriori (solid black line). Three right panels: phase-folded light curve to the period of the planet (black points correspond to one-sixth of the transit duration for each planet). The three panels are associated with the b, c, and d planets, respectively. The dark area to the left on the far right plot is due to the transit of the other exoplanets that happened close to the transit of TOI-663 d.</figcaption></figure>

$r_2$ to guarantee full exploration of physically plausible values in the (p,b) plane. We assume circular orbits. The TESS PDC-SAP light curves are detrended from contamination, so we set the TESS dilution factor to one. In addition, we added in quadrature a jitter term $\sigma_{TESS}$ to the TESS photometric uncertainties, which may be underestimated due to additional systematics in the space-based photometry. Figure 6 depicts the best-fit transit model for the two TESS sectors. In the same figure, we also showed the phase-folded light curve for each planet.

## 4.2. Radial velocity analysis

For the RV modeling using *juliet*, we tested a three-planet and four-planet model including or not a linear long-term trend and GPs and we compared the evidence to select the one with the highest probability. We used as prior for each candidate the period and time-of-transit center estimated from the TESS photometry analysis and assume circular orbits. After fitting all of our datasets with a three-planet model including a trend seen in the HARPS dataset, we saw that the model did not fit properly our set of RV measurements for the different instruments. As a model with only three planets and a trend does not seem to match our data, we tested two different models: adding a GP or considering the possibility of having a fourth non-transiting exoplanet in the system. We looked for another periodic signal in the RV datasets. We found a potential peak at longer period around 25.8 days in the HARPS dataset compatible at the ~1.4-sigma level with the estimated rotational period derived from the log $R'_{HK}$ value (see Sect. 2). We fitted a four-planet model to our RV datasets (considering this new planet is not transiting and with a period range between 8 and 30 days) and found that the evidence of this model (ln Z = 293.33 ± 0.66) was slightly lower than the evidence for the model with 3 exoplanets, a long-term trend and a GP (ln Z = 294.60 ± 0.67), although the model comparison does not favor one model over the other (Kass & Raftery 1995). For the joint fit of photometry and RVs we adopted the model with three exoplanets, a long-term trend and a GP (Matern kernel with non-informative loguniform priors), as the signal of a possible fourth exoplanet was not convincing and the GP allows a better description of the residuals of the three-planet model. The long-term trend was estimated at $0.0229_{-0.0085}^{+0.0091}$ m day$^{-1}$ in the joint fit presented in the next section.

## 4.3. Joint fit of all data

We performed a joint analysis using *juliet* of the TESS photometry complemented by some ground-based photometry from ExTrA, LCOGT, and MuSCAT and velocity data from HARPS, ESPRESSO, MAROON-X, and IRD to obtain the most precise

parameters of the TOI-663 system. Each individual ground-based light curve presented in Sect. 3 was analyzed for the transit time variation (TTV) analysis (see Sect. 4.5). From this, we decided to only select a few of these transits for the global analysis to improve on the radius measurement. To choose which transits to use, we kept a first selection of transits for which the posterior of the time-of-transit center was more constrained after the fit (a nonuniform posterior on $t_0$), and then we only selected the nights with the smallest systematics visible in the data. We selected all transits available for TOI-663.03 as it was the least constrained by the TESS photometry.

After a verification that the depth of the transits were consistent in the different wavelength ranges, we decided to fit only one set of $r_1$ and $r_2$ parameters for each planet. To detrend the photometry data, we used time-dependent GPs. We chose the approximate Matern kernel proposed by Foreman-Mackey et al. (2017) because there are no obvious quasi-periodic oscillations in the light curves. We do not have a sufficient precision on the RV data to constrain the eccentricity of the different planets in this system so we decided to fix them at 0 to decrease the number of free parameters. Multi-planet systems tend to have low eccentricities but this decision is explained in more detail in the analysis in Sect. 4.4. We seek the posterior distribution of 127 free parameters in total, and due to the large number of parameters, we perform this fit using *dynesty*. We constrained our priors using previous *juliet* results before the final run for the joint modeling to optimize the computation time due to the large parameter space. Because we used rather broad priors for the planetary parameters in the separate fits of TESS and HARPS data, nested sampling is an efficient method for exploring the parameter space, and because most of the planetary parameters are specific to a given data type, they would not change significantly in a joint fit.

Table 3 shows the posteriors of the planetary parameters derived using the stellar parameters presented in Sect. 2, and Table C.1 shows the posteriors of the jitter parameters. Figures 7–9 show the results of the photometry of our joint fit to the data for the three exoplanets. In Figs. 10 and 11, we present the RV component of the joint fit and the phase-folded data. The corner plots of the RV semi-amplitude for the three planets are presented in Fig. B.1.

## 4.4. Stability of the system

From our global fit, we find that the three planets currently have period ratios $P_c/P_b=1.807$ and $P_d/P_c=1.513$. We investigated the short-term dynamics in this three-planet system. Our
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Modeling of the TESS data. Two left panels: TESS photometry time series from Sectors 9 and 35 along with the maximum a posteriori (solid black line). Three right panels: phase-folded light curve to the period of the planet (black points correspond to one-sixth of the transit duration for each planet). The three panels are associated with the b, c, and d planets, respectively. The dark area to the left on the far right plot is due to the transit of the other exoplanets that happened close to the transit of TOI-663 d.</figcaption></figure>

$r_2$ to guarantee full exploration of physically plausible values in the (p,b) plane. We assume circular orbits. The TESS PDC-SAP light curves are detrended from contamination, so we set the TESS dilution factor to one. In addition, we added in quadrature a jitter term $\sigma_{TESS}$ to the TESS photometric uncertainties, which may be underestimated due to additional systematics in the space-based photometry. Figure 6 depicts the best-fit transit model for the two TESS sectors. In the same figure, we also showed the phase-folded light curve for each planet.

## 4.2. Radial velocity analysis

For the RV modeling using *juliet*, we tested a three-planet and four-planet model including or not a linear long-term trend and GPs and we compared the evidence to select the one with the highest probability. We used as prior for each candidate the period and time-of-transit center estimated from the TESS photometry analysis and assume circular orbits. After fitting all of our datasets with a three-planet model including a trend seen in the HARPS dataset, we saw that the model did not fit properly our set of RV measurements for the different instruments. As a model with only three planets and a trend does not seem to match our data, we tested two different models: adding a GP or considering the possibility of having a fourth non-transiting exoplanet in the system. We looked for another periodic signal in the RV datasets. We found a potential peak at longer period around 25.8 days in the HARPS dataset compatible at the ~1.4-sigma level with the estimated rotational period derived from the log $R'_{HK}$ value (see Sect. 2). We fitted a four-planet model to our RV datasets (considering this new planet is not transiting and with a period range between 8 and 30 days) and found that the evidence of this model (ln Z = 293.33 ± 0.66) was slightly lower than the evidence for the model with 3 exoplanets, a long-term trend and a GP (ln Z = 294.60 ± 0.67), although the model comparison does not favor one model over the other (Kass & Raftery 1995). For the joint fit of photometry and RVs we adopted the model with three exoplanets, a long-term trend and a GP (Matern kernel with non-informative loguniform priors), as the signal of a possible fourth exoplanet was not convincing and the GP allows a better description of the residuals of the three-planet model. The long-term trend was estimated at $0.0229_{-0.0085}^{+0.0091}$ m day$^{-1}$ in the joint fit presented in the next section.

## 4.3. Joint fit of all data

We performed a joint analysis using *juliet* of the TESS photometry complemented by some ground-based photometry from ExTrA, LCOGT, and MuSCAT and velocity data from HARPS, ESPRESSO, MAROON-X, and IRD to obtain the most precise

parameters of the TOI-663 system. Each individual ground-based light curve presented in Sect. 3 was analyzed for the transit time variation (TTV) analysis (see Sect. 4.5). From this, we decided to only select a few of these transits for the global analysis to improve on the radius measurement. To choose which transits to use, we kept a first selection of transits for which the posterior of the time-of-transit center was more constrained after the fit (a nonuniform posterior on $t_0$), and then we only selected the nights with the smallest systematics visible in the data. We selected all transits available for TOI-663.03 as it was the least constrained by the TESS photometry.

After a verification that the depth of the transits were consistent in the different wavelength ranges, we decided to fit only one set of $r_1$ and $r_2$ parameters for each planet. To detrend the photometry data, we used time-dependent GPs. We chose the approximate Matern kernel proposed by Foreman-Mackey et al. (2017) because there are no obvious quasi-periodic oscillations in the light curves. We do not have a sufficient precision on the RV data to constrain the eccentricity of the different planets in this system so we decided to fix them at 0 to decrease the number of free parameters. Multi-planet systems tend to have low eccentricities but this decision is explained in more detail in the analysis in Sect. 4.4. We seek the posterior distribution of 127 free parameters in total, and due to the large number of parameters, we perform this fit using *dynesty*. We constrained our priors using previous *juliet* results before the final run for the joint modeling to optimize the computation time due to the large parameter space. Because we used rather broad priors for the planetary parameters in the separate fits of TESS and HARPS data, nested sampling is an efficient method for exploring the parameter space, and because most of the planetary parameters are specific to a given data type, they would not change significantly in a joint fit.

Table 3 shows the posteriors of the planetary parameters derived using the stellar parameters presented in Sect. 2, and Table C.1 shows the posteriors of the jitter parameters. Figures 7–9 show the results of the photometry of our joint fit to the data for the three exoplanets. In Figs. 10 and 11, we present the RV component of the joint fit and the phase-folded data. The corner plots of the RV semi-amplitude for the three planets are presented in Fig. B.1.

## 4.4. Stability of the system

From our global fit, we find that the three planets currently have period ratios $P_c/P_b=1.807$ and $P_d/P_c=1.513$. We investigated the short-term dynamics in this three-planet system. Our
<table><caption>**Table 3.** Posterior stellar and planetary parameters obtained from our joint photometric and RV juliet analysis for TOI-663.</caption><thead><tr><th>Parameter name</th><th>Posterior estimate for TOI-663 b</th><th>Posterior estimate for TOI-663 c</th><th>Posterior estimate for TOI-663 d</th></tr></thead><tbody><tr><td colspan="4"><strong>Posterior parameters</strong></td></tr><tr><td>Period (days)</td><td>2.5989020<sup>+0.0000032</sup><sub>-0.0000030</sub></td><td>4.6955507<sup>+0.0000095</sup><sub>-0.0000094</sub></td><td>7.1027268<sup>+0.0000297</sup><sub>-0.0000279</sub></td></tr><tr><td>t<sub>0</sub> (BJD<sub>TDB</sub>)</td><td>2 458 544.38918<sup>+0.00068</sup><sub>-0.00071</sub></td><td>2 458 543.40593<sup>+0.00105</sup><sub>-0.00107</sub></td><td>2 458 547.11249<sup>+0.00230</sup><sub>-0.00215</sub></td></tr><tr><td>r<sub>1</sub></td><td>0.458<sup>+0.083</sup><sub>-0.081</sub></td><td>0.483<sup>+0.070</sup><sub>-0.086</sub></td><td>0.445<sup>+0.101</sup><sub>-0.077</sub></td></tr><tr><td>r<sub>2</sub></td><td>0.0407<sup>+0.0012</sup><sub>-0.0012</sub></td><td>0.0404<sup>+0.0013</sup><sub>-0.0013</sub></td><td>0.0345<sup>+0.0021</sup><sub>-0.0021</sub></td></tr><tr><td>K (m s<sup>-1</sup>)</td><td>3.25<sup>+0.45</sup><sub>-0.47</sub></td><td>2.18<sup>+0.58</sup><sub>-0.57</sub></td><td>&lt;2.7 at 99%</td></tr><tr><td colspan="4"><strong>Derived transit and RVs parameters</strong></td></tr><tr><td>p = R<sub>p</sub>/R<sub>★</sub></td><td>0.0407<sup>+0.0012</sup><sub>-0.0012</sub></td><td>0.0404<sup>+0.0013</sup><sub>-0.0013</sub></td><td>0.0345<sup>+0.0021</sup><sub>-0.0021</sub></td></tr><tr><td>b = (a/R<sub>★</sub>) cos(i<sub>p</sub>)</td><td>0.19<sup>+0.12</sup><sub>-0.12</sub></td><td>0.23<sup>+0.11</sup><sub>-0.13</sub></td><td>0.17<sup>+0.15</sup><sub>-0.12</sub></td></tr><tr><td>i<sub>p</sub> (deg)</td><td>89.14<sup>+0.57</sup><sub>-0.60</sub></td><td>89.30<sup>+0.41</sup><sub>-0.35</sub></td><td>89.60<sup>+0.27</sup><sub>-0.36</sub></td></tr><tr><td colspan="4"><strong>Derived physical parameters<sup>(*)</sup></strong></td></tr><tr><td>M<sub>p</sub>(M<sub>⊕</sub>)</td><td>4.45<sup>+0.63</sup><sub>-0.66</sub></td><td>3.65<sup>+0.97</sup><sub>-0.97</sub></td><td>&lt;5.2 at 99%</td></tr><tr><td>R<sub>p</sub>(R<sub>⊕</sub>)</td><td>2.27<sup>+0.10</sup><sub>-0.09</sub></td><td>2.26<sup>+0.10</sup><sub>-0.10</sub></td><td>1.92<sup>+0.13</sup><sub>-0.13</sub></td></tr><tr><td>ρ<sub>p</sub>(g cm<sup>-3</sup>)</td><td>2.07<sup>+0.42</sup><sub>-0.38</sub></td><td>1.73<sup>+0.54</sup><sub>-0.49</sub></td><td>&lt;4.0 at 99%</td></tr><tr><td>a<sub>p</sub>(au)</td><td>0.0295 ± 0.0011</td><td>0.0437 ± 0.0016</td><td>0.0576<sup>+0.0021</sup><sub>-0.0022</sub></td></tr><tr><td>T<sub>eq</sub>(K)<sup>(†)</sup></td><td>674 ± 30</td><td>553 ± 24</td><td>482 ± 21</td></tr><tr><td>S(S<sub>⊕</sub>)</td><td>49.1<sup>+9.3</sup><sub>-8.1</sub></td><td>22.3<sup>+4.2</sup><sub>-3.7</sub></td><td>12.9<sup>+2.4</sup><sub>-2.1</sub></td></tr></tbody></table>

**Notes.** The posterior estimate corresponds to the median value. Error bars denote the 68% posterior credible intervals. (*)We sample from a normal distribution for the stellar mass, stellar radius, and stellar temperature that is based on the results from Sect. 2. (†)Equilibrium temperature was calculated assuming 0.3 Bond albedo and the semimajor axis distance considering an uniform temperature across the entire planet's surface.

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Phase-folded light curves for TOI-663 b for the TESS (17 transits), ExTrA (three transits with two telescopes each), and LCOGT (four transits) photometry. The black line is the maximum a posteriori model, the black points are data binned to one-sixth of the transit duration.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Phase-folded light curves for TOI-663 c for the TESS (eight transits), ExTrA (two transits with two telescopes each), MuSCAT (one transit with three filters), and LCOGT (two transits) photometry. The black line is the maximum a posteriori model, the black points are data binned to one-sixth of the transit duration.</figcaption></figure>
<table><caption>**Table 3.** Posterior stellar and planetary parameters obtained from our joint photometric and RV juliet analysis for TOI-663.</caption><thead><tr><th>Parameter name</th><th>Posterior estimate for TOI-663 b</th><th>Posterior estimate for TOI-663 c</th><th>Posterior estimate for TOI-663 d</th></tr></thead><tbody><tr><td colspan="4"><strong>Posterior parameters</strong></td></tr><tr><td>Period (days)</td><td>2.5989020<sup>+0.0000032</sup><sub>-0.0000030</sub></td><td>4.6955507<sup>+0.0000095</sup><sub>-0.0000094</sub></td><td>7.1027268<sup>+0.0000297</sup><sub>-0.0000279</sub></td></tr><tr><td>t<sub>0</sub> (BJD<sub>TDB</sub>)</td><td>2 458 544.38918<sup>+0.00068</sup><sub>-0.00071</sub></td><td>2 458 543.40593<sup>+0.00105</sup><sub>-0.00107</sub></td><td>2 458 547.11249<sup>+0.00230</sup><sub>-0.00215</sub></td></tr><tr><td>r<sub>1</sub></td><td>0.458<sup>+0.083</sup><sub>-0.081</sub></td><td>0.483<sup>+0.070</sup><sub>-0.086</sub></td><td>0.445<sup>+0.101</sup><sub>-0.077</sub></td></tr><tr><td>r<sub>2</sub></td><td>0.0407<sup>+0.0012</sup><sub>-0.0012</sub></td><td>0.0404<sup>+0.0013</sup><sub>-0.0013</sub></td><td>0.0345<sup>+0.0021</sup><sub>-0.0021</sub></td></tr><tr><td>K (m s<sup>-1</sup>)</td><td>3.25<sup>+0.45</sup><sub>-0.47</sub></td><td>2.18<sup>+0.58</sup><sub>-0.57</sub></td><td>&lt;2.7 at 99%</td></tr><tr><td colspan="4"><strong>Derived transit and RVs parameters</strong></td></tr><tr><td>p = R<sub>p</sub>/R<sub>★</sub></td><td>0.0407<sup>+0.0012</sup><sub>-0.0012</sub></td><td>0.0404<sup>+0.0013</sup><sub>-0.0013</sub></td><td>0.0345<sup>+0.0021</sup><sub>-0.0021</sub></td></tr><tr><td>b = (a/R<sub>★</sub>) cos(i<sub>p</sub>)</td><td>0.19<sup>+0.12</sup><sub>-0.12</sub></td><td>0.23<sup>+0.11</sup><sub>-0.13</sub></td><td>0.17<sup>+0.15</sup><sub>-0.12</sub></td></tr><tr><td>i<sub>p</sub> (deg)</td><td>89.14<sup>+0.57</sup><sub>-0.60</sub></td><td>89.30<sup>+0.41</sup><sub>-0.35</sub></td><td>89.60<sup>+0.27</sup><sub>-0.36</sub></td></tr><tr><td colspan="4"><strong>Derived physical parameters<sup>(*)</sup></strong></td></tr><tr><td>M<sub>p</sub>(M<sub>⊕</sub>)</td><td>4.45<sup>+0.63</sup><sub>-0.66</sub></td><td>3.65<sup>+0.97</sup><sub>-0.97</sub></td><td>&lt;5.2 at 99%</td></tr><tr><td>R<sub>p</sub>(R<sub>⊕</sub>)</td><td>2.27<sup>+0.10</sup><sub>-0.09</sub></td><td>2.26<sup>+0.10</sup><sub>-0.10</sub></td><td>1.92<sup>+0.13</sup><sub>-0.13</sub></td></tr><tr><td>ρ<sub>p</sub>(g cm<sup>-3</sup>)</td><td>2.07<sup>+0.42</sup><sub>-0.38</sub></td><td>1.73<sup>+0.54</sup><sub>-0.49</sub></td><td>&lt;4.0 at 99%</td></tr><tr><td>a<sub>p</sub>(au)</td><td>0.0295 ± 0.0011</td><td>0.0437 ± 0.0016</td><td>0.0576<sup>+0.0021</sup><sub>-0.0022</sub></td></tr><tr><td>T<sub>eq</sub>(K)<sup>(†)</sup></td><td>674 ± 30</td><td>553 ± 24</td><td>482 ± 21</td></tr><tr><td>S(S<sub>⊕</sub>)</td><td>49.1<sup>+9.3</sup><sub>-8.1</sub></td><td>22.3<sup>+4.2</sup><sub>-3.7</sub></td><td>12.9<sup>+2.4</sup><sub>-2.1</sub></td></tr></tbody></table>

**Notes.** The posterior estimate corresponds to the median value. Error bars denote the 68% posterior credible intervals. (*)We sample from a normal distribution for the stellar mass, stellar radius, and stellar temperature that is based on the results from Sect. 2. (†)Equilibrium temperature was calculated assuming 0.3 Bond albedo and the semimajor axis distance considering an uniform temperature across the entire planet's surface.

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Phase-folded light curves for TOI-663 b for the TESS (17 transits), ExTrA (three transits with two telescopes each), and LCOGT (four transits) photometry. The black line is the maximum a posteriori model, the black points are data binned to one-sixth of the transit duration.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Phase-folded light curves for TOI-663 c for the TESS (eight transits), ExTrA (two transits with two telescopes each), MuSCAT (one transit with three filters), and LCOGT (two transits) photometry. The black line is the maximum a posteriori model, the black points are data binned to one-sixth of the transit duration.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Phase-folded light curves for TOI-663 d for the TESS photometry (seven transits) and ExTrA photometry (two transits, one with one telescope and the other with two telescopes). The black line is the maximum a posteriori model, the black points are data binned to one-sixth of the transit duration. The dark area to the left on the far right plot is due to the transit of the other exoplanets that happened close to the transit of TOI-663 d.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 10.** RV measurements from HARPS, ESPRESSO, MAROON-X, and IRD for TOI-663. The data points with errors are shown in colors; the black line and gray band correspond respectively to the median and 1σ of 1000 randomly chosen posterior samples from the joint fit.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 11.** RVs phase-folded to the periods of the planets. The black line is the best-fit model from the joint fit.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Phase-folded light curves for TOI-663 d for the TESS photometry (seven transits) and ExTrA photometry (two transits, one with one telescope and the other with two telescopes). The black line is the maximum a posteriori model, the black points are data binned to one-sixth of the transit duration. The dark area to the left on the far right plot is due to the transit of the other exoplanets that happened close to the transit of TOI-663 d.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 10.** RV measurements from HARPS, ESPRESSO, MAROON-X, and IRD for TOI-663. The data points with errors are shown in colors; the black line and gray band correspond respectively to the median and 1σ of 1000 randomly chosen posterior samples from the joint fit.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 11.** RVs phase-folded to the periods of the planets. The black line is the best-fit model from the joint fit.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>Fig. 12. Chaos map around our best-fit solution. The map has a resolution of 121 × 121 pixels, which represent as many system configurations that were numerically integrated to derive a NAFF chaos estimation as defined in Eq. (1). We report the results in the form of a color code, where the bluer regions of the map refer to systems with smaller NAFF, and hence weaker chaos. The nearby principal MMR are indicated, and the vertical line marks the orbital period of planet c as estimated from our global fit.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>Fig. 13. TTVs. The observation minus the calculated linear ephemeris (O-C) diagram derived using all photometric datasets (with a transit) available through TFOP. Top: TOI-663.01. Middle: TOI-663.02. Bottom: TOI-663.03.</figcaption></figure>

objective is both to get a first hint on the range of orbital eccentricities compatible with the system stability, and to look for potential mean-motion resonances (MMRs). As such, we built a two-dimensional chaos map. It consists in a 121 × 121 grid of simulations in the parameter space exploring the orbital eccentricity of the middle planet $e_c$ on one axis and the orbital period of that same planet $P_c$ on the other axis (Fig. 12). All the other dynamical parameters were initially fixed to their estimation resulting from the global fit, hence with the hypothesis that $e_b = 0 = e_d$ and $\omega_i = 0$ for i = b,c,d. Each subsequent box in the grid defines a unique set of system orbital parameters that we used as initial conditions for numerical integrations of the equations of motion. We used the IAS15 adaptive timestep integrator to perform the numerical integrations (Rein & Spiegel 2015). IAS15 is part of the python package REBOUND<fnref n="7" /> (Rein & Liu 2012). We also included a correction for general relativistic effects following Anderson et al. (1975) and implemented in REBOUNDx (Tamayo et al. 2020). Each set of initial conditions was integrated over 10 kyr, during which the mean longitudes of the planets were regularly recorded in output. After the numerical integrations, we used the mean longitude time series to derive the numerical analysis of fundamental frequencies (NAFF) fast chaos indicator (Laskar 1990, 1993). This quantity informs about the strength of chaos in a planetary system via the variations of the orbital mean-motions, according to the following formula:

$$ \text{NAFF} = \max_j \left[ \log_{10} \frac{|n_{j,2} - n_{j,1}|}{n_{j,0}} \right]. \qquad (1) $$

In this equation, $n_1$ and $n_2$ depict the orbital mean-motion computed over the first and second half of the integration, respectively. $n_0$ is the initial mean motion. The mean-motion difference between the two integration halves is computed on each planet j, and the maximum value defines the NAFF of the planetary system. We refer to Stalport et al. (2022) for a more detailed synthesis of this process, and for a discussion of the link between NAFF and the orbital stability.

In Fig. 12, we report the results as a color code, where the bluer regions indicate smaller NAFF, and hence, weaker chaos. Orbital eccentricities $e_c$ larger than 0.1 are strongly disfavored, given the resulting strong chaos. While this exact threshold should be taken with caution since it results from a partial exploration of the parameter space, it provides us with insights into the range of eccentricity allowed by orbital stability. Only small values are not destructive for the system. This justifies the assumption of null eccentricities to favor the convergence of our global fit. Additionally, two vertical features are detected. They overlap with the 9:5 MMR between planets b and c, and the 3:2 MMR between planets c and d. A vertical line is drawn at the orbital period of planet c as derived from our global fit, which shows that the system does not fall in any of those features. We do not find the planets in any low-order MMR.

In the previous section, the possibility of a fourth planet in the TOI-663 system is mentioned. This additional candidate would have a significantly longer orbital period (more than three times the period of planet d), hence likely not affecting the orbital stability of the inner system.

## 4.5. TTV analysis

Using juliet, we calculated transit timing variations (TTVs) for all photometric datasets. Rather than fitting a period P and a time of transit center $t_0$, juliet searches for individual transit times. We combined multiple transits happening at the same time in order to reduce the uncertainties on the estimation of the time of transit center. The results of the analysis showing the difference between the observed transit times and the calculated linear ephemeris from TESS transits and for the different planets are presented in Fig. 13. The data show no significant variation at

<footnote n="7">⁷ https://rebound.readthedocs.io/en/latest/</footnote>
<figure><img src="image_13.png" /><figcaption>Fig. 12. Chaos map around our best-fit solution. The map has a resolution of 121 × 121 pixels, which represent as many system configurations that were numerically integrated to derive a NAFF chaos estimation as defined in Eq. (1). We report the results in the form of a color code, where the bluer regions of the map refer to systems with smaller NAFF, and hence weaker chaos. The nearby principal MMR are indicated, and the vertical line marks the orbital period of planet c as estimated from our global fit.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>Fig. 13. TTVs. The observation minus the calculated linear ephemeris (O-C) diagram derived using all photometric datasets (with a transit) available through TFOP. Top: TOI-663.01. Middle: TOI-663.02. Bottom: TOI-663.03.</figcaption></figure>

objective is both to get a first hint on the range of orbital eccentricities compatible with the system stability, and to look for potential mean-motion resonances (MMRs). As such, we built a two-dimensional chaos map. It consists in a 121 × 121 grid of simulations in the parameter space exploring the orbital eccentricity of the middle planet $e_c$ on one axis and the orbital period of that same planet $P_c$ on the other axis (Fig. 12). All the other dynamical parameters were initially fixed to their estimation resulting from the global fit, hence with the hypothesis that $e_b = 0 = e_d$ and $\omega_i = 0$ for i = b,c,d. Each subsequent box in the grid defines a unique set of system orbital parameters that we used as initial conditions for numerical integrations of the equations of motion. We used the IAS15 adaptive timestep integrator to perform the numerical integrations (Rein & Spiegel 2015). IAS15 is part of the python package REBOUND<fnref n="7" /> (Rein & Liu 2012). We also included a correction for general relativistic effects following Anderson et al. (1975) and implemented in REBOUNDx (Tamayo et al. 2020). Each set of initial conditions was integrated over 10 kyr, during which the mean longitudes of the planets were regularly recorded in output. After the numerical integrations, we used the mean longitude time series to derive the numerical analysis of fundamental frequencies (NAFF) fast chaos indicator (Laskar 1990, 1993). This quantity informs about the strength of chaos in a planetary system via the variations of the orbital mean-motions, according to the following formula:

$$ \text{NAFF} = \max_j \left[ \log_{10} \frac{|n_{j,2} - n_{j,1}|}{n_{j,0}} \right]. \qquad (1) $$

In this equation, $n_1$ and $n_2$ depict the orbital mean-motion computed over the first and second half of the integration, respectively. $n_0$ is the initial mean motion. The mean-motion difference between the two integration halves is computed on each planet j, and the maximum value defines the NAFF of the planetary system. We refer to Stalport et al. (2022) for a more detailed synthesis of this process, and for a discussion of the link between NAFF and the orbital stability.

In Fig. 12, we report the results as a color code, where the bluer regions indicate smaller NAFF, and hence, weaker chaos. Orbital eccentricities $e_c$ larger than 0.1 are strongly disfavored, given the resulting strong chaos. While this exact threshold should be taken with caution since it results from a partial exploration of the parameter space, it provides us with insights into the range of eccentricity allowed by orbital stability. Only small values are not destructive for the system. This justifies the assumption of null eccentricities to favor the convergence of our global fit. Additionally, two vertical features are detected. They overlap with the 9:5 MMR between planets b and c, and the 3:2 MMR between planets c and d. A vertical line is drawn at the orbital period of planet c as derived from our global fit, which shows that the system does not fall in any of those features. We do not find the planets in any low-order MMR.

In the previous section, the possibility of a fourth planet in the TOI-663 system is mentioned. This additional candidate would have a significantly longer orbital period (more than three times the period of planet d), hence likely not affecting the orbital stability of the inner system.

## 4.5. TTV analysis

Using juliet, we calculated transit timing variations (TTVs) for all photometric datasets. Rather than fitting a period P and a time of transit center $t_0$, juliet searches for individual transit times. We combined multiple transits happening at the same time in order to reduce the uncertainties on the estimation of the time of transit center. The results of the analysis showing the difference between the observed transit times and the calculated linear ephemeris from TESS transits and for the different planets are presented in Fig. 13. The data show no significant variation at

<footnote n="7">⁷ https://rebound.readthedocs.io/en/latest/</footnote>
<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Mass-radius diagram of the detected exoplanets (Otegi et al. 2020) and the TOI-663 planetary system. The green, gray, and dark blue lines show the mass-radius relation for an Earth-like, a pure-iron composition, and a 50% water content with and Earth-like interior, respectively. The blue region shows the M-R relation of pure water for the equilibrium temperature range of the planets in the TOI-663 system.</figcaption></figure>

the level of 4, 5, and 8 minutes RMS for TOI-663 b, TOI-663 c, and TOI-663 d, respectively, consistent with the uncertainties on the measured values.

**5. Discussion and conclusion**

In this paper we present the validation and characterization of three mini-Neptunes transiting the M dwarf TOI-663. The exoplanets were detected by the TESS mission, then confirmed and characterized from ground-based transit follow-up observations and precise RV measurements.

**5.1. Mass-radius diagram and internal structure**

TOI-663 b, c, and d, with respective radii of $2.27 \pm 0.10 R_⊕$, $2.26 \pm 0.10 R_⊕$, and $1.92 \pm 0.13 R_⊕$ and masses of $4.45 \pm 0.65 M_⊕$, and $3.65 \pm 0.97 M_⊕$, $1.9 \pm 1.3 M_⊕$, lie at the edge of the population of mini-Neptunes. However, the determination of the mass of planet d is at the ~1.7 sigma level, which is not sufficient to precisely characterize this exoplanet. We can only provide an upper limit for planet d of <5.2 $M_⊕$ at 99%. We plan to observe TOI-663 in the coming semesters using a combination of HARPS and the Near Infra Red Planet Searcher (NIRPS; Bouchy et al. 2017) to characterize the long-term trend found in the HARPS dataset, improve the mass measurements, and better characterize possible additional signals in this system.

The three detected planets of the TOI-663 system are low-mass mini-Neptunes. Their radii are considerably larger than rocky analogs (Fig. 14), and hence their radii must be dominated by volatiles, for example water.

To estimate the range of compositions for TOI-663 b, TOI-663 c, and TOI-663 d, we employed a layered interior model consisting of up to four layers: a H-He atmosphere, a water layer, a silicate mantle, and an iron core (see Dorn et al. 2017). We solved the standard structure equations for each layer and used a nested sampling algorithm to construct structure models that reproduce the measured masses and radii of the planets,

to quantify the degeneracy between various interior parameters, and to produce posterior probability distributions (Buchner et al. 2014) assuming a solar Mg/Si ratio (Lodders 2020) and an age of 5 Gyr. For the water layer, we used the AQUA equation of state (Haldemann et al. 2020). In the atmosphere layer, we adopted the semi-gray atmosphere model of Guillot (2010). This model takes the stellar irradiation and the intrinsic luminosity of the planet into account; the latter was estimated using the scaling relation from Mordasini (2020).

In a first scenario, we allowed the water mass fraction to be up to 50% underneath a H-He layer. The results are summarized in the upper half of Table 4. Within these constraints, we find that the three TOI-663 planets have a similar composition of roughly 30% core, 40% mantle, and 30% water layer mass fraction, with a small H-He atmosphere of $log(M_{atm}/M_p) = -4$ to $-5$.

In order to estimate the maximum water mass fraction in the planets, we also consider internal structures with no H-He atmospheres ($M_{atm} = 0$) and without limiting the total water mass fraction. Additionally, we followed Dorn & Lichtenberg (2021) in assuming an isothermal temperature profile until 0.1 bar, followed by an adiabatic profile. The results for this case are summarized in the lower half of Table 4. We find that the water mass fraction is about 20% higher than in the constrained case (i.e., approximately 50%). Accordingly, the core and mantle mass fractions are reduced by ~10% (i.e., $M_{core}/M_p \sim 20\%$ and $M_{mantle}/M_p \sim 30\%$).

Overall, the TOI-663 planets are expected to be water-rich and may have small H-He atmospheres. It should be noted, however, that while disfavored by the nested sampling algorithm, an Earth-like composition with a small H-He atmosphere (<1%) can also fit the observed mass and radius of the TOI-663 planetary system. Moreover, pollution of the H-He envelope could lead to a contraction of the atmosphere and hence a higher atmospheric mass fraction (Lozovsky et al. 2018). In addition, we have estimated here the amount of surface water and neglected any water dissolved deep in the interior, which could result in higher water mass fractions (Dorn & Lichtenberg 2021). Also, the host star elemental mass fractions can differ from solar
<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Mass-radius diagram of the detected exoplanets (Otegi et al. 2020) and the TOI-663 planetary system. The green, gray, and dark blue lines show the mass-radius relation for an Earth-like, a pure-iron composition, and a 50% water content with and Earth-like interior, respectively. The blue region shows the M-R relation of pure water for the equilibrium temperature range of the planets in the TOI-663 system.</figcaption></figure>

the level of 4, 5, and 8 minutes RMS for TOI-663 b, TOI-663 c, and TOI-663 d, respectively, consistent with the uncertainties on the measured values.

**5. Discussion and conclusion**

In this paper we present the validation and characterization of three mini-Neptunes transiting the M dwarf TOI-663. The exoplanets were detected by the TESS mission, then confirmed and characterized from ground-based transit follow-up observations and precise RV measurements.

**5.1. Mass-radius diagram and internal structure**

TOI-663 b, c, and d, with respective radii of $2.27 \pm 0.10 R_⊕$, $2.26 \pm 0.10 R_⊕$, and $1.92 \pm 0.13 R_⊕$ and masses of $4.45 \pm 0.65 M_⊕$, and $3.65 \pm 0.97 M_⊕$, $1.9 \pm 1.3 M_⊕$, lie at the edge of the population of mini-Neptunes. However, the determination of the mass of planet d is at the ~1.7 sigma level, which is not sufficient to precisely characterize this exoplanet. We can only provide an upper limit for planet d of <5.2 $M_⊕$ at 99%. We plan to observe TOI-663 in the coming semesters using a combination of HARPS and the Near Infra Red Planet Searcher (NIRPS; Bouchy et al. 2017) to characterize the long-term trend found in the HARPS dataset, improve the mass measurements, and better characterize possible additional signals in this system.

The three detected planets of the TOI-663 system are low-mass mini-Neptunes. Their radii are considerably larger than rocky analogs (Fig. 14), and hence their radii must be dominated by volatiles, for example water.

To estimate the range of compositions for TOI-663 b, TOI-663 c, and TOI-663 d, we employed a layered interior model consisting of up to four layers: a H-He atmosphere, a water layer, a silicate mantle, and an iron core (see Dorn et al. 2017). We solved the standard structure equations for each layer and used a nested sampling algorithm to construct structure models that reproduce the measured masses and radii of the planets,

to quantify the degeneracy between various interior parameters, and to produce posterior probability distributions (Buchner et al. 2014) assuming a solar Mg/Si ratio (Lodders 2020) and an age of 5 Gyr. For the water layer, we used the AQUA equation of state (Haldemann et al. 2020). In the atmosphere layer, we adopted the semi-gray atmosphere model of Guillot (2010). This model takes the stellar irradiation and the intrinsic luminosity of the planet into account; the latter was estimated using the scaling relation from Mordasini (2020).

In a first scenario, we allowed the water mass fraction to be up to 50% underneath a H-He layer. The results are summarized in the upper half of Table 4. Within these constraints, we find that the three TOI-663 planets have a similar composition of roughly 30% core, 40% mantle, and 30% water layer mass fraction, with a small H-He atmosphere of $log(M_{atm}/M_p) = -4$ to $-5$.

In order to estimate the maximum water mass fraction in the planets, we also consider internal structures with no H-He atmospheres ($M_{atm} = 0$) and without limiting the total water mass fraction. Additionally, we followed Dorn & Lichtenberg (2021) in assuming an isothermal temperature profile until 0.1 bar, followed by an adiabatic profile. The results for this case are summarized in the lower half of Table 4. We find that the water mass fraction is about 20% higher than in the constrained case (i.e., approximately 50%). Accordingly, the core and mantle mass fractions are reduced by ~10% (i.e., $M_{core}/M_p \sim 20\%$ and $M_{mantle}/M_p \sim 30\%$).

Overall, the TOI-663 planets are expected to be water-rich and may have small H-He atmospheres. It should be noted, however, that while disfavored by the nested sampling algorithm, an Earth-like composition with a small H-He atmosphere (<1%) can also fit the observed mass and radius of the TOI-663 planetary system. Moreover, pollution of the H-He envelope could lead to a contraction of the atmosphere and hence a higher atmospheric mass fraction (Lozovsky et al. 2018). In addition, we have estimated here the amount of surface water and neglected any water dissolved deep in the interior, which could result in higher water mass fractions (Dorn & Lichtenberg 2021). Also, the host star elemental mass fractions can differ from solar
<table><caption>**Table 4.** Nested sampling results for a maximum water mass fraction of 50% (top) and assuming no H-He atmosphere and no limit on the water mass fraction (bottom).</caption><thead><tr><th>Planet</th><th><i>M</i><sub>core</sub>/<i>M</i><sub>p</sub></th><th><i>M</i><sub>mantle</sub>/<i>M</i><sub>p</sub></th><th><i>M</i><sub>water</sub>/<i>M</i><sub>p</sub></th><th>log(<i>M</i><sub>atm</sub>/<i>M</i><sub>p</sub>)</th></tr></thead><tbody><tr><td>TOI-663 b</td><td>0.30 ± 0.19</td><td>0.42 ± 0.23</td><td>0.28 ± 0.13</td><td>-4.84 ± 1.81</td></tr><tr><td>TOI-663 c</td><td>0.30 ± 0.19</td><td>0.42 ± 0.23</td><td>0.28 ± 0.13</td><td>-4.09 ± 1.34</td></tr><tr><td>TOI-663 d</td><td>0.39 ± 0.23</td><td>0.39 ± 0.24</td><td>0.22 ± 0.12</td><td>-5.00 ± 1.43</td></tr><tr><td>TOI-663 b</td><td>0.16 ± 0.10</td><td>0.28 ± 0.15</td><td>0.56 ± 0.09</td><td>N/A</td></tr><tr><td>TOI-663 c</td><td>0.12 ± 0.08</td><td>0.26 ± 0.17</td><td>0.62 ± 0.13</td><td>N/A</td></tr><tr><td>TOI-663 d</td><td>0.18 ± 0.13</td><td>0.35 ± 0.22</td><td>0.47 ± 0.13</td><td>N/A</td></tr></tbody></table>

abundances, which can, for example, change the composition of the planet's mantle. Furthermore, the system's age is unknown, which could affect the characterization of the TOI-663 planetary system. Nonetheless, while these uncertainties could change the exact values of the nested sampling algorithm, it is clear that the composition of the TOI-663 planets is water-rich.

Interestingly, Luque & Pallé (2022) proposed a population of planets around M dwarfs that are water-rich worlds. On one hand, the planets of TOI-663 may indeed be part of water-rich worlds, according to our analysis. On the other hand, the planets do not follow the density trend of the proposed water-rich population as obtained from Zeng et al. (2019), which, however, neglects density variations due to differing equilibrium temperatures or ages. Furthermore, Rogers et al. (2023) show that the water-rich worlds can be equally well described by compositions that include H-He envelopes. In general, additional constraints are needed to further reduce the existing degeneracy between hydrogen- and water-dominated envelopes, for example, spectral data from JWST.

From a formation point of view, significant water budgets on the TOI-663 planets are feasible. High water budgets on the TOI-663 planets may be inherited from a formation beyond the snow line (Mordasini et al. 2009) as well as the production of endogenic water (Kimura & Ikoma 2020). How much endogenic water can be produced by redox reactions between a primordial H-He envelope in contact with an oxidized mantle for such small sub-Neptunes as the TOI-663 planets is still poorly known, but recent studies indicate that this is an important aspect to consider (Misener et al. 2023).

Despite evidence that M dwarfs tend to form rockier planets, the mass-radius diagram of known exoplanets with precise masses (<25%) and radii (<8%), shown in Fig. 14, shows that sub-Neptunes around M dwarfs appear to be less dense than those around other stars. The three exoplanets of the TOI-663 system more or less conform to this observation.

## 5.2. Position around the radius valley

Using data from the *Kepler* and K2 missions, Cloutier & Menou (2020) discovered that the radius valley separating rocky and gaseous planets has a negative slope with insolation around low-mass stars, as opposed to the positive slope found around solar-type stars (Martinez et al. 2019). This supports models of the direct formation of gas-poor terrestrial planets whose atmospheres were not removed through photo-evaporation (Owen & Wu 2017; Jin & Mordasini 2018; Lopez & Rice 2018), core-powered mass loss (Ginzburg et al. 2018; Gupta & Schlichting 2019, 2020), or impact erosion (Schlichting et al. 2015; Wyatt et al. 2020). In contrast to Cloutier & Menou (2020), Van Eylen et al. (2021) show the opposite trend of the slope of the radius valley with orbital period, which is similar to the

positive slope found around solar-type stars. To resolve the discrepancy, they state that a larger sample of well-studied planets orbiting M dwarf stars is required. Unfortunately, the three TOI-663 exoplanets are located above the radius valley and do not help us constrain either of the two models.

## 5.3. Similarity of multi-planet systems

M dwarfs tend to host multi-planet systems, which can provide much stronger tests for planetary formation theories than single-planet systems (Fang & Margot 2012; Steffen & Hwang 2015; Ballard & Johnson 2016). Weiss et al. (2018) have demonstrated a correlation between the size of planets in a multi-planet system and a regular spacing between them, described as the “peas in the pod” theory. This concept has been employed to restrict planet formation models. According to Otegi et al. (2022), who examined multi-planet systems observed by K2 and TESS for which both the radii and the masses were measured precisely, M-dwarf multi-planet systems exhibit more similarity in terms of density than in mass, indicating a possible connection in their formation mechanism. The densities of the three exoplanets of the TOI-663 system are consistent with each other within 1 sigma and therefore compatible with the Otegi et al. (2022) observation.

## 5.4. Potential for atmospheric characterization

Regarding atmospheric characterization, we calculated the transmission spectroscopic metric (TSM) in the *J* magnitude using the equation from Kempton et al. (2018) and obtained the following values: TSM = $76_{-12}^{+16}$ for TOI-663 b, TSM = $75_{-17}^{+29}$ for TOI-663 c, and TSM = $78_{-35}^{+120}$ for TOI-663 d. This value takes into account the equilibrium temperature of the exoplanet computed at its semimajor axis distance. Unfortunately, they all fall below the criteria of Kempton et al. (2018) for being the best candidates for transmission spectroscopy.

*Acknowledgements.* We are grateful to the ESO/La Silla staff for their continuous support. We acknowledge funding from the European Research Council under the ERC Grant Agreement n. 337591-ExTRA. This paper includes data collected by the TESS mission. Funding for the TESS mission is provided by the NASA Explorer Program. We acknowledge the use of public TESS Alert data from the 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. We thank the Swiss National Science Foundation (SNSF) and the Geneva University for their continuous support to our planet search programs. This work has been in particular carried out in the frame of the National Centre for Competence in Research 'PlanetS' supported by the Swiss National Science Foundation (SNSF). This work made use of tfpplotter by J. Lillo-Box (publicly available in www.github.com/jlillo/tpfplotter), which also made use of the python packages astropy, lightkurve, matplotlib and numpy. This research was carried out in part at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics
<table><caption>**Table 4.** Nested sampling results for a maximum water mass fraction of 50% (top) and assuming no H-He atmosphere and no limit on the water mass fraction (bottom).</caption><thead><tr><th>Planet</th><th><i>M</i><sub>core</sub>/<i>M</i><sub>p</sub></th><th><i>M</i><sub>mantle</sub>/<i>M</i><sub>p</sub></th><th><i>M</i><sub>water</sub>/<i>M</i><sub>p</sub></th><th>log(<i>M</i><sub>atm</sub>/<i>M</i><sub>p</sub>)</th></tr></thead><tbody><tr><td>TOI-663 b</td><td>0.30 ± 0.19</td><td>0.42 ± 0.23</td><td>0.28 ± 0.13</td><td>-4.84 ± 1.81</td></tr><tr><td>TOI-663 c</td><td>0.30 ± 0.19</td><td>0.42 ± 0.23</td><td>0.28 ± 0.13</td><td>-4.09 ± 1.34</td></tr><tr><td>TOI-663 d</td><td>0.39 ± 0.23</td><td>0.39 ± 0.24</td><td>0.22 ± 0.12</td><td>-5.00 ± 1.43</td></tr><tr><td>TOI-663 b</td><td>0.16 ± 0.10</td><td>0.28 ± 0.15</td><td>0.56 ± 0.09</td><td>N/A</td></tr><tr><td>TOI-663 c</td><td>0.12 ± 0.08</td><td>0.26 ± 0.17</td><td>0.62 ± 0.13</td><td>N/A</td></tr><tr><td>TOI-663 d</td><td>0.18 ± 0.13</td><td>0.35 ± 0.22</td><td>0.47 ± 0.13</td><td>N/A</td></tr></tbody></table>

abundances, which can, for example, change the composition of the planet's mantle. Furthermore, the system's age is unknown, which could affect the characterization of the TOI-663 planetary system. Nonetheless, while these uncertainties could change the exact values of the nested sampling algorithm, it is clear that the composition of the TOI-663 planets is water-rich.

Interestingly, Luque & Pallé (2022) proposed a population of planets around M dwarfs that are water-rich worlds. On one hand, the planets of TOI-663 may indeed be part of water-rich worlds, according to our analysis. On the other hand, the planets do not follow the density trend of the proposed water-rich population as obtained from Zeng et al. (2019), which, however, neglects density variations due to differing equilibrium temperatures or ages. Furthermore, Rogers et al. (2023) show that the water-rich worlds can be equally well described by compositions that include H-He envelopes. In general, additional constraints are needed to further reduce the existing degeneracy between hydrogen- and water-dominated envelopes, for example, spectral data from JWST.

From a formation point of view, significant water budgets on the TOI-663 planets are feasible. High water budgets on the TOI-663 planets may be inherited from a formation beyond the snow line (Mordasini et al. 2009) as well as the production of endogenic water (Kimura & Ikoma 2020). How much endogenic water can be produced by redox reactions between a primordial H-He envelope in contact with an oxidized mantle for such small sub-Neptunes as the TOI-663 planets is still poorly known, but recent studies indicate that this is an important aspect to consider (Misener et al. 2023).

Despite evidence that M dwarfs tend to form rockier planets, the mass–radius diagram of known exoplanets with precise masses (<25%) and radii (<8%), shown in Fig. 14, shows that sub-Neptunes around M dwarfs appear to be less dense than those around other stars. The three exoplanets of the TOI-663 system more or less conform to this observation.

## 5.2. Position around the radius valley

Using data from the *Kepler* and K2 missions, Cloutier & Menou (2020) discovered that the radius valley separating rocky and gaseous planets has a negative slope with insolation around low-mass stars, as opposed to the positive slope found around solar-type stars (Martinez et al. 2019). This supports models of the direct formation of gas-poor terrestrial planets whose atmospheres were not removed through photo-evaporation (Owen & Wu 2017; Jin & Mordasini 2018; Lopez & Rice 2018), core-powered mass loss (Ginzburg et al. 2018; Gupta & Schlichting 2019, 2020), or impact erosion (Schlichting et al. 2015; Wyatt et al. 2020). In contrast to Cloutier & Menou (2020), Van Eylen et al. (2021) show the opposite trend of the slope of the radius valley with orbital period, which is similar to the

positive slope found around solar-type stars. To resolve the discrepancy, they state that a larger sample of well-studied planets orbiting M dwarf stars is required. Unfortunately, the three TOI-663 exoplanets are located above the radius valley and do not help us constrain either of the two models.

## 5.3. Similarity of multi-planet systems

M dwarfs tend to host multi-planet systems, which can provide much stronger tests for planetary formation theories than single-planet systems (Fang & Margot 2012; Steffen & Hwang 2015; Ballard & Johnson 2016). Weiss et al. (2018) have demonstrated a correlation between the size of planets in a multi-planet system and a regular spacing between them, described as the “peas in the pod” theory. This concept has been employed to restrict planet formation models. According to Otegi et al. (2022), who examined multi-planet systems observed by K2 and TESS for which both the radii and the masses were measured precisely, M-dwarf multi-planet systems exhibit more similarity in terms of density than in mass, indicating a possible connection in their formation mechanism. The densities of the three exoplanets of the TOI-663 system are consistent with each other within 1 sigma and therefore compatible with the Otegi et al. (2022) observation.

## 5.4. Potential for atmospheric characterization

Regarding atmospheric characterization, we calculated the transmission spectroscopic metric (TSM) in the *J* magnitude using the equation from Kempton et al. (2018) and obtained the following values: TSM = $76_{-12}^{+16}$ for TOI-663 b, TSM = $75_{-17}^{+29}$ for TOI-663 c, and TSM = $78_{-35}^{+120}$ for TOI-663 d. This value takes into account the equilibrium temperature of the exoplanet computed at its semimajor axis distance. Unfortunately, they all fall below the criteria of Kempton et al. (2018) for being the best candidates for transmission spectroscopy.

*Acknowledgements.* We are grateful to the ESO/La Silla staff for their continuous support. We acknowledge funding from the European Research Council under the ERC Grant Agreement n. 337591-ExTRA. This paper includes data collected by the TESS mission. Funding for the TESS mission is provided by the NASA Explorer Program. We acknowledge the use of public TESS Alert data from the 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. We thank the Swiss National Science Foundation (SNSF) and the Geneva University for their continuous support to our planet search programs. This work has been in particular carried out in the frame of the National Centre for Competence in Research ‘PlanetS’ supported by the Swiss National Science Foundation (SNSF). This work made use of tfpplotter by J. Lillo-Box (publicly available in www.github.com/jlillo/tpfplotter), which also made use of the python packages astropy, lightkurve, matplotlib and numpy. This research was carried out in part at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics
and Space Administration (NASA). Some of the observations in this paper made use of the High-Resolution Imaging instrument Zorro and were obtained under Gemini LLP Proposal Number: GN/S-2021A-LP-105. Zorro was funded by the NASA Exoplanet Exploration Program and built at the NASA Ames Research Center by Steve B. Howell, Nic Scott, Elliott P. Horch, and Emmett Quigley. Zorro was mounted on the Gemini South telescope of the international Gemini Observatory, a program of NSF's OIR Lab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation on behalf of the Gemini partnership: the National Science Foundation (United States), National Research Council (Canada), Agencia Nacional de Investigación y Desarrollo (Chile), Ministerio de Ciencia, Tecnología e Innovación (Argentina), Ministério da Ciência, Tecnologia, Inovações e Comunicações (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea). Some observations in this paper were obtained with MAROON-X, a visitor instrument on the Gemini North telescope, located within the Maunakea Science Reserve and adjacent to the summit of Maunakea. We are grateful for the privilege of observing the Universe from a place that is unique in both its astronomical quality and its cultural significance. This work is based upon observations carried out at the Observatorio Astronómico Nacional on the Sierra de San Pedro Mártí (OAN-SPM), Baja California, México. SAINT-EX observations and team were supported by the Swiss National Science Foundation (PP00P2-163967 and PP00P2-190080), the Centre for Space and Habitability (CSH) of the University of Bern, the National Centre for Competence in Research PlanetS, supported by the Swiss National Science Foundation (SNSF), and UNAM PAPIIT-IG101321. This work makes use of observations from the LCOGT network. Part of the LCOGT telescope time was granted by NOIRLab through the Mid-Scale Innovations Program (MSIP). MSIP is funded by NSF. This research has made use of the Exoplanet Follow-up Observation Program (ExoFOP; DOI: 10.26134/ExoFOP5) website, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. Funding for the TESS mission is provided by NASA's Science Mission Directorate. This work is partly supported by JSPS KAKENHI Grant Numbers JP17H04574, JP18H05439, JP 18H05442, JP19K14783, JP20K14521, JP21H00035, JP21K13955, JP21K20376, Grant-in-Aid for JSPS Fellows, Grant Number JP20J21872, JST CREST Grant Number JPMJCR1761, and Astrobiology Center SATELLITE Research project AB022006. This article is based on observations made with the MuSCAT2 instrument, developed by ABC, at Teleoscopio Carlos Sánchez operated on the island of Tenerife by the IAC in the Spanish Observatorio del Teide. M.T. is supported by JSPS KAKENHI grant Nos.18H05442. C.D. acknowledges support from the Swiss National Science Foundation under grant TMSGI2_211313. This material is based upon work supported by NASA through the NASA Hubble Fellowship grant HST-HF2-51519.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. J.R.M. and B.L.C.M. acknowledge continuous grants from the brazilian funding agencies CAPES (Finance Code 001) and CNPq. The material is based upon work supported by NASA under award number 80GSFC21M0002.

## References

Allard, F., Homeier, D., & Freytag, B. 2012, Phil. Trans. R. Soc. London Ser. A, 370, 2765

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

Anderson, J. D., Esposito, P. B., Martin, W., Thornton, C. L., & Muhleman, D. O. 1975, ApJ, 200, 221

Astellidillo-Defru, N., Delfosse, X., Bonfils, X., et al. 2017a, A&A, 600, A13

Astellidillo-Defru, N., Forveille, T., Bonfils, X., et al. 2017b, A&A, 602, A88

Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Mantelet, G., & Andrae, R. 2018, AJ, 156, 58

Ballard, S., & Johnson, J. A. 2016, ApJ, 816, 66

Barros, S. C. C., Almenara, J. M., Demangeon, O., et al. 2015, MNRAS, 454, 4267

Batalha, N. E., Lewis, N. K., Line, M. R., Valenti, J., & Stevenson, K. 2018, ApJ, 856, L34

Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 018002

Bonfils, X., Almenara, J. M., Jocou, L., et al. 2015, SPIE Conf. Ser., 9605, 96051L

Bouchy, F., Doyon, R., Artigau, É., et al. 2017, The Messenger, 169, 21

Brown, T. M., Baliber, N., Bianco, F. B., et al. 2013, PASP, 125, 1031

Buchner, J., Georgakakis, A., Nandra, K., et al. 2014, A&A, 564, A125

Canto Martins, B. L., Gomes, R. L., Messias, Y. S., et al. 2020, ApJS, 250, 20

Chen, Y., Girardi, L., Bressan, A., et al. 2014, MNRAS, 444, 2525

Ciardi, D. R., Beichman, C. A., Horch, E. P., & Howell, S. B. 2015, ApJ, 805, 16

Cloutier, R., & Menou, K. 2020, AJ, 159, 211

Cointepas, M., Almenara, J. M., Bonfils, X., et al. 2021, A&A, 650, A145

Collins, K. 2019, AAS Meeting Abs., 233, 140.05

Collins, K. A., Kielkopf, J. F., Stassun, K. G., & Hessman, F. V. 2017, AJ, 153, 77

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246

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

Delisle, J. B. 2017, A&A, 605, A96

Demory, B. O., Pozuelos, F. J., Gómez Maqueo Chew, Y., et al. 2020, A&A, 642, A49

Díaz, R. F., Almenara, J. M., Santerne, A., et al. 2014, MNRAS, 441, 983

Dorn, C., & Lichtenberg, T. 2021, ApJ, 922, L4

Dorn, C., Venturini, J., Khan, A., et al. 2017, A&A, 597, A37

Dotter, A., Chaboyer, B., Jevremović, D., et al. 2008, ApJS, 178, 89

Espinoza, N. 2018, Res. Notes Am. Astron. Soc., 2, 209

Espinoza, N., & Jordán, A. 2015, MNRAS, 450, 1879

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

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

Fang, J., & Margot, J.-L. 2012, ApJ, 761, 92

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

Fukui, A., Narita, N., Tristram, P. J., et al. 2011, PASJ, 63, 287

Fulton, B. J., Petigura, E. A., Howard, A. W., et al. 2017, AJ, 154, 109

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

Furlan, E., Ciardi, D. R., Everett, M. E., et al. 2017, AJ, 153, 71

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

Gillon, M., Triaud, A. H. M. J., Demory, B.-O., et al. 2017, Nature, 542, 456

Ginzburg, S., Schlichting, H. E., & Sari, R. 2018, MNRAS, 476, 759

Gomes da Silva, J., Santos, N. C., & Bonfils, X. 2011, ASP Conf. Ser., 448, 1117

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

Guillot, T. 2010, A&A, 520, A27

Gupta, A., & Schlichting, H. E. 2019, MNRAS, 487, 24

Gupta, A., & Schlichting, H. E. 2020, MNRAS, 493, 792

Haldemann, J., Alibert, Y., Mordasini, C., & Benz, W. 2020, A&A, 643, A105

Hardegree-Ullman, K. K., Zink, J. K., Christiansen, J. L., et al. 2020, ApJS, 247, 28

Hirano, T., Kuzuhara, M., Kotani, T., et al. 2020, PASJ, 72, 93

Howell, S. B., Everett, M. E., Sherry, W., Horch, E., & Ciardi, D. R. 2011, AJ, 142, 19

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, 9913E

Jenkins, J. M., Tenenbaum, P., Seader, S., et al. 2020, Kepler Data Processing Handbook: Transiting Planet Search (Berlin: Springer)

Jensen, E. 2013, Astrophysics Source Code Library [record ascl:1306.007]

Jin, S., & Mordasini, C. 2018, ApJ, 853, 163

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

Kimura, T., & Ikoma, M. 2020, MNRAS, 496, 3755

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

Kotani, T., Tamura, M., Nishikawa, J., et al. 2018, SPIE Conf. Ser., 10702, 1070211

Kreidberg, L. 2015, PASP, 127, 1161

Kürster, M., Endl, M., Rouesnel, F., et al. 2003, ESA Spe. Pub., 539, 485

Laskar, J. 1990, Icarus, 88, 266

Laskar, J. 1993, Physica D, 67, 257

Lee, E. J., & Connors, N. J. 2021, ApJ, 908, 32

Li, J., Tenenbaum, P., Twicken, J. D., et al. 2019, PASP, 131, 024506

Lillo-Box, J., Figueira, P., Leleu, A., et al. 2020, A&A, 642, A121

Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2

Lodders, K. 2020, Solar Elemental Abundances (Oxford: Oxford University Press)

Lopez, E. D., & Rice, K. 2018, MNRAS, 479, 5303

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

Luger, R., Sestovic, M., Kruse, E., et al. 2017, Nat. Astron., 1, 0129

Luque, R., & Pallé, E. 2022, Science, 377, 1211

Maŕz Apellániz, J., & Weiler, M. 2018, A&A, 619, A180

Mann, A. W., Feiden, G. A., Gaidos, E., Boyajian, T., & von Braun, K. 2015, ApJ, 804, 64

Mann, A. W., Dupuy, T., Kraus, A. L., et al. 2019, ApJ, 871, 63

Martinez, C. F., Cunha, K., Ghezzi, L., & Smith, V. V. 2019, ApJ, 875, 29

Mayo, A. W., Vanderburg, A., Latham, D. W., et al. 2018, AJ, 155, 136
and Space Administration (NASA). Some of the observations in this paper made use of the High-Resolution Imaging instrument Zorro and were obtained under Gemini LLP Proposal Number: GN/S-2021A-LP-105. Zorro was funded by the NASA Exoplanet Exploration Program and built at the NASA Ames Research Center by Steve B. Howell, Nic Scott, Elliott P. Horch, and Emmett Quigley. Zorro was mounted on the Gemini South telescope of the international Gemini Observatory, a program of NSF's OIR Lab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation on behalf of the Gemini partnership: the National Science Foundation (United States), National Research Council (Canada), Agencia Nacional de Investigación y Desarrollo (Chile), Ministerio de Ciencia, Tecnología e Innovación (Argentina), Ministério da Ciência, Tecnologia, Inovações e Comunicações (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea). Some observations in this paper were obtained with MAROON-X, a visitor instrument on the Gemini North telescope, located within the Maunakea Science Reserve and adjacent to the summit of Maunakea. We are grateful for the privilege of observing the Universe from a place that is unique in both its astronomical quality and its cultural significance. This work is based upon observations carried out at the Observatorio Astronómico Nacional on the Sierra de San Pedro Mártí (OAN-SPM), Baja California, México. SAINT-EX observations and team were supported by the Swiss National Science Foundation (PP00P2-163967 and PP00P2-190080), the Centre for Space and Habitability (CSH) of the University of Bern, the National Centre for Competence in Research PlanetS, supported by the Swiss National Science Foundation (SNSF), and UNAM PAPIIT-IG101321. This work makes use of observations from the LCOGT network. Part of the LCOGT telescope time was granted by NOIRLab through the Mid-Scale Innovations Program (MSIP). MSIP is funded by NSF. This research has made use of the Exoplanet Follow-up Observation Program (ExoFOP; DOI: 10.26134/ExoFOP5) website, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. Funding for the TESS mission is provided by NASA's Science Mission Directorate. This work is partly supported by JSPS KAKENHI Grant Numbers JP17H04574, JP18H05439, JP 18H05442, JP19K14783, JP20K14521, JP21H00035, JP21K13955, JP21K20376, Grant-in-Aid for JSPS Fellows, Grant Number JP20J21872, JST CREST Grant Number JPMJCR1761, and Astrobiology Center SATELLITE Research project AB022006. This article is based on observations made with the MuSCAT2 instrument, developed by ABC, at Tele scopio Carlos Sánchez operated on the island of Tenerife by the IAC in the Spanish Observatorio del Teide. M.T. is supported by JSPS KAKENHI grant Nos.18H05442. C.D. acknowledges support from the Swiss National Science Foundation under grant TMSG12_211313. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE 1746045. G.S. acknowledges support provided by NASA through the NASA Hubble Fellowship grant HST-HF2-51519.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. J.R.M. and B.L.C.M. acknowledge continuous grants from the brazilian funding agencies CAPES (Finance Code 001) and CNPq. The material is based upon work supported by NASA under award number 80GSFC21M0002.

## References

Allard, F., Homeier, D., & Freytag, B. 2012, Phil. Trans. R. Soc. London Ser. A, 370, 2765

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

Anderson, J. D., Esposito, P. B., Martin, W., Thornton, C. L., & Muhleman, D. O. 1975, ApJ, 200, 221

Astellidillo-Defru, N., Delfosse, X., Bonfils, X., et al. 2017a, A&A, 600, A13

Astellidillo-Defru, N., Forveille, T., Bonfils, X., et al. 2017b, A&A, 602, A88

Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Mantelet, G., & Andrae, R. 2018, AJ, 156, 58

Ballard, S., & Johnson, J. A. 2016, ApJ, 816, 66

Barros, S. C. C., Almenara, J. M., Demangeon, O., et al. 2015, MNRAS, 454, 4267

Batalha, N. E., Lewis, N. K., Line, M. R., Valenti, J., & Stevenson, K. 2018, ApJ, 856, L34

Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 018002

Bonfils, X., Almenara, J. M., Jocou, L., et al. 2015, SPIE Conf. Ser., 9605, 96051L

Bouchy, F., Doyon, R., Artigau, É., et al. 2017, The Messenger, 169, 21

Brown, T. M., Baliber, N., Bianco, F. B., et al. 2013, PASP, 125, 1031

Buchner, J., Georgakakis, A., Nandra, K., et al. 2014, A&A, 564, A125

Canto Martins, B. L., Gomes, R. L., Messias, Y. S., et al. 2020, ApJS, 250, 20

Chen, Y., Girardi, L., Bressan, A., et al. 2014, MNRAS, 444, 2525

Ciardi, D. R., Beichman, C. A., Horch, E. P., & Howell, S. B. 2015, ApJ, 805, 16

Cloutier, R., & Menou, K. 2020, AJ, 159, 211

Cointepas, M., Almenara, J. M., Bonfils, X., et al. 2021, A&A, 650, A145

Collins, K. 2019, AAS Meeting Abs., 233, 140.05

Collins, K. A., Kielkopf, J. F., Stassun, K. G., & Hessman, F. V. 2017, AJ, 153, 77

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246

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

Delisle, J. B. 2017, A&A, 605, A96

Demory, B. O., Pozuelos, F. J., Gómez Maqueo Chew, Y., et al. 2020, A&A, 642, A49

Díaz, R. F., Almenara, J. M., Santerne, A., et al. 2014, MNRAS, 441, 983

Dorn, C., & Lichtenberg, T. 2021, ApJ, 922, L4

Dorn, C., Venturini, J., Khan, A., et al. 2017, A&A, 597, A37

Dotter, A., Chaboyer, B., Jevremović, D., et al. 2008, ApJS, 178, 89

Espinoza, N. 2018, Res. Notes Am. Astron. Soc., 2, 209

Espinoza, N., & Jordán, A. 2015, MNRAS, 450, 1879

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

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

Fang, J., & Margot, J.-L. 2012, ApJ, 761, 92

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

Fukui, A., Narita, N., Tristram, P. J., et al. 2011, PASJ, 63, 287

Fulton, B. J., Petigura, E. A., Howard, A. W., et al. 2017, AJ, 154, 109

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

Furlan, E., Ciardi, D. R., Everett, M. E., et al. 2017, AJ, 153, 71

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

Gillon, M., Triaud, A. H. M. J., Demory, B.-O., et al. 2017, Nature, 542, 456

Ginzburg, S., Schlichting, H. E., & Sari, R. 2018, MNRAS, 476, 759

Gomes da Silva, J., Santos, N. C., & Bonfils, X. 2011, ASP Conf. Ser., 448, 1117

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

Guillot, T. 2010, A&A, 520, A27

Gupta, A., & Schlichting, H. E. 2019, MNRAS, 487, 24

Gupta, A., & Schlichting, H. E. 2020, MNRAS, 493, 792

Haldemann, J., Alibert, Y., Mordasini, C., & Benz, W. 2020, A&A, 643, A105

Hardegree-Ullman, K. K., Zink, J. K., Christiansen, J. L., et al. 2020, ApJS, 247, 28

Hirano, T., Kuzuhara, M., Kotani, T., et al. 2020, PASJ, 72, 93

Howell, S. B., Everett, M. E., Sherry, W., Horch, E., & Ciardi, D. R. 2011, AJ, 142, 19

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, 9913E

Jenkins, J. M., Tenenbaum, P., Seader, S., et al. 2020, Kepler Data Processing Handbook: Transiting Planet Search (Berlin: Springer)

Jensen, E. 2013, Astrophysics Source Code Library [record ascl:1306.007]

Jin, S., & Mordasini, C. 2018, ApJ, 853, 163

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

Kimura, T., & Ikoma, M. 2020, MNRAS, 496, 3755

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

Kotani, T., Tamura, M., Nishikawa, J., et al. 2018, SPIE Conf. Ser., 10702, 1070211

Kreidberg, L. 2015, PASP, 127, 1161

Kürster, M., Endl, M., Rouesnel, F., et al. 2003, ESA Spe. Pub., 539, 485

Laskar, J. 1990, Icarus, 88, 266

Laskar, J. 1993, Physica D, 67, 257

Lee, E. J., & Connors, N. J. 2021, ApJ, 908, 32

Li, J., Tenenbaum, P., Twicken, J. D., et al. 2019, PASP, 131, 024506

Lillo-Box, J., Figueira, P., Leleu, A., et al. 2020, A&A, 642, A121

Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2

Lodders, K. 2020, Solar Elemental Abundances (Oxford: Oxford University Press)

Lopez, E. D., & Rice, K. 2018, MNRAS, 479, 5303

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

Luger, R., Sestovic, M., Kruse, E., et al. 2017, Nat. Astron., 1, 0129

Luque, R., & Pallé, E. 2022, Science, 377, 1211

Maŕz Apellániz, J., & Weiler, M. 2018, A&A, 619, A180

Mann, A. W., Feiden, G. A., Gaidos, E., Boyajian, T., & von Braun, K. 2015, ApJ, 804, 64

Mann, A. W., Dupuy, T., Kraus, A. L., et al. 2019, ApJ, 871, 63

Martinez, C. F., Cunha, K., Ghezzi, L., & Smith, V. V. 2019, ApJ, 875, 29

Mayo, A. W., Vanderburg, A., Latham, D. W., et al. 2018, AJ, 155, 136
Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20

McCully, C., Volgenau, N. H., Harbeck, D.-R., et al. 2018, SPIE Conf. Ser., 10707, 107070K

Millholland, S., Wang, S., & Laughlin, G. 2017, ApJ, 849, L33

Misener, W., Schlichting, H. E., & Young, E. D. 2023, MNRAS, 524, 981

Mordasini, C. 2020, A&A, 638, A52

Mordasini, C., Alibert, Y., & Benz, W. 2009, A&A, 501, 1139

Mugrauer, M., & Michel, K.-U. 2020, Astron. Nachr., 341, 996

Mugrauer, M., & Michel, K.-U. 2021, Astron. Nachr., 342, 840

Narita, N., Fukui, A., Kusakabe, N., et al. 2015, J. Astron. Teles. Instrum. Syst., 1, 045001

Narita, N., Fukui, A., Kusakabe, N., et al. 2019, J. Astron. Teles. Instrum. Syst., 5, 015001

Otegi, J. F., Bouchy, F., & Helled, R. 2020, A&A, 634, A43

Otegi, J. F., Helled, R., & Bouchy, F. 2022, A&A, 658, A107

Owen, J. E., & Wu, Y. 2017, ApJ, 847, 29

Parviainen, H., Palle, E., Zapatero-Osorio, M. R., et al. 2020, A&A, 633, A28

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

Pepe, F., Cristiani, S., Rebolo, R., et al. 2021, A&A, 645, A96

Pojmanski, G. 2002, Acta Astron., 52, 397

Reefe, M., Alfaro, O., Foster, S., et al. 2022, J. Astron. Teles. Instrum. Syst., 8, 027002

Rein, H., & Liu, S. F. 2012, A&A, 537, A128

Rein, H., & Spiegel, D. S. 2015, MNRAS, 446, 1424

Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telesc. Instrum. Syst., 1, 014003

Rogers, J. G., Schlichting, H. E., & Owen, J. E. 2023, ApJ, 947, L19

Schlichting, H. E., Sari, R., & Yalinewich, A. 2015, Icarus, 247, 81

Scott, N. J., Howell, S. B., Gnilka, C. L., et al. 2021, Front. Astron. Space Sci., 8, 138

Seifahrt, A., Stürmer, J., Bean, J. L., & Schwab, C. 2018, SPIE Conf. Ser., 10702, 107026D

Seifahrt, A., Bean, J. L., Kasper, D., et al. 2022, SPIE Conf. Ser., 12184, 121841G

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

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

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

Stalport, M., Delisle, J. B., Udry, S., et al. 2022, A&A, 664, A53

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

Steffen, J. H., & Hwang, J. A. 2015, MNRAS, 448, 1956

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

Tamayo, D., Rein, H., Shi, P., & Hernandez, D. M. 2020, MNRAS, 491, 2885

Tamura, M., Suto, H., Nishikawa, J., et al. 2012, SPIE Conf. Ser., 8446, 84461T

Trifonov, T., Caballero, J. A., Morales, J. C., et al. 2021, Science, 371, 1038

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

Van Eylen, V., Astudillo-Defru, N., Bonfils, X., et al. 2021, MNRAS, 507, 2154

Weiss, L. M., Marcy, G. W., Petigura, E. A., et al. 2018, AJ, 155, 48

Wizinowich, P., Acton, D. S., Shelton, C., et al. 2000, PASP, 112, 315

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

Wyatt, M. C., Kral, Q., & Sinclair, C. A. 2020, MNRAS, 491, 782

Yee, S. W., Petigura, E. A., & von Braun, K. 2017, ApJ, 836, 77

Zechmeister, M., Reiners, A., Amado, P. J., et al. 2018, A&A, 609, A12

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

Ziegler, C., Tokovinin, A., Briceño, C., et al. 2020, AJ, 159, 19

⁷ Institute for Computational Science, University of Zurich, Winterthurerstr. 190, 8057 Zurich, Switzerland

⁸ Institute for Particle Physics and Astrophysics, ETH Zürich, Wolfgang-Pauli-Strasse 27, 8093 Zurich, Switzerland

⁹ Departamento de Física, Universidade Federal do Rio Grande do Norte, 59072-970 Natal, RN, Brazil

¹⁰ Instituto de Astrofísica de Canarias (IAC), 38200 La Laguna, Tenerife, Spain

¹<fnref n="1" /> Dept. Astrofísica, Universidad de La Laguna (ULL), 38206 La Laguna, Tenerife, Spain

¹<fnref n="2" /> Department of Physics and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA

¹<fnref n="3" /> Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA

¹<fnref n="4" /> Department of Aeronautics and Astronautics, MIT, 77 Massachusetts Avenue, Cambridge, MA 02139, USA

¹<fnref n="5" /> NASA Sagan Fellow; Department of Astrophysical Sciences, Princeton University, NJ 08544, USA

¹<fnref n="6" /> NASA Ames Research Center, Moffett Field, CA 94035, USA

¹⁷ Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove drive, Pasadena, CA 91109, USA

¹⁸ U.S. Naval Observatory, 3450 Massachusetts Avenue NW, Washington, DC 20392, USA

¹⁹ Villa '39 Observatory, Landers, CA 92285, USA

²⁰ Department of Physics and Astronomy, Union College, 807 Union St., Schenectady, NY 12308, USA

²<fnref n="1" /> SUPA Physics and Astronomy, University of St. Andrews, Fife, KY16 9SS Scotland, UK

²<fnref n="2" /> Space sciences, Technologies and Astrophysics Research (STAR) Institute, University of Liège, Allée du 6 Août 19C, 4000 Liège, Belgium

²<fnref n="3" /> Perth Exoplanet Survey Telescope, Perth, Western Australia, Australia

²<fnref n="4" /> Komaba Institute for Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan

²<fnref n="5" /> Okayama Observatory, Kyoto University, 3037-5 Honjo, Kamogatacho, Asakuchi, Okayama 719-0232, Japan

²<fnref n="6" /> Department of Multi-Disciplinary Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo, Japan

²⁷ Department of Physical Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan

²⁸ Astrobiology Center, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan

²⁹ National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan

³⁰ Department of Astronomical Science, The Graduate University for Advanced Studies, SOKENDAI, 2-21-1, Osawa, Mitaka, Tokyo, 181-8588, Japan

³<fnref n="1" /> Department of Astronomy, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan

³<fnref n="2" /> Department of Astronomy and Astrophysics, The University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA

³<fnref n="3" /> Landessternwarte, Zentrum für Astronomie der Universität Heidelberg, Königstuhl 12, 69117 Heidelberg, Germany

³<fnref n="4" /> Center for Space and Habitability, University of Bern, Gesellschaftsstrasse 6, 3012 Bern, Switzerland

³<fnref n="5" /> Universidad Nacional Autónoma de Mexico, Instituto de Astronomía, AP 70-264, CDMX 04510, Mexico

³<fnref n="6" /> Universidad Nacional Autónoma de Mexico, Instituto de Astronomía, AP 106, Ensenada 22800, BC, Mexico

³⁷ SETI Institute, Mountain View, CA 94043, USA

³⁸ Department of Astronomy, University of Maryland, College Park, MD 20742, USA

<footnote n="1">¹ Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France e-mail: marion.cointepas@unige.ch</footnote>

<footnote n="2">² Observatoire de Genève, Département d'Astronomie, Université de Genève, Chemin Pegasi 51b, 1290 Versoix, Switzerland</footnote>

<footnote n="3">³ Departamento de Matemática y Física Aplicadas, Universidad Católica de la Santísima Concepción, Alonso de Rivera 2850, Concepción, Chile</footnote>

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

<footnote n="5">⁵ George Mason University, 4400 University Drive, Fairfax, VA, 22030 USA</footnote>

<footnote n="6">⁶ Department of Astronomy and Tsinghua Centre for Astrophysics, Tsinghua University, Beijing 100084, PR China</footnote>
Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20

McCully, C., Volgenau, N. H., Harbeck, D.-R., et al. 2018, SPIE Conf. Ser., 10707, 107070K

Millholland, S., Wang, S., & Laughlin, G. 2017, ApJ, 849, L33

Misener, W., Schlichting, H. E., & Young, E. D. 2023, MNRAS, 524, 981

Mordasini, C. 2020, A&A, 638, A52

Mordasini, C., Alibert, Y., & Benz, W. 2009, A&A, 501, 1139

Mugrauer, M., & Michel, K.-U. 2020, Astron. Nachr., 341, 996

Mugrauer, M., & Michel, K.-U. 2021, Astron. Nachr., 342, 840

Narita, N., Fukui, A., Kusakabe, N., et al. 2015, J. Astron. Teles. Instrum. Syst., 1, 045001

Narita, N., Fukui, A., Kusakabe, N., et al. 2019, J. Astron. Teles. Instrum. Syst., 5, 015001

Otegi, J. F., Bouchy, F., & Helled, R. 2020, A&A, 634, A43

Otegi, J. F., Helled, R., & Bouchy, F. 2022, A&A, 658, A107

Owen, J. E., & Wu, Y. 2017, ApJ, 847, 29

Parviainen, H., Palle, E., Zapatero-Osorio, M. R., et al. 2020, A&A, 633, A28

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

Pepe, F., Cristiani, S., Rebolo, R., et al. 2021, A&A, 645, A96

Pojmanski, G. 2002, Acta Astron., 52, 397

Reefe, M., Alfaro, O., Foster, S., et al. 2022, J. Astron. Teles. Instrum. Syst., 8, 027002

Rein, H., & Liu, S. F. 2012, A&A, 537, A128

Rein, H., & Spiegel, D. S. 2015, MNRAS, 446, 1424

Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telesc. Instrum. Syst., 1, 014003

Rogers, J. G., Schlichting, H. E., & Owen, J. E. 2023, ApJ, 947, L19

Schlichting, H. E., Sari, R., & Yalinewich, A. 2015, Icarus, 247, 81

Scott, N. J., Howell, S. B., Gnilka, C. L., et al. 2021, Front. Astron. Space Sci., 8, 138

Seifahrt, A., Stürmer, J., Bean, J. L., & Schwab, C. 2018, SPIE Conf. Ser., 10702, 107026D

Seifahrt, A., Bean, J. L., Kasper, D., et al. 2022, SPIE Conf. Ser., 12184, 121841G

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

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

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

Stalport, M., Delisle, J. B., Udry, S., et al. 2022, A&A, 664, A53

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

Steffen, J. H., & Hwang, J. A. 2015, MNRAS, 448, 1956

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

Tamayo, D., Rein, H., Shi, P., & Hernandez, D. M. 2020, MNRAS, 491, 2885

Tamura, M., Suto, H., Nishikawa, J., et al. 2012, SPIE Conf. Ser., 8446, 84461T

Trifonov, T., Caballero, J. A., Morales, J. C., et al. 2021, Science, 371, 1038

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

Van Eylen, V., Astudillo-Defru, N., Bonfils, X., et al. 2021, MNRAS, 507, 2154

Weiss, L. M., Marcy, G. W., Petigura, E. A., et al. 2018, AJ, 155, 48

Wizinowich, P., Acton, D. S., Shelton, C., et al. 2000, PASP, 112, 315

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

Wyatt, M. C., Kral, Q., & Sinclair, C. A. 2020, MNRAS, 491, 782

Yee, S. W., Petigura, E. A., & von Braun, K. 2017, ApJ, 836, 77

Zechmeister, M., Reiners, A., Amado, P. J., et al. 2018, A&A, 609, A12

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

Ziegler, C., Tokovinin, A., Briceño, C., et al. 2020, AJ, 159, 19

⁷ Institute for Computational Science, University of Zurich, Winterthurerstr. 190, 8057 Zurich, Switzerland

⁸ Institute for Particle Physics and Astrophysics, ETH Zürich, Wolfgang-Pauli-Strasse 27, 8093 Zurich, Switzerland

⁹ Departamento de Física, Universidade Federal do Rio Grande do Norte, 59072-970 Natal, RN, Brazil

¹⁰ Instituto de Astrofísica de Canarias (IAC), 38200 La Laguna, Tenerife, Spain

¹<fnref n="1" /> Dept. Astrofísica, Universidad de La Laguna (ULL), 38206 La Laguna, Tenerife, Spain

¹<fnref n="2" /> Department of Physics and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA

¹<fnref n="3" /> Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA

¹<fnref n="4" /> Department of Aeronautics and Astronautics, MIT, 77 Massachusetts Avenue, Cambridge, MA 02139, USA

¹<fnref n="5" /> NASA Sagan Fellow; Department of Astrophysical Sciences, Princeton University, NJ 08544, USA

¹<fnref n="6" /> NASA Ames Research Center, Moffett Field, CA 94035, USA

¹⁷ Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove drive, Pasadena, CA 91109, USA

¹⁸ U.S. Naval Observatory, 3450 Massachusetts Avenue NW, Washington, DC 20392, USA

¹⁹ Villa '39 Observatory, Landers, CA 92285, USA

²⁰ Department of Physics and Astronomy, Union College, 807 Union St., Schenectady, NY 12308, USA

²<fnref n="1" /> SUPA Physics and Astronomy, University of St. Andrews, Fife, KY16 9SS Scotland, UK

²<fnref n="2" /> Space sciences, Technologies and Astrophysics Research (STAR) Institute, University of Liège, Allée du 6 Août 19C, 4000 Liège, Belgium

²<fnref n="3" /> Perth Exoplanet Survey Telescope, Perth, Western Australia, Australia

²<fnref n="4" /> Komaba Institute for Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan

²<fnref n="5" /> Okayama Observatory, Kyoto University, 3037-5 Honjo, Kamogatacho, Asakuchi, Okayama 719-0232, Japan

²<fnref n="6" /> Department of Multi-Disciplinary Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo, Japan

²⁷ Department of Physical Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan

²⁸ Astrobiology Center, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan

²⁹ National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan

³⁰ Department of Astronomical Science, The Graduate University for Advanced Studies, SOKENDAI, 2-21-1, Osawa, Mitaka, Tokyo, 181-8588, Japan

³<fnref n="1" /> Department of Astronomy, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan

³<fnref n="2" /> Department of Astronomy and Astrophysics, The University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA

³<fnref n="3" /> Landessternwarte, Zentrum für Astronomie der Universität Heidelberg, Königstuhl 12, 69117 Heidelberg, Germany

³<fnref n="4" /> Center for Space and Habitability, University of Bern, Gesellschaftsstrasse 6, 3012 Bern, Switzerland

³<fnref n="5" /> Universidad Nacional Autónoma de Mexico, Instituto de Astronomía, AP 70-264, CDMX 04510, Mexico

³<fnref n="6" /> Universidad Nacional Autónoma de Mexico, Instituto de Astronomía, AP 106, Ensenada 22800, BC, Mexico

³⁷ SETI Institute, Mountain View, CA 94043, USA

³⁸ Department of Astronomy, University of Maryland, College Park, MD 20742, USA

<footnote n="1">¹ Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France e-mail: marion.cointepas@unige.ch</footnote>

<footnote n="2">² Observatoire de Genève, Département d'Astronomie, Université de Genève, Chemin Pegasi 51b, 1290 Versoix, Switzerland</footnote>

<footnote n="3">³ Departamento de Matemática y Física Aplicadas, Universidad Católica de la Santísima Concepción, Alonso de Rivera 2850, Concepción, Chile</footnote>

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

<footnote n="5">⁵ George Mason University, 4400 University Drive, Fairfax, VA, 22030 USA</footnote>

<footnote n="6">⁶ Department of Astronomy and Tsinghua Centre for Astrophysics, Tsinghua University, Beijing 100084, PR China</footnote>
Appendix A: Radial velocity measurements

<table><caption>Table A.1. RV measurements for the different spectrographs used in this analysis.</caption><thead><tr><th>Time (BJD - 2400000)</th><th>RV (km/s)</th><th>RV err (km/s)</th><th>Instrument</th></tr></thead><tbody><tr><td>58910.723746</td><td>34.324020</td><td>0.004000</td><td>HARPS</td></tr><tr><td>58911.771067</td><td>34.318880</td><td>0.004440</td><td>HARPS</td></tr><tr><td>58912.716950</td><td>34.329700</td><td>0.005490</td><td>HARPS</td></tr><tr><td>58915.705378</td><td>34.323300</td><td>0.004940</td><td>HARPS</td></tr><tr><td>58916.751278</td><td>34.326390</td><td>0.003840</td><td>HARPS</td></tr><tr><td>58917.771100</td><td>34.334450</td><td>0.004290</td><td>HARPS</td></tr><tr><td>58918.715956</td><td>34.327300</td><td>0.004800</td><td>HARPS</td></tr><tr><td>58924.721557</td><td>34.323460</td><td>0.004360</td><td>HARPS</td></tr><tr><td>58925.718218</td><td>34.324960</td><td>0.005640</td><td>HARPS</td></tr><tr><td>58926.725780</td><td>34.323560</td><td>0.004520</td><td>HARPS</td></tr><tr><td>58927.736396</td><td>34.322920</td><td>0.004490</td><td>HARPS</td></tr><tr><td>58928.709603</td><td>34.321810</td><td>0.004390</td><td>HARPS</td></tr><tr><td>58929.717392</td><td>34.314330</td><td>0.007000</td><td>HARPS</td></tr><tr><td>58930.707725</td><td>34.312610</td><td>0.008150</td><td>HARPS</td></tr><tr><td>58931.737397</td><td>34.322390</td><td>0.005710</td><td>HARPS</td></tr><tr><td>59171.848568</td><td>34.323730</td><td>0.004830</td><td>HARPS</td></tr><tr><td>59176.826197</td><td>34.334770</td><td>0.005100</td><td>HARPS</td></tr><tr><td>59196.785378</td><td>34.321760</td><td>0.004010</td><td>HARPS</td></tr><tr><td>59197.791841</td><td>34.322360</td><td>0.004700</td><td>HARPS</td></tr><tr><td>59198.795178</td><td>34.332740</td><td>0.004190</td><td>HARPS</td></tr><tr><td>59220.777917</td><td>34.336650</td><td>0.004230</td><td>HARPS</td></tr><tr><td>59221.761564</td><td>34.337840</td><td>0.003190</td><td>HARPS</td></tr><tr><td>59223.747478</td><td>34.335450</td><td>0.003470</td><td>HARPS</td></tr><tr><td>59224.732789</td><td>34.328460</td><td>0.003270</td><td>HARPS</td></tr><tr><td>59262.759243</td><td>34.322320</td><td>0.003830</td><td>HARPS</td></tr><tr><td>59263.749363</td><td>34.340640</td><td>0.004060</td><td>HARPS</td></tr><tr><td>59264.712802</td><td>34.333680</td><td>0.003770</td><td>HARPS</td></tr><tr><td>59265.741670</td><td>34.324630</td><td>0.004580</td><td>HARPS</td></tr><tr><td>59266.698914</td><td>34.329960</td><td>0.004020</td><td>HARPS</td></tr><tr><td>59267.741355</td><td>34.332680</td><td>0.004740</td><td>HARPS</td></tr><tr><td>59268.788980</td><td>34.328840</td><td>0.004230</td><td>HARPS</td></tr><tr><td>59269.672306</td><td>34.338130</td><td>0.004470</td><td>HARPS</td></tr><tr><td>59286.794959</td><td>34.329480</td><td>0.004600</td><td>HARPS</td></tr><tr><td>59287.761950</td><td>34.322910</td><td>0.004330</td><td>HARPS</td></tr><tr><td>59288.761531</td><td>34.320690</td><td>0.006520</td><td>HARPS</td></tr><tr><td>59290.717159</td><td>34.319910</td><td>0.004200</td><td>HARPS</td></tr><tr><td>59366.521206</td><td>34.331450</td><td>0.005130</td><td>HARPS</td></tr><tr><td>59367.552764</td><td>34.326720</td><td>0.005090</td><td>HARPS</td></tr><tr><td>59369.521491</td><td>34.331190</td><td>0.005140</td><td>HARPS</td></tr><tr><td>59372.544956</td><td>34.338960</td><td>0.004120</td><td>HARPS</td></tr><tr><td>59674.553825</td><td>34.274761</td><td>0.001426</td><td>ESPRESSO</td></tr><tr><td>59675.576856</td><td>34.263909</td><td>0.001523</td><td>ESPRESSO</td></tr><tr><td>59676.670671</td><td>34.268350</td><td>0.002146</td><td>ESPRESSO</td></tr><tr><td>59677.512074</td><td>34.266289</td><td>0.003088</td><td>ESPRESSO</td></tr><tr><td>59678.659022</td><td>34.263947</td><td>0.002010</td><td>ESPRESSO</td></tr><tr><td>59679.625607</td><td>34.268052</td><td>0.001628</td><td>ESPRESSO</td></tr><tr><td>59680.566038</td><td>34.259502</td><td>0.001265</td><td>ESPRESSO</td></tr><tr><td>59682.524050</td><td>34.272004</td><td>0.001389</td><td>ESPRESSO</td></tr><tr><td>59685.570917</td><td>34.265878</td><td>0.002843</td><td>ESPRESSO</td></tr><tr><td>59687.548051</td><td>34.277600</td><td>0.001226</td><td>ESPRESSO</td></tr><tr><td>59267.876851</td><td>33.052791</td><td>0.002311</td><td>MAROON-X blue</td></tr><tr><td>59268.007143</td><td>33.048796</td><td>0.003156</td><td>MAROON-X blue</td></tr><tr><td>59269.852498</td><td>33.056575</td><td>0.002331</td><td>MAROON-X blue</td></tr><tr><td>59269.981134</td><td>33.052285</td><td>0.002319</td><td>MAROON-X blue</td></tr><tr><td>59276.915064</td><td>33.064205</td><td>0.002425</td><td>MAROON-X blue</td></tr><tr><td>59277.880398</td><td>33.057519</td><td>0.002386</td><td>MAROON-X blue</td></tr><tr><td>59278.032393</td><td>33.057208</td><td>0.002369</td><td>MAROON-X blue</td></tr><tr><td>59323.793794</td><td>33.059731</td><td>0.001462</td><td>MAROON-X blue</td></tr><tr><td>59327.768365</td><td>33.057727</td><td>0.001732</td><td>MAROON-X blue</td></tr><tr><td>59332.910732</td><td>33.063903</td><td>0.002816</td><td>MAROON-X blue</td></tr><tr><td>59333.769718</td><td>33.064511</td><td>0.001553</td><td>MAROON-X blue</td></tr><tr><td>59333.894284</td><td>33.067126</td><td>0.002444</td><td>MAROON-X blue</td></tr><tr><td>59334.771413</td><td>33.059801</td><td>0.001489</td><td>MAROON-X blue</td></tr><tr><td>59334.858713</td><td>33.059432</td><td>0.001687</td><td>MAROON-X blue</td></tr><tr><td>59267.876851</td><td>33.048031</td><td>0.001835</td><td>MAROON-X red</td></tr><tr><td>59268.007143</td><td>33.050574</td><td>0.002479</td><td>MAROON-X red</td></tr><tr><td>59269.852498</td><td>33.055762</td><td>0.001778</td><td>MAROON-X red</td></tr><tr><td>59269.981134</td><td>33.052448</td><td>0.001904</td><td>MAROON-X red</td></tr><tr><td>59276.915064</td><td>33.062313</td><td>0.001993</td><td>MAROON-X red</td></tr><tr><td>59277.880398</td><td>33.059240</td><td>0.001935</td><td>MAROON-X red</td></tr><tr><td>59278.032393</td><td>33.055978</td><td>0.001772</td><td>MAROON-X red</td></tr><tr><td>59323.793794</td><td>33.059169</td><td>0.001161</td><td>MAROON-X red</td></tr><tr><td>59327.768365</td><td>33.059375</td><td>0.001376</td><td>MAROON-X red</td></tr><tr><td>59332.910732</td><td>33.062870</td><td>0.001697</td><td>MAROON-X red</td></tr><tr><td>59333.769718</td><td>33.065211</td><td>0.001285</td><td>MAROON-X red</td></tr><tr><td>59333.894284</td><td>33.064958</td><td>0.001632</td><td>MAROON-X red</td></tr><tr><td>59334.771413</td><td>33.061625</td><td>0.001232</td><td>MAROON-X red</td></tr><tr><td>59334.858713</td><td>33.061741</td><td>0.001275</td><td>MAROON-X red</td></tr><tr><td>59248.003442</td><td>-0.00625</td><td>0.00376</td><td>IRD</td></tr><tr><td>59275.916993</td><td>-0.00858</td><td>0.00484</td><td>IRD</td></tr><tr><td>59275.975987</td><td>-0.00166</td><td>0.00446</td><td>IRD</td></tr><tr><td>59338.810860</td><td>0.00222</td><td>0.00353</td><td>IRD</td></tr><tr><td>59373.773176</td><td>0.00134</td><td>0.00364</td><td>IRD</td></tr></tbody></table>
Appendix A: Radial velocity measurements

<table><caption>Table A.1. RV measurements for the different spectrographs used in this analysis.</caption><thead><tr><th>Time (BJD - 2400000)</th><th>RV (km/s)</th><th>RV err (km/s)</th><th>Instrument</th></tr></thead><tbody><tr><td>58910.723746</td><td>34.324020</td><td>0.004000</td><td>HARPS</td></tr><tr><td>58911.771067</td><td>34.318880</td><td>0.004440</td><td>HARPS</td></tr><tr><td>58912.716950</td><td>34.329700</td><td>0.005490</td><td>HARPS</td></tr><tr><td>58915.705378</td><td>34.323300</td><td>0.004940</td><td>HARPS</td></tr><tr><td>58916.751278</td><td>34.326390</td><td>0.003840</td><td>HARPS</td></tr><tr><td>58917.771100</td><td>34.334450</td><td>0.004290</td><td>HARPS</td></tr><tr><td>58918.715956</td><td>34.327300</td><td>0.004800</td><td>HARPS</td></tr><tr><td>58924.721557</td><td>34.323460</td><td>0.004360</td><td>HARPS</td></tr><tr><td>58925.718218</td><td>34.324960</td><td>0.005640</td><td>HARPS</td></tr><tr><td>58926.725780</td><td>34.323560</td><td>0.004520</td><td>HARPS</td></tr><tr><td>58927.736396</td><td>34.322920</td><td>0.004490</td><td>HARPS</td></tr><tr><td>58928.709603</td><td>34.321810</td><td>0.004390</td><td>HARPS</td></tr><tr><td>58929.717392</td><td>34.314330</td><td>0.007000</td><td>HARPS</td></tr><tr><td>58930.707725</td><td>34.312610</td><td>0.008150</td><td>HARPS</td></tr><tr><td>58931.737397</td><td>34.322390</td><td>0.005710</td><td>HARPS</td></tr><tr><td>59171.848568</td><td>34.323730</td><td>0.004830</td><td>HARPS</td></tr><tr><td>59176.826197</td><td>34.334770</td><td>0.005100</td><td>HARPS</td></tr><tr><td>59196.785378</td><td>34.321760</td><td>0.004010</td><td>HARPS</td></tr><tr><td>59197.791841</td><td>34.322360</td><td>0.004700</td><td>HARPS</td></tr><tr><td>59198.795178</td><td>34.332740</td><td>0.004190</td><td>HARPS</td></tr><tr><td>59220.777917</td><td>34.336650</td><td>0.004230</td><td>HARPS</td></tr><tr><td>59221.761564</td><td>34.337840</td><td>0.003190</td><td>HARPS</td></tr><tr><td>59223.747478</td><td>34.335450</td><td>0.003470</td><td>HARPS</td></tr><tr><td>59224.732789</td><td>34.328460</td><td>0.003270</td><td>HARPS</td></tr><tr><td>59262.759243</td><td>34.322320</td><td>0.003830</td><td>HARPS</td></tr><tr><td>59263.749363</td><td>34.340640</td><td>0.004060</td><td>HARPS</td></tr><tr><td>59264.712802</td><td>34.333680</td><td>0.003770</td><td>HARPS</td></tr><tr><td>59265.741670</td><td>34.324630</td><td>0.004580</td><td>HARPS</td></tr><tr><td>59266.698914</td><td>34.329960</td><td>0.004020</td><td>HARPS</td></tr><tr><td>59267.741355</td><td>34.332680</td><td>0.004740</td><td>HARPS</td></tr><tr><td>59268.788980</td><td>34.328840</td><td>0.004230</td><td>HARPS</td></tr><tr><td>59269.672306</td><td>34.338130</td><td>0.004470</td><td>HARPS</td></tr><tr><td>59286.794959</td><td>34.329480</td><td>0.004600</td><td>HARPS</td></tr><tr><td>59287.761950</td><td>34.322910</td><td>0.004330</td><td>HARPS</td></tr><tr><td>59288.761531</td><td>34.320690</td><td>0.006520</td><td>HARPS</td></tr><tr><td>59290.717159</td><td>34.319910</td><td>0.004200</td><td>HARPS</td></tr><tr><td>59366.521206</td><td>34.331450</td><td>0.005130</td><td>HARPS</td></tr><tr><td>59367.552764</td><td>34.326720</td><td>0.005090</td><td>HARPS</td></tr><tr><td>59369.521491</td><td>34.331190</td><td>0.005140</td><td>HARPS</td></tr><tr><td>59372.544956</td><td>34.338960</td><td>0.004120</td><td>HARPS</td></tr><tr><td>59674.553825</td><td>34.274761</td><td>0.001426</td><td>ESPRESSO</td></tr><tr><td>59675.576856</td><td>34.263909</td><td>0.001523</td><td>ESPRESSO</td></tr><tr><td>59676.670671</td><td>34.268350</td><td>0.002146</td><td>ESPRESSO</td></tr><tr><td>59677.512074</td><td>34.266289</td><td>0.003088</td><td>ESPRESSO</td></tr><tr><td>59678.659022</td><td>34.263947</td><td>0.002010</td><td>ESPRESSO</td></tr><tr><td>59679.625607</td><td>34.268052</td><td>0.001628</td><td>ESPRESSO</td></tr><tr><td>59680.566038</td><td>34.259502</td><td>0.001265</td><td>ESPRESSO</td></tr><tr><td>59682.524050</td><td>34.272004</td><td>0.001389</td><td>ESPRESSO</td></tr><tr><td>59685.570917</td><td>34.265878</td><td>0.002843</td><td>ESPRESSO</td></tr><tr><td>59687.548051</td><td>34.277600</td><td>0.001226</td><td>ESPRESSO</td></tr><tr><td>59267.876851</td><td>33.052791</td><td>0.002311</td><td>MAROON-X blue</td></tr><tr><td>59268.007143</td><td>33.048796</td><td>0.003156</td><td>MAROON-X blue</td></tr><tr><td>59269.852498</td><td>33.056575</td><td>0.002331</td><td>MAROON-X blue</td></tr><tr><td>59269.981134</td><td>33.052285</td><td>0.002319</td><td>MAROON-X blue</td></tr><tr><td>59276.915064</td><td>33.064205</td><td>0.002425</td><td>MAROON-X blue</td></tr><tr><td>59277.880398</td><td>33.057519</td><td>0.002386</td><td>MAROON-X blue</td></tr><tr><td>59278.032393</td><td>33.057208</td><td>0.002369</td><td>MAROON-X blue</td></tr><tr><td>59323.793794</td><td>33.059731</td><td>0.001462</td><td>MAROON-X blue</td></tr><tr><td>59327.768365</td><td>33.057727</td><td>0.001732</td><td>MAROON-X blue</td></tr><tr><td>59332.910732</td><td>33.063903</td><td>0.002816</td><td>MAROON-X blue</td></tr><tr><td>59333.769718</td><td>33.064511</td><td>0.001553</td><td>MAROON-X blue</td></tr><tr><td>59333.894284</td><td>33.067126</td><td>0.002444</td><td>MAROON-X blue</td></tr><tr><td>59334.771413</td><td>33.059801</td><td>0.001489</td><td>MAROON-X blue</td></tr><tr><td>59334.858713</td><td>33.059432</td><td>0.001687</td><td>MAROON-X blue</td></tr><tr><td>59267.876851</td><td>33.048031</td><td>0.001835</td><td>MAROON-X red</td></tr><tr><td>59268.007143</td><td>33.050574</td><td>0.002479</td><td>MAROON-X red</td></tr><tr><td>59269.852498</td><td>33.055762</td><td>0.001778</td><td>MAROON-X red</td></tr><tr><td>59269.981134</td><td>33.052448</td><td>0.001904</td><td>MAROON-X red</td></tr><tr><td>59276.915064</td><td>33.062313</td><td>0.001993</td><td>MAROON-X red</td></tr><tr><td>59277.880398</td><td>33.059240</td><td>0.001935</td><td>MAROON-X red</td></tr><tr><td>59278.032393</td><td>33.055978</td><td>0.001772</td><td>MAROON-X red</td></tr><tr><td>59323.793794</td><td>33.059169</td><td>0.001161</td><td>MAROON-X red</td></tr><tr><td>59327.768365</td><td>33.059375</td><td>0.001376</td><td>MAROON-X red</td></tr><tr><td>59332.910732</td><td>33.062870</td><td>0.001697</td><td>MAROON-X red</td></tr><tr><td>59333.769718</td><td>33.065211</td><td>0.001285</td><td>MAROON-X red</td></tr><tr><td>59333.894284</td><td>33.064958</td><td>0.001632</td><td>MAROON-X red</td></tr><tr><td>59334.771413</td><td>33.061625</td><td>0.001232</td><td>MAROON-X red</td></tr><tr><td>59334.858713</td><td>33.061741</td><td>0.001275</td><td>MAROON-X red</td></tr><tr><td>59248.003442</td><td>-0.00625</td><td>0.00376</td><td>IRD</td></tr><tr><td>59275.916993</td><td>-0.00858</td><td>0.00484</td><td>IRD</td></tr><tr><td>59275.975987</td><td>-0.00166</td><td>0.00446</td><td>IRD</td></tr><tr><td>59338.810860</td><td>0.00222</td><td>0.00353</td><td>IRD</td></tr><tr><td>59373.773176</td><td>0.00134</td><td>0.00364</td><td>IRD</td></tr></tbody></table>
Appendix B: Corner plot of some parameters from the joint fit of all data

<figure><img src="image_16.png" /><figcaption>**Fig. B.1.** Corner plots showing the posterior distributions of the period (days), time of transit center (BJD<sub>TDB</sub>), and RV semi-amplitude (km/s) for the three planets from the joint modeling of photometry and RV datasets.</figcaption></figure>
## Appendix B: Corner plot of some parameters from the joint fit of all data

<figure><img src="image_16.png" /><figcaption>**Fig. B.1.** Corner plots showing the posterior distributions of the period (days), time of transit center ($BJD_{TDB}$), and RV semi-amplitude (km/s) for the three planets from the joint modeling of photometry and RV datasets.</figcaption></figure>
## Appendix C: Jitter posterior parameters obtained from our joint analysis

<table><caption>**Table C.1.** Jitter posterior parameters obtained from our joint photometric and RV *juliet* analysis for TOI-663.</caption><thead><tr><th>Parameter name</th><th>Posterior estimate</th></tr></thead><tbody><tr><td colspan="2"><strong>Posterior parameters for RV instruments</strong></td></tr><tr><td>σ<sub>w,HARPS</sub> (m/s)</td><td>0.29<sup>+1.68</sup><sub>-0.26</sub></td></tr><tr><td>σ<sub>w,ESPRESSO</sub> (m/s)</td><td>0.13<sup>+0.65</sup><sub>-0.11</sub></td></tr><tr><td>σ<sub>w,MAROON-Xblue</sub> (m/s)</td><td>0.11<sup>+0.46</sup><sub>-0.09</sub></td></tr><tr><td>σ<sub>w,MAROON-Xred</sub> (m/s)</td><td>0.10<sup>+0.36</sup><sub>-0.26</sub></td></tr><tr><td>σ<sub>w,IRD</sub> (m/s)</td><td>0.25<sup>+1.98</sup><sub>-0.23</sub></td></tr><tr><td colspan="2"><strong>Posterior parameters for photometry instruments</strong></td></tr><tr><td>σ<sub>w,TESS</sub> (ppm)</td><td>8.5<sup>+28.9</sup><sub>-6.5</sub></td></tr><tr><td>σ<sub>w,LCOn1</sub> (ppm)</td><td>213.7<sup>+411.5</sup><sub>-206.5</sub></td></tr><tr><td>σ<sub>w,LCOn2</sub> (ppm)</td><td>948.4<sup>+38.0</sup><sub>-73.4</sub></td></tr><tr><td>σ<sub>w,LCOn3</sub> (ppm)</td><td>901.4<sup>+70.3</sup><sub>-120.6</sub></td></tr><tr><td>σ<sub>w,LCOn4</sub> (ppm)</td><td>885.1<sup>+82.6</sup><sub>-142.3</sub></td></tr><tr><td>σ<sub>w,LCOn5</sub> (ppm)</td><td>792.0<sup>+91.1</sup><sub>-97.4</sub></td></tr><tr><td>σ<sub>w,LCOn6</sub> (ppm)</td><td>842.2<sup>+104.9</sup><sub>-149.4</sub></td></tr><tr><td>σ<sub>w,MuSCATg</sub> (ppm)</td><td>594.9<sup>+333.2</sup><sub>-577.6</sub></td></tr><tr><td>σ<sub>w,MuSCATtr</sub> (ppm)</td><td>720.1<sup>+204.5</sup><sub>-659.1</sub></td></tr><tr><td>σ<sub>w,MuSCATz</sub> (ppm)</td><td>920.2<sup>+56.3</sup><sub>-90.8</sub></td></tr><tr><td>σ<sub>w,ExTrA2n1</sub> (ppm)</td><td>27.9<sup>+209.6</sup><sub>-24.8</sub></td></tr><tr><td>σ<sub>w,ExTrA3n1</sub> (ppm)</td><td>27.4<sup>+218.9</sup><sub>-24.3</sub></td></tr><tr><td>σ<sub>w,ExTrA2n2</sub> (ppm)</td><td>28.2<sup>+230.9</sup><sub>-25.2</sub></td></tr><tr><td>σ<sub>w,ExTrA3n2</sub> (ppm)</td><td>28.8<sup>+213.5</sup><sub>-25.7</sub></td></tr><tr><td>σ<sub>w,ExTrA2n3</sub> (ppm)</td><td>24.5<sup>+249.4</sup><sub>-21.7</sub></td></tr><tr><td>σ<sub>w,ExTrA3n3</sub> (ppm)</td><td>28.6<sup>+227.4</sup><sub>-25.5</sub></td></tr><tr><td>σ<sub>w,ExTrA3n4</sub> (ppm)</td><td>26.7<sup>+246.0</sup><sub>-23.6</sub></td></tr><tr><td>σ<sub>w,ExTrA2n5</sub> (ppm)</td><td>28.9<sup>+258.4</sup><sub>-25.9</sub></td></tr><tr><td>σ<sub>w,ExTrA3n5</sub> (ppm)</td><td>30.0<sup>+285.3</sup><sub>-27.0</sub></td></tr></tbody></table>

The posterior estimate corresponds to the median value. Error bars denote the 68% posterior credible intervals.

<footnote>(*)To be read as ExTrA, telescope, night number.</footnote>
## Appendix C: Jitter posterior parameters obtained from our joint analysis

<table><caption>Table C.1. Jitter posterior parameters obtained from our joint photometric and RV juliet analysis for TOI-663.</caption><thead><tr><th>Parameter name</th><th>Posterior estimate</th></tr></thead><tbody><tr><td colspan="2"><strong>Posterior parameters for RV instruments</strong></td></tr><tr><td>σ<sub>w,HARPS</sub> (m/s)</td><td>0.29<sup>+1.68</sup><sub>-0.26</sub></td></tr><tr><td>σ<sub>w,ESPRESSO</sub> (m/s)</td><td>0.13<sup>+0.65</sup><sub>-0.11</sub></td></tr><tr><td>σ<sub>w,MAROON-Xblue</sub> (m/s)</td><td>0.11<sup>+0.46</sup><sub>-0.09</sub></td></tr><tr><td>σ<sub>w,MAROON-Xred</sub> (m/s)</td><td>0.10<sup>+0.36</sup><sub>-0.26</sub></td></tr><tr><td>σ<sub>w,IRD</sub> (m/s)</td><td>0.25<sup>+1.98</sup><sub>-0.23</sub></td></tr><tr><td colspan="2"><strong>Posterior parameters for photometry instruments</strong></td></tr><tr><td>σ<sub>w,TESS</sub> (ppm)</td><td>8.5<sup>+28.9</sup><sub>-6.5</sub></td></tr><tr><td>σ<sub>w,LCOn1</sub> (ppm)</td><td>213.7<sup>+411.5</sup><sub>-206.5</sub></td></tr><tr><td>σ<sub>w,LCOn2</sub> (ppm)</td><td>948.4<sup>+38.0</sup><sub>-73.4</sub></td></tr><tr><td>σ<sub>w,LCOn3</sub> (ppm)</td><td>901.4<sup>+70.3</sup><sub>-120.6</sub></td></tr><tr><td>σ<sub>w,LCOn4</sub> (ppm)</td><td>885.1<sup>+82.6</sup><sub>-142.3</sub></td></tr><tr><td>σ<sub>w,LCOn5</sub> (ppm)</td><td>792.0<sup>+91.1</sup><sub>-97.4</sub></td></tr><tr><td>σ<sub>w,LCOn6</sub> (ppm)</td><td>842.2<sup>+104.9</sup><sub>-149.4</sub></td></tr><tr><td>σ<sub>w,MuSCATg</sub> (ppm)</td><td>594.9<sup>+333.2</sup><sub>-577.6</sub></td></tr><tr><td>σ<sub>w,MuSCATtr</sub> (ppm)</td><td>720.1<sup>+204.5</sup><sub>-659.1</sub></td></tr><tr><td>σ<sub>w,MuSCATz</sub> (ppm)</td><td>920.2<sup>+56.3</sup><sub>-90.8</sub></td></tr><tr><td>σ<sub>w,ExTrA2n1</sub> (ppm)</td><td>27.9<sup>+209.6</sup><sub>-24.8</sub></td></tr><tr><td>σ<sub>w,ExTrA3n1</sub> (ppm)</td><td>27.4<sup>+218.9</sup><sub>-24.3</sub></td></tr><tr><td>σ<sub>w,ExTrA2n2</sub> (ppm)</td><td>28.2<sup>+230.9</sup><sub>-25.2</sub></td></tr><tr><td>σ<sub>w,ExTrA3n2</sub> (ppm)</td><td>28.8<sup>+213.5</sup><sub>-25.7</sub></td></tr><tr><td>σ<sub>w,ExTrA2n3</sub> (ppm)</td><td>24.5<sup>+249.4</sup><sub>-21.7</sub></td></tr><tr><td>σ<sub>w,ExTrA3n3</sub> (ppm)</td><td>28.6<sup>+227.4</sup><sub>-25.5</sub></td></tr><tr><td>σ<sub>w,ExTrA3n4</sub> (ppm)</td><td>26.7<sup>+246.0</sup><sub>-23.6</sub></td></tr><tr><td>σ<sub>w,ExTrA2n5</sub> (ppm)</td><td>28.9<sup>+258.4</sup><sub>-25.9</sub></td></tr><tr><td>σ<sub>w,ExTrA3n5</sub> (ppm)</td><td>30.0<sup>+285.3</sup><sub>-27.0</sub></td></tr></tbody></table>

The posterior estimate corresponds to the median value. Error bars denote the 68% posterior credible intervals.

<footnote>(*)To be read as ExTrA, telescope, night number.</footnote>