# The GAPS Programme at TNG. XLII. A characterisation study of the multi-planet system around the 400 Myr-old star HD 63433 (TOI-1726)

M. Damasso, D. Locci, S. Benatti, A. Maggio, D. Nardiello, M. Baratella, K.

Biazzo, A. S. Bonomo, S. Desidera, V. d'Orazi, et al.

► To cite this version:

M. Damasso, D. Locci, S. Benatti, A. Maggio, D. Nardiello, et al., The GAPS Programme at TNG. XLII. A characterisation study of the multi-planet system around the 400 Myr-old star HD 63433 (TOI-1726). Astronomy & Astrophysics, 2023, 672, 10.1051/0004-6361/202245391. insu-04479094

HAL Id: insu-04479094

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

Submitted on 3 Mar 2024

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

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

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

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

## XLII. A characterisation study of the multi-planet system around the 400 Myr-old star HD 63433 (TOI-1726)\*

M. Damasso<fnref n="1" />, D. Locci<fnref n="2" />, S. Benatti<fnref n="2" />, A. Maggio<fnref n="2" />, D. Nardiello<fnref n="3" />,<fnref n="4" />, M. Baratella<fnref n="5" />, K. Biazzo<fnref n="6" />, A. S. Bonomo<fnref n="1" />, S. Desidera<fnref n="3" />, V. D'Orazi<fnref n="3" />, M. Mallonn<fnref n="5" />, A. F. Lanza<fnref n="7" />, A. Sozzetti<fnref n="1" />, F. Marzari<fnref n="3" />, F. Borsa<fnref n="8" />, J. Maldonado<fnref n="2" />, L. Mancini<fnref n="9" />,<fnref n="1" />,<fnref n="10" />, E. Poretti<fnref n="11" />, G. Scandariato<fnref n="7" />, A. Bignamini<fnref n="12" />, L. Borsato<fnref n="3" />, R. Capuzzo Dolcetta<fnref n="13" />, M. Cecconi<fnref n="11" />, R. Claudi<fnref n="3" />, R. Cosentino<fnref n="11" />, E. Covino<fnref n="14" />, A. Fiorenzano<fnref n="11" />, A. Harutyunyan<fnref n="11" />, A. W. Mann<fnref n="15" />, G. Micela<fnref n="2" />, E. Molinari<fnref n="16" />, M. Molinaro<fnref n="12" />, I. Pagano<fnref n="7" />, M. Pedani<fnref n="11" />, M. Pinamonti<fnref n="1" />, G. Piotto<fnref n="17" />, and H. Stoev<fnref n="11" />

Received 7 November 2022 / Accepted 31 January 2023

## ABSTRACT

**Context.** The GAPS collaboration is carrying out a spectroscopic and photometric follow-up of a sample of young stars with planets (age $\le600$ Myr) to characterise planetary systems at the early stages of their evolution.

**Aims.** For more than 2 yr, we monitored with the HARPS-N spectrograph the 400 Myr-old star HD 63433, which hosts two close-in (orbital periods $P_b \sim 7.1$ and $P_c \sim 20.5$ days) sub-Neptunes detected by the TESS space telescope, and it was announced in 2020. Using radial velocities and additional TESS photometry, we aim to provide the first measurement of their masses, improve the measure of their size and orbital parameters, and study the evolution of the atmospheric mass-loss rate due to photoevaporation.

**Methods.** We tested state-of-the-art analysis techniques and different models to mitigate the dominant signals due to stellar activity that are detected in the radial velocity time series. We used a hydro-based analytical description of the atmospheric mass-loss rate, coupled with a core-envelope model and stellar evolutionary tracks, to study the past and future evolution of the planetary masses and radii.

**Results.** We derived new measurements of the planetary orbital periods and radii ($P_b = 7.10794 \pm 0.000009$ days, $r_b = 2.02_{-0.05}^{+0.06} R_⊕$; $P_c = 20.54379 \pm 0.00002$ days, $r_c = 2.44 \pm 0.07 R_⊕$), and determined mass upper limits ($m_b \le 11 M_⊕$; $m_c \le 31 M_⊕$; 95% confidence level), with evidence at a $2.1–2.7\sigma$ significance level that HD 63433 c might be a dense mini-Neptune with a Neptune-like mass. For a grid of test masses below our derived dynamical upper limits, we found that HD 63433 b has very likely lost any gaseous H-He envelope, supporting HST-based observations that are indicative of there being no ongoing atmospheric will keep evaporating over the next ~5 Gyr if its current mass is $m_c \le 15 M_⊕$, while it should be hydrodynamically stable for higher masses.

**Key words.** stars: individual: HD63433 – planetary systems – techniques: photometric – techniques: radial velocities

<footnote n="1">¹ INAF – Osservatorio Astrofisico di Torino, Via Osservatorio 20, 10025 Pino Torinese, Italy
e-mail: mario.damasso@inaf.it</footnote>

<footnote n="2">² INAF – Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134, Palermo, Italy</footnote>

<footnote n="3">³ INAF – Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 35122, Padova, Italy</footnote>

<footnote n="4">⁴ Aix-Marseille Univ, CNRS, CNES, LAM, Marseille, France</footnote>

<footnote n="5">⁵ Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany</footnote>

<footnote n="6">⁶ INAF – Osservatorio Astronomico di Roma, Via Frascati 33, 00078 Monte Porzio Catone (Roma), Italy</footnote>

<footnote n="7">⁷ INAF – Osservatorio Astrofisico di Catania, Via S. Sofia 78, 95123 Catania, Italy</footnote>

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

<footnote n="9">⁹ Department of Physics, University of Rome “Tor Vergata”, Via della Ricerca Scientifica 1, 00133, Rome, Italy</footnote>

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

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

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

<footnote n="13">¹³ Dipartimento di Fisica – Università di Roma La Sapienza, P.le A. Moro 5, 00185 Roma, Italy</footnote>

<footnote n="14">¹⁴ INAF – Osservatorio Astronomico di Capodimonte, Salita Moiariello 16, 80131, Napoli, Italy</footnote>

<footnote n="15">¹⁵ Department of Physics and Astronomy, The University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA</footnote>

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

<footnote n="17">¹⁷ Dipartimento di Fisica e Astronomia “G. Galilei” – Università degli Studi di Padova, Vicolo dell'Osservatorio 3, 35122 Padova, Italy</footnote>

<footnote>\* Based on observations made with the Italian Telescopio Nazionale Galileo (TNG) operated by the Fundación Galileo Galilei (FGG) of the Istituto Nazionale di Astrofisica (INAF) at the Observatorio del Roque de los Muchachos (La Palma, Canary Islands, Spain).</footnote>
# 1. Introduction

Thanks to the constantly increasing number of characterised extrasolar planets, there is a growing number of studies related to the exoplanet demographics, defined by Gaudi et al. (2021) as a research area mainly dedicated to study planets' distribution as a function of the physical parameters that may influence planet formation and evolution, exploring a broad range of the parameter space. Among the many physical parameters involved, a particularly crucial one is the age of the exoplanet systems. Drawing a relationship between the observed properties of the systems and their ages is a necessary step to take in order to build a general picture of how planets form and evolve with time, and to constrain the timescales of the main processes at work within the first hundreds of millions years from planet formation. That is why it is crucial to detect planets at different stages of their evolution (infant: age <20 Myr; young: 20< age <100 Myr; intermediate-age: up to 800 Myr), and determine their orbital and main physical parameters. Among several key aspects, characterising young and close-in super-Earths and sub-Neptunes spotted at different ages offers the opportunity of theoretically studying how the atmospheric mass-loss rate evolves with time, in response to the high stellar high-energy irradiation, after dissipation of the protoplanetary disk.

Currently, the statistical sample of exoplanet systems younger than ~800 Myr is still small if compared to the total number of discoveries achieved so far. This prevents drawing results from comparative analyses with more evolved systems. Detecting infant and young planets is particularly challenging, especially for blind searches carried out with the radial velocity (RV) method because of the very high level of scatter (up to several hundreds of m s⁻¹) mostly due to stellar magnetic activity, which heavily hampers the detection of periodic signals produced by planetary companions. The photometric transit method proved to be a more fruitful detection technique in this case, in particular thanks to the observations of the *Kepler/K2* and the *Transiting Exoplanet Survey Satellite (TESS)* space telescopes. They yielded discoveries of planetary-sized companions, which triggered RV follow-up campaigns to measure their masses and bulk densities. Among these, the long-term Italian project Global Architecture of Planetary Systems (GAPS; Covino et al. 2013; Poretti et al. 2016) is investing a good fraction of the allocated observing time to follow-up stars with an age between ~2–600 Myr using the High Accuracy Radial velocity Planet Searcher (HARPS-N) spectrograph (Cosentino et al. 2012) located in the Northern Hemisphere at Telescopio Nazionale Galileo (TNG). After an initial effort in revising previously announced and debated planet detections (e.g. Carleo et al. 2018, 2020; Damasso et al. 2020), the main focus of the GAPS young planets sub-programme shifted to the measurement of the fundamental properties of transiting planets. GAPS has observed and characterised several systems so far (e.g. Carleo et al. 2021; Suárez Mascareño et al. 2021; Nardiello et al. 2022), and further characterisation studies have been carried out in parallel in the Southern Hemisphere with the HARPS spectrograph (e.g. Benatti et al. 2021; Desidera et al. 2023). The measurement of the planetary dynamical masses remains crucial to investigate the evolution of exo-atmospheres even when it is only possible to provide upper limits, and the current escape rates can be assessed through atmospheric evaporation models (e.g. Benatti et al. 2021; Carleo et al. 2021; Poppenhaeger et al. 2021; Maggio et al. 2022).

In this paper, we present a characterisation study of the two-planet system detected by the TESS space telescope (Ricker

et al. 2016) around the bright and 400-Myr old Sun-like star HD 63433, a member of the Ursa Major moving group, also known as TESS object of interest (TOI) 1726. The system was first validated and partly characterised by Mann et al. (2020). In particular, they revealed that the orbit of the innermost planet is prograde, and emphasised that HD 63433 represents a benchmark system for investigating the planet's evolution after the first hundreds millions years. The first results for this system led the GAPS collaboration to select HD 63433 as a promising target for an intense RV follow-up for further characterisation. The architecture of HD 63433 was also investigated by Dai et al. (2020), who measured the sky-projected spin-orbit angle of planet c, consistent with an aligned and prograde orbit. Results from Mann et al. (2020) and Dai et al. (2020) indicate that the system apparently did not experience catastrophic events capable of heavily influencing the dynamical evolution. The presence of a cold Jupiter in the system was examined by Hirsch et al. (2021), who analysed sparse RVs with a long time baseline, and found evidence only for a signal clearly ascribable to stellar activity.

The main purpose of the GAPS follow-up campaign of HD 63433 was measuring the masses of the two known planets. Accomplishing this goal becomes even more compelling to interpret the recent observational results of Zhang et al. (2022). They found evidence of Lyα absorption during a transit of HD 63433 c, suggesting ongoing atmospheric evaporation from this planet, while they suggest that HD 63433 b has already lost its primordial atmosphere: these observations have to be necessarily reconciled with mass measurements within the framework of current theoretical models. The study by Zhang et al. (2022) emphasises a key aspect: the two planets represent a great opportunity for a comparative study within the same system, to examine how their different orbital configurations led to the apparently diverse atmospheric mass-loss timescales.

In this work, we focus on the planetary mass measurements through the exploitation of RVs calculated from HARPS-N spectra, and we investigate the time evolution of the planetary atmospheres through evaporation models and updated stellar evolutionary tracks. The paper is organised as follows. In Sect. 2, we describe the photometric and spectroscopic datasets, and in Sect. 3 we present our measured fundamental stellar parameters. In Sect. 4, we present a frequency content analysis of the spectroscopic dataset, as a preparatory step to a more sophisticated analysis discussed in Sect. 5. In Sect. 6, we use the mass-radius diagram for a general comment on the possible planetary structures constrained by our results, especially in the context of known young exoplanets. Our results for planets b and c are then used for investigating the atmospheric evolutionary history in Sect. 7. Conclusions are summarised in Sect. 8.

# 2. Observations and data reduction

## 2.1. TESS photometry

TESS observed HD 63433 in Sector 20 (from 24 December 2019 to 21 January 2020) during Cycle 2 (GO22038, PI: Roettenbacher; GO22203, PI: Ge; GO22032, PI: Metcalfe), and again in Sectors 44–47 (from 12 October 2021 to 28 January 2022) during Cycle 4 (GO4104, PI: Huber; GO4242, PI: Mayo; GO4191, PI: Burt; GO4060, PI: ; GO4039, PI: Mann). In our study, we used the short-cadence (2 min) light curves. We did not make use of the official Presearch Data Conditioning Simple Aperture Photometry (PDCSAP, Smith et al. 2012; Stumpe et al. 2012, 2014) light curves, because they are affected by systematics due to overcorrections and/or injection of spurious signals. We corrected
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** TESS light curve of HD 63433 for Sector 20 (upper panel), and Sectors from 44 to 47 (lower panel). Planetary transits and bad-quality points have been removed from the dataset. The red curve represents the best-fit GP model.</figcaption></figure>

the Simple Aperture Photometry (SAP) data by using cotrend-
ing basis vectors obtained with the tools developed by Nardiello
et al. (2020). The extracted light curve for all Sectors is shown
in Fig. 1. In order to model and remove the clearly visible vari-
ability due to the stellar activity, we used the package PyORBIT
(Malavolta 2016; Malavolta et al. 2016, 2018). First, we masked
all the transits from the light curve, then we modelled the stel-
lar activity by using Gaussian rocess (GP) regression performed
with the package celerite2 (Foreman-Mackey et al. 2017;
Foreman-Mackey 2018). We interpolated each excluded point
within the masked transits using the best-fit GP solution, and we
corrected the light curve variability by dividing the observed flux
for the best-fit model. The final result is a flattened light curve
with all the transits preserved, which we used in the combined
light curve-RV analysis (Sect. 5). We already verified in previous
works (see, e.g. Nardiello 2020; Nardiello et al. 2021) that this
procedure for detrending the light curve does not alter the transit
signal, by affecting its shape and depth. At the time of writing,
according to the Web TESS Viewing Tool HD 63433 won't be
observed by TESS anymore up to September 2023 (Cycle 5).

## 2.2. STELLA photometry

HD 63433 was observed with the STELLA telescope and its
wide-field imager WiFSIP (Strassmeier et al. 2004) between 29
September 2020 and 17 March 2021 in Johnson V and Cousin
I bands (82 and 81 epochs, respectively). The follow-up was
intended to monitor the photometric evolution in response to
a changing activity during one of the spectroscopic observing
seasons. Each nightly observing block consisted of seven expo-
sures with an exposure time of 1 second for each band. To avoid

saturation of this bright target (V<sub>T</sub> ~ 7 mag), the telescope was
defocussed to get about 5 arcsec for the full width at half max-
imum (FWHM) of the object point spread function (PSF). The
data reduction followed that detailed description in Mallonn &
Strassmeier (2016). We averaged the nightly exposures for each
filter. The light curves are shown in Fig. C.1 (left column). The
scatter and median uncertainty for the V and I-band light curve
are 0.009 mag and 0.01 mag, respectively.

2.3. HARPS-N spectroscopic data

The spectroscopic follow-up of HD 63433 was carried out with
the HARPS-N spectrograph (Cosentino et al. 2012) mounted at
the Telescopio Nazionale Galileo (TNG) on the island of La
Palma (Canary Islands, Spain). The observations consist of 103
spectra collected between 26 February 2020 and 17 April 2022
(781 days), with a typical exposure time of 600 s, and a median
S/N of 203 measured at a wavelength of ~550 nm. The spectra
have been reduced with the standard Data Reduction Software
(DRS) pipeline (version 3.7.1) through the YABI workflow inter-
face (Hunter et al. 2012), which is maintained by the Italian
centre for Astronomical Archive (IA2)<fnref n="1" />. The RVs and activ-
ity diagnostics such as the FWHM and bisector velocity span
(BIS) have been derived from the DRS cross-correlation function
(CCF), which was calculated by adopting a reference mask for a
star of spectral type G2, and a half-window width of 200 km s⁻<fnref n="1" />.
The RVs have a median internal error σ<sub>RV, DRS</sub> = 0.81 m s⁻<fnref n="1" />, and
an rms of 23.3 m s⁻<fnref n="1" />.

As a sanity check, we measured the RVs also using the Tem-
plate Enhanced Radial velocity Reanalysis Application (TERRA)
pipeline (v1.8; Anglada-Escudé & Butler 2012), which, for
young and active stars, is in principle more effective in provid-
ing RVs less affected by stellar activity than those calculated
from the CCF. TERRA dataset is characterised by a median inter-
nal error σ<sub>RV, TERRA</sub> = 1.25 m s<sup>-1</sup>, and an rms of 24.1 m s<sup>-1</sup>. In
Sect. 5, we investigate the results obtained using both DRS and
TERRA RVs.

We calculated two additional chromospheric activity diag-
nostics from the HARPS-N spectra, namely the log R'_{HK} from the
CaII H&K lines (provided by the DRS), and the index based on
the H-alpha line calculated with the tool ACTIN (v1.3.9; Gomes
da Silva et al. 2018). RVs and activity diagnostics are listed in
Tables D.1 and D.2.

We note that there is very little overlap between the begin-
ning of the last season of the RVs and the final part of the TESS
light curve of Sector 47. The lack of simultaneous photometric
monitoring with the RVs makes it impossible to use the light
curve as an ancillary dataset to filter out the activity term in
the RVs. That is especially true for a star such as HD 63433,
which shows rapidly changing photospheric active regions in the
light curve that will generate RV variations independent of any
planet. The light curve observed by STELLA is characterised by
sparse sampling, and overlaps only partially with the second sea-
son covered by HARPS-N. Therefore, even the STELLA dataset
is not of practical use to effectively correct stellar activity in
the RVs.

**3. Fundamental stellar parameters**

HD 63433 is a dwarf star that belongs to the Ursa Major mov-
ing group, with an estimated age of roughly 400 Myr (see Mann
et al. 2020). Given its relatively young age, the spectroscopic
analysis might be hampered by the intense magnetic fields that

<footnote n="1">¹ https://ia2.inaf.it</footnote>
alter the structure of the stellar photosphere, in particular the upper layers. This alteration affects the formation of spectral lines in these layers for which abundances show a trend with optical depths (see Baratella et al. 2020b for further details). As a consequence, the standard spectroscopic approach, that is the equivalent width (EW) method, based on the use of iron (Fe) lines, could fail. In particular, when deriving the microturbulence velocity ($\xi$) by imposing that weak and strong iron lines have the same abundance, $\xi$ could be over-estimated, which leads to an under-estimation of the iron abundance ([Fe/H]).

Following the same strategy and the same line list as in Baratella et al. (2020a), we applied a new method that consists of using a combination of Fe and titanium (Ti) lines to derive $T_{\text{eff}}$ (through excitation equilibrium), and using only Ti lines to derive $\log g$ (through ionisation equilibrium) and $\xi$ (by zeroing the trend between individual abundances and strength, or EW, of the lines). The code MOOG (Sneden 1973) was used to analyse the HARPS-N co-added spectrum. We measured line EWs with the ARES v2 code (Sousa et al. 2015) and discarded lines with errors larger than 10% and with EW > 120 mÅ. The 1D LTE model atmospheres linearly interpolated from the ATLAS9 grid of Castelli & Kurucz (2003), with new opacities (ODFNEW), were adopted.

From the spectral type of the star (G2V), we estimated an input $T_{\text{eff}}$ of 5660 K from the relation by Pecaut & Mamajek (2013). Using different colour indexes (the reddening is negligible as the star is 22 pc distant from the Sun) and different calibrated relations (Casagrande et al. 2010, Mucciarelli et al. 2021), the photometric estimates of $T_{\text{eff}}$ range between 5589 K (from $B_p - R_p$), 5619 K (from $V - K$) and 5762 K (from $J - K$). The derived surface gravity based on the Gaia DR3 parallax (Gaia Collaboration 2022) is between 4.50 and 4.54 dex (error of the order of 0.05), depending on the photometric $T_{\text{eff}}$ considered. Finally, using the Dutra-Ferreira et al. (2016) relation, we estimated an initial $\xi$ of 0.095 ± 0.05 km s⁻¹.

The final values obtained with the new spectroscopic approach are: $T_{\text{eff}}=5700 \pm 75$ K, $\log g=4.54 \pm 0.05$ and $\xi=1.03 \pm 0.10$ km s⁻¹. The iron abundance is [Fe/H] I=0.02±0.06, [Fe/H] II=0.05±0.05, while the titanium abundance is [Ti/H] I=0.04±0.07, [Ti/H] II=0.05±0.04. The uncertainties reported are the quadratic sum of the scatter due to the EW measurements and the contribution of the parameter uncertainties to the final abundances. Our results are in perfect agreement with the initial guesses and with the recent analysis of Mann et al. (2020).

Fixing $T_{\text{eff}}$, $\log g$, $\xi$, and [Fe/H] I to the values found above, we measured the stellar projected rotational velocity ($v \sin i_*$) using the same MOOG code and applying the spectral synthesis of two regions around 6200 and 6700 Å. We adopted the same grid of model atmosphere and, after fixing the microturbulence velocity to the value of 3.2 km s⁻¹ from the relationship by Brewer et al. (2016), we find a $v \sin i_*$ of 7.2 ± 0.7 km s⁻¹, consistent with the result by Mann et al. (2020).

Finally, we also derived the lithium abundance $\log A(\text{Li})_{\text{NLTE}}$ from the measured lithium EW (=84.5±2.0 mÅ) and considering our stellar parameters previously derived together with the NLTE corrections by Lind et al. (2009). The values of the lithium abundance is 2.56 ± 0.06 and the position of the target in a $\log A(\text{Li})_{\text{NLTE}}-T_{\text{eff}}$ diagram is compatible with the UMa group, as expected (see Mann et al. 2020).

We calculated the stellar radius and mass with the EXOFASTv2 tool (Eastman et al. 2019), by fitting the stellar Spectral Energy Distribution (SED) and providing the stellar luminosity derived from the SED as input to the MIST stellar evolutionary

<figure><img src="image_3.png" /><figcaption>Fig. 2. Spectral energy distribution of the host star HD 63433 with the best-fit model overplotted (solid line). Red and blue points correspond to the observed and predicted values, respectively.</figcaption></figure>

tracks (Dotter 2016). For the SED we considered the *Tycho B* and *V* magnitudes (Høg et al. 2000), the 2MASS near-IR *J*, *H* and *K* magnitudes (Cutri et al. 2003), and the WISE mid-IR W1, W2, W3 and W4 magnitudes (Cutri et al. 2021). We imposed Gaussian priors on (i) the stellar effective temperature and metallicity from our analysis of the HARPS-N spectra; (ii) the parallax 44.6848 ± 0.0228 mas from the *Gaia* EDR3 (Gaia Collaboration 2016, 2021) and (iii) the stellar age 413 ± 23 Myr from the most updated value for the Ursa Major association (Jones et al. 2015). We used uninformative priors for all the other parameters. The best fit of the SED is shown in Fig. 2. We found $R_\star = 0.897 \pm 0.019 R_\odot$ and $M_\star = 0.994 \pm 0.028 M_\odot$, in excellent agreement with the previous determination by Mann et al. (2020).

From modelling the activity of the TESS light curve (Fig. 1) with a GP regression, we obtained a rotation period P = 6.48 ± 0.08 days. This result is in agreement with the measurement by Mann et al. (2020), and it is confirmed by looking at the periodogram of the V-band STELLA light curve calculated with the generalised Lomb-Scargle (GLS, Zechmeister & Kürster 2009) tool (Fig. C.1), which shows the main and significant peak at 6.4 days. The fundamental stellar properties of HD 63433 are summarised in Table 1.

## 4. Frequency content analysis of radial velocities and activity diagnostics

We analysed the frequency content of RVs and activity indicators to investigate the presence of sinusoidal periodic signals. The time series of the RVs (both DRS and TERRA), and the corresponding GLS periodograms are shown in the panels of the first two columns of Fig. 3. The main peak is statistically significant (bootstrap false alarm probability FAP <1%) and located at the second harmonic of the stellar rotation period. The periodograms of the residuals, after removing the best-fit sinusoid calculated by GLS (third column of Fig. 3), show the main peak at the first harmonic of the rotation period, with low significance in the case of TERRA RVs. No other significant peaks are detected, especially at the orbital frequencies of the planets, indicating that the RVs scatter is dominated by stellar activity.

The time series and GLS periodograms of the spectroscopic activity diagnostics are shown in Fig. 4. The log $R'_{\text{HK}}$ index and BIS contain a quite significant signal related to the stellar rotation period. The main peaks in the periodograms occur at a frequency corresponding to $P_{\text{rot},*}$ and to the second harmonic of
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Time series of the RVs calculated from HARPS-N spectra with the DRS and TERRA pipelines, and their GLS periodograms. First column: RV timeseries. The mean value has been subtracted from the original DRS data. Second column: The corresponding GLS periodograms of the two RV dataset. FAP levels are of 1% and 10% are indicated by yellow and red horizontal lines, respectively. They are calculated through a bootstrap analysis. Third column: GLS periodograms of the pre-whitened RV data.</figcaption></figure>

<table><caption>**Table 1.** Fundamental parameters of HD 63433 (TOI-1726).</caption><thead><tr><th>Parameter</th><th>Value</th><th>Ref.</th></tr></thead><tbody><tr><td><i>B</i><sub>T</sub> (mag)</td><td>7.749 ± 0.016</td><td>[1]</td></tr><tr><td><i>V</i><sub>T</sub> (mag)</td><td>6.987 ± 0.010</td><td>[1]</td></tr><tr><td><i>J</i> (mag)</td><td>5.624 ± 0.043</td><td>[2]</td></tr><tr><td><i>H</i> (mag)</td><td>5.359 ± 0.026</td><td>[2]</td></tr><tr><td><i>K</i><sub>S</sub> (mag)</td><td>5.258 ± 0.016</td><td>[2]</td></tr><tr><td>WISE1 (mag)</td><td>5.246 ± 0.178</td><td>[3]</td></tr><tr><td>WISE2 (mag)</td><td>5.129 ± 0.087</td><td>[3]</td></tr><tr><td>WISE3 (mag)</td><td>5.297 ± 0.016</td><td>[3]</td></tr><tr><td>WISE4 (mag)</td><td>5.163 ± 0.031</td><td>[3]</td></tr><tr><td>parallax, ϖ (mas)</td><td>44.685 ± 0.023</td><td>[4]</td></tr><tr><td><i>T</i><sub>eff</sub> (K)</td><td>5700 ± 75</td><td>[5]</td></tr><tr><td>log g (dex)</td><td>4.54 ± 0.05</td><td>[5]</td></tr><tr><td><i>v</i><sub>t</sub> (km s<sup>-1</sup>)</td><td>1.03 ± 0.10</td><td>[5]</td></tr><tr><td>[Fe/H] I (dex)</td><td>0.02 ± 0.06</td><td>[5]</td></tr><tr><td>A(Li)</td><td>2.56 ± 0.06</td><td>[5]</td></tr><tr><td><i>v</i> sin i<sub>★</sub> (km s<sup>-1</sup>)</td><td>7.2 ± 0.7</td><td>[5]</td></tr><tr><td><i>P</i><sub>rot,★</sub> (days)</td><td>6.48 ± 0.08</td><td>[5]</td></tr><tr><td>Mass, <i>M</i><sub>★</sub> (<i>M</i><sub>☉</sub>)</td><td>0.994<sup>+0.027</sup><sub>-0.028</sub></td><td>[5]</td></tr><tr><td>Radius, <i>R</i><sub>★</sub> (<i>R</i><sub>☉</sub>)</td><td>0.897 ± 0.019</td><td>[5]</td></tr><tr><td>Density, ρ<sub>★</sub> (ρ<sub>☉</sub>)</td><td>1.376 ± 0.078</td><td>[5]</td></tr><tr><td><i>L</i><sub>★</sub> (<i>L</i><sub>☉</sub>)</td><td>0.765<sup>+0.053</sup><sub>-0.050</sub></td><td>[5]</td></tr><tr><td><i>L</i><sub>★,X</sub> (erg s<sup>-1</sup>)</td><td>7.5 × 10<sup>28</sup></td><td>[6]</td></tr><tr><td>Age (Myr)</td><td>414 ± 23</td><td>[7]</td></tr></tbody></table>

**References.** [1] Tycho catalogue, Høg et al. (2000); [2] 2MASS, Cutri et al. (2003); [3] AllWISE, Cutri et al. (2021); [4] Gaia eDR3, Gaia Collaboration (2016, 2021, 2022); [5] This work; [6] Zhang et al. (2022); [7] Jones et al. (2015).

$P_{\text{rot},★}$, respectively. A non-significant (FAP = 41%) peak at a frequency compatible with the rotational frequency is also detected in the residuals of the H-alpha index, after removing a long-term curvature clearly visible in the time series, for which we cannot assess a periodicity, if actually present. The periodogram of the FWHM time series has the main peak located at ~162 days (FAP = 0.3%), whose origin is not clear. After a pre-whitening, the periodogram of the FWHM residuals shows a peak at 2.15 days with FAP = 51%. Given the very low statistical significance of this signal, which corresponds to the second harmonic of $P_{\text{rot},★}$, in the following analysis we do not use the FWHM to correct the RVs for the stellar activity term.

We note that the RV and BIS data are significantly correlated. The Spearman's rank correlation coefficient is $ρ = -0.88$, calculated using the `r_correlate` function in `IDL`.

## 5. Joint analysis of transits and RVs

We modelled the RV time series jointly with the detrended and flattened TESS light curve, which includes only the part of the data with the transits, and a sufficiently long out-of-transit baseline for a proper modelling of the planet ingress and egress. A key point in our analysis is the modelling of the stellar activity term which is the dominant signal in the RV time series. We took special care of it by testing different models, mostly based on the GP regression analysis. A schematic description of the test models is provided in Table 2, with more details given in Appendix A.

For all the tested models, we explored the full parameter space using the publicly available Monte Carlo (MC) nested sampler and Bayesian inference tool MultiNest v3.10 (e.g. Feroz et al. 2019), through the `pyMultiNest` wrapper (Buchner et al. 2014). Our MC set-up included 300 live points, a sampling
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Time series (**left column**) and periodograms (**right column**) of spectroscopic activity diagnostics. The periodogram of the H-alpha index refers to pre-whitened data, after removing the long-term curvature visible in the time series. FAPs are calculated through a bootstrap analysis.</figcaption></figure>

efficiency of 0.5, and a Bayesian tolerance of 0.3. We used the code `batman` (Kreidberg 2015) for modelling the photometric transits. In a few cases, when the same dataset has been used, we performed a Bayesian model selection by comparing the values of the Bayesian evidence $ln Z$ calculated by MultiNest, and using the empirical scale indicated in Table 1 of Feroz et al. (2011) to assess their relative statistical significance<fnref n="2" />. We assigned the same a-priori probability to each model.

Concerning the light curve, we modelled the limb darkening with a quadratic law, and fitted the coefficients $LD_{c1}$ and $LD_{c2}$ using the formalism and uniform priors given by Kipping (2013, see Eqs. (15) and (16) therein). We also introduced constant jitters $\sigma_{jit,TESS}$ added in quadrature to the nominal photometric uncertainties, one jitter term for Sector 20, and one for Sectors 44–47, to take into account the change in the TESS performance over nearly 2 yr. The complete list of priors used for all the test models is provided in Table A.1.

The main goals of our analysis are the improvement of the orbital parameters for HD 63433 b and HD 63433 c, and the measure of the fundamental planetary parameters radius, mass, and average density. Based on the statistical results and mass measurements summarised in Table 2, we notice that:

- the GP quasi-periodic and quasi-periodic with cosine (QPC) kernels perform equally well (models M1 and M8), and better than the double simple harmonic oscillator (dSHO) kernel (model M7), when modelling only the RV time series. It

is not statistically advantageous treating each observing season separately (model M2). The stellar rotation period is recovered with high precision (GP hyper-parameter $\theta = 6.390^{+0.007}_{-0.005}$ days for M1), even with an adopted uniform and large prior $U(0,10)$ days;

- results for all the test models show that the planetary orbits can be assumed circular. The data do not allow for significant non-zero eccentricities (see Table 3);

- we cannot constrain the mass $m_b$ of HD 63433 b, but only derive upper limits, with values that change depending on the model. We notice that the orbital period $P_b$ and the stellar rotation period $P_{rot,*}$ are not dissimilar, and this could make more difficult the identification of a planetary Doppler signal with a small semi-amplitude, despite the very precise transit ephemeris recovered;

- most of the models allow us to retrieve the mass $m_c$ of HD 63433 c with a statistical significance (defined as the ratio $m_c/\sigma_{m_c}^-$) in the range 2.1–2.7$\sigma$ (with the exception of models M6, M7, and M9, for which the significance is lower);

- using the more complex multi-dimensional GP (quasi-periodic kernel), we find lower upper limits for the masses, with no improvement in their precision.

Despite the fact that we followed several pathways to model the stellar activity in the RV dataset, the data do not allow for a significant and accurate measurement of the mass $m_c$, but we can conclude with high confidence that the mass of HD 63433 c is less than twice the mass of Neptune. Nonetheless, we note that in a few cases the median values of Gaussian-like posterior distributions are very similar, and suggest that $m_c$ could be actually consistent with the mass of Neptune. Examples of posteriors for $m_b$ and $m_c$ are shown in Fig. 5, with reference to Table 2. We also note that, in a few cases, masses of young planets have been

<footnote n="2">² According to that scale for interpreting model probabilities, ($\log Z_1 - \log Z_2$) < 1 indicates inconclusive evidence in favour of model 1 over model 2; if ($\log Z_1 - \log Z_2$) ~ 1, there is weak evidence in favour of model 1; if ($\log Z_1 - \log Z_2$) ~ 2.5, there is moderate evidence in favour of model 1; if ($\log Z_1 - \log Z_2$) ≥ 5, model 1 is strongly favoured over model 2.</footnote>
<table><caption>**Table 2.** Summary of the different RV (DRS) + transit photometry + activity diagnostics (when specified) models tested in this work.</caption><thead><tr><th>Model ID</th><th>Description</th><th>Bayesian evidence <sup>(a)</sup><br>(ln Z)</th><th>Planetary mass <sup>(b)</sup><br>M<sub>⊕</sub></th><th>Planetary density<br>g cm<sup>-3</sup></th></tr></thead><tbody><tr><td>M1</td><td>GP quasi-periodic kernel applied to RVs</td><td>103927.9 ± 0.6</td><td>m<sub>b</sub> = 4.1<sup>+3.7</sup><sub>-2.7</sub> (10.6)<br>m<sub>c</sub> = 16.6<sup>+8.1</sup><sub>-7.9</sub> (29.4)</td><td>ρ<sub>b</sub> = 2.7<sup>+2.5</sup><sub>-1.9</sub> (7.3)<br>ρ<sub>c</sub> = 6.0<sup>+3.1</sup><sub>-2.8</sub> (11.0)</td></tr><tr><td>M2</td><td>GP quasi-periodic kernel applied to RVs, with a different set<br>of hyper-parameters (h,w,λ) for each observing season</td><td>103925.5 ± 0.6</td><td>m<sub>b</sub> = 3.5<sup>+3.7</sup><sub>-2.5</sub> (10.5)<br>m<sub>c</sub> = 20.5<sup>+8.7</sup><sub>-9.4</sub> (34.1)</td><td>ρ<sub>b</sub> = 2.3<sup>+2.6</sup><sub>-1.7</sub> (7.0)<br>ρ<sub>c</sub> = 7.7<sup>+3.3</sup><sub>-3.6</sub> (13.2)</td></tr><tr><td>M3</td><td>GP quasi-periodic kernel trained <sup>(c)</sup> on BIS</td><td>103417.5 ± 0.6</td><td>m<sub>b</sub> = 2.9<sup>+3.5</sup><sub>-1.7</sub> (9.4)<br>m<sub>c</sub> = 19.7<sup>+6.8</sup><sub>-7.3</sub> (31.3)</td><td>ρ<sub>b</sub> = 1.9<sup>+2.3</sup><sub>-1.3</sub> (6.2)<br>ρ<sub>c</sub> = 7.4 ± 2.8 (12.0)</td></tr><tr><td>M4</td><td>GP quasi-periodic kernel trained <sup>(c)</sup> on log R'<sub>HK</sub></td><td>104237.8 ± 0.6</td><td>m<sub>b</sub> = 2.4<sup>+3.0</sup><sub>-1.8</sub> (8.1)<br>m<sub>c</sub> = 18.9 ± 7.0 (30.4)</td><td>ρ<sub>b</sub> = 1.6<sup>+1.9</sup><sub>-1.1</sub> (5.4)<br>ρ<sub>c</sub> = 7.1<sup>+2.8</sup><sub>-2.6</sub> (11.7)</td></tr><tr><td>M5</td><td>Multi-dimensional GP framework <sup>(d)</sup><br>using log R'<sub>HK</sub> and BIS, and the quasi-periodic kernel; circular orbits</td><td>836.6 ± 0.3</td><td>m<sub>b</sub> = 1.2<sup>+1.4</sup><sub>-0.9</sub> (3.9)<br>m<sub>c</sub> = 7.0<sup>+3.2</sup><sub>-3.0</sub> (12.1)</td><td>ρ<sub>b</sub> = 0.8<sup>+1.0</sup><sub>-0.9</sub> (2.6)<br>ρ<sub>c</sub> = 2.6<sup>+0.9</sup><sub>-1.1</sub> (4.6)</td></tr><tr><td>M6</td><td>Multi-dimensional GP framework <sup>(d)</sup><br>using log R'<sub>HK</sub> and the quasi-periodic kernel; circular orbits</td><td>565.6 ± 0.5</td><td>m<sub>b</sub> = 1.6<sup>+2.2</sup><sub>-1.2</sub> (5.6)<br>m<sub>c</sub> = 8.7<sup>+4.6</sup><sub>-1.6</sub> (17.1)</td><td>ρ<sub>b</sub> = 1.1<sup>+1.5</sup><sub>-0.8</sub> (3.8)<br>ρ<sub>c</sub> = 3.3<sup>+2.0</sup><sub>-1.7</sub> (6.6)</td></tr><tr><td>M7</td><td>GP dSHO kernel applied to RVs</td><td>103923.6 ± 0.5</td><td>m<sub>b</sub> = 2.9<sup>+3.0</sup><sub>-2.0</sub> (8.4)<br>m<sub>c</sub> = 11.8<sup>+7.7</sup><sub>-6.8</sub> (25.0)</td><td>ρ<sub>b</sub> = 1.9<sup>+2.1</sup><sub>-1.3</sub> (5.8)<br>ρ<sub>c</sub> = 4.3<sup>+2.9</sup><sub>-2.5</sub> (9.1)</td></tr><tr><td>M8</td><td>GP QPC kernel applied to RVs</td><td>103928.1 ± 0.6</td><td>m<sub>b</sub> = 4.3<sup>+4.2</sup><sub>-2.0</sub> (11.4)<br>m<sub>c</sub> = 16.6<sup>+8.7</sup><sub>-8.3</sub> (30.5)</td><td>ρ<sub>b</sub> = 2.8<sup>+2.8</sup><sub>-1.9</sub> (7.5)<br>ρ<sub>c</sub> = 6.2<sup>+3.3</sup><sub>-3.1</sub> (11.7)</td></tr><tr><td>M9</td><td>Activity corrected through a linear regression fit between RVs and BIS</td><td>103991.2 ± 0.5</td><td>m<sub>b</sub> = 2.1<sup>+2.1</sup><sub>-1.4</sub> (6.2)<br>m<sub>c</sub> = 8.4 ± 4.9 (16.7)</td><td>ρ<sub>b</sub> = 1.4<sup>+1.4</sup><sub>-1.0</sub> (4.2)<br>ρ<sub>c</sub> = 3.1 ± 1.8 (6.3)</td></tr></tbody></table>

**Notes.** All include two Keplerians for planets b and c, except for models M5 and M6, for which we adopted circular orbits. <sup>(a)</sup>We note that the analysed datasets related to each model are not always the same, therefore model comparison based on the values of the Bayesian evidence ln Z can be performed only when the same datasets are involved. <sup>(b)</sup>The 95th percentile is provided in parenthesis. <sup>(c)</sup>With “trained” we mean that the same GP kernel is used to model the RVs and the activity diagnostic. In the case of a quasi-periodic kernel, three hyper-parameters over four, namely theta, w, and λ, are shared by the two datasets, while two distinct amplitudes h are used. <sup>(d)</sup>To speed up the analysis, we did not fit the photometric transits, and we assumed circular orbits (the results from other models show that the eccentricities are consistent with zero). The priors for the planet ephemeris are defined based on the results of the M4 model. We adopted planet radii calculated from M4 to derive the densities.

claimed with a precision of the order of 3–3.5σ. For example, Desidera et al. (2023) found that TOI-179 b (further discussed in Sect. 6) has a mass of $24.1_{-7.7}^{+7.1} M_⊕$ (3.1σ precision). The host TOI-179 has a similar mass, and an activity-induced RV scatter very similar to that measured for HD 63433 (~24 m s⁻¹), while TOI-179 b has an orbital period nearly five times shorter than HD 63433 c, making a more precise RV detection of the planet potentially easier. K2-100 b (Barragán et al. 2019b), is a 750 Myr-old planet with an orbital period of 1.7-d, and a measured mass of 21.8±6.2 $M_⊕$ (3.5σ precision). Similarly, mass measurements have been claimed for TOI-560 b (480 Myr; $m = 10.2_{-3.1}^{+3.4} M_⊕$, i.e. 3.3σ precision; Barragán et al. 2022b), and for Kepler-411 d (200 Myr; $m = 15.2 \pm 5.1$, i.e. 3σ precision, measured through transit timing variations; Sun et al. 2019). Since we determined the mass of HD 63433 c with a significance of 2.7σ at best, we do not consider this result significant enough for claiming a mass measurement, but indeed it deserves further attention to be confirmed with additional data.

Concerning the transit ephemeris and other specific parameters of the transit model, we did not find significant differences among the results obtained for all the tested models (except for the models M5 and M6, where the light curve is not fitted). We provide in Table 3 the results of model M4, selected among the solutions available, which are all equivalent concerning the transit ephemeris. We underline that the election of one of the test models to our reference model is a tricky issue in this case, because there is not a representation of the activity term in the RVs which we can significantly define the best. Therefore, the choice of model M4 should be considered mainly for illustrative purposes. We found that for model M4 the fitted RV jitter $\sigma_{jitt, HARPS-N}$ is nearly half that of model M1, for

which we got $\sigma_{jitt, HARPS-N} = 12.1 \pm 1.8 \text{ m s}^{-1}$, suggesting that, at some level, the activity modelling benefits from using the time series of the log R'HK index to constrain the GP hyper-parameters. Nonetheless, we note that even the RV uncorrelated jitter is lower for model M5 ($\sigma_{jitt, HARPS-N} = 3.8_{-2.1}^{+1.8} \text{ m s}^{-1}$), but the masses of planet c are different: for model M4, it is similar to that of Neptune at a ~2.7σ level, while for model M5 we find 7 $M_⊕$ with ~2σ significance. We show the best-fit solution (model M4) for transits and spectroscopic orbits in Fig. 6.

With more data from new TESS Sectors, we improved the precision of the transit ephemeris and other transit-related parameters, like the planetary radii, with respect to the results of Mann et al. (2020), confirming that the planets have co-planar orbits. Nevertheless, some caution is required about the precision of our determined planetary radii, which, for each planet, is calculated by fitting a light curve obtained by combining all the available transits together, and it is likely optimistic. The unspotted level of HD 63433 is unknown and this introduces a systematic error on the planetary radii. Considering that the flux in the TESS passband is modulated at the level of about 2%, we estimate a systematic error at the level of 1% on the planetary radii from transits close to light maxima and minima, respectively. This comes from the fact that the square of the ratio between the planetary and the stellar radii is proportional to the relative flux deficit at the centre of a transit. Therefore, a change in the out-of-transit flux by about 2% between light maxima and minima will affect the planetary radius derived from the corresponding transits by about 1%. This error estimate is indeed a lower limit because starspots uniformly distributed in longitude do not contribute to the amplitude of the photometric modulation, but affect the measurement of the radius ratio.
<table><caption>**Table 3.** Best-fit values of the free parameters of model M4 (RVs calculated by the DRS).</caption><thead><tr><th>Parameter</th><th>Best-fit value <sup>(a)</sup></th></tr></thead><tbody><tr><td>Fitted</td><td></td></tr><tr><td>RV stellar activity GP term:</td><td></td></tr><tr><td><i>h</i> (m s<sup>-1</sup>)</td><td>26.0<sup>+3.2</sup><sub>-3.0</sub></td></tr><tr><td><i>λ</i> (days)</td><td>45.7<sup>+14.0</sup><sub>-9.9</sub></td></tr><tr><td><i>w</i></td><td>0.23 ± 0.02</td></tr><tr><td><i>θ</i> (days)</td><td>6.39 ± 0.01</td></tr><tr><td>Planet-related parameters:</td><td></td></tr><tr><td><i>K</i><sub>b</sub> (m s<sup>-1</sup>)</td><td>0.8<sup>+1.0</sup><sub>-0.6</sub> (2.7)</td></tr><tr><td>Orbital period, <i>P</i><sub>b</sub> (days)</td><td>7.10794 ± 0.000009</td></tr><tr><td><i>T</i><sub>conj, b</sub> (BJD-2450000)</td><td>9584.5991 ± 0.0004</td></tr><tr><td>√<i>e</i><sub>b</sub> cos <i>ω</i><sub>★, b</sub></td><td>-0.029<sup>+0.335</sup><sub>-0.303</sub></td></tr><tr><td>√<i>e</i><sub>b</sub> sin <i>ω</i><sub>★, b</sub></td><td>-0.133<sup>+0.092</sup><sub>-0.138</sub></td></tr><tr><td><i>K</i><sub>c</sub> (m s<sup>-1</sup>)</td><td>4.5 ± 1.7 (7.3)</td></tr><tr><td>Orbital period, <i>P</i><sub>c</sub> (days)</td><td>20.54379 ± 0.00002</td></tr><tr><td><i>T</i><sub>conj, c</sub> (BJD-2450000)</td><td>9583.6367 ± 0.0004</td></tr><tr><td>√<i>e</i><sub>c</sub> cos <i>ω</i><sub>★, c</sub></td><td>-0.277<sup>+0.373</sup><sub>-0.255</sub></td></tr><tr><td>√<i>e</i><sub>c</sub> sin <i>ω</i><sub>★, c</sub></td><td>-0.060<sup>+0.128</sup><sub>-0.140</sub></td></tr><tr><td>Acceleration, <i>γ</i> (m s<sup>-1</sup> day<sup>-1</sup>)</td><td>0.004 ± 0.017</td></tr><tr><td><i>R</i><sub>b</sub>/<i>R</i><sub>★</sub></td><td>0.02063<sup>+0.00052</sup><sub>-0.00028</sub></td></tr><tr><td>Inclination, <i>i</i><sub>b</sub> (deg)</td><td>88.96<sup>+0.70</sup><sub>-0.82</sub></td></tr><tr><td><i>R</i><sub>c</sub>/<i>R</i><sub>★</sub></td><td>0.0249<sup>+0.0004</sup><sub>-0.0003</sub></td></tr><tr><td>Inclination, <i>i</i><sub>c</sub> (deg)</td><td>89.41<sup>+0.33</sup><sub>-0.26</sub></td></tr><tr><td>RV-related parameters:</td><td></td></tr><tr><td><i>σ</i><sub>jit, HARPS-N</sub> (m s<sup>-1</sup>)</td><td>6.6<sup>+2.9</sup><sub>-4.5</sub></td></tr><tr><td><i>γ</i><sub>HARPS-N</sub> (m s<sup>-1</sup>)</td><td>-15806.1<sup>+4.5</sup><sub>-4.6</sub></td></tr><tr><td>Light curve-related parameters:</td><td></td></tr><tr><td><i>σ</i><sub>jit, TESS</sub> sect. 20</td><td>0.000341 ± 0.000004</td></tr><tr><td><i>σ</i><sub>jit, TESS</sub> sect. 44–47</td><td>0.000262 ± 0.000002</td></tr><tr><td><i>LD</i><sub>c1</sub></td><td>0.50<sup>+0.13</sup><sub>-0.14</sub></td></tr><tr><td><i>LD</i><sub>c2</sub></td><td>-0.05<sup>+0.22</sup><sub>-0.17</sub></td></tr><tr><td>Derived</td><td></td></tr><tr><td>Eccentricity, <i>e</i><sub>b</sub></td><td>0.09<sup>+0.14</sup><sub>-0.07</sub> (0.37)</td></tr><tr><td><i>a</i><sub>b</sub>/<i>R</i><sub>★</sub><sup>(b)</sup></td><td>16.8<sup>+0.4</sup><sub>-1.0</sub></td></tr><tr><td><i>a</i><sub>b</sub> (au)</td><td>0.0722 ± 0.0007</td></tr><tr><td>Impact param., <i>b</i><sub>b</sub></td><td>0.32<sup>+0.28</sup><sub>-0.21</sub></td></tr><tr><td>Transit duration, <i>T</i><sub>b, 1,4</sub> (days)</td><td>0.129<sup>+0.005</sup><sub>-0.011</sub></td></tr><tr><td>Radius <sup>(c)</sup>, <i>r</i><sub>b</sub> [<i>R</i><sub>⊕</sub>]</td><td>2.02<sup>+0.06</sup><sub>-0.05</sub></td></tr><tr><td>Mass, <i>m</i><sub>b</sub> (<i>M</i><sub>⊕</sub>)</td><td>2.4<sup>+3.0</sup><sub>-1.8</sub> (8.1)</td></tr><tr><td>Mean density, <i>ρ</i><sub>b</sub> (g cm<sup>-3</sup>)</td><td>1.6<sup>+2.0</sup><sub>-1.2</sub> (5.4)</td></tr><tr><td><i>T</i><sub>eq, b</sub> (K)</td><td>969 ± 17</td></tr><tr><td>Insolation, <i>S</i><sub>b</sub> (<i>S</i><sub>⊕</sub>)</td><td>146 ± 10</td></tr><tr><td>Eccentricity, <i>e</i><sub>c</sub></td><td>0.11<sup>+0.18</sup><sub>-0.07</sub> (0.43)</td></tr><tr><td><i>a</i><sub>c</sub>/<i>R</i><sub>★</sub><sup>(b)</sup></td><td>36.4<sup>+1.0</sup><sub>-1.5</sub></td></tr><tr><td><i>a</i><sub>c</sub> (au)</td><td>0.147 ± 0.001</td></tr></tbody></table>

**Table 3.** continued.

<table><tr><td>Impact param., <i>b</i><sub>c</sub></td><td>0.34<sup>+0.17</sup><sub>-0.19</sub></td></tr><tr><td>Transit duration, <i>T</i><sub>c, 1,4</sub> (days)</td><td>0.173<sup>+0.006</sup><sub>-0.008</sub></td></tr><tr><td>Radius <sup>(c)</sup>, <i>r</i><sub>c</sub> (<i>M</i><sub>⊕</sub>)</td><td>2.44 ± 0.07</td></tr><tr><td>Mass, <i>m</i><sub>c</sub> (<i>M</i><sub>⊕</sub>)</td><td>18.9<sup>+7.0</sup><sub>-6.9</sub> (30.3)</td></tr><tr><td>Mean density, <i>ρ</i><sub>c</sub> (g cm<sup>-3</sup>)</td><td>7.1<sup>+2.8</sup><sub>-2.6</sub> (11.7)</td></tr><tr><td><i>T</i><sub>eq, c</sub> (K)</td><td>680 ± 12</td></tr><tr><td>Insolation, <i>S</i><sub>c</sub> (<i>S</i><sub>⊕</sub>)</td><td>36 ± 3</td></tr></table>

**Notes.** <sup>(a)</sup>The uncertainties are given as the 16th and 84th percentiles of the posterior distributions. For some of the parameters, we provide the 95th percentile in parenthesis. <sup>(b)</sup>In place of *a/R*<sub>★</sub>, we used the stellar density *ρ*<sub>★</sub> as a free parameter with Gaussian prior *N*(1.376,0.078) *ρ*<sub>⊕</sub>, from which we derived the *a*<sub>p</sub>/<i>R</i><sub>★</sub> ratios at each step of the MC sampling. <sup>(c)</sup>A systematic error of at least 1% of the best-fit radius should be added, to take into account effects due to stellar activity (Sect. 5).

An additional systematic effect is produced by the occultations of starspots during transits, but its amplitude is probably lower by about 50% than the effect produced by the changes in the out-of-transit light reference level, if we assume that spots in this G dwarf have a temperature deficit of about 1000 K with respect to the unperturbed photosphere and a filling factor of about 2% in the latitude bands occulted by the planets (e.g. Ballerini et al. 2012). We did not find a significant RV acceleration ($\dot{\gamma}$ = 0.004 ± 0.017 m s<sup>-1</sup>d<sup>-1</sup>).

Limits on outer Jovian planets around HD 63433 can be placed by the HIPPARCOS-Gaia proper motion anomaly technique (Brandt 2021; Kervella et al. 2022). For example, the sensitivity limits based on this approach would be compatible with the presence of a ≥1.0 $M_{Jup}$ companion in the “sweet spot” separation range 3–10 au (Kervella et al. 2022). Interestingly, we note that the reported S/N of the proper motion anomaly for HD 63433 is in the range 2–2.5 (with a growing trend from Gaia DR2 to DR3). While low, it is intriguing, as in this regime of orbital separations Jupiter- and super-Jupiter-mass companions might still be missed by existing RV datasets (see e.g. Fig. 9 of Hirsch et al. 2021), particularly if they were to lie on significantly non-coplanar orbits. Given that the star is nearby, there are good prospects for improved sensitivity to such companions based on Gaia DR4 results, which are expected to be released in late 2025.

We did not measure different nor more precise results for the planetary masses when using the RVs extracted with TERRA. As an example, for the case of models M1, M4, and M6 we got $m_b = 3.8_{-2.6}^{+3.9}$ and $m_c = 14.8_{-7.7}^{+8.0} M_⊕$, $m_b = 2.6_{-1.8}^{+3.2}$ and $m_c = 15.1_{-7.1}^{+7.5} M_⊕$, and $m_b = 1.8_{-1.3}^{+2.3}$ and $m_c = 5.9_{-3.8}^{+4.6} M_⊕$ respectively.

## 6. Considerations on planetary structure based on the mass-radius diagram

The more conservative conclusion that we can draw from the analysis described in Sect. 5 is that we determined upper limits for the planet masses, $m_b \le 11 M_⊕$ and $m_c \le 31 M_⊕$. The mass of HD 63433 b remains undetected (and we cannot exclude the possibility that this result is mainly influenced by the choice of a specific model). The mass of HD 63433 c remains (very) uncertain, but the clue that it might be Neptune-like would imply a bulk density that, if confirmed, would be interestingly high for a planet with the size of a mini-Neptune. For this reason, in this
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Posterior distributions for the planetary masses for some of the models tested in this work (Table 2). Vertical dashed lines indicate the 95% percentiles.</figcaption></figure>

Section we comment on this and other possibilities for the planets' physical structures which are compatible with our results, by inspecting the location of the two planets in a mass-radius diagram for exoplanets. Specifically, our considerations concern the solutions which correspond to models M4 and M5 of Table 2).

We show in the upper panel of Fig. 7 a mass-radius diagram which includes well-characterised planets (with masses and radii measured with a precision lower than 30% and 10%, respectively). Assuming the lower value for $m_c$ (model M5), HD 63433 c is located in a well populated region of the diagram, corresponding to water worlds, i.e. planets containing significant amounts of H$_2$O-dominated fluid/ice in addition to rock and gas. Instead, if the estimate of $m_c$ from model M4 is the one more representative of the real mass of HD 63433 c, the planet would occupy a position on the diagram which is less populated, but not empty, suggestive of a structure with a rocky core

and a water layer < 50% by mass. The lower panel of Fig. 7 shows a mass-radius diagram which includes only planets with age <900 Myr, with no selection based on the precision of their masses and radii. We note that, assuming the higher mass (density) estimate from model M4, HD 63433 c would have a similar structure of TOI-179 b (Desidera et al. 2023; Vines et al. 2023) and Kepler-411 b (Sun et al. 2019). TOI-179 is a K2V star with an age similar to that of the system HD 63433 (400±100 Myr; Desidera et al. 2023). TOI-179 b is a dense mini-Neptune that moves on an eccentric orbit ($e = 0.34^{+0.07}_{-0.09}$) with a shorter orbital semi-major axis ($a = 0.0481 \pm 0.0004$ au) and higher equilibrium temperature ($T_{eq} \sim 990$ K, in the case of a circular orbit with same semi-major axis $a$) than HD 63433 c. Desidera et al. (2023) show that the typical composition of TOI-179 b can be assumed 75% rock + 25% water, and conclude that the planet has likely lost most of its primordial atmosphere, and it is stable against hydrodynamic evaporation. Kepler-411 b is a member of a fairly compact four-planet system orbiting a ~200 Myr old K2V star with a period similar to that of TOI-179 b.

Based on our derived upper limits to the mass of HD 63433 b, we can conclude that its composition could be compatible with that of a water world, and we cannot exclude that it is similar to that of the outermost companion, especially if we consider the results of the M5 model. As for the case of TOI-179 b, noticing that the two planets have similar equilibrium temperatures, HD 63433 b could have lost any primordial gaseous H-He envelope (as we discuss in detail in the next Section), and reached its definitive position on the mass-radius diagram.

## 7. Planetary atmosphere mass-loss due to photoevaporation

We evaluated the mass-loss rate of the planetary atmosphere of planet b and c using the hydro-based approximation developed by Kubyshkina et al. (2018a,b), coupled with the planetary core-envelope model by Lopez & Fortney (2014) and the MESA Stellar Tracks (MIST; Choi et al. 2016). For the stellar X-ray emission at different ages, we adopted the analytic description by Penz et al. (2008), anchored to the current value of the X-ray luminosity, $L_{\star,\text{X}} = 7.5 \times 10^{28}$ erg s$^{-1}$ in the band 5–100 Å, derived from XMM-Newton spectra (Zhang et al. 2022). The stellar EUV luminosity (100–920 Å) was computed by means of the new scaling law by Sanz-Forcada et al. (2022). More details on our modelling of atmospheric evaporation are provided in Appendix B. Following Maggio et al. (2022), we expect small but non negligible differences in the evaporation efficiency and timescales if a different X-ray to EUV scaling law is assumed (e.g. King & Wheatley 2021, Johnstone et al. 2021), but a detailed comparison of results is not warranted at this stage, given the uncertainties on the planetary masses.

We performed several simulations of the past and future evolution of the planetary atmospheres, assuming different possible values for the planetary masses at the current age. The grids of test masses are defined by taking into account the upper limits presented in Sect. 5. For planet b, we explored the mass range 2–12 $M_\oplus$, while for planet c we limited our analysis to masses in the range 5–15 $M_\oplus$, because we found that for $m_c > 15 M_\oplus$ the planet is stable against photoevaporation.

First, we considered the cases with the planetary parameters adopted by Zhang et al. (2022) in their 3D hydrodynamic modelling of the hydrogen escape and of the possible Ly$\alpha$ and He I 10833 Å absorption features. In particular, they assumed $m_b = 5.5 M_\oplus$ and $R_b = 2.15 R_\oplus$ for planet b, and $m_c = 7.3 M_\oplus$ and
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** *First two rows:* TESS light curve and best-fit transit models (**left panels**), and spectroscopic orbits (HARPS-N DRS RVs; **right panels**) related to planets HD 63433 b and HD 63433 c. Here, we show the solution obtained with model M4. The **upper right plot** shows the RV residuals, after removing the Doppler signal due to HD 63433 c and the stellar activity term from the original data, phased to the orbital period of 7.01794 d: it appears clear that the spectroscopic orbit corresponding to HD 63433 b remains not characterised. The phase-folded plot for HD 63433 c shows the RV residuals after removing the stellar activity term from the original data. *Third row:* Zoomed view of the stellar activity term in the HARPS-N RVs, after removing the Doppler signal due to HD 63433 c from the original RVs, as modelled by a GP quasi-periodic kernel. In model M5, the GP quasi-periodic kernel was trained on the time series of the log $R'_{HK}$ chromospheric activity diagnostic.</figcaption></figure>

$R_c = 2.67 R_⊕$ for planet c. Using these values, we predict a current atmospheric mass fraction of ~0.5% for planet b and ~1.4% for planet c, which are in good agreement with the corresponding values of 0.6% and 2% derived by Zhang et al. (2022). We also checked that our X-ray to EUV luminosity scaling provides consistent fluxes at 1 au distance, at the current age: we derived a value of 95 erg s⁻¹ cm⁻², to be compared with the estimate of 91 ± 27 erg s⁻¹ cm⁻² by Zhang et al. (2022), who employed a nominal stellar spectrum from 5 Å to 5 µm. Finally, we compared the predictions of the current mass-loss rates, which resulted in a factor of 1.7 higher for planet b and 2.1 lower for planet c, with respect to the values calculated by Zhang et al. (2022). This relatively small difference is not surprising, given the quite different modelling approaches. However, we stress that these planetary properties are just instantaneous snapshots, which do not take into account the past evolutionary history.

In the following, we report the results of our backward and forward in time simulations, which include the time evolution of the XUV irradiation and of the planetary structure in response to the stellar behaviour. We are interested in investigating how the planetary masses, radii, and atmospheric mass-loss rates change with time due to photoevaporation. The results of the simulations are shown in Figs. 8–9 and in Table 4.

**Planet b.** First, we computed atmospheric mass fractions and mass-loss rates at the present age for the grid of planet masses, and we simulated the evolution in the future. In Table 4 we report the e-folding mass-loss timescales. These evaporation

timescales are very short (1–100 Myr) in any case. Hence, we confirm the conclusion of Zhang et al. (2022), based on observations, that the inner planet has probably lost its atmosphere entirely. Then we explored the possible backward evolutionary pathways, in order to assess the initial planetary masses and radii at the age of 10 Myr (see Appendix B for details on the modelling). In practice, we found that it is always possible to find initial conditions such that the atmosphere is depleted within about 400 Myr. For each assumed planetary mass at the current age, we were able to determine what could be the highest initial mass at 10 Myr. Based on our results for the mass-loss timescales discussed above, we adopt the boundary condition that the atmospheric evaporation is completed at the current age of 414 Myr. These evolutionary histories are shown in Fig. 8. For example, if planet b has $m_b = 12 M_⊕$ at the present age, it could have started with a mass $m_{b,t=10 \text{Myr}} \sim 12.02 M_⊕$ and an atmospheric fraction $f_{\text{atm},t=10 \text{Myr}} = 0.14\%$ at 10 Myr.

For planetary masses $\le 7 M_⊕$ we found a remarkable change in evaporation efficiency. This is due to the relative change of the Jeans escape parameter, $\Lambda$, with respect to the control parameter $e^\Sigma$ in the hydro-based formulation by Kubyshkina et al. (2018b, see Appendix B for details). While for relatively high planetary masses the atmospheric escape is mainly driven by the photoevaporation, for lower planetary masses the mass-loss rate is increasingly larger because the atmospheric escape is more due to the thermal energy and the lower gravity of the planet. Hence, in order to reach the current age with an atmospheric mass
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Upper panel: Mass-radius diagram for exoplanets selected from the TEPCAT sample, available at https://www.astro.keele.ac.uk/jkt/tepcat/ (updated to 7 October 2022; Southworth 2011). Black dots represent planets with mass and radius measured with a relative precision lower than 30% and 10%, respectively. Planets of the HD 63433 system are indicated by reddish triangles (planet b), which denote mass upper limits, and squares (planet c), with reference to the corresponding models listed in Table 2. Theoretical curves for some planet compositions are overplotted, as calculated by Zeng et al. (2019) assuming 1 milli-bar surface pressure. Models for a planet with an H<sub>2</sub>O-gaseous atmosphere (50% and 100% H<sub>2</sub>O by mass, cyan and blue curves, respectively), and for a planet with different percentages of an H<sub>2</sub> gaseous envelope over a 50% water-rich layer (magenta curves) are calculated for an isothermal fluid/steam envelope equilibrium temperature of 700 K, similar to that of HD 63433 c (680±12 K). Lower panel: Mass-radius diagram for planets with age < 900 Myr. Triangles are used to identify planets for which only mass upper limits are available. No selection is made based on the precision of mass and radius measurements. We used Desidera et al. (2023) as a reference for the mass and radius of TOI-179 b.</figcaption></figure>

fraction $f_{\text{atm},t=414\text{Myr}} = 0$, the initial configuration of the planet is characterised by a massive and inflated atmospheric envelope. We found that if the current mass of planet b is $7 M_{\oplus}$, $m_{b,t=10\text{Myr}}$ could have been up to about $40 M_{\oplus}$, $f_{\text{atm},t=10\text{Myr}} = 82\%$, and the corresponding initial radius was $R_{b,t=10\text{Myr}} \sim 14.25 R_{\oplus}$. We remark that these are maximum mass values, such that complete evaporation occurs in about 400 Myr. Lower initial masses are also possible, with shorter evaporation timescales and identical atmosphere-depleted status at the present age.

We would like to remind readers that for the backward in time modelling we are considering as a boundary condition that the current mass-loss rate of planet b is zero. That explains the difference with the non-null mass-loss rates indicated in Table 4. In particular, for the cases of a planet with 2 and 4.5 $M_{\oplus}$, the boundary condition implies the runaway values reached in the late stages before the present age.

A further caveat about the simulations for planets with masses $\lesssim 7 M_{\oplus}$ is that a significant decrease in the planetary mass at early ages implies an outward migration of the planet, with a rate depending on the amount of angular momentum driven away from the system by the planetary wind flow (Fujita et al. 2022). An increase in the star-planet distance yields a decrease in XUV irradiation and equilibrium temperature, hence a decrease in mass-loss rate. Exploring this highly non-linear regime is beyond the scope of the present paper.

**Planet c.** In this case, the evolution is more simple, with no predicted change of regime (Fig. 9), due to the larger distance of the planet from the host star, a lower equilibrium temperature, and lower high-energy irradiation. In most of the cases that we have examined, i.e. for masses in the range 7.5–15 $M_{\oplus}$ at the current age, the planet mass and radius are scarcely affected by atmospheric evaporation, although the atmospheric mass
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Evolutionary history of planet HD 63433 b. (a) Mass-loss rate vs. time for different values of the planetary mass at the present age. (b) Fraction of planetary mass lost vs. time. Panels *c* and *d* show the evolution of the planetary mass and radius, respectively. Dashed lines indicate backward time evolution from the current age, solid lines forward time evolution. Note that all the simulations were constructed with the boundary condition that the planet completely lost the atmosphere at the present age of 414 Myr.</figcaption></figure>

fraction can change significantly during the system lifetime. The e-folding mass-loss timescales are always > 5 Gyr, even for the case of a planet with $M_p = 7.3 M_⊕$ considered by Zhang et al. (2022). The mass-loss timescale of 0.9 Gyr predicted by these authors was significantly lower since it was computed assuming a constant XUV irradiation and mass-loss rate. Instead, in our evolutionary models, the XUV luminosity drops to about 1/3 the current value already at 1.3 Gyr, and it keeps decreasing at later times. Only if a current mass of 5 $M_⊕$ is considered, the mass and radius change substantially starting from an age of 10 Myr.

## 8. Conclusions

In this paper, we analysed new photometric and spectroscopic data of the 400 Myr-old star HD 63433 (TOI-1726), which hosts

a multi-planet system. The presence of two planets in this relatively young system offers the opportunity to examine differences in planetary evolution pathways over the first hundreds of millions years after the system formation, under the influence of the same host star. We provided a characterisation of the system, with special focus on the measurement of the planet radii and masses, and on the investigation of the past and future planetary atmosphere evolution.

Thanks to new TESS photometry and follow-up with HARPS-N, we could improve the planetary ephemeris with respect to the study of Mann et al. (2020), with statistical evidence in favour of circular orbits for the two planets. We revised the measure and improve the precision of the planetary radii ($r_b = 2.02_{-0.05}^{+0.06} R_⊕$ and $r_c = 2.44 \pm 0.07 R_⊕$), and provide the first dynamic constraints on the planetary masses. After testing several approaches to mitigate the dominant stellar activity
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Evolutionary history of planet HD 63433 c. *Panel a:* mass-loss rate vs. time for different values of the planetary mass at the present age. *Panel b:* fraction of planetary mass lost vs. time. *Panels c and d:* time evolution of the planetary mass and radius, respectively. Dashed lines indicate backward time evolution from the current age, solid lines forward time evolution.</figcaption></figure>

contribution in the RV time series, we can only provide upper limits for the mass $m_b$ of the innermost planet b, conservatively concluding that $m_b \le 11 M_\oplus$ at 95% of confidence. As for the larger planet c, not all the models converge to similar mass estimates, but results from a few test models suggest a Neptune-like mass with a significance of 2.1–2.7$\sigma$, and that HD 63433 c could be a dense mini-Neptune, with a few known counterparts in the mass-radius diagram. The more conservative conclusion is that $m_c \le 31 M_\oplus$ at 95% of confidence.

Based on the constraints to the masses derived from our analysis, we explored the evolution of the atmospheres of both planets, taking into account the decay of stellar activity and XUV irradiation with time. We estimate that the current mass-loss timescale of planet b is sufficiently short to justify the assumption that its atmosphere is already evaporated. This finding is supported by the results of Zhang et al. (2022), who reported no evidence of Ly$\alpha$ or He I absorption due to photoevaporation in dedicated photometric observations with HST.

However, the past history of the planet is quite uncertain, because the evolution of the mass-loss rate is strongly dependent on the assumed mass and current atmospheric mass fraction. For masses $\gtrsim 7 M_\oplus$, the atmospheric mass fraction was already small at early ages (~10 Myr), hence the planetary mass and radius changed little in time. Instead, a planet with a lower mass at the present age could have started as an inflated body with a few tens of Earth masses, with a structure dominated by a heavy atmosphere. In any case, the current planetary structure is essentially depleted of a gaseous envelope.

For HD 63433 c, our models indicate that the atmospheric evaporation will keep going over the next ~4.5 Gyr, because the decay of the XUV irradiation determines a steady decrease of the mass-loss rate with time. Our prediction for the mass-loss timescale is different from the shorter value determined by Zhang et al. (2022), who did not take into account the evolution of the stellar activity. Another effect included in our modelling is the time evolution of the stellar bolometric luminosity and
<table><caption>**Table 4.** Results of the planetary mass-loss simulations from the current age to 5 Gyr (Sect. 7).</caption><thead><tr><th>Current<br>planet mass<br>(M<sub>⊕</sub>)</th><th>Mass-loss rate<br>at 414 Myr<br>(g s<sup>-1</sup>)</th><th>Core mass<br>(M<sub>⊕</sub>)</th><th>Core radius<br>(R<sub>⊕</sub>)</th><th>Mass-loss<br>timescale<br>(Myr)</th></tr></thead><tbody><tr><td colspan="5"><strong>Planet b</strong></td></tr><tr><td>12</td><td>3.7 &times; 10<sup>9</sup></td><td>11.99</td><td>1.93</td><td>50</td></tr><tr><td>10.5</td><td>6.1 &times; 10<sup>9</sup></td><td>10.46</td><td>1.86</td><td>103</td></tr><tr><td>9</td><td>1.1 &times; 10<sup>10</sup></td><td>8.99</td><td>1.79</td><td>113</td></tr><tr><td>7</td><td>2.6 &times; 10<sup>10</sup></td><td>6.99</td><td>1.67</td><td>76</td></tr><tr><td>5.5 <sup>(a)</sup></td><td>1.1 &times; 10<sup>11</sup></td><td>5.47</td><td>1.56</td><td>43</td></tr><tr><td>4.5</td><td>1.3 &times; 10<sup>11</sup></td><td>4.48</td><td>1.49</td><td>21</td></tr><tr><td>2</td><td>2.5 &times; 10<sup>12</sup></td><td>1.99</td><td>1.19</td><td>1</td></tr><tr><td colspan="5"><strong>Planet c</strong></td></tr><tr><td>15</td><td>0.0</td><td>14.93</td><td>2.05</td><td>—</td></tr><tr><td>13</td><td>1.1 &times; 10<sup>9</sup></td><td>12.93</td><td>1.97</td><td>&gt; 5 &times; 10<sup>3</sup></td></tr><tr><td>11</td><td>1.8 &times; 10<sup>9</sup></td><td>10.93</td><td>1.89</td><td>&gt; 5 &times; 10<sup>3</sup></td></tr><tr><td>10</td><td>2.4 &times; 10<sup>9</sup></td><td>9.93</td><td>1.84</td><td>&gt; 5 &times; 10<sup>3</sup></td></tr><tr><td>7.5</td><td>5.8 &times; 10<sup>9</sup></td><td>7.43</td><td>1.70</td><td>&gt; 5 &times; 10<sup>3</sup></td></tr><tr><td>7.3 <sup>(a)</sup></td><td>1.0 &times; 10<sup>10</sup></td><td>7.2</td><td>1.69</td><td>&gt; 5 &times; 10<sup>3</sup></td></tr><tr><td>5</td><td>1.9 &times; 10<sup>10</sup></td><td>4.94</td><td>1.53</td><td>1.6 &times; 10<sup>3</sup></td></tr></tbody></table>

**Notes.** We assumed $r_b = 2.02$ and $r_c = 2.44 R_⊕$ as the current values of the planetary radii. <sup>(a)</sup>Mass estimated by Zhang et al. (2022). They assumed $r_b = 2.15 R_⊕$ and $r_c = 2.67 R_⊕$ as the current values of the planetary radii.

planetary equilibrium temperature. Neglecting this effect, the atmosphere of a $13 M_⊕$ planet would stop evaporating within 0.6 Gyr, while we found that such a planet keeps losing its atmosphere with a timescale longer than 5 Gyr. Summarising, if the current mass of HD 63433 c is $15 M_⊕$, the atmosphere should be stable against evaporation; for masses in the range 7–13 $M_⊕$ the planet mass and radius do not change appreciably in time, while for the case with $m_c = 5 M_⊕$ we predict a decrease of 0.8% in mass and 25% in radius over the next ~4.5 Gyr.

We acknowledge that HD 63433 has also been followed-up with the CARMENES spectrograph (Mallorquín Díaz et al. 2023). The GAPS and CARMENES teams have coordinated the submission of two studies, which have been intentionally carried out in an independent way.

**Acknowledgements.** This work has been supported by the PRIN-INAFF 2019 "Planetary systems at young ages (PLATEA)" and ASI-INAFF agreement no. 2018-16-HH.0. A.M. and D.L. acknowledge partial financial support from the ASI-INAFF agreement no. 2018-16-HH.0 (THE StellaR PAth project), and from the ARIEL ASI-INAFF agreement no. 2021-5-HH.0. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

## References

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

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

Baratella, M., D'Orazi, V., Biazzo, K., et al. 2020a, A&A, 640, A123

Baratella, M., D'Orazi, V., Carraro, G., et al. 2020b, A&A, 634, A34

Barragán, O., Gandolfi, D., & Antoniciello, G. 2019a, MNRAS, 482, 1017

Barragán, O., Aigrain, S., Kubyshkina, D., et al. 2019b, MNRAS, 490, 698

Barragán, O., Aigrain, S., Rajpaul, V. M., & Zicher, N. 2022a, MNRAS, 509, 866

Barragán, O., Armstrong, D. J., Gandolfi, D., et al. 2022b, MNRAS, 514, 1606

Benatti, S., Damasso, M., Borsa, F., et al. 2021, A&A, 650, A66

Brandt, T. D. 2021, ApJS, 254, 42

Brewer, J. M., Fischer, D. A., Valenti, J. A., & Piskunov, N. 2016, ApJS, 225, 32

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

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

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

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

Casagrande, L., Ramírez, I., Meléndez, J., Bessell, M., & Asplund, M. 2010, A&A, 512, A54

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

Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102

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

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

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. 2021, VizieR Online Data Catalog: II/328

Dai, F., Roy, A., Fulton, B., et al. 2020, AJ, 160, 193

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

Desidera, S., Damasso, M., Gratton, R., et al. 2023, A&A, in press, https://doi.org/10.1051/0004-6361/202244611

Dotter, A. 2016, ApJS, 222, 8

Dutra-Ferreira, L., Pasquini, L., Smiljanic, R., Porto de Mello, G. F., & Steffen, M. 2016, A&A, 585, A75

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

Feroz, F., Balan, S. T., & Hobson, M. P. 2011, MNRAS, 415, 3462

Foreman-Mackey, D. 2018, RNAAS, 2, 31

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

Feroz, F., Hobson, M. P., Cameron, E., & Pettitt, A. N. 2019, Open J. Astrophys., 2, 10

Fossati, L., Erkaev, N. V., Lammer, H., et al. 2017, A&A, 598, A90

Fujita, N., Hori, Y., & Sasaki, T. 2022, ApJ, 928, 105

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

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

Gaia Collaboration 2022, VizieR Online Data Catalog: I/355

Gaudí, B. S., Meyer, M., & Christiansen, J. 2021, in ExoFrontiers; Big Questions in Exoplanetary Science, ed. N. Madhusudhan, 2

Georgieva, I. Y., Persson, C. M., Barragán, O., et al. 2021, MNRAS, 505, 4684

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

Haywood, R. D., Collier Cameron, A., Queloz, D., et al. 2014, MNRAS, 443, 2517

Hirsch, L. A., Rosenthal, L., Fulton, B. J., et al. 2021, AJ, 161, 134

Hög, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, A27

Hunter, A., Macgregor, A. B., Szabo, T., Wellington, C., & Bellgard, M. I. 2012, Source Code for Biology and Medicine, 7, 1

Johnstone, C. P., Bartel, M., & Güdel, M. 2021, A&A, 649, A96

Jones, J., White, R. J., Boyajian, T., et al. 2015, ApJ, 813, 58

Kervella, P., Arenou, F., & Thévenin, F. 2022, A&A, 657, A7

King, G. W., & Wheatley, P. J. 2021, MNRAS, 501, A28

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

Kreidberg, L. 2015, PASP, 127, 1161

Kubyshkina, D., Fossati, L., Erkaev, N. V., et al. 2018a, ApJ, 866, L18

Kubyshkina, D., Fossati, L., Erkaev, N. V., et al. 2018b, A&A, 619, A151

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

Lopez, E. D., & Fortney, J. J. 2014, ApJ, 792, 1

Lopez, E. D., Fortney, J. J., & Miller, N. 2012, ApJ, 761, 59

Maggio, A., Locci, D., Pillitteri, I., et al. 2022, ApJ, 925, 172

Malvolta, L. 2016, Astrophysics Source Code Library, [record ascl:1612.008]

Malvolta, L., Nascimbeni, V., Piotto, G., et al. 2016, A&A, 588, A118

Malvolta, L., Mayo, A. W., Louden, T., et al. 2018, AJ, 155, 107

Mallonn, M., & Strassmeier, K. G. 2016, A&A, 590, A100

Mallorquín Díaz, M., Béjar, V. J. S., Lodieu, N., et al. 2023, A&A, 671, A19

Mann, A. W., Johnson, M. C., Vanderbilt, A., et al. 2020, AJ, 160, 179

Mucciarelli, A., Bellazzini, M., & Massari, D. 2021, A&A, 653, A90

Nardiello, D. 2020, MNRAS, 498, 5972

Nardiello, D., Piotto, G., Deleuil, M., et al. 2020, MNRAS, 495, 4924

Nardiello, D., Deleuil, M., Mantovan, G., et al. 2021, MNRAS, 505, 3767

Nardiello, D., Malvolta, L., Desidera, S., et al. 2022, A&A, 664, A163
Pecaut, M. J., & Mamajek, E. E. 2013, ApJS, 208, 9

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

Perger, M., Anglada-Escudé, G., Ribas, I., et al. 2021, A&A, 645, A58

Poppenhaeger, K., Ketzer, L., & Mallonn, M. 2021, MNRAS, 500, 4560

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

Rajpaul, V., Aigrain, S., Osborne, M. A., Reece, S., & Roberts, S. 2015, MNRAS, 452, 2269

Ricker, G. R., Vanderspek, R., Winn, J., et al. 2016, SPIE Conf. Ser., 9904, 99042B

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

Sanz-Forcada, J., López-Puertas, M., Nortmann, L., & Lampón, M. 2022, 21st Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, https://backoffice.inviteo.com/upload/compte153/Base/inscriptions_projets/fichier/108345-posterjsanzcs21.pdf

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

Sneden, C. 1973, ApJ, 184, 839

Sousa, S. G., Santos, N. C., Adibekyan, V., Delgado-Mena, E., & Israeliian, G. 2015, A&A, 577, A67

Southworth, J. 2011, MNRAS, 417, 2166

Strassmeier, K. G., Granzer, T., Weber, M., et al. 2004, Astron. Nachr., 325, 527

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

Suárez Mascareño, A., Damasso, M., Lodieu, N., et al. 2021, Nat. Astron., 6, 232

Sun, L., Ioannidis, P., Gu, S., et al. 2019, A&A, 624, A15

Vines, J. I., Jenkins, J. S., Berdiñas, Z., et al. 2023, MNRAS, 518, 2627

Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577

Zeng, L., Jacobsen, S. B., Sasselov, D. D., et al. 2019, PNAS, 116, 9723

Zhang, M., Knutson, H. A., Wang, L., et al. 2022, AJ, 163, 68
# Appendix A: Description of the GP models and model priors

Hereafter we provide details about the different GP kernels tested in this work to filter out the stellar activity contribution to the observed RV variability.

**GP quasi-periodic (QP) kernel.** An element of the QP covariance matrix (e.g. Haywood et al. 2014) used in our work is defined as follows:

$$k_{QP}(t, t') = h^2 \exp \left[ -\frac{(t-t')^2}{2\lambda_{QP}^2} - \frac{\sin^2\left(\frac{\pi(t-t')}{\theta}\right)}{2w^2} \right] + (\sigma_{RV}^2(t) + \sigma_{jit}^2) \delta_{t,t'} \quad (A.1)$$

Here, $t$ and $t'$ represent two different epochs of observations, $\sigma_{RV}$ is the radial velocity uncertainty, and $\delta_{t,t'}$ is the Kronecker delta. Our analysis takes into account other sources of uncorrelated noise – instrumental and/or astrophysical – by including a constant jitter term $\sigma_{jit}$ which is added in quadrature to the formal uncertainties $\sigma_{RV}$. The GP hyper-parameters are $h$, which denotes the scale amplitude of the correlated signal; $\theta$, which represents the periodic timescale of the correlated signal, and corresponds to the stellar rotation period; $w$, which describes the "weight" of the rotation period harmonic content within a complete stellar rotation (i.e. a low value of $w$ indicates that the periodic variations contain a significant contribution from the harmonics of the rotation periods); and $\lambda_{QP}$, which represents the decay timescale of the correlations, and is related to the temporal evolution of the magnetically active regions responsible for the correlated signal observed in the RVs.

**GP quasi-periodic with cosine (QPC) kernel.** The QPC kernel has been introduced by Perger et al. (2021). The covariance matrix element implemented in our work is described by the equation

$$k_{QPC}(t, t') = \exp\left(-2\frac{(t-t')^2}{\lambda_{QPC}^2}\right) \cdot \left[h_1^2 \exp\left(-\frac{1}{2w^2}\sin^2\left(\frac{\pi(t-t')}{\theta}\right)\right) + h_2^2 \cos\left(\frac{4\pi(t-t')}{\theta}\right)\right] + (\sigma_{RV}^2(t) + \sigma_{jit}^2) \delta_{t,t'} \quad (A.2)$$

Again, $t$ and $t'$ represent two different epochs of observations; $h_1$ and $h_2$ are scale amplitudes; $\theta$ still represents the periodic time-scale of the modelled signal, and corresponds to the stellar rotation period; $w$ still describes the weight of the rotation period harmonic content within a complete stellar rotation; $\lambda_{QPC}$ is defined as $2 \cdot \lambda_{QP}$, better representing the average lifetime of the activity-related features responsible for the stellar correlated signal in the RVs according to Perger et al. (2021); $\sigma_{RV}(t)$ and $\sigma_{jit}$ are the radial velocity uncertainty at epoch $t$ and the uncorrelated jitter, respectively, and $\delta_{t,t'}$ is the Kronecker delta.

**GP rotational or double simple harmonic oscillator (dSHO) kernel.** The “rotational” kernel is defined by a mixture of two stochastically driven, damped simple harmonic oscillators (SHOs) with undamped periods of $P_{*,rot}$ and $P_{*,rot}/2$. This can be obtained by combining two SHOTerm kernels included in the package celerite (Foreman-Mackey et al. 2017)<fnref n="3" />. The power

spectral density corresponding to this kernel is

$$S(\omega) = \sqrt{\frac{2}{\pi}} \frac{S_1 \omega_1^4}{(\omega^2 - \omega_1^2)^2 + 2\omega_1^2 \omega^2} + \sqrt{\frac{2}{\pi}} \frac{S_2 \omega_2^4}{(\omega^2 - \omega_2^2)^2 + 2\omega_2^2 \omega^2/Q^2}, \quad (A.3)$$

where

$$S_1 = \frac{A^2}{\omega_1 Q_1 (1+f)}, \quad (A.4)$$

$$S_2 = \frac{A^2}{\omega_2 Q_2 (1+f)} \cdot f, \quad (A.5)$$

$$\omega_1 = \frac{4\pi Q_1}{P_{\text{rot}} \sqrt{4Q_1^2 - 1}}, \quad (A.6)$$

$$\omega_2 = \frac{8\pi Q_2}{P_{\text{rot}} \sqrt{4Q_2^2 - 1}}, \quad (A.7)$$

$$Q_1 = \frac{1}{2} + Q_0 + \Delta Q, \quad (A.8)$$

$$Q_2 = \frac{1}{2} + Q_0. \quad (A.9)$$

The parameters in [A.4-A.9], where the subscripts 1 and 2 refer to the primary ($P_{*,rot}$) and secondary ($P_{*,rot}/2$) modes, represent the inputs to the SHOTerm kernels. However, instead of using them directly, we adopt a different parametrization using the following variables as free hyper-parameters in the MC analysis, from which $S_i$, $Q_i$, and $\omega_i$ are derived through Eq. [A.4-A.9]: the variability amplitude $A$, the stellar rotation period $P_{*,rot}$, the quality factor $Q_0$, the difference $\Delta Q$ between the quality factors of the first and second terms, and the fractional amplitude $f$ of the secondary mode relative to the primary.

**Multi-dimensional GP framework.** This framework has been first introduced by Rajpaul et al. (2015). A python version of the code is implemented in the package pyaneti (Barragán et al. 2019a)<fnref n="4" />. In our work, we included the module pyaneti in our own code to make it work with the nested sampler MultiNest. In brief, the multi-dimensional GP framework uses ancillary data sensitive to stellar activity in a combined fit with the RVs (in our case, we used the activity diagnostics BIS and $\log R'_{HK}$), with each dataset represented in the time domain by linear combinations of the same GP-drawn function $G(t)$ and its derivative $\dot{G}(t)$. For all the details of the multi-dimensional GP implementation see Barragán et al. (2022a), especially Sect. 3.1.2 and Eq. (8) therein. For our test, we adopted a QP kernel, as described above, and we fixed to zero the coefficient of $\dot{G}(t)$ for the case of $\log R'_{HK}$.

The prior distributions for the free parameters of the models tested in this work are provided in Table A.1.

<footnote n="3">³ https://github.com/dfm/celerite/blob/main/celerite/terms.py</footnote>

<footnote n="4">⁴ https://github.com/oscaribv/pyaneti</footnote>
<table><caption>Table A.1. Priors used for the models listed in Table 2.</caption><thead><tr><th>Parameter</th><th>Prior<sup>b</sup></th></tr></thead><tbody><tr><td colspan="2"><strong>Stellar Activity - QP kernel</strong></td></tr><tr><td>h [m s<sup>-1</sup>]</td><td>U(0, 500)</td></tr><tr><td>θ [d]</td><td>U(0, 10)</td></tr><tr><td>λ<sub>QP</sub> [d]</td><td>U(0, 1000)</td></tr><tr><td>w</td><td>U(0, 1)</td></tr><tr><td colspan="2"><strong>Stellar Activity - QPC kernel</strong></td></tr><tr><td>h<sub>1</sub> [m s<sup>-1</sup>]</td><td>U(0, 500)</td></tr><tr><td>h<sub>2</sub> [m s<sup>-1</sup>]</td><td>U(0, 500)</td></tr><tr><td>θ [d]</td><td>U(0, 10)</td></tr><tr><td>λ<sub>QPC</sub> [d]</td><td>U(0, 1000)</td></tr><tr><td>w</td><td>U(0, 1)</td></tr><tr><td colspan="2"><strong>Stellar Activity - dSHO kernel</strong></td></tr><tr><td>log A</td><td>U(0.05,10)</td></tr><tr><td>θ [d]</td><td>U(0,10)</td></tr><tr><td>log Q<sub>0</sub></td><td>U(-10,10)</td></tr><tr><td>log ΔQ</td><td>U(-10,10)</td></tr><tr><td>f</td><td>U(0,10)</td></tr><tr><td colspan="2"><strong>Stellar Activity - Multi-dim. QP kernel<sup>b</sup></strong></td></tr><tr><td>A<sub>RV</sub>, B<sub>RV</sub> [km s<sup>-1</sup>]</td><td>U(-0.1,0.1), U(-0.1,0.1)</td></tr><tr><td>A<sub>BIS</sub>, B<sub>BIS</sub> [km s<sup>-1</sup>]</td><td>U(-0.1,0.1), U(-0.1,0.1)</td></tr><tr><td>A<sub>R'HK</sub>, B<sub>R'HK</sub> [dex]</td><td>U(-0.1,0.1), fixed to 0</td></tr><tr><td colspan="2"><strong>Keplerian of planet b</strong></td></tr><tr><td>K<sub>b</sub> [m s<sup>-1</sup>]</td><td>U(0,30)</td></tr><tr><td>P<sub>b</sub> [d]</td><td>U(6.9,7.3)</td></tr><tr><td>T<sub>conj,b</sub> [BJD-2 450 000]</td><td>U(9584.4,9584.7)</td></tr><tr><td>√e<sub>b</sub> cos ω<sub>*,b</sub></td><td>U(-1,1)</td></tr><tr><td>√e<sub>b</sub> sin ω<sub>*,b</sub></td><td>U(-1,1)</td></tr><tr><td colspan="2"><strong>Keplerian of planet c</strong></td></tr><tr><td>K<sub>c</sub> [m s<sup>-1</sup>]</td><td>U(0,30)</td></tr><tr><td>P<sub>c</sub> [d]</td><td>U(19.5,21.5)</td></tr><tr><td>T<sub>conj,c</sub> [BJD-2 450 000]</td><td>U(9583.5,9583.8)</td></tr><tr><td>√e<sub>c</sub> cos ω<sub>*,c</sub></td><td>U(-1,1)</td></tr><tr><td>√e<sub>c</sub> sin ω<sub>*,c</sub></td><td>U(-1,1)</td></tr><tr><td>RV acceleration, ˙γ [m s<sup>-1</sup> d<sup>-1</sup>]</td><td>U(-1,1)</td></tr><tr><td colspan="2"><strong>Transit parameters for planet b</strong></td></tr><tr><td>R<sub>b</sub>/R<sub>*</sub></td><td>U(0,1)</td></tr><tr><td>orbital plane inclination, i<sub>b</sub> [deg]</td><td>U(80,90)</td></tr><tr><td colspan="2"><strong>Transit parameters for planet c</strong></td></tr><tr><td>R<sub>c</sub>/R<sub>*</sub></td><td>U(0,1)</td></tr><tr><td>orbital plane inclination, i<sub>c</sub> [deg]</td><td>U(80,90)</td></tr><tr><td colspan="2"><strong>Instrument-related</strong></td></tr><tr><td>γ<sub>RV</sub> [m s<sup>-1</sup>]</td><td>DRS: U(-16000,-15500)<br>TERRA: U(-100,100)</td></tr><tr><td>σ<sub>jit,RV</sub> [m s<sup>-1</sup>]</td><td>U(0,50)</td></tr><tr><td>σ<sub>jit, TESS sect. 20</sub></td><td>U(0,300e-5)</td></tr><tr><td>σ<sub>jit, TESS sect. 44-47</sub></td><td>U(0,300e-5)</td></tr><tr><td>Limb darkening param.<sup>c</sup> q<sub>1</sub></td><td>U(0,1)</td></tr><tr><td>Limb darkening param.<sup>c</sup> q<sub>2</sub></td><td>U(0,1)</td></tr></tbody></table>

Notes. (a) $U(a, b)$ denotes a prior drawn from an uninformative distribution in the range (a,b). $N(a, b)$ denotes a prior drawn from a Gaussian distribution with mean $a$ and sigma $b$. (b) We follow the notation of Eq. (8) in Barragán et al. (2022a) (c) Following Kipping (2013), $q_1$ and $q_2$ are related to the limb darkening coefficients LD$_{c1}$ and LD$_{c2}$ by the relations $LD_{c1} = 2q_2 \sqrt{q_1}$ and $LD_{c2} = \sqrt{q_1}(1 - 2q_2)$.

## Appendix B: Details on atmospheric evaporation modelling

The hydro-based approximation developed by Kubyshkina et al. (2018a,b) provides an analytic expression of the planetary mass-loss rate as a function of planetary mass and radius, star-planet distance, and stellar high-energy flux. We employed this approximation to model the time evolution of the planet characteristics, taking into account both the stellar evolutionary track and the long-term change of the XUV irradiation.

For the X-ray evolution as a function of time, we used the description by Penz et al. (2008), based on open cluster and field G-type stars. It describes the saturation and decay phases of stellar activity with simple analytic power laws. Maggio et al. (2022) showed that they represent a good approximation of the more complex behaviour modelled by Johnstone et al. (2021) for late-type stars. In order to evaluate the stellar irradiation in an extreme ultraviolet band we adopted a scaling law between EUV (100–920 Å) and X-ray (5–100 Å) luminosities, calibrated on stars observed in both bands, derived by Sanz-Forcada et al. (2022), namely

$$ \log L_{\text{EUV}} = (0.793 \pm 0.058) \log L_{\chi} + (6.53 \pm 1.61) \quad (\text{B.1}) $$

which is an updated version of the better-known and widely used relationship by Sanz-Forcada et al. (2011). The predicted evolution of the X-ray and EUV stellar luminosities is shown in Fig. C.4.

We also took into account the evolution of the planetary radius according to Lopez & Fortney (2014). This relation was developed for H-He dominated atmospheres, and provides the envelope radius $R_{\text{env}}$, as a function of the planetary mass, the atmospheric mass fraction $f_{\text{atm}}$, the bolometric flux received, and the age of the system. It takes into account also the cooling and contraction of the envelope as a consequence of its thermal evolution, based on the grid of models presented in Lopez et al. (2012).

The track followed by HD 63433 in the theoretical temperature-luminosity diagram is shown in Fig. C.3. It was computed for a star with a mass 0.994$M_{\odot}$ and metallicity $[\text{Fe/H}] = +0.02$ (Sect. 3) using the web-based interpolator<fnref n="5" /> of the MESA Isochrones and Stellar Tracks (MIST, Choi et al. 2016).

The Lopez & Fortney (2014) relation can be inverted to estimate $f_{\text{atm}}$ at the current age, starting with the known (measured) planetary radius and with a core radius, $R_{\text{core}}$, computed as in Benatti et al. (2021), assuming that all the planetary mass is concentrated in the core. Different assumptions (or future measurements) of the planetary mass, yield different $R_{\text{core}}$ and $f_{\text{atm}}$ values.

We then proceeded with the time evolution. For each time step of the simulation, we computed the mass-loss rate and we updated $f_{\text{atm}}$ and the planetary mass, thus obtaining a new value of $R_{\text{env}}$ with the relation by Lopez & Fortney (2014). The latter quantity added to the core radius (assumed constant) provides the updated planetary radius. For each starting planetary mass, we let the system evolve from the current age (414 Myr for HD 63433) until 5 Gyr. In order to check if the atmosphere is hydrodynamically stable or not, we use the Jeans escape parameter

$$ \Lambda = \frac{G m_{\text{H}} M_{\text{p}}}{k_{\text{B}} T_{\text{eq}} R_{\text{p}}} \quad (\text{B.2}) $$

<footnote n="5">⁵ https://waps.cfa.harvard.edu/MIST/ interp tracks.html</footnote>
where G is the gravitational constant, $m_H$ is the hydrogen mass, and $k_B$ is the Boltzman constant, while $M_p$ and $R_p$ are the planet mass and radius, and $T_{eq}$ is the planet equilibrium temperature, computed as

$$T_{\text{eq}} = T_{\text{eff}} \left[ f_p (1 - A_B) \right]^{1/4} \left( \frac{R_*}{2d} \right)^{1/2} \qquad (B.3)$$

where $T_{\text{eff}}$ and $R_*$ are the stellar effective temperature and radius, respectively, $A_B$ is the Bond albedo, $d$ the star-planet distance, $f_p$ is a parameter that takes into account if the planet is tidally locked or not. We assumed $A_B = 0.5$, and $f_p = \frac{2}{3}$, because tidal locking is assumed as in (Zhang et al. 2022).

Note that $\Lambda$ depends on the stellar $T_{\text{eff}}$ through $T_{\text{eq}}$, and $\Lambda$ enters also in the analytic expression for the atmospheric mass-loss rate derived by Kubyshkina et al. (2018b), which is also a function of the XUV flux, $F_{XUV}$, at the planetary orbital position. According to Fossati et al. (2017), photoevaporation of planetary atmospheres occurs when the Jeans escape parameter is smaller than a critical value ($\Lambda < 80$), but two different evaporation regimes are possible, depending on how $\Lambda$ compares with a control parameter $e^\Sigma$, determined by $d$ and by time-dependent values of $R_p$ and $F_{XUV}$ (Kubyshkina et al. 2018b): if $\Lambda > e^\Sigma$, the evaporation is controlled by the XUV flux, while for smaller values of $\Lambda$ the atmospheric escape is mainly driven by the planetary intrinsic thermal energy and low gravity.

In this work, assuming an X-ray luminosity of $5 \times 10^{28}$ erg/s at the current age of the system, and allowing for the evolution of the bolometric luminosity and the effective temperature of the star as described above (Fig. C.3), we investigated the hydrodynamic atmospheric evolution of the two planets of HD 63433 in the future. Furthermore, we studied the system's evolution back in time, as described by Georgieva et al. (2021). According to the aforementioned scenario and assuming a core which does not change in size or mass, we created a synthetic population of planets with different initial masses, and we assigned to them the core radius and core mass fixed by the total mass assumed for the planet at its present age. This framework leaves $f_{\text{atm}}$ to dictate the total initial mass, while the total radius is again based on the analytic fit by Lopez & Fortney (2014). Then we searched for the planet in the synthetic set which ends up with mass, radius and $f_{\text{atm}}$ most similar to the values at the present age, and we selected its predicted evolutionary history. We started all the simulations in the past at an age of 10 Myr, when we suppose that the circumstellar disk is already disappeared and all planets are in their final, stable orbit. These simulations always end at an age of 414 Myr.
**Appendix C: Additional plots**

<figure><img src="image_11.png" /><figcaption>**Fig. C.1.** Time series and GLS periodograms of the light curves collected by STELLA in V and I bands. The horizontal dashed line denotes the FAP level of 10%, the vertical dashed line identifies the stellar rotation period.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. C.2.** RVs as a function of the CCF-BIS, showing the significant correlation between the two dataset.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Fig. C.3.** Evolutionary track of HD 63433 up to an age of 5 Gyr in the effective temperature-bolometric luminosity plane. The current location of the star on the track is marked by a black cross.</figcaption></figure>
<figure><img src="image_14.png" /></figure>

**Fig. C.4.** Time evolution of X-ray (5–100 Å), EUV (100–920 Å), and total XUV luminosity of HD 63433, according to Penz et al. (2008) and the X-ray/EUV scaling by Sanz-Forcada et al. (2022). The grey area is the observed $1\sigma$ spread around the median (dark grey line) of the X-ray luminosity distributions for dG stars in the Hyades and Pleiades open clusters. Uncertainties on the age and X-ray luminosity of HD 63433 are within the size of the square symbol.
## Appendix D: Dataset

<table><caption>Table D.1. RVs derived from HARPS-N spectra and used in this work.</caption><thead><tr><th>Time [BJD]</th><th><i>RV</i><sub>DRS</sub> [m s<sup>-1</sup>]</th><th><i>σ</i><sub>RV, DRS</sub> [m s<sup>-1</sup>]</th><th><i>RV</i><sub>TERRA</sub> [m s<sup>-1</sup>]</th><th><i>σ</i><sub>RV, TERRA</sub> [m s<sup>-1</sup>]</th></tr></thead><tbody><tr><td>2458906.346145</td><td>-15775.70</td><td>0.74</td><td>32.74</td><td>1.24</td></tr><tr><td>2458908.346263</td><td>-15834.58</td><td>0.70</td><td>-29.73</td><td>1.10</td></tr><tr><td>2458909.336032</td><td>-15819.72</td><td>0.64</td><td>-12.85</td><td>0.94</td></tr><tr><td>2458910.333356</td><td>-15767.51</td><td>0.75</td><td>41.23</td><td>1.15</td></tr><tr><td>2458911.344394</td><td>-15826.59</td><td>1.40</td><td>-22.43</td><td>1.78</td></tr><tr><td>2458920.541754</td><td>-15825.27</td><td>0.91</td><td>-20.51</td><td>1.45</td></tr><tr><td>2458921.474998</td><td>-15810.30</td><td>0.74</td><td>-4.86</td><td>1.09</td></tr><tr><td>2458923.506504</td><td>-15815.76</td><td>0.72</td><td>-10.45</td><td>1.26</td></tr><tr><td>2458924.532470</td><td>-15788.20</td><td>0.90</td><td>16.66</td><td>1.18</td></tr><tr><td>2458974.419605</td><td>-15765.59</td><td>0.71</td><td>42.03</td><td>1.25</td></tr><tr><td>2458976.369195</td><td>-15811.03</td><td>0.62</td><td>-5.36</td><td>0.91</td></tr><tr><td>2458977.382239</td><td>-15760.20</td><td>0.60</td><td>47.91</td><td>1.22</td></tr><tr><td>2458978.378076</td><td>-15850.47</td><td>0.83</td><td>-47.17</td><td>1.36</td></tr><tr><td>2458979.374190</td><td>-15800.85</td><td>0.91</td><td>4.99</td><td>1.24</td></tr><tr><td>2458982.404317</td><td>-15844.09</td><td>0.77</td><td>-39.61</td><td>1.39</td></tr><tr><td>2458983.379788</td><td>-15761.20</td><td>1.56</td><td>47.97</td><td>1.68</td></tr><tr><td>2458985.373704</td><td>-15814.69</td><td>0.81</td><td>-8.34</td><td>1.04</td></tr><tr><td>2458986.372046</td><td>-15788.50</td><td>0.76</td><td>18.13</td><td>1.10</td></tr><tr><td>2458988.370885</td><td>-15818.16</td><td>0.98</td><td>-13.09</td><td>1.41</td></tr><tr><td>2458989.359219</td><td>-15800.55</td><td>1.11</td><td>4.54</td><td>1.51</td></tr><tr><td>2458990.358976</td><td>-15805.23</td><td>1.01</td><td>2.18</td><td>1.39</td></tr><tr><td>2458996.373074</td><td>-15779.65</td><td>0.87</td><td>28.57</td><td>1.18</td></tr><tr><td>2459125.758727</td><td>-15824.03</td><td>2.01</td><td>-20.72</td><td>2.01</td></tr><tr><td>2459126.739655</td><td>-15805.39</td><td>1.02</td><td>1.09</td><td>1.34</td></tr><tr><td>2459127.755228</td><td>-15819.27</td><td>0.84</td><td>-12.24</td><td>1.25</td></tr><tr><td>2459130.744963</td><td>-15793.85</td><td>1.08</td><td>16.01</td><td>1.33</td></tr><tr><td>2459153.732409</td><td>-15821.53</td><td>1.33</td><td>-14.29</td><td>0.71</td></tr><tr><td>2459154.676285</td><td>-15842.67</td><td>1.74</td><td>-32.27</td><td>2.31</td></tr><tr><td>2459169.686433</td><td>-15830.39</td><td>1.23</td><td>-23.63</td><td>0.68</td></tr><tr><td>2459170.680175</td><td>-15788.51</td><td>0.61</td><td>20.72</td><td>0.54</td></tr><tr><td>2459172.580845</td><td>-15798.06</td><td>0.74</td><td>4.53</td><td>0.43</td></tr><tr><td>2459197.786154</td><td>-15832.84</td><td>0.78</td><td>-26.09</td><td>0.41</td></tr><tr><td>2459212.568090</td><td>-15807.42</td><td>0.72</td><td>2.61</td><td>0.96</td></tr><tr><td>2459229.676377</td><td>-15822.43</td><td>1.52</td><td>-13.77</td><td>1.53</td></tr><tr><td>2459235.587963</td><td>-15831.56</td><td>0.51</td><td>-24.12</td><td>0.83</td></tr><tr><td>2459236.520232</td><td>-15811.93</td><td>0.71</td><td>-5.25</td><td>0.92</td></tr><tr><td>2459236.606006</td><td>-15808.53</td><td>1.12</td><td>-0.93</td><td>1.54</td></tr><tr><td>2459237.585692</td><td>-15832.90</td><td>0.65</td><td>-26.83</td><td>0.89</td></tr><tr><td>2459239.569178</td><td>-15807.09</td><td>0.65</td><td>-0.11</td><td>0.91</td></tr><tr><td>2459240.567944</td><td>-15809.60</td><td>0.55</td><td>-2.16</td><td>0.76</td></tr><tr><td>2459244.604205</td><td>-15797.66</td><td>0.47</td><td>10.58</td><td>0.73</td></tr><tr><td>2459245.554757</td><td>-15762.05</td><td>0.57</td><td>47.67</td><td>1.13</td></tr><tr><td>2459246.374595</td><td>-15848.49</td><td>1.50</td><td>-41.22</td><td>1.91</td></tr><tr><td>2459246.503411</td><td>-15850.74</td><td>0.42</td><td>-43.41</td><td>1.06</td></tr><tr><td>2459247.608735</td><td>-15787.92</td><td>0.53</td><td>19.56</td><td>0.93</td></tr><tr><td>2459253.606306</td><td>-15788.81</td><td>0.61</td><td>19.21</td><td>0.92</td></tr><tr><td>2459254.521233</td><td>-15812.67</td><td>0.51</td><td>-4.90</td><td>0.79</td></tr><tr><td>2459255.554371</td><td>-15814.67</td><td>0.56</td><td>-6.64</td><td>0.78</td></tr><tr><td>2459256.531086</td><td>-15800.52</td><td>0.70</td><td>7.63</td><td>0.97</td></tr><tr><td>2459258.595699</td><td>-15780.58</td><td>0.69</td><td>28.10</td><td>0.97</td></tr><tr><td>2459272.584028</td><td>-15797.32</td><td>0.71</td><td>11.59</td><td>0.79</td></tr><tr><td>2459275.463016</td><td>-15773.16</td><td>0.77</td><td>36.47</td><td>1.18</td></tr><tr><td>2459276.476808</td><td>-15829.40</td><td>0.62</td><td>-21.49</td><td>0.93</td></tr><tr><td>2459277.498226</td><td>-15783.24</td><td>0.60</td><td>26.27</td><td>0.95</td></tr><tr><td>2459288.524057</td><td>-15797.25</td><td>0.93</td><td>10.73</td><td>1.10</td></tr><tr><td>2459290.355022</td><td>-15779.53</td><td>0.53</td><td>29.68</td><td>0.87</td></tr><tr><td>2459293.348057</td><td>-15832.87</td><td>0.52</td><td>-25.32</td><td>0.80</td></tr><tr><td>2459299.513734</td><td>-15826.13</td><td>0.86</td><td>-19.04</td><td>1.12</td></tr><tr><td>2459303.465955</td><td>-15787.00</td><td>0.85</td><td>21.79</td><td>1.06</td></tr><tr><td>2459304.439774</td><td>-15803.78</td><td>1.01</td><td>3.31</td><td>1.11</td></tr><tr><td>2459305.455985</td><td>-15811.91</td><td>1.17</td><td>-3.59</td><td>1.28</td></tr><tr><td>2459307.428793</td><td>-15806.26</td><td>0.70</td><td>-1.68</td><td>0.83</td></tr><tr><td>2459327.354815</td><td>-15847.69</td><td>0.57</td><td>-42.36</td><td>1.05</td></tr><tr><td>2459328.357969</td><td>-15774.60</td><td>0.72</td><td>34.63</td><td>1.28</td></tr><tr><td>2459330.442524</td><td>-15773.70</td><td>0.93</td><td>33.45</td><td>1.12</td></tr><tr><td>2459338.380324</td><td>-15826.70</td><td>0.75</td><td>-19.87</td><td>0.85</td></tr><tr><td>2459343.387805</td><td>-15793.89</td><td>1.67</td><td>13.73</td><td>1.84</td></tr><tr><td>2459344.386754</td><td>-15812.97</td><td>1.06</td><td>-3.13</td><td>1.24</td></tr><tr><td>2459345.409001</td><td>-15795.02</td><td>1.13</td><td>11.16</td><td>1.32</td></tr><tr><td>2459346.388649</td><td>-15857.68</td><td>1.00</td><td>-53.50</td><td>1.46</td></tr><tr><td>2459352.368704</td><td>-15837.85</td><td>1.61</td><td>-34.23</td><td>2.29</td></tr><tr><td>2459358.388673</td><td>-15805.17</td><td>1.10</td><td>2.74</td><td>1.31</td></tr><tr><td>2459541.678592</td><td>-15789.95</td><td>0.62</td><td>19.35</td><td>0.45</td></tr><tr><td>2459560.749308</td><td>-15758.89</td><td>0.62</td><td>50.50</td><td>0.48</td></tr></tbody></table>
Table D.1. continued.

<table><thead><tr><th>Time [BJD]</th><th><i>RV</i><sub>DRS</sub> [m s<sup>-1</sup>]</th><th><i>σ</i><sub>RV, DRS</sub> [m s<sup>-1</sup>]</th><th><i>RV</i><sub>TERRA</sub> [m s<sup>-1</sup>]</th><th><i>σ</i><sub>RV, TERRA</sub> [m s<sup>-1</sup>]</th></tr></thead><tbody><tr><td>2459565.539049</td><td>-15814.38</td><td>0.67</td><td>-7.12</td><td>0.50</td></tr><tr><td>2459589.663616</td><td>-15811.75</td><td>0.95</td><td>-7.32</td><td>1.35</td></tr><tr><td>2459601.673776</td><td>-15779.54</td><td>0.77</td><td>28.53</td><td>1.31</td></tr><tr><td>2459614.563832</td><td>-15791.46</td><td>1.23</td><td>14.31</td><td>1.87</td></tr><tr><td>2459616.597751</td><td>-15815.08</td><td>1.16</td><td>-9.99</td><td>1.36</td></tr><tr><td>2459617.395140</td><td>-15775.45</td><td>1.65</td><td>29.21</td><td>1.85</td></tr><tr><td>2459618.567952</td><td>-15829.05</td><td>1.92</td><td>-24.25</td><td>2.09</td></tr><tr><td>2459619.461850</td><td>-15816.31</td><td>1.04</td><td>-11.76</td><td>1.70</td></tr><tr><td>2459620.613258</td><td>-15778.00</td><td>2.69</td><td>24.84</td><td>2.92</td></tr><tr><td>2459628.408285</td><td>-15805.65</td><td>0.84</td><td>0.95</td><td>1.68</td></tr><tr><td>2459629.419883</td><td>-15824.68</td><td>0.70</td><td>-20.32</td><td>1.52</td></tr><tr><td>2459642.494913</td><td>-15799.08</td><td>0.75</td><td>4.61</td><td>1.28</td></tr><tr><td>2459644.571088</td><td>-15836.18</td><td>0.77</td><td>-31.63</td><td>1.47</td></tr><tr><td>2459646.370440</td><td>-15794.65</td><td>0.58</td><td>11.40</td><td>1.19</td></tr><tr><td>2459647.458259</td><td>-15786.00</td><td>0.60</td><td>20.26</td><td>1.27</td></tr><tr><td>2459648.542710</td><td>-15835.27</td><td>0.61</td><td>-30.67</td><td>1.46</td></tr><tr><td>2459655.483703</td><td>-15788.79</td><td>1.38</td><td>4.85</td><td>1.62</td></tr><tr><td>2459656.538090</td><td>-15800.28</td><td>1.55</td><td>-5.22</td><td>1.99</td></tr><tr><td>2459658.467214</td><td>-15809.20</td><td>1.55</td><td>-16.37</td><td>1.96</td></tr><tr><td>2459659.454730</td><td>-15769.76</td><td>0.55</td><td>24.23</td><td>1.59</td></tr><tr><td>2459660.465796</td><td>-15788.29</td><td>1.85</td><td>8.20</td><td>2.30</td></tr><tr><td>2459661.505656</td><td>-15827.78</td><td>1.29</td><td>-32.00</td><td>2.00</td></tr><tr><td>2459662.478287</td><td>-15793.48</td><td>1.24</td><td>1.29</td><td>1.86</td></tr><tr><td>2459682.409509</td><td>-15781.05</td><td>0.66</td><td>12.83</td><td>1.47</td></tr><tr><td>2459683.445090</td><td>-15852.01</td><td>1.07</td><td>-57.51</td><td>1.85</td></tr><tr><td>2459684.444817</td><td>-15792.74</td><td>1.35</td><td>1.81</td><td>2.67</td></tr><tr><td>2459685.465805</td><td>-15788.30</td><td>2.26</td><td>9.12</td><td>2.94</td></tr><tr><td>2459686.470475</td><td>-15819.33</td><td>0.91</td><td>-24.76</td><td>1.72</td></tr><tr><td>2459687.367088</td><td>-15796.11</td><td>2.52</td><td>0.00</td><td>2.86</td></tr></tbody></table>
<table><caption>Table D.2. Activity diagnostics derived from HARPS-N spectra and used in this work.</caption><thead><tr><th>Time [BJD]</th><th>BIS<sub>DRS</sub> [m s<sup>-1</sup>]</th><th>σ<sub>BIS</sub> [m s<sup>-1</sup>]</th><th>FWHM<sub>DRS</sub> [m s<sup>-1</sup>]</th><th>σ<sub>FWHM</sub> [m s<sup>-1</sup>]</th><th>log R'<sub>HK</sub> [dex]</th><th>σ<sub>log R'<sub>HK</sub></sub> [dex]</th><th>Hα</th><th>σ<sub>Hα</sub></th></tr></thead><tbody><tr><td>2458906.346145</td><td>-2.54</td><td>1.48</td><td>11463.30</td><td>1.48</td><td>-4.378</td><td>0.0010</td><td>0.1580</td><td>0.0003</td></tr><tr><td>2458908.346263</td><td>39.46</td><td>1.41</td><td>11399.39</td><td>1.41</td><td>-4.407</td><td>0.0009</td><td>0.1545</td><td>0.0003</td></tr><tr><td>2458909.336032</td><td>33.79</td><td>1.29</td><td>11409.17</td><td>1.29</td><td>-4.423</td><td>0.0009</td><td>0.1513</td><td>0.0002</td></tr><tr><td>2458910.333356</td><td>-68.54</td><td>1.51</td><td>11474.73</td><td>1.51</td><td>-4.400</td><td>0.0010</td><td>0.1611</td><td>0.0003</td></tr><tr><td>2458911.344394</td><td>19.68</td><td>2.79</td><td>11446.55</td><td>2.79</td><td>-4.394</td><td>0.0021</td><td>0.1553</td><td>0.0004</td></tr><tr><td>2458920.541754</td><td>59.38</td><td>1.82</td><td>11492.41</td><td>1.82</td><td>-4.387</td><td>0.0012</td><td>0.1554</td><td>0.0003</td></tr><tr><td>2458921.474998</td><td>-8.13</td><td>1.48</td><td>11411.82</td><td>1.48</td><td>-4.425</td><td>0.0013</td><td>0.1525</td><td>0.0002</td></tr><tr><td>2458923.506504</td><td>62.76</td><td>1.44</td><td>11455.82</td><td>1.44</td><td>-4.396</td><td>0.0009</td><td>0.1588</td><td>0.0002</td></tr><tr><td>2458924.532470</td><td>-33.60</td><td>1.80</td><td>11440.11</td><td>1.80</td><td>-4.391</td><td>0.0011</td><td>0.1540</td><td>0.0002</td></tr><tr><td>2458974.419605</td><td>-44.05</td><td>1.42</td><td>11444.22</td><td>1.42</td><td>-4.388</td><td>0.0009</td><td>0.1563</td><td>0.0002</td></tr><tr><td>2458976.369195</td><td>7.00</td><td>1.24</td><td>11331.89</td><td>1.24</td><td>-4.405</td><td>0.0008</td><td>0.1594</td><td>0.0002</td></tr><tr><td>2458977.382239</td><td>-50.08</td><td>1.21</td><td>11491.63</td><td>1.21</td><td>-4.394</td><td>0.0008</td><td>0.1600</td><td>0.0002</td></tr><tr><td>2458978.378076</td><td>70.27</td><td>1.66</td><td>11410.23</td><td>1.66</td><td>-4.411</td><td>0.0011</td><td>0.1516</td><td>0.0002</td></tr><tr><td>2458979.374190</td><td>-17.86</td><td>1.82</td><td>11366.34</td><td>1.82</td><td>-4.421</td><td>0.0013</td><td>0.1508</td><td>0.0003</td></tr><tr><td>2458982.404317</td><td>73.69</td><td>1.54</td><td>11392.66</td><td>1.54</td><td>-4.407</td><td>0.0011</td><td>0.1539</td><td>0.0003</td></tr><tr><td>2458983.379788</td><td>-54.20</td><td>3.12</td><td>11399.07</td><td>3.12</td><td>-4.398</td><td>0.0024</td><td>0.1544</td><td>0.0005</td></tr><tr><td>2458985.373704</td><td>0.29</td><td>1.63</td><td>11319.15</td><td>1.63</td><td>-4.420</td><td>0.0011</td><td>0.1505</td><td>0.0002</td></tr><tr><td>2458986.372046</td><td>-27.89</td><td>1.51</td><td>11398.43</td><td>1.51</td><td>-4.420</td><td>0.0011</td><td>0.1520</td><td>0.0002</td></tr><tr><td>2458988.370885</td><td>32.11</td><td>1.97</td><td>11466.71</td><td>1.97</td><td>-4.385</td><td>0.0013</td><td>0.1626</td><td>0.0003</td></tr><tr><td>2458989.359219</td><td>-23.65</td><td>2.22</td><td>11369.92</td><td>2.22</td><td>-4.394</td><td>0.0015</td><td>0.1562</td><td>0.0003</td></tr><tr><td>2458990.358976</td><td>-22.23</td><td>2.02</td><td>11519.07</td><td>2.02</td><td>-4.393</td><td>0.0014</td><td>0.1568</td><td>0.0003</td></tr><tr><td>2458996.373074</td><td>-44.24</td><td>1.74</td><td>11447.19</td><td>1.74</td><td>-4.401</td><td>0.0012</td><td>0.1545</td><td>0.0003</td></tr><tr><td>2459125.758727</td><td>0.01</td><td>4.01</td><td>11172.08</td><td>4.01</td><td>-4.405</td><td>0.0025</td><td>0.1498</td><td>0.0004</td></tr><tr><td>2459126.739655</td><td>-24.39</td><td>2.04</td><td>11362.44</td><td>2.04</td><td>-4.408</td><td>0.0014</td><td>0.1513</td><td>0.0003</td></tr><tr><td>2459127.755228</td><td>29.60</td><td>1.67</td><td>11442.62</td><td>1.67</td><td>-4.404</td><td>0.0011</td><td>0.1517</td><td>0.0002</td></tr><tr><td>2459130.744963</td><td>4.05</td><td>2.16</td><td>11476.10</td><td>2.16</td><td>-4.385</td><td>0.0014</td><td>0.1636</td><td>0.0003</td></tr><tr><td>2459153.732409</td><td>19.33</td><td>2.66</td><td>11461.25</td><td>2.66</td><td>-4.412</td><td>0.0020</td><td>0.1519</td><td>0.0004</td></tr><tr><td>2459154.676285</td><td>28.76</td><td>3.49</td><td>11314.73</td><td>3.49</td><td>-4.428</td><td>0.0030</td><td>0.1489</td><td>0.0005</td></tr><tr><td>2459169.686433</td><td>29.10</td><td>2.46</td><td>11377.57</td><td>2.46</td><td>-4.386</td><td>0.0017</td><td>0.1636</td><td>0.0004</td></tr><tr><td>2459170.680175</td><td>-21.88</td><td>1.23</td><td>11449.45</td><td>1.23</td><td>-4.389</td><td>0.0008</td><td>0.1575</td><td>0.0002</td></tr><tr><td>2459172.580845</td><td>-12.34</td><td>1.48</td><td>11180.50</td><td>1.48</td><td>-4.397</td><td>0.0008</td><td>0.1535</td><td>0.0002</td></tr><tr><td>2459197.786154</td><td>58.33</td><td>1.56</td><td>11193.44</td><td>1.56</td><td>-4.389</td><td>0.0008</td><td>0.1608</td><td>0.0002</td></tr><tr><td>2459212.568090</td><td>6.42</td><td>1.44</td><td>11301.57</td><td>1.44</td><td>-4.408</td><td>0.0009</td><td>0.1539</td><td>0.0002</td></tr><tr><td>2459229.676377</td><td>-2.09</td><td>3.04</td><td>11411.19</td><td>3.04</td><td>-4.401</td><td>0.0024</td><td>0.1551</td><td>0.0005</td></tr><tr><td>2459235.587963</td><td>36.10</td><td>1.03</td><td>11417.11</td><td>1.03</td><td>-4.411</td><td>0.0007</td><td>0.1534</td><td>0.0002</td></tr><tr><td>2459236.520232</td><td>-1.22</td><td>1.43</td><td>11406.18</td><td>1.43</td><td>-4.412</td><td>0.0010</td><td>0.1581</td><td>0.0003</td></tr><tr><td>2459236.606006</td><td>-0.54</td><td>2.23</td><td>11405.12</td><td>2.23</td><td>-4.418</td><td>0.0016</td><td>0.1520</td><td>0.0003</td></tr><tr><td>2459237.585692</td><td>31.91</td><td>1.29</td><td>11450.56</td><td>1.29</td><td>-4.404</td><td>0.0008</td><td>0.1542</td><td>0.0002</td></tr><tr><td>2459239.569178</td><td>18.19</td><td>1.30</td><td>11557.06</td><td>1.30</td><td>-4.379</td><td>0.0008</td><td>0.1644</td><td>0.0002</td></tr><tr><td>2459240.567944</td><td>-10.94</td><td>1.10</td><td>11396.16</td><td>1.10</td><td>-4.387</td><td>0.0007</td><td>0.1580</td><td>0.0002</td></tr><tr><td>2459244.604205</td><td>0.37</td><td>0.94</td><td>11352.52</td><td>0.94</td><td>-4.409</td><td>0.0006</td><td>0.1616</td><td>0.0002</td></tr><tr><td>2459245.554757</td><td>-52.17</td><td>1.14</td><td>11501.32</td><td>1.14</td><td>-4.386</td><td>0.0007</td><td>0.1604</td><td>0.0002</td></tr><tr><td>2459246.374595</td><td>79.91</td><td>2.99</td><td>11471.85</td><td>2.99</td><td>-4.385</td><td>0.0023</td><td>0.1588</td><td>0.0005</td></tr><tr><td>2459246.503411</td><td>75.60</td><td>0.85</td><td>11452.11</td><td>0.85</td><td>-4.387</td><td>0.0005</td><td>0.1607</td><td>0.0002</td></tr><tr><td>2459247.608735</td><td>-15.50</td><td>1.05</td><td>11432.80</td><td>1.05</td><td>-4.398</td><td>0.0007</td><td>0.1570</td><td>0.0002</td></tr><tr><td>2459253.606306</td><td>-22.40</td><td>1.23</td><td>11335.92</td><td>1.23</td><td>-4.391</td><td>0.0008</td><td>0.1607</td><td>0.0002</td></tr><tr><td>2459254.521233</td><td>42.08</td><td>1.03</td><td>11508.35</td><td>1.03</td><td>-4.404</td><td>0.0006</td><td>0.1550</td><td>0.0002</td></tr><tr><td>2459255.554371</td><td>5.38</td><td>1.11</td><td>11349.79</td><td>1.11</td><td>-4.416</td><td>0.0007</td><td>0.1576</td><td>0.0002</td></tr><tr><td>2459256.531086</td><td>-15.02</td><td>1.39</td><td>11458.13</td><td>1.39</td><td>-4.415</td><td>0.0009</td><td>0.1523</td><td>0.0002</td></tr><tr><td>2459258.595699</td><td>-28.94</td><td>1.39</td><td>11467.53</td><td>1.39</td><td>-4.390</td><td>0.0009</td><td>0.1555</td><td>0.0002</td></tr><tr><td>2459272.584028</td><td>14.92</td><td>1.41</td><td>11336.17</td><td>1.41</td><td>-4.396</td><td>0.0009</td><td>0.1565</td><td>0.0002</td></tr><tr><td>2459275.463016</td><td>-70.04</td><td>1.54</td><td>11467.22</td><td>1.54</td><td>-4.404</td><td>0.0010</td><td>0.1559</td><td>0.0003</td></tr><tr><td>2459276.476808</td><td>48.88</td><td>1.24</td><td>11405.41</td><td>1.24</td><td>-4.405</td><td>0.0008</td><td>0.1570</td><td>0.0003</td></tr><tr><td>2459277.498226</td><td>-23.93</td><td>1.20</td><td>11436.29</td><td>1.20</td><td>-4.395</td><td>0.0008</td><td>0.1575</td><td>0.0002</td></tr><tr><td>2459288.524057</td><td>-18.01</td><td>1.86</td><td>11415.00</td><td>1.86</td><td>-4.412</td><td>0.0013</td><td>0.1584</td><td>0.0003</td></tr><tr><td>2459290.355022</td><td>-9.75</td><td>1.07</td><td>11456.17</td><td>1.07</td><td>-4.402</td><td>0.0007</td><td>0.1568</td><td>0.0002</td></tr><tr><td>2459293.348057</td><td>26.70</td><td>1.04</td><td>11430.86</td><td>1.04</td><td>-4.410</td><td>0.0007</td><td>0.1604</td><td>0.0002</td></tr><tr><td>2459299.513734</td><td>19.96</td><td>1.73</td><td>11442.26</td><td>1.73</td><td>-4.396</td><td>0.0011</td><td>0.1554</td><td>0.0002</td></tr><tr><td>2459303.465955</td><td>-8.15</td><td>1.70</td><td>11439.02</td><td>1.70</td><td>-4.384</td><td>0.0011</td><td>0.1578</td><td>0.0002</td></tr><tr><td>2459304.439774</td><td>-4.83</td><td>2.02</td><td>11393.90</td><td>2.02</td><td>-4.398</td><td>0.0014</td><td>0.1567</td><td>0.0003</td></tr><tr><td>2459305.455985</td><td>10.20</td><td>2.34</td><td>11448.21</td><td>2.34</td><td>-4.397</td><td>0.0016</td><td>0.1556</td><td>0.0003</td></tr><tr><td>2459307.428793</td><td>-10.12</td><td>1.40</td><td>11195.83</td><td>1.40</td><td>-4.412</td><td>0.0007</td><td>0.1561</td><td>0.0002</td></tr><tr><td>2459327.354815</td><td>57.11</td><td>1.14</td><td>11404.00</td><td>1.14</td><td>-4.427</td><td>0.0008</td><td>0.1560</td><td>0.0002</td></tr><tr><td>2459328.357969</td><td>-63.05</td><td>1.45</td><td>11394.44</td><td>1.45</td><td>-4.416</td><td>0.0010</td><td>0.1564</td><td>0.0002</td></tr><tr><td>2459330.442524</td><td>-26.86</td><td>1.86</td><td>11408.83</td><td>1.86</td><td>-4.393</td><td>0.0012</td><td>0.1581</td><td>0.0002</td></tr><tr><td>2459338.380324</td><td>44.82</td><td>1.49</td><td>11362.00</td><td>1.49</td><td>-4.419</td><td>0.0010</td><td>0.1499</td><td>0.0002</td></tr><tr><td>2459343.387805</td><td>27.82</td><td>3.35</td><td>11440.00</td><td>3.35</td><td>-4.388</td><td>0.0025</td><td>0.1555</td><td>0.0004</td></tr><tr><td>2459344.386754</td><td>8.38</td><td>2.11</td><td>11237.21</td><td>2.11</td><td>-4.398</td><td>0.0012</td><td>0.1530</td><td>0.0002</td></tr><tr><td>2459345.409001</td><td>-9.05</td><td>2.25</td><td>11377.46</td><td>2.25</td><td>-4.412</td><td>0.0016</td><td>0.1516</td><td>0.0003</td></tr><tr><td>2459346.388649</td><td>78.09</td><td>1.99</td><td>11350.37</td><td>1.99</td><td>-4.433</td><td>0.0014</td><td>0.1527</td><td>0.0002</td></tr><tr><td>2459352.368704</td><td>59.49</td><td>3.23</td><td>11448.40</td><td>3.23</td><td>-4.434</td><td>0.0026</td><td>0.1474</td><td>0.0004</td></tr><tr><td>2459358.388673</td><td>-11.48</td><td>2.21</td><td>11429.58</td><td>2.21</td><td>-4.418</td><td>0.0016</td><td>0.1512</td><td>0.0003</td></tr><tr><td>2459541.678592</td><td>-37.56</td><td>1.23</td><td>11464.09</td><td>1.23</td><td>-4.397</td><td>0.0009</td><td>0.1542</td><td>0.0002</td></tr><tr><td>2459560.749308</td><td>-87.49</td><td>1.23</td><td>11443.54</td><td>1.23</td><td>-4.386</td><td>0.0008</td><td>0.1520</td><td>0.0002</td></tr></tbody></table>
<table><caption>Table D.2. continued.</caption><thead><tr><th>Time [BJD]</th><th>BIS<sub>DRS</sub> [m s<sup>-1</sup>]</th><th>σ<sub>BIS</sub> [m s<sup>-1</sup>]</th><th>FWHM<sub>DRS</sub> [m s<sup>-1</sup>]</th><th>σ<sub>FWHM</sub> [m s<sup>-1</sup>]</th><th>log R'<sub>HK</sub> [dex]</th><th>σ<sub>log R'<sub>HK</sub></sub> [dex]</th><th>Hα</th><th>σ<sub>Hα</sub></th></tr></thead><tbody><tr><td>2459565.539049</td><td>-12.47</td><td>1.35</td><td>11401.77</td><td>1.35</td><td>-4.402</td><td>0.0010</td><td>0.1502</td><td>0.0002</td></tr><tr><td>2459589.663616</td><td>7.66</td><td>1.90</td><td>11457.76</td><td>1.90</td><td>-4.380</td><td>0.0015</td><td>0.1534</td><td>0.0003</td></tr><tr><td>2459601.673776</td><td>-10.02</td><td>1.53</td><td>11457.29</td><td>1.53</td><td>-4.368</td><td>0.0010</td><td>0.1527</td><td>0.0002</td></tr><tr><td>2459614.563832</td><td>-9.15</td><td>2.45</td><td>11441.36</td><td>2.45</td><td>-4.382</td><td>0.0020</td><td>0.1585</td><td>0.0005</td></tr><tr><td>2459616.597751</td><td>26.39</td><td>2.32</td><td>11378.78</td><td>2.32</td><td>-4.379</td><td>0.0018</td><td>0.1519</td><td>0.0003</td></tr><tr><td>2459617.395140</td><td>-54.23</td><td>3.30</td><td>11417.37</td><td>3.30</td><td>-4.385</td><td>0.0029</td><td>0.1511</td><td>0.0004</td></tr><tr><td>2459618.567952</td><td>34.20</td><td>3.84</td><td>11400.53</td><td>3.84</td><td>-4.393</td><td>0.0037</td><td>0.1507</td><td>0.0006</td></tr><tr><td>2459619.461850</td><td>16.30</td><td>2.08</td><td>11342.25</td><td>2.08</td><td>-4.391</td><td>0.0016</td><td>0.1498</td><td>0.0003</td></tr><tr><td>2459620.613258</td><td>-22.19</td><td>5.38</td><td>11401.43</td><td>5.38</td><td>-4.392</td><td>0.0059</td><td>0.1512</td><td>0.0007</td></tr><tr><td>2459628.408285</td><td>14.19</td><td>1.68</td><td>11467.33</td><td>1.68</td><td>-4.376</td><td>0.0012</td><td>0.1579</td><td>0.0003</td></tr><tr><td>2459629.419883</td><td>28.93</td><td>1.40</td><td>11412.13</td><td>1.40</td><td>-4.395</td><td>0.0010</td><td>0.1496</td><td>0.0003</td></tr><tr><td>2459642.494913</td><td>-37.57</td><td>1.51</td><td>11360.45</td><td>1.51</td><td>-4.407</td><td>0.0011</td><td>0.1458</td><td>0.0002</td></tr><tr><td>2459644.571088</td><td>30.47</td><td>1.55</td><td>11436.43</td><td>1.55</td><td>-4.415</td><td>0.0011</td><td>0.1434</td><td>0.0002</td></tr><tr><td>2459646.370440</td><td>-10.34</td><td>1.16</td><td>11371.16</td><td>1.16</td><td>-4.404</td><td>0.0008</td><td>0.1465</td><td>0.0002</td></tr><tr><td>2459647.458259</td><td>-24.90</td><td>1.19</td><td>11492.72</td><td>1.19</td><td>-4.388</td><td>0.0008</td><td>0.1504</td><td>0.0002</td></tr><tr><td>2459648.542710</td><td>40.87</td><td>1.21</td><td>11317.05</td><td>1.21</td><td>-4.407</td><td>0.0009</td><td>0.1509</td><td>0.0002</td></tr><tr><td>2459655.483703</td><td>-26.49</td><td>2.77</td><td>11155.91</td><td>2.77</td><td>-4.400</td><td>0.0022</td><td>0.1518</td><td>0.0004</td></tr><tr><td>2459656.538090</td><td>11.97</td><td>3.11</td><td>11185.47</td><td>3.11</td><td>-4.397</td><td>0.0027</td><td>0.1473</td><td>0.0004</td></tr><tr><td>2459658.467214</td><td>-14.85</td><td>3.10</td><td>11110.32</td><td>3.10</td><td>-4.430</td><td>0.0027</td><td>0.1416</td><td>0.0004</td></tr><tr><td>2459659.454730</td><td>-37.78</td><td>1.10</td><td>11166.31</td><td>1.10</td><td>-4.411</td><td>0.0007</td><td>0.1516</td><td>0.0002</td></tr><tr><td>2459660.465796</td><td>12.08</td><td>3.70</td><td>11282.75</td><td>3.70</td><td>-4.385</td><td>0.0029</td><td>0.1511</td><td>0.0006</td></tr><tr><td>2459661.505656</td><td>49.54</td><td>2.59</td><td>11176.03</td><td>2.59</td><td>-4.404</td><td>0.0021</td><td>0.1471</td><td>0.0004</td></tr><tr><td>2459662.478287</td><td>5.25</td><td>2.49</td><td>11169.92</td><td>2.49</td><td>-4.407</td><td>0.0019</td><td>0.1447</td><td>0.0003</td></tr><tr><td>2459682.409509</td><td>-33.98</td><td>1.32</td><td>11252.66</td><td>1.32</td><td>-4.390</td><td>0.0008</td><td>0.1521</td><td>0.0002</td></tr><tr><td>2459683.445090</td><td>78.31</td><td>2.14</td><td>11135.60</td><td>2.14</td><td>-4.410</td><td>0.0016</td><td>0.1491</td><td>0.0003</td></tr><tr><td>2459684.444817</td><td>-20.93</td><td>2.70</td><td>11158.18</td><td>2.70</td><td>-4.422</td><td>0.0022</td><td>0.1469</td><td>0.0004</td></tr><tr><td>2459685.465805</td><td>-20.70</td><td>4.52</td><td>11169.35</td><td>4.52</td><td>-4.415</td><td>0.0046</td><td>0.1426</td><td>0.0006</td></tr><tr><td>2459686.470475</td><td>34.51</td><td>1.82</td><td>11199.68</td><td>1.82</td><td>-4.402</td><td>0.0013</td><td>0.1501</td><td>0.0003</td></tr><tr><td>2459687.367088</td><td>-11.12</td><td>5.03</td><td>11179.72</td><td>5.03</td><td>-4.425</td><td>0.0052</td><td>0.1491</td><td>0.0007</td></tr></tbody></table>