Optical and near-infrared stellar activity
characterization of the early M dwarf Gl 205 with
SOPHIE and SPIRou

P. Cortés-Zuleta, I. Boisse, B. Klein, E. Martioli, P. Cristofari, A.

Antoniadis-Karnavas, J.-F. Donati, X. Delfosse, C. Cadieux, N. Heidari, et al.

► To cite this version:

P. Cortés-Zuleta, I. Boisse, B. Klein, E. Martioli, P. Cristofari, et al. Optical and near-infrared stellar activity characterization of the early M dwarf Gl 205 with SOPHIE and SPIRou. Astronomy & Astrophysics - A&A, 2023, 673, pp.A14. 10.1051/0004-6361/202245131. hal-04252193

HAL Id: hal-04252193

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

Submitted on 26 Oct 2023

**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 CC BY 4.0 - Attribution - International License
# Optical and near-infrared stellar activity characterization of the early M dwarf Gl 205 with SOPHIE and SPIRou

P. Cortés-Zuleta<sup>1,2</sup> ◦, I. Boisse<sup>1</sup>, B. Klein<sup>3</sup> ◦, E. Martioli<sup>4,12</sup> ◦, P. I. Cristofari<sup>5</sup> ◦, A. Antoniadis-Karnavas<sup>6,7</sup>, J.-F. Donati<sup>5</sup> ◦, X. Delfosse<sup>8</sup>, C. Cadieux<sup>9</sup> ◦, N. Heidari<sup>15,16,1</sup> ◦, É. Artigau<sup>9,13</sup> ◦, S. Bellotti<sup>5,11</sup> ◦, X. Bonfils<sup>8</sup>, A. Carmona<sup>8</sup> ◦, N. J. Cook<sup>9</sup> ◦, R. F. Díaz<sup>10</sup> ◦, R. Doyon<sup>9,13</sup>, P. Fouqué<sup>5</sup> ◦, C. Moutou<sup>5</sup> ◦, P. Petit<sup>5</sup> ◦, T. Vandal<sup>9</sup>, L. Acuña<sup>1</sup> ◦, L. Arnold<sup>19</sup> ◦, N. Astudillo-Defru<sup>22</sup> ◦, V. Bourrier<sup>20</sup>, F. Bouchy<sup>20</sup>, R. Cloutier<sup>17</sup> ◦, S. Dalal<sup>23</sup> ◦, M. Deleuil<sup>1</sup>, O. D. S. Demangeon<sup>6,7</sup>, X. Dumusque<sup>20</sup> ◦, T. Forveille<sup>8</sup> ◦, J. Gomes da Silva<sup>6</sup> ◦, N. Hara<sup>20</sup>, G. Hébrard<sup>12,14</sup>, S. Hoyer<sup>1</sup> ◦, G. Hussain<sup>11</sup>, F. Kiefer<sup>12,21</sup> ◦, J. Morin<sup>18</sup> ◦, A. Santerne<sup>1</sup> ◦, N. C. Santos<sup>6,7</sup>, D. Segransan<sup>20</sup>, M. Stalport<sup>20</sup>, and S. Udry<sup>20</sup> ◦

*(Affiliations can be found after the references)*

Received 4 October 2022 / Accepted 19 January 2023

## ABSTRACT

**Context.** The stellar activity of M dwarfs is the main limiting factor in the discovery and characterization of the exoplanets orbiting them, because it induces quasi-periodic radial velocity (RV) variations.

**Aims.** We aim to characterize the magnetic field and stellar activity of the early, moderately active M dwarf Gl 205 in the optical and near-infrared (NIR) domains.

**Methods.** We obtained high-precision quasi-simultaneous spectra in the optical and NIR with the SOPHIE spectrograph and SPIRou spectropolarimeter between 2019 and 2022. We computed the RVs from both instruments and the SPIRou Stokes V profiles. We used Zeeman-Doppler imaging (ZDI) to map the large-scale magnetic field over the time span of the observations. We studied the temporal behavior of optical and NIR RVs and activity indicators with the Lomb-Scargle periodogram and a quasi-periodic Gaussian process regression (GPR). In the NIR, we studied the equivalent width of Al I, Ti I, K I, Fe I, and He I. We modeled the activity-induced RV jitter using a multi-dimensional GPR with activity indicators as ancillary time series.

**Results.** The optical and NIR RVs show similar scatter but NIR shows a more complex temporal evolution. We observe an evolution of the magnetic field topology from a poloidal dipolar field in 2019 to a dominantly toroidal field in 2022. We measured a stellar rotation period of $P_{\text{rot}} = 34.4 \pm 0.5$ days in the longitudinal magnetic field. Using ZDI, we measure the amount of latitudinal differential rotation (DR) shearing the stellar surface, yielding rotation periods of $P_{\text{eq}} = 32.0 \pm 1.8$ days at the stellar equator and $P_{\text{pol}} = 45.5 \pm 0.3$ days at the poles. We observed inconsistencies in the periodicities of the activity indicators that could be explained by these DR values. The multi-dimensional GP modeling yields an RMS of the RV residuals down to the noise level of $3 \text{ m s}^{-1}$ for both instruments while using $H\alpha$ and the BIS in the optical and the full width at half maximum (FWHM) in the NIR as ancillary time series.

**Conclusions.** The RV variations observed in Gl 205 are due to stellar activity, with a complex evolution and different expressions in the optical and NIR revealed thanks to an extensive follow-up. Spectropolarimetry remains the best technique to constrain the stellar rotation period over standard activity indicators, particularly for moderately active M dwarfs.

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

## 1. Introduction

M dwarfs are the most abundant stars in the Milky Way (Henry et al. 2006; Reylé et al. 2021) and have become key targets for exoplanetary surveys (e.g., Morley et al. 2017). Their low mass favors the detection of the planets orbiting them as the Doppler variations are larger than those seen for Sun-like stars for a given planetary mass and equilibrium temperature. Dedicated surveys to discover and characterize exoplanets around M dwarfs are currently being carried out using transit photometry (e.g., MEarth: Nutzman & Charbonneau 2008; TRAPPIST: Gillon et al. 2011; SPECULOOS: Delrez et al. 2018; Tierras: Garcia-Mejia et al. 2020) and Doppler spectroscopy (e.g., HPF: Mahadevan et al. 2012; HARPS: Bonfils et al. 2013; HARPS-N: Covino et al. 2013; IRD: Kotani et al. 2014; CARMENES: Quirrenbach et al. 2018; MAROON-X: Seifahrt et al. 2018; SOPHIE: Hobson et al. 2018; SPIRou: Donati et al. 2020).

Transit and radial velocity (RV) surveys have shown that M dwarfs can host Earth-like planets (Bonfils et al. 2013; Kopparapu 2013; Dressing & Charbonneau 2013, 2015; Sabotta et al. 2021; Pinamonti et al. 2022). However, despite efforts to detect and characterize Earth-like planets around M dwarfs, the presence of stellar activity remains a significant limiting factor in such endeavors. Cool stars, such as M dwarfs, are known to host magnetic fields that cause the phenomena known as stellar activity (Reiners 2012).

Stellar magnetic activity produces quasi-periodic RV signals with amplitudes of a few meters per second that can easily lead to a false-positive exoplanet detection (e.g., Queloz et al. 2001; Desidera et al. 2004; Huélamo et al. 2008; Carolo et al. 2014; Bortle et al. 2021). Features over the stellar surface, particularly dark spots and bright plages, can survive and evolve on timescales of a few stellar rotation cycles and can therefore
modulate the RVs at the period of the stellar rotation (Boisse
et al. 2011; Scandariato et al. 2017).

Several techniques have been developed during the last
decade to mitigate the effects of stellar activity in order to
improve the detection of exoplanets around active stars. When
the stellar rotation period is known, one can model the RV vari-
ations induced by spots by fitting sinusoidal signals and their
harmonics (Boisse et al. 2011) or by simulations of the active
regions (Dumusque et al. 2014a). Another approach is to use
activity indicators built from spectral information. These can be
proxies of the shape of the cross-correlation function (CCF), for
example the full width at half maximum (FWHM) and the bisec-
tor inverse slope (BIS; Queloz et al. 2001; Boisse et al. 2011), or
of the chromospheric emission in spectral lines sensitive to activ-
ity, such as the Ca II H&K and Hα (Boisse et al. 2009). A third
approach, among others, is to model the stellar surface structures
generating RV variations (e.g., Hébrard et al. 2016; Klein et al.
2021, 2022).

The activity-induced behavior of the spectral lines in M dwarf spectra has been studied only for a handful of lines, such as Hα, Ca II H&K, Ca II IRT, Na I D, and He I, with the Hα and CaII H&K lines the most extensively explored (e.g., Gomes da Silva et al. 2011; Newton et al. 2017; Fuhrmeister et al. 2019; Lafarga et al. 2021). The majority of the well-known activity tracers are located at optical wavelengths. Recently, some near-infrared (NIR) lines were studied, for example the He I triplet (Schöfer et al. 2019; Fuhrmeister et al. 2019, 2020) and the K I line (Fuhrmeister et al. 2022; Terrien et al. 2022). The He I triplet in absorption is not detected over the range of the whole M dwarf spectral type and disappears for stars later than M5. When in emission, it could be related to flaring events (Fuhrmeister et al. 2019). Moreover, Fuhrmeister et al. (2020) found that the variability in the He I triplet may be correlated with the Hα variability only for active M dwarfs. In the particular case of the active M dwarf AU Mic, the He I flux correlated well with the RVs (Klein et al. 2021). Regarding the K I line emission, it has been found that it is rarely correlated or anti-correlated with Hα (Fuhrmeister et al. 2022). Moreover, Terrien et al. (2022) found clear signals of Zeeman broadening in this line that are modulated by the stellar rotation in Gl 699 and Teegarden's star data.

More recently, the use of data-driven techniques, namely Gaussian processes regression (GPR; e.g., Haywood et al. 2014; Rajpaul et al. 2015; Jones et al. 2017; Gilbertson et al. 2020; Klein et al. 2021; Barragán et al. 2022a; Delisle et al. 2022) or principal component analysis (PCA; e.g., Davis et al. 2017; Cretignier et al. 2022), has shown good results in modeling the contribution of stellar activity in RVs. Nowadays, the use of GPs has become the standard procedure for modeling and filtering stellar activity in RV and transits exoplanet searches, because the GPR can easily fit the quasi-periodic activity signal. This technique is particularly successful when the information from activity indicators is used when fitting the GPs on the RV time series along with the Keplerian signal (Rajpaul et al. 2015; Barragán et al. 2019a; Suárez Mascareño et al. 2020; Faria et al. 2022; Zicher et al. 2022).

Spectropolarimetry data have been widely used to measure
and constrain the properties of the large-scale magnetic field of
M dwarfs (Donati et al. 2008; Morin et al. 2008, 2010). In partic-
ular, it has been shown that the longitudinal magnetic field Bℓ is
a reliable magnetic activity tracer for M dwarfs and therefore can
be used to determine the stellar rotation period (e.g., Morin et al.
2008, 2010; Folsom et al. 2016; Hébrard et al. 2016; Martioli
et al. 2022).

Another way to mitigate stellar activity is to observe in NIR
wavelengths, because for M dwarfs, the spot-induced RV jit-
ter decreases as a function of wavelength (Martín et al. 2006;
Desort et al. 2007; Reiners et al. 2010; Mahmud et al. 2011). The
flux contrast between the dark spots and the stellar surface is
smaller at NIR wavelengths than in the optical, generating less
activity-induced RV jitter. However, this effect strongly depends
on the spectral type and the spot configuration on the stellar
surface (Reiners et al. 2010; Andersen & Korhonen 2015). On
the contrary, the RV jitter due to the Zeeman effect increases at
longer wavelengths as the induced variation is proportional to
the wavelength and the magnetic field strength (Hébrard et al.
2014). In cases of low-temperature contrast, the expected gain of
the NIR is hampered by the increasing Zeeman effect (Reiners
et al. 2013; Klein & Donati 2020). Simultaneous observations at
the optical and NIR can help to disentangle the activity contri-
bution from the two wavelength domains (Reiners et al. 2013;
Robertson et al. 2020). Moreover, multiwavelength observations
have been crucial to reject the planetary nature of RV signals
(e.g., TW Hydrae: Huélamo et al. 2008; AD Leo: Carleo et al.
2020; Carmona et al. 2023)

Gl 205 is an early (M1.5), relatively nearby (5.7 pc), slow-
rotating M1.5 dwarf with moderate levels of activity. It has
a mass of 0.55 ± 0.03 M⊙ and a radius of 0.56 ± 0.03 R⊙
(Schweitzer et al. 2019). More stellar parameters are listed in
Table 1. Previously, Gl 205 was monitored with HARPS (Bonfils
et al. 2013), and with HARPS-Pol and NARVAL (Hébrard et al.
2016). Using the Hα and Ca II H&K indices, Bonfils et al. (2013)
found the stellar rotational period to be close to 33 days, and
using spectropolarimetric data, Hébrard et al. (2016) determined
a rotation period of P<sub>rot</sub> = 33.63 ± 0.37 days. Photometry data in
the V band showed a periodic signal of 33 d and a possible long-
trend magnetic cycle of ~1500 days (Kiraga & Stepien 2007).
The discrepancies between the periodicities of several activity
indicators could be explained by differential rotation shearing
the stellar surface and this may induce a difference of ~10 days
in the rotation period between the equator and the pole (Hébrard
et al. 2016).

In this work, we analyze the magnetic field and stellar
activity of the early-M dwarf Gl 205, which was intensively
monitored over a two-year period with the SOPHIE optical spec-
trograph and with the SPIRou NIR spectropolarimeter. This
paper is structured as follows. In Sect. 2, we describe the obser-
vations and the reduction of the SOPHIE and SPIRou data.
Section 3 is dedicated to the stellar characterization of Gl 205
and Sect. 4 to the analysis of TESS photometry. We describe the
magnetic field properties using the SPIRou spectropolarimetric
data in Sect. 5. In Sect. 6 we compare the optical and NIR RVs
and in Sect. 7 we describe and analyze the activity indicators. We
filter the RV variations due to activity using a multi-dimensional
GP framework in Sect. 8. In Sect. 9 we discuss our RV detection
limit of Keplerian signals. In Sect. 10 we discuss our results and
present the principal conclusions of this work.

2. Observations and data reduction

2.1. SOPHIE

SOPHIE is a high-resolution fiber-fed, cross-dispersed échelle
spectrograph mounted on the 1.93m telescope at the Observa-
toire de Haute-Provence (OHP; Perruchot et al. 2008; Bouchy
et al. 2013). It covers a wavelength domain from 3870 to
6940 Å across 39 spectral orders.
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** SOPHIE and SPIRou data of Gl 205. Left: RV time series of Gl 205 from SOPHIE in blue and SPIRou in red with an offset equal to the mean RV of each instrument, RV<sub>SOPHIE</sub> = 8727.8 m s<sup>-1</sup> and RV<sub>SPIRou</sub> = 9345.5 m s<sup>-1</sup>. The scatter in the time series is 4.7 m s<sup>-1</sup> in the SPIRou data set and 4.8 m s<sup>-1</sup> in SOPHIE. The labels above the time series depict the four subsets of observations described in Sect. 6. Right: GLS periodogram of the full RV time series, combining SOPHIE and SPIRou observations. The highest peak is marked with an arrow at 34.3 days.</figcaption></figure>

<table><caption>**Table 1.** Stellar parameters of Gl 205.</caption><thead><tr><th>Parameter</th><th>Units</th><th>Value</th><th>References</th></tr></thead><tbody><tr><td>RA</td><td>h m s</td><td>05 31 27.395785</td><td>1</td></tr><tr><td>Dec</td><td>d m s</td><td>-03 40 38.024004</td><td>1</td></tr><tr><td>Parallax</td><td>mas</td><td>175.31 ± 0.02</td><td>1</td></tr><tr><td>Distance</td><td>pc</td><td>5.7040 ± 0.0006</td><td>1</td></tr><tr><td>Spectral type</td><td></td><td>M1.5V</td><td>2</td></tr><tr><td>V</td><td>mag</td><td>7.968</td><td>3</td></tr><tr><td>B</td><td>mag</td><td>9.443</td><td>3</td></tr><tr><td>J</td><td>mag</td><td>4.83</td><td>4</td></tr><tr><td>K</td><td>mag</td><td>3.90</td><td>4</td></tr><tr><td>Mass</td><td>M<sub>☉</sub></td><td>0.549 ± 0.029</td><td>5</td></tr><tr><td>Radius</td><td>R<sub>☉</sub></td><td>0.556 ± 0.033</td><td>5</td></tr><tr><td>L<sub>☉</sub></td><td>L<sub>☉</sub></td><td>0.061 ± 0.006</td><td>6</td></tr><tr><td>Inclination</td><td>°</td><td>60 ± 10</td><td>7</td></tr><tr><td>v sin i</td><td>km s<sup>-1</sup></td><td>0.7 ± 0.1</td><td>This work</td></tr><tr><td>P<sub>rot</sub></td><td>days</td><td>34.4 ± 0.5</td><td>This work</td></tr><tr><td>log g</td><td>dex</td><td>4.70 ± 0.1</td><td>This work <sup>(a)</sup></td></tr><tr><td>T<sub>eff</sub>, SOPHIE</td><td>K</td><td>3878 ± 81</td><td>This work <sup>(b)</sup></td></tr><tr><td>T<sub>eff</sub>, SPIRou</td><td>K</td><td>3770 ± 30</td><td>This work <sup>(a)</sup></td></tr><tr><td>[Fe/H]</td><td>dex</td><td>0.21 ± 0.06</td><td>This work <sup>(b)</sup></td></tr><tr><td>[M/H]</td><td>dex</td><td>0.43 ± 0.1</td><td>This work <sup>(a)</sup></td></tr><tr><td>[α/Fe]</td><td>dex</td><td>-0.08 ± 0.04</td><td>This work <sup>(a)</sup></td></tr><tr><td>Age</td><td>Gyr</td><td>3.6<sup>+5.0</sup><sub>-1.7</sub></td><td>This work <sup>(a)</sup></td></tr></tbody></table>

**Notes.** <sup>(a)</sup>Derived from SPIRou spectra. <sup>(b)</sup>Derived from SOPHIE spectra.

**References.** (1) Gaia Collaboration (2020); (2) Kiraga & Stepien (2007); (3) Koen et al. (2010); (4) Ducati (2002); (5) Schweitzer et al. (2019); (6) Baraffe et al. (2015); (7) Hébrard et al. (2016).

Observations of Gl 205 were carried out between November 2019 and January 2022 as a part of the Subprogram 3 (SP3) of the SOPHIE exoplanet consortium, which is dedicated to hunting for exoplanets around M dwarfs. So far, the main results of this program include the detection of the exoplanets Gl 96b (Hobson et al. 2018), Gl 378b (Hobson et al. 2019), and Gl 411b (Díaz et al. 2019), and the independent confirmation of Gl 617Ab (Hobson et al. 2018).

The observations were gathered using the high-resolution (HR) mode, in which the spectrograph reaches a resolving power of $\lambda/\Delta\lambda \sim 75\,000$. In order to measure the instrumental drift, simultaneous calibrations with a Fabry-Pérot (FP) étalon were performed. In total, we gathered 74 spectra of the star, with a median exposure time of 930 s and a median signal-to-noise ratio (S/N) per pixel at 550 nm of 95. After removing the observations with airmass >1.6, S/N < 80, or those affected by moonlight pollution, a total of 62 observations remained.

We used the SOPHIE Data Reduction Software (DRS, Bouchy et al. 2009a) to reduce and extract the spectra. The SOPHIE observations are corrected for the charge transfer inefficiency (CTI) effect of the CCD following Bouchy et al. (2009a) and Hobson et al. (2018). The DRS automatically computes the RVs using the CCF technique, which is obtained by cross-correlation of the spectra with an empirical M2V mask. The DRS also uses the CCF to deliver stellar activity indicators, such as the depth of the CCF defined as CCF contrast, CCF FWHM, and BIS.

However, in the case of M dwarf spectra, the CCF method is not an optimal approach because of the large number of absorption lines. To use most of the Doppler information available, we used a template-matching algorithm to obtain high-precision RVs (NAIRA, Astudillo-Defru et al. 2015, 2017). First, all the available spectra are normalized by the blaze functions and scaled to unity by their median. Second, the spectra are shifted using the RVs computed by the DRS. Third, these spectra are co-added to build a high-S/N stellar template. Finally, the maximum-likelihood RV is the minimum of the Chi-square profile obtained by shifting the stellar template over an array of RVs.

A list of standard stars was used to correct the long-term variations of the RV zero-point in SOPHIE, an effect described by Courcol et al. (2015). To do so, we systematically monitored a group of G-type stars: HD 185144, HD 9407, HD 89269A, and a group of M dwarfs: Gl 514, Gl 15A, and Gl 686 (Hobson et al. 2018) on each night when this was possible.

The final RVs from NAIRA used in this work with their error bars are listed in the Table A.1, along with the activity indicators: CCF FWHM, CCF contrast, BIS, the S index, and $H\alpha$ line (see Sect. 7). Figure 1 shows the SOPHIE RV time series whose average error bars are 1.9 m s<sup>-1</sup> and the scatter is 4.8 m s<sup>-1</sup>.
## 2.2. SPIRou

The Spectro-Polarimetre InfraRouge (SPIRou, Donati et al. 2020) is a high-resolution NIR spectropolarimeter and velocimeter mounted at the Canada-France-Hawaii Telescope (CFHT) in Hawaii. With a nominal spectral range from 9800 to 23 500 Å, it covers the Y, J, H, and K bands of the infrared spectrum at a spectral resolving power of $\lambda/\Delta\lambda \sim 70\,000$. The observations of Gl 205 are part of the Planet Search program (WP1) of the SPIRou Legacy Survey (SLS, Donati et al. 2020), whose main goal is to perform a systematic RV monitoring of nearby M dwarfs.

SPIRou operates as a spectropolarimeter as well as a spectrograph. A spectropolarimetric sequence consists of four subexposures, each one with a different rotation angle of the half-wave Fresnel rhombs in the polarimeter. We used the spectropolarimetric mode in order to obtain a set of Stokes I (unpolarized) and Stokes V (circularly polarized) profiles spectra per spectropolarimetric sequence of four subexposures (see Sect. 5).

The star was observed from September 2019 to January 2022, collecting a total of 156 sequences of four subexposures. The median exposure time of the observations is 61 s per rhomb position (244 s for a complete sequence) and the median S/N per pixel at 1670 nm is 290. We removed 11 sequences from the analysis, each with $S/N < 150$ or airmass $>1.7$.

The data were reduced using the SPIRou data-reduction software APERO<fnref n="1" /> v0.7.194 (Cook et al. 2022). APERO performs an automatic reduction of the 4096×4096 pixels raw images, including correction for detector effects, and identification and removal of bad pixels and cosmic rays. The data are calibrated by performing flat and blaze corrections, and are then optimally extracted (Horne 1986) from both science channels (fibers A and B that carry orthogonal polarimetric states of the incoming light) and from the calibration channel (fiber C).

The pixel-to-wavelength calibration is done using an UNe hollow cathode lamp and a Fabry-Pérot etalon, following Hobson et al. (2021), in order to obtain the wavelength at the observatory rest frame. APERO then utilizes the barycorrpy<fnref n="2" /> (Kanodia & Wright 2018) Python code to compute the Barycentric Earth RV (BERV) and barycentric Julian date (BJD) of each exposure. APERO performs a telluric and night emission correction in two steps. The first step is to obtain an atmospheric absorption model built with TAPAS (Transmissions of the Atmosphere for Astronomical data; Bertaux et al. 2014), which is applied to pre-clean the science frames. This model only leaves percentage-level residuals in deep (>50%) lines of H₂O and dry absorption molecules (e.g., CH₄, O₂, CO₂, N₂O, and O₃). This procedure of building a TAPAS model is also done for a set of rapid rotating hot stars observed under different atmospheric conditions in order to build a library of telluric residual models. This grid of telluric residual models has three degrees of freedom (optical depths of H₂O and dry components, and a constant). The second step is to subtract this telluric residual model from the pre-cleaned data to obtain telluric-corrected spectra.

In a standard procedure, APERO automatically computes the RVs by cross-correlation of the telluric-corrected spectra with a given binary mask of stellar absorption lines. However, in this work, we made use of the RVs computed by the line-by-line (LBL) method based on the Bouchy et al. (2001) framework and optimized for SPIRou data. The LBL method is fully described by Artigau et al. (2022). This algorithm exploits the RV content

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** TESS light curves of Gl 205. Top: sector 6 data. Middle: zoom onto the region of the two flares identified in Sector 6 using the stella package. The red line depicts the best-fit flare model. Bottom: sector 32 data.</figcaption></figure>

per line on the spectra to obtain one single RV measurement. Usually, an M dwarf observed by SPIRou will have about 16 000 individual spectral lines. A finite-mixture model approach deals with the high-sigma outliers in the RVs of the lines that come from cosmic rays or errors in the correction for tellurics. After outlier removal, the final RV is the mean of a Gaussian distribution containing the individual RVs of all the lines, and its uncertainty is derived following Bouchy et al. (2001). This method is also applied to the simultaneous calibration fiber to correct for the instrumental drift. The LBL RVs are corrected for the instrumental drift and for the long-term zero point using a GPR with data of the most observed stars in the SPIRou Legacy Survey. The details of this procedure will be described in Vandal et al. (in prep). The log of the observations and the RV measurements from the LBL method are listed in Table A.2. Figure 1 shows the SPIRou RV time series whose average error bar is 1.9 m s⁻<fnref n="1" /> and scatter is 4.4 m s⁻<fnref n="1" />.

## 2.3. Transiting Exoplanet Survey Satellite

Currently, the Transiting Exoplanet Survey Satellite (TESS; Ricker 2014) is in its extended mission, after successfully completing an all-sky survey of bright stars during its primary mission. TESS observed Gl 205 with two-minute cadence during Sector 6, between December 12, 2018, and January 6, 2019, and again in Sector 32, from November 19 to December 16, 2020. We obtained the Presearch Data Conditioning (PDC) flux time series processed by the TESS Science Processing Operations Center (SPOC) from the Mikulski Archive for Space Telescopes (MAST)<fnref n="3" />. The light curves are shown in Fig. 2. We used the quality flags given by the pipeline to remove bad regions of the

<footnote n="1">¹ https://github.com/njcuk9999/apero-drs</footnote>

<footnote n="2">² https://github.com/shbhuk/barycorrpy</footnote>

<footnote n="3">³ mast.stsci.edu</footnote>
light curves, and therefore only kept data points with a quality flag equal to zero.

## 3. Stellar characterization

We use the high-resolution spectra from SOPHIE and SPIRou independently to derive the atmospheric stellar parameters of Gl 205. For the SOPHIE optical part, we made use of the ODUSSEAS (Antoniadis-Karnavas et al. 2020) code to compute the effective temperature, $T_{\text{eff}}$, and the metallicity, [Fe/H]. For the SPIRou NIR spectra, we follow Cristofari et al. (2022b) to derive $T_{\text{eff}}$, surface gravity (log g), overall metallicity ([M/H]), and alpha-enhancement ([α/H]). The $T_{\text{eff}}$ measured from the optical and NIR spectra are in good agreement. The values of stellar parameters derived for Gl 205 are listed in Table 1. In this section, we describe both techniques in detail.

### 3.1. SOPHIE spectra

A detailed description of the machine-learning tool ODUSSEAS can be found in Antoniadis-Karnavas et al. (2020). The method is based on measuring the pseudo equivalent widths (pEWs) of absorption lines and blended lines in the range between 5300 Å and 6900 Å. The line list consists of 4104 absorption features, the same as that used by Neves et al. (2014).

ODUSSEAS receives 1D spectra and their resolutions as input. The tool contains a supervised machine-learning algorithm based on the scikit-learn Python package (Pedregosa et al. 2011) in order to determine the $T_{\text{eff}}$ and [Fe/H] of M dwarf stars.

Applied to new spectra, ODUSSEAS measures the pEWs of their lines and compares them to the models generated from the reference HARPS spectra sample, which are convolved to the respective resolution of the new spectra to be analyzed. The reference data set is built with spectra taken from the HARPS M dwarf sample. In the case of the SOPHIE spectrum of the star Gl 205, the HARPS reference spectra are convolved from their original resolution of 115 000 to the SOPHIE resolution of 75 000.

For the analysis of Gl 205, the new reference data set of the upgraded ODUSSEAS version has been applied and includes spectra from 47 M dwarfs. The references for training and testing the models are the pEWs of the 47 HARPS spectra, used together with interferometry-based $T_{\text{eff}}$ (Rabus et al. 2019; Khata et al. 2021) and [Fe/H] derived by applying the method by Neves et al. (2012) using updated values of parallaxes from Gaia DR3.

The resulting stellar parameters of the star are calculated from the mean values of 100 determinations obtained by randomly shuffling and splitting the training (80% of the reference sample, i.e., 37 stars) and testing groups (remaining 20%, i.e., 10 stars) each time. We report parameter uncertainties derived by quadratically adding the dispersion of the resulting stellar parameters and the mean absolute errors of the machine-learning models at this resolution. We obtained an $T_{\text{eff}}$ = 3878 ± 81 K and [Fe/H] = 0.21 ± 0.06 dex.

### 3.2. SPIRou spectra

We estimate the $T_{\text{eff}}$, log g, [M/H], and [α/H] from a high-resolution template spectrum built from the over 500 spectra from subexposures acquired with SPIRou. The process relies on the direct comparison of the template spectrum with a grid of synthetic spectra computed from MARCS model atmospheres (Cristofari et al. 2022b,a). The comparison is performed on

carefully selected spectral windows containing about 20 atomic lines and 40 molecular lines.

Prior to the comparison, synthetic spectra are broadened to account for instrumental effects, and the local continuum of the models is adjusted on windows built around the selected lines. Comparing the template spectrum to a grid of synthetic spectra with various $T_{\text{eff}}$, log g, [M/H], and [α/H] results in the computation of a 3D $\chi^2$ grid on which we fit a 3D second-degree polynomial to retrieve a minimum. With this model, we derive $T_{\text{eff}}$ = 3770 ± 30 K, log g = 4.70 ± 0.1 dex, [M/H] = 0.43 ± 0.1 dex, and [α/Fe] = -0.08 ± 0.04.

### 3.3. Age and galactic population

We computed the age of Gl 205 with the stardate<fnref n="4" /> (Angus et al. 2019) Python package, which combines isochrones fitting with gyrochronology. The inputs of this code are the stellar parameters derived in Sects. 3.2 and 3.1: $T_{\text{eff}}$, log g, [Fe/H], and the parallax from Table 1. As we obtained two estimations of $T_{\text{eff}}$, we tested both values to see if we obtain similar results. Using the $T_{\text{eff}}$ from SOPHIE spectra we measure a stellar age of 3.7$^{+5.0}_{-1.6}$ Gyr. On the other hand, using the SPIRou $T_{\text{eff}}$ we estimate the age at 3.6$^{+5.0}_{-1.7}$ Gyr. Both estimations are in complete agreement.

The rotation period of Gl 205 is higher than that of stars of similar mass in the 4 Gyr-old open cluster M67 sequence (Dungee et al. 2022). Calibrating the gyrochronology relationship with this sequence results in an age determination of 5.2 ± 0.7 (Fouqué et al. 2023), which is within the error bars of our independent estimate.

In order to identify the galactic population (thin disk, thick disk, or halo) to which Gl 205 belongs, we followed Reddy et al. (2006) to obtain the probabilities of belonging to the three populations based on the galactic velocities of the star. For Gl 205, the galactic velocities are (U, V, W) = (49, 13, 31) km s⁻¹, computed using the position, proper motion, and parallax from Gaia DR3 (Gaia Collaboration 2020). We find a 99% probability that Gl 205 belongs to the thin disk of the Milky Way.

The ages of the stars in the thin disk span a wide range. However, it has been stated that most of them have ages of less than 5 Gyr, reaching up to 14 Gyr (Reddy et al. 2006; Haywood 2008; Holmberg et al. 2009). Our estimation of the age of Gl 205 is in agreement with these results. Regarding the metallicity, Allende Prieto et al. (2004) found that the mean metallicity of the stars in the thin disk is $\langle$[Fe/H]$\rangle$ = -0.09 ± 0.19. We estimated a high metallicity for Gl 205 of 0.21 ± 0.06, which is located within $2\sigma$ of the mean value.

### 3.4. Comparison with the literature

We compared our results with previous studies of Gl 205. Schweitzer et al. (2019) derived the photospheric parameters $T_{\text{eff}}$, log g, and [Fe/H] from CARMENES VIS spectra and the stellar radius using the Stefan-Boltzmann's law. Using the log g and the stellar radius, the authors derived the stellar mass; their results of stellar mass and radius are listed in Table 1. These authors found $T_{\text{eff}}$ = 3891 ± 51 K, log g = 4.64 ± 0.07 dex, and [Fe/H] = 0.23 ± 0.16 dex. Our results of [Fe/H] and $T_{\text{eff}}$ computed from SOPHIE spectra and the log g from SPIRou spectra are in particular good agreement within $1\sigma$ with the values of CARMENES. The $T_{\text{eff}}$ derived using SPIRou spectra is within $3\sigma$ of the CARMENES result.

<footnote n="4">⁴ https://github.com/RuthAngus/stardate</footnote>
Neves et al. (2014) obtained [Fe/H] and $T_{eff}$ for Gl 205 from HARPS spectra with values of [Fe/H] = 0.19 ± 0.09 and $T_{eff}$ = 3670 ± 110 K. Our estimation of the metallicity from SOPHIE spectra is in good agreement with the one from HARPS and our $T_{eff}$ value lies within $3\sigma$. However, the HARPS metallicity is in better agreement with our SPIRou $T_{eff}$.

Using spectra of Gl 205 from the APOGEE survey, Souto et al. (2022) determined $T_{eff}$ = 3820 ± 110 K, $log g$ = 4.67 ± 0.2 dex, and [Fe/H] = 0.29 ± 0.10 dex. Their result for $T_{eff}$ is in particularly good agreement with our value derived from SPIRou spectra. Maldonado et al. (2015) used HARPS and HARPS-N spectra to derive stellar parameters of early-M dwarfs, finding a metallicity of Gl 205 of [Fe/H] = -0.03 ± 0.19 dex. Our estimation of [Fe/H] is within $2\sigma$ of their value.

# 4. Light-curve analysis

In this section, we describe the analysis of the TESS light curves of Gl 205 to identify stellar flares and search for exoplanet transits. As each TESS sector has a duration of ~27 days, the rotation period of Gl 205 (see Table 1) is not covered in one single sector and therefore determination of the rotation period from the TESS data is not possible. Moreover, there is a time gap of 684 d between the sector 6 and sector 32.

## 4.1. Flares identification

For the identification of stellar flares, we used the Python package `stella`<fnref n="5" /> (Feinstein et al. 2020a). This open-source code uses convolutional neural networks (CNN) to identify flare events in the two-minute cadence TESS light curves and delivers the probability of such an event for a given time. The reported flare probability is the average prediction of the ten models available in `stella`. After identification, `stella` uses an empirical flare model (Walkowicz et al. 2011; Davenport et al. 2014) to obtain the best-fit parameters through a $\chi^2$ fit. The flare model includes a sharp Gaussian rise and an exponential decay (Feinstein et al. 2020b). Candidate events with probability higher than 50% of being a flare were considered for modeling. In the modeling process, the part of the light curves that includes the flare is detrended to account for stellar variability.

We applied the algorithm in the available TESS sectors of Gl 205 using the available trained CNNs from Feinstein et al. (2020b). We identified two flare events during sector 6 and none during sector 32 (see Fig. 2). Table 2 reports the time of the flare's peak, its amplitude, the equivalent duration (ED) – which is measured as the area of the flare event –, the rise and fall parameters from the flare model, and the probability.

The second flare identified seems more prominent than the first, reaching a higher flux amplitude. However, both events have the same equivalent duration of about 10 hours, meaning that the energies of these events are similar but cover different timescales; while the first flare lasted 2.5 h, the second flare lasted 0.5 h.

Günther et al. (2020) performed a study of stellar flares in the first data release of TESS, finding that mid to late M dwarfs show the highest fraction of flaring stars. However, the authors warned that this may be influenced by the TESS target selection. Only 10% of the early M dwarfs in their sample are flaring stars. Moreover, fast rotators ($P < 5$ days) may have higher flare rates than slow rotators. Our findings in Gl 205 agree with the results

<table><caption>**Table 2.** Parameters of the flare events identified in the TESS light curves of Gl 205.</caption><thead><tr><th>t<sub>peak</sub><br>BJD</th><th>Amplitude<br>(rel. flux)</th><th>ED<br>(h)</th><th>Rise<br>–</th><th>Fall<br>–</th><th>Probability<br>(%)</th></tr></thead><tbody><tr><td>2458471.056</td><td>0.0018</td><td>9.97</td><td>0.0002</td><td>0.025</td><td>67</td></tr><tr><td>2458484.149</td><td>0.0041</td><td>10.01</td><td>0.0002</td><td>0.004</td><td>99</td></tr></tbody></table>

Notes. Time of the flare's peak $t_{peak}$, equivalent duration of the flare event ED, Gaussian rise of the flare, and its exponential decay.

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Detrended TESS light curves of Gl 205. The red line depicts the applied detrending model to remove the stellar variability using the wōtan package. Below each panel are the residuals of the model with a standard deviation of 210 ppm.</figcaption></figure>

of Günther et al. (2020), as we observed a low rate of flaring for this star, which is expected for early and slow-rotating M dwarfs.

## 4.2. Planet transit search

The TESS data-validation report of Gl 205 identifies a planet candidate with an orbital period of 22.2 days. However, the transit-event candidate proposed by the automatic pipeline is clearly affected by the flare of 2458471.056 BJD. We conduct our own analysis in order to identify possible transits in the light curves.

First, we removed outliers using a $3\sigma$ clipping procedure and excluded the data affected by the flares. As the light curves are affected by stellar variability, we used the Python package `wōtan`<fnref n="6" /> (Hippke et al. 2019) to remove the trends. This code includes several methods to perform light-curve detrending. We applied a method based on a time-windowed sliding filter with an iterative robust location estimator following the results of Hippke et al. (2019). For the detrending model, we excluded the edges of the light curves, because these zones are usually affected by strong systematic errors. Figure 3 shows the detrending model of the light curves and the residuals, which have a standard deviation of 210 ppm.

After detrending the light curves, we used the transit least squares<fnref n="7" /> algorithm (Hippke & Heller 2019) to search for periodic transit events. The algorithm searches for transit features in an automatically built grid of orbital periods and transit

<footnote n="5">⁵ https://github.com/afeinstein20/stella</footnote>

<footnote n="6">⁶ https://github.com/hippke/wotan</footnote>

<footnote n="7">⁷ https://github.com/hippke/tls</footnote>
<figure><img src="image_5.png" /><figcaption>Fig. 4. Time series of the SPIRou longitudinal magnetic field $B_\ell$. The GP model using a quasi-periodic kernel is shown in red. The residuals of this model have a scatter of 0.97 G and are shown in the bottom panel. At the top of the figure, we mark the four subsets of data used in Sect. 5.2.</figcaption></figure>

durations. The grid of periods depends on the time span of the data and the transit durations come from an empirical relation detailed in Hippke & Heller (2019).

We first applied the algorithm in the whole data set, including sector 6 and sector 32. The grid of periods was automatically set between 0.6 and 365 days. The maximum signal-detection efficiency (SDE) is found at an orbital period of 11.86 d but it has low significance (SDE = 8) and the periodogram is highly degenerate. The S/N of this stacked transit signal is 2.5 with only three events in the data set: one in sector 6 and two in sector 32. To confirm or rule-out this period we applied the algorithm in sector 32 independently, because two transits were identified in this sector. We find no transits in this sector and therefore discarded transit events in the two available TESS sectors of Gl 205.

## 5. Magnetic field analysis

The SPIRou spectropolarimetric products, in particular the Stokes I and V profiles, were obtained using the Libre-ESPRIT pipeline described in Donati et al. (1997). We applied least-squares deconvolution (LSD; Donati et al. 1997) to compute the average Stokes I (unpolarized) and Stokes V (circularly polarized) line profiles for all our SPIRou observations. We used a mask of atomic lines spanning the SPIRou spectral domain and computed from a ATLAS9 local thermodynamical equilibrium (LTE) atmosphere model (Kurucz 1993), assuming an effective temperature of 3750 K and a surface gravity of log g = 5.0. We only selected lines with a relative absorption of greater than 3% (from the unpolarized continuum) to avoid an overrepresentation of weak lines in our final line list. Lines affected by strong tellurics, that is, of relative absorption deeper than 20% within ±30 km s⁻¹ from the line center, are masked out in the LSD process. The extracted Stokes I and V profiles feature a mean central wavelength of 1700 nm, an effective Landé factor of 1.25, and a relative depth of 12% with respect to the continuum.

### 5.1. Stellar rotation period

The disk-integrated longitudinal magnetic field $B_\ell$ was computed using the Stokes V and Stokes I profiles following the method of Donati et al. (1997):

$$B_{\ell}(\text{G}) = \frac{-2.14 \times 10^{11}}{\lambda_0 \cdot g_{\text{eff}} \cdot c} \frac{\int uV(u)du}{\int (I_c - I)du}, \quad (1)$$

where I and V are the unpolarized and circularly polarized LSD profiles, $I_c$ is the continuum level, u is the velocity in km s⁻¹, $\lambda_0$ is the mean wavelength, $g_{\text{eff}}$ is the effective Landé factor, and c the speed of light in km s⁻¹.

The $B_\ell$ time series of Gl 205 is listed in Table A.2. The mean value of the $B_\ell$ time series is 1.6 G, the standard deviation is 2.9 G, and the mean of the error bars is 1.0 G. We expect that the longitudinal magnetic field is modulated by the rotational period of the star, because the large-scale magnetic topology at the surface of the star is expected to evolve on a longer timescale. To measure the stellar rotation period, we first computed the generalized Lomb-Scargle (GLS) periodogram (Lomb 1976; Scargle 1982; Zechmeister & Kürster 2009) implemented in the astropy (Astropy Collaboration 2013, 2018) Python package. The GLS periodogram of the $B_\ell$ time series shows a strong peak of periodicity at 32.7 days with a false-alarm probability (FAP) of below 1% (see Fig. 10).

However, it has been demonstrated that stellar activity shows a quasi-periodic behavior (e.g., Haywood et al. 2014; Angus et al. 2018) rather than a single sinusoidal signal, as is the case of the periodicity seen in a GLS periodogram. Moreover, the magnetic field of Gl 205 is known to evolve on a timescale of a few rotation cycles (Hébrard et al. 2016). We therefore take advantage of the flexibility of the GPs to constrain the rotational period of the star using the $B_\ell$ time series. For this purpose, it is important to highlight that the three well-defined observational seasons of the $B_\ell$ (see Fig. 4) are long enough to cover a few stellar rotation cycles.

The kernel of the GPR is set to be the quasi-periodic as defined in Roberts et al. (2012):

$$k_{i,j} = A^2 \exp\left[-\frac{(x_i - x_j)^2}{2l^2} - \frac{1}{\beta^2} \sin^2\left(\frac{\pi(x_i - x_j)}{P_{\text{rot}}}\right)\right] + \sigma^2 \delta_{i,j}, \quad (2)$$

where $x_i$ and $x_j$ are two observation dates, A is the amplitude of the covariance, l is the decay time, β is the smoothing factor, $P_{\text{rot}}$ is the stellar rotation period, and σ is the uncorrelated white noise also known as the jitter term.

The posterior distributions of the hyperparameters were sampled from a Markov-chain Monte Carlo (MCMC) routine using the package emcee (Foreman-Mackey et al. 2013). We set up 50 walkers and 5000 steps after a burn-in phase of 500 steps. The priors used and the final posterior distributions are listed in
<table><caption>**Table 3.** Priors and best-fit hyperparameters of the GP model using a quasi-periodic kernel in the SPIRou longitudinal magnetic field data $B_t$.</caption><thead><tr><th>Parameter</th><th>Units</th><th>Prior</th><th>Posterior</th></tr></thead><tbody><tr><td>Mean value μ</td><td>G</td><td>U(-∞, +∞)</td><td>1.3 ± 0.9</td></tr><tr><td>jitter σ</td><td>G</td><td>U(0, +∞)</td><td>0.5 ± 0.2</td></tr><tr><td>Amplitude α</td><td>G</td><td>U(0, +∞)</td><td>3.1<sup>+0.6</sup><sub>-0.4</sub></td></tr><tr><td>decay time l</td><td>days</td><td>U(30, 300)</td><td>62<sup>+15</sup><sub>-12</sub></td></tr><tr><td>Smoothing factor β</td><td></td><td>U(0.25, 1.25)</td><td>0.6 ± 0.1</td></tr><tr><td>Rotation period P<sub>rot</sub></td><td>days</td><td>U(1, 100)</td><td>34.4 ± 0.5</td></tr><tr><td>RMS of residuals</td><td>G</td><td></td><td>0.97</td></tr><tr><td>χ<sup>2</sup></td><td></td><td></td><td>0.92</td></tr></tbody></table>

Notes. The symbol $U(a, b)$ defines a uniform prior with a and b being the minimum and maximum limits, respectively.

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** GLS periodogram of the residuals of the GPR applied in the longitudinal magnetic field $B_t$.</figcaption></figure>

Table 3 and Fig. B.1 shows the corner plot of the posterior distributions. The best-fitting model of the GPR is illustrated in Fig. 4 and has a reduced $\chi^2$ of 0.92. As a sanity check, we plot the GLS periodogram of the GP model residuals in Fig. 5, where we see that there is no periodic signal left over.

We measured a stellar rotation period of $P_{rot} = 34.4 \pm 0.5$ days and a decay time of $l = 62^{+15}_{-12}$ days, which implies that the active features could evolve relatively quickly, that is, on a timescale of about two rotation cycles, as the GP decay time has been shown to be a good indicator of the average timescale of evolution of the active features (Nicholson & Aigrain 2022). With this new estimation of the rotation period, we can derive the rotational velocity $v \sin i$. Assuming a stellar radius of $0.556 \pm 0.033 M_\odot$ and inclination of $60^\circ \pm 10^\circ$ (see Table 1), we obtained a $v \sin i$ of $0.7 \pm 0.1$ km s$^{-1}$. (Hébrard et al. 2016) obtained a $v \sin i$ of $1.0 \pm 0.5$ km s$^{-1}$ from their ZDI analysis, which is in agreement with our results.

## 5.2. Zeeman-Doppler imaging

We use Zeeman-Doppler imaging (Semel 1989; Brown et al. 1991; Donati & Brown 1997) to invert the observed Stokes V LSD profiles into distributions of the large-scale magnetic field at the surface of Gl 205. As described in Donati et al. (2006), ZDI decomposes the large-scale field vector into its poloidal and toroidal components, both expressed as weighted sums of spherical harmonics. For a given field distribution, local Stokes V profiles are computed for each resolution element of the visible hemisphere of the star using analytical expressions from the

Unno-Rachkovski's solution to the radiative transfer equation in a plane-parallel Milne-Eddington atmosphere (Unno 1956). These profiles are then (i) shifted to the local projected rotational velocity, (ii) weighted according to stellar inclination and the limb-darkening law (assumed linear with a coefficient of 0.3, Claret & Bloemen 2011), and then (iii) combined into global profiles at the times of the observations. ZDI uses a conjugate gradient algorithm to iteratively compare the synthetic profiles to the observed Stokes V LSD profiles down to a given reduced $\chi^2$. The degeneracy between multiple magnetic maps is lifted by imposing a maximum-entropy regularization condition to the fit, assuming that the map with the minimum amount of information is the most reliable (Skilling & Bryan 1984).

The intrinsic evolution of the magnetic field is not yet included in our ZDI code despite notable progress over the last few years (Yu et al. 2019; Finociety & Donati 2022). To prevent the code from focusing on the intrinsic evolution of the field rather than on its rotational modulation, we divide our observations into four subsets of data, namely $S_1$, $S_2$, $S_3$, and $S_4$, containing respectively 36, 31, 19, and 31 observations. The time windows spanned by these four data sets are indicated in Table 4 and depicted in Fig. 4. We then apply ZDI independently to each of these data sets. We note that these four seasons were defined so that the longitudinal field remains roughly periodic on timescales longer than a rotation period, ensuring that the magnetic topology does not dramatically evolve during each season and that the variation of the Stokes V profiles primarily reflects the stellar modulation rather than the intrinsic evolution of the field.

Our ZDI reconstruction includes the modeling of latitudinal differential rotation (DR) shearing the large-scale field at the surface of Gl 205 (Donati et al. 2000; Petit et al. 2002). The stellar rotation rate $\Omega$ is assumed to vary as a function of the colatitude $\theta$ such that

$$ \Omega(\theta) = \Omega_{\text{eq}} - (\cos \theta)^2 \, d\Omega, \quad (3) $$

where $\Omega_{\text{eq}}$ and $d\Omega$ stand for the rotation rate at the stellar equator and the difference in rotation rate between the equator and the pole, respectively. The best DR parameters explaining the observations are estimated using the method described in Donati et al. (2000). For a wide range of DR parameters, we apply ZDI to a fixed level of entropy. The best parameters with error bars are estimated by fitting a 2D paraboloid to the resulting $\chi^2$ distribution in the $(\Omega_{\text{eq}}, d\Omega)$ space.

The best-fitting maps of the large-scale field vector for the four subsets of data are shown in Fig. 6, and the associated magnetic properties are listed in Table 4. The best fits to the observed Stokes V LSD profiles are shown in Appendix C for the four seasons. Data sets $S_1$ to $S_4$ have been fitted to reduced $\chi^2$ of 1.0, 1.08, 1.1, and 1.05, starting from reduced $\chi^2$ of 2.4, 1.5, 4.1, and 1.4. As expected from the time series of longitudinal magnetic field shown in Fig. 4, we observe a significant fluctuation in the field topology throughout our observations. In Seasons $S_1$ and $S_3$, the field is found to be a dipole of 13–15 G (consistent with the topology found in Hébrard et al. 2016), and to be tilted at $36^\circ$ and $56^\circ$ to the rotation axis towards phases 0.78 and 0.84, respectively. A weaker and more complex field is found in Season $S_2$. Though dominantly poloidal, the magnetic energy budget is much more spread into the dipolar (45%), quadripolar (20%), and octupolar (31%) than in Seasons $S_1$ and $S_3$ (both dipolar at more than 90%). Between Seasons $S_3$ and
<table><caption>**Table 4.** Main results of the spectropolarimetric analysis of each of subset of data of the Gl 205 SPIRou observations with ZDI.</caption><thead><tr><th>Name</th><th>Start date</th><th>End date</th><th>ΔT<br>(days)</th><th>⟨B⟩<br>(G)</th><th>f<sub>pol</sub><br>(%)</th><th>f<sub>axi</sub><br>(%)</th><th>Ω<sub>eq</sub><br>(rad day<sup>-1</sup>)</th><th>dΩ<br>(rad day<sup>-1</sup>)</th><th>P<sub>eq</sub><br>(days)</th><th>P<sub>pol</sub><br>(days)</th></tr></thead><tbody><tr><td>S<sub>1</sub></td><td>11-09-19</td><td>12-12-19</td><td>92</td><td>13.2</td><td>90</td><td>67</td><td>0.205 ± 0.010</td><td>0.046 ± 0.021</td><td>30.6 ± 1.4</td><td>39.5 ± 5.8</td></tr><tr><td>S<sub>2</sub></td><td>26-08-20</td><td>08-10-20</td><td>43</td><td>8.9</td><td>72</td><td>45</td><td>0.19 ± 0.01</td><td>0.050 ± 0.046</td><td>33.1 ± 1.8</td><td>44.8 ± 15.2</td></tr><tr><td>S<sub>3</sub></td><td>13-08-21</td><td>24-09-21</td><td>42</td><td>14.8</td><td>77</td><td>46</td><td>0.207 ± 0.008</td><td>0.077 ± 0.020</td><td>30.3 ± 1.2</td><td>48.4 ± 8.1</td></tr><tr><td>S<sub>4</sub></td><td>16-11-21</td><td>30-01-22</td><td>75</td><td>12.3</td><td>35</td><td>77</td><td>0.190 ± 0.005</td><td>0.061 ± 0.011</td><td>33.0 ± 0.8</td><td>48.5 ± 4.6</td></tr></tbody></table>

**Notes.** Total time span of each subset of data ΔT, average large-scale magnetic strength ⟨B⟩, fraction of poloidal field fpol, fraction of axisymmetric field faxi, best-fitting DR parameters Ωeq and dΩ, equatorial and polar rotation periods, Peq and Ppol.

<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Radial, azimuthal and meridional components of the large-scale magnetic field at the surface of Gl 205 for the four subsets of data defined in Sect. 5.2. From top to bottom, the rows show the results for Seasons S₁ (September–December 2019), S₂ (August–October 2020), S₃ (August–September 2021), and S₄ (November 2021–January 2022). Each panel is a flattened polar view of the stellar surface. The circles indicate the equator (black solid line) and the –30, 30 and 60° parallels (dashed lines) and the ticks around the star mark the spectropolarimetric observations in unit of stellar rotation phase. The color bar depicts the magnetic flux in Gauss.</figcaption></figure>

S₄, the field quickly evolves, that is, on a two-month timescale, from a poloidal dipole to a dominantly toroidal topology<fnref n="8" />.

The best-fitting DR parameters and associated equatorial and polar stellar rotation periods are listed in Table 4, whereas the 2D χ² maps associated to the fit in the (Ωeq, dΩ) space are shown

extended without dramatically affecting the resulting magnetic topologies, which confirms that the season lengths are relatively well suited to the ZDI analysis.

<footnote n="8">⁸ We note that removing a small fraction of randomly selected Stokes V profiles in each season leads to fully consistent magnetic properties, which confirms that the field topology remains globally consistent throughout each season. In contrast, seasons S₂, S₃, and S₄ cannot be</footnote>
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Distribution of reduced $\chi^2$ as a function of $\Omega_{eq}$ and $d\Omega$ extracted from the Stokes V LSD profiles in Season S<sub>1</sub> to S<sub>4</sub> defined in Table 4. In each map, the 1 and 3 $\sigma$ contours are indicated by the black solid lines.</figcaption></figure>

in Fig. 7. Consistent latitudinal DRs of 0.05–0.08 rad d⁻¹ are detected in Seasons S<sub>1</sub>, S<sub>3</sub>, and S<sub>4</sub>, whereas solid-body rotation cannot be firmly excluded in Season S<sub>2</sub>. These values of latitudinal DRs are of the same order of magnitude as those of other early-M dwarfs with slightly faster rotation (Donati et al. 2008).

We note that more complex magnetic topologies, such as those seen in Seasons S<sub>2</sub> and S<sub>4</sub>, are associated with larger equatorial rotation periods than simpler dipolar fields. This could reflect the fact that the field is dominated by smaller-scale

structures mostly located at mid-to-high latitudes and, thus, only partially covering the stellar equator. However, the accuracy of our DR parameters is strongly limited by the relatively low latitudinal precision of our reconstruction, which is due to the very low $v \sin i$ of the star. As a consequence, no conclusions can be drawn as to the potential evolution of the DR throughout the observations. We report an inverse-variance-weighted average of the DR values from Season S<sub>1</sub> to S<sub>4</sub> of $P_{\text{eq}} = 32.0 \pm 1.2$ d for the rotation period at the stellar equator, and $P_{\text{pol}} = 45.5 \pm 0.3$ d for the rotation period at the poles.

# 6. Optical and NIR RVs

## 6.1. Periodogram analysis

The almost simultaneous acquisition of SOPHIE and SPIRou RVs allows us to quantify the stellar activity jitter in the optical and NIR domains. We computed the GLS periodogram to look for periodicities in the RVs. The RV time series including the SOPHIE and SPIRou data is shown in Fig. 1 along with its GLS periodogram. We identified a highly significant (FAP $\ll 0.1\%$) period at 34.3 days and a second peak, also with low FAP, at 37.9 days. There is also a third peak with low FAP at ~0.5 day, which is an alias of the periodicity at 1 day due to the daily sampling.

The SOPHIE and SPIRou data sets do not show the same periodicities. In the SOPHIE RVs, there is no significant periodicity found by the GLS periodogram, however the strongest peak occurs close to the $P_{\text{rot}}/2$ harmonic at ~17 days (see the top-right panel of Fig. 9). The SPIRou RVs, on the contrary, show a clear peak of periodicity at 34.4 days (see the top-left panel of Fig. 10) which is the expected stellar rotation period. In terms of the scatter of the RVs, the data sets have a root mean square (RMS) of 4.7 m s⁻¹ for SPIRou and 4.8 m s⁻¹ for SOPHIE.

## 6.2. Global modeling

In order to obtain a first approximation of the magnitude of the RV jitter, we applied a GP with a quasi-periodic kernel in the SOPHIE and SPIRou RV time series, following Sect. 5.1. The MCMC procedure to obtain the posterior distribution is set up with 50 walkers and 5000 steps after a burn-in phase of 500 steps. The priors applied in the hyperparameters are listed in Table 5. As we measured the stellar rotation period from the spectropolarimetric data, we applied a more constraining prior to the GP period of the RVs.

The best-fit models of the GP regressions are shown in Fig. 8 and the final posterior distributions of the GP hyperparameters are listed in Table 5. The periods found in the SOPHIE and SPIRou RVs of $34.2_{-0.9}^{+2.8}$ days and $39.2 \pm 2.4$ days, respectively, are consistent at a $3\sigma$ level with the period of the longitudinal magnetic field of $34.4 \pm 0.5$ d (see Sect. 5.1). The decay time found for the SOPHIE RVs is larger than for SPIRou and is consistent with more than two stellar rotations ($l = 96_{-47}^{+41}$ days). However, in the case of SPIRou, the decay time of $38_{-2}^{+4}$ days is comparable with the rotation period of $39.2 \pm 2.4$ days. This effect is problematic for the quasi-periodicity of the signal. When there is a short evolution timescale, $l \le P_{\text{rot}}$, the periodicity of the signal is meaningless (see e.g., Rajpaul et al. 2015; Barragán et al. 2022a).

We examined the effect of applying a different prior in the decay time for the SPIRou RVs, keeping it uniform but limiting it to between 70 and 150 days to ensure at least two rotation cycles. As a result, we obtained consistent posterior distributions
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Time series of the SOPHIE and SPIRou RVs. *Top*: the best-fitting GP models using a quasi-periodic kernel in the SOPHIE and SPIRou data are shown in blue and red, respectively. The light-blue and light-red colors depict the $3\sigma$ level of the uncertainties. The residuals of this model have a scatter of $2.6 \text{ m s}^{-1}$ for SOPHIE, and $2.4 \text{ m s}^{-1}$ for SPIRou and are shown below the time series. *Middle and bottom*: zoom onto the SOPHIE and SPIRou RVs divided into the seasons defined in Table 6 with their best-fitting GP models.</figcaption></figure>

– within the error bars – with the values from Table 5, except for the decay time. With a different prior, we obtained a decay time of $73_{-2}^{+5}$ days. However, the log-likelihood of the first model (log-$L = -374$) is slightly higher than the model with a new prior in the decay time (log-$L = -384$), meaning that the first model shows a better fit to the data.

## 6.3. Seasonal analysis

The RVs observation time span is clearly divided in three seasons (see Fig. 1), with the first epoch from 58738 BJD to 58922 BJD (September 2019 to March 2020), the second epoch from 59088 BJD to 59297 BJD (August 2020 to March 2021), and the third epoch from 59440 BJD to 59607 BJD (August 2021 to January 2022). These subsets are longer than the subsets defined in Sect. 5.2 because in order to look for periodic variability, a longer time-span is required than for the ZDI analysis. The details of the RV subsets, such as the start and ending dates as well as the number of data points included, are listed in Table 6. These RV subsets, which are depicted in Fig. 1, are

<table><caption>**Table 5.** Priors and best-fit hyperparameters of the GP model using a quasi-periodic kernel in the SOPHIE and SPIRou RVs.</caption><thead><tr><th>Parameter</th><th>Units</th><th>Prior</th><th>SOPHIE<br>posterior</th><th>SPIRou<br>posterior</th></tr></thead><tbody><tr><td>Mean value μ</td><td>m s<sup>-1</sup></td><td>U(-∞, +∞)</td><td>0.6 ± 1.4</td><td>0.5 ± 1.7</td></tr><tr><td>RV jitter σ</td><td>m s<sup>-1</sup></td><td>U(0, +∞)</td><td>2.5 ± 0.8</td><td>2.0 ± 0.3</td></tr><tr><td>Amplitude α</td><td>m s<sup>-1</sup></td><td>U(0, +∞)</td><td>4.1 ± 1.1</td><td>5.4 ± 1.2</td></tr><tr><td>Decay time l</td><td>days</td><td>U(30, 300)</td><td>96<sup>+41</sup><sub>-47</sub></td><td>38<sup>+4</sup><sub>-2</sub></td></tr><tr><td>Smoothing factor β</td><td></td><td>U(0.25, 1.25)</td><td>0.5 ± 0.2</td><td>0.9 ± 0.3</td></tr><tr><td>Rotation period P<sub>rot</sub></td><td>days</td><td>U(30, 50)</td><td>34.2<sup>+2.8</sup><sub>-0.9</sub></td><td>39.2 ± 2.4</td></tr><tr><td>RMS of residuals</td><td>m s<sup>-1</sup></td><td></td><td>2.6</td><td>2.4</td></tr><tr><td>χ<sup>2</sup></td><td></td><td></td><td>2.1</td><td>1.8</td></tr></tbody></table>

**Notes.** The symbol $U(a, b)$ defines a uniform prior with $a$ and $b$ being the minimum and maximum limits, respectively.

defined as $S_1^{RV}$, $S_2^{RV}$, and $S_3^{RV}$ in order to differentiate them from the seasons defined in Sect. 5.2.
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** SOPHIE RVs and optical activity indicator. Left: time series; middle left: data folded at the $P_{\text{rot}}$ of 34.4 days; middle right: correlation with the RVs where the red line depicts the linear correlation between both quantities; right: GLS periodogram where the horizontal dashed black line indicates the FAP level at 1% and the vertical color lines mark the $P_{\text{rot}}$ at 34.4 days in yellow, the first harmonic at $P_{\text{rot}}/2$ in light blue, and the time range of differential rotation in light green.</figcaption></figure>

Even though the seasons $S_1^{\text{RV}}$ and $S_2^{\text{RV}}$ overlap to some extent with $S_1$ and $S_2$, we cannot directly compare $S_3$ with the RV data because only three SOPHIE observations are within this season. Therefore, we defined a fourth season of RV observations, $S_4^{\text{RV}}$, from 59500 BJD to 59607 BJD (October 2021 to January 2022) that overlaps with $S_4$ and allows us to measure the differences between seasons $S_3$ and $S_4$ of the $B_\ell$ in the RVs data sets.

To gain insight into the variability of the RV signal, we define the amplitude $A$ as the peak-to-valley difference and we list it alongside the RMS of the model residuals and the reduced $\chi^2$ in Table 6 for all the data sets and each RV season.

It is clearly seen in Fig. 8 that the best-fitting GP model for the SOPHIE RVs is more stable in time than that for the SPIRou RVs, and is similar to a double sinusoidal model at $P_{\text{rot}}$ and its first harmonic. As seen in the middle and bottom panels of Fig. 8, the signal in the SOPHIE RVs remains approximately consistent between $S_1^{\text{RV}}$ and $S_2^{\text{RV}}$, with similar shapes and amplitudes of $10.8 \pm 2.9$ and $9.8 \pm 3.0 \, \text{m s}^{-1}$, respectively. The goodness of the model is similar for the seasons $S_1^{\text{RV}}$ through $S_4^{\text{RV}}$, with an RMS of between $2.3$ and $2.8 \, \text{m s}^{-1}$. However, during $S_3^{\text{RV}}$ and $S_4^{\text{RV}}$, the

reduced $\chi^2$ becomes slightly poorer than for the other seasons. This may be due to the presence of outliers in the subdata set.

On the contrary, the SPIRou RV signal is not consistent from one season to the next (see bottom panels of Fig. 8), and exhibits a high level of variability. The amplitude of the RV signal rises slightly between $S_1^{\text{RV}}$ and $S_4^{\text{RV}}$, from $13.7 \pm 2.3$ to $15.2 \pm 2.2 \, \text{m s}^{-1}$, but this rise is not significant and the difference is within the error bars. As in the SOPHIE RVs, there are no important differences between the RMSs or reduced $\chi^2$ of each season, but the goodness of the model fit deteriorates for $S_3^{\text{RV}}$. As in the case of the SOPHIE RVs, this may be due to the presence of outliers; if these outliers are removed, the RMS improves to $2.4 \, \text{m s}^{-1}$.

At this point, we cannot tell whether or not there are significant differences between $S_3^{\text{RV}}$ and $S_4^{\text{RV}}$ of the SOPHIE and SPIRou RVs, and we require this information in order to measure the possible impact of changes in the magnetic field topology on the RVs (see Sect. 5.2). As listed in Table 5, there is little discrepancy between these two seasons in the amplitudes, RMSs, and reduced $\chi^2$. Despite the fact that the $B_\ell$ shows a high level of variability at these epochs, the scatter of the RVs does not
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Same as Fig. 9 but for SPIRou data.</figcaption></figure>

increase. This may mean that the variability in the longitudinal
magnetic field is not directly expressed as high RV dispersion.

**7. Activity indicators**

**7.1. Optical activity indicators**

Using the SOPHIE data, we computed standard activity indica-
tors that quantify distortions in the CCF profile, such as the BIS,

the FWHM, and the contrast. The BIS is defined as the velocity
span between the top (55% < CCF depth < 80%) and the bottom
(20% < CCF depth < 40%) of the CCF. The FWHM and the con-
trast are direct measurements of the width and the depth of the
CCF, respectively.

The time series of the BIS, FWHM, and contrast of the CCF
from the SOPHIE data are shown in Fig. 9, with their corre-
sponding GLS periodograms. No significant periodicity close to
<table><caption>**Table 6.** Results of the RV analysis in the observational seasons of the SOPHIE and SPIRou.</caption><thead><tr><th rowspan="2">Name</th><th rowspan="2">Start</th><th rowspan="2">End</th><th rowspan="2">N</th><th colspan="3">SOPHIE</th><th colspan="3">SPIRou</th></tr><tr><th>A<br>(m s<sup>-1</sup>)</th><th>RMS<br>(m s<sup>-1</sup>)</th><th>χ<sup>2</sup><br>–</th><th>A<br>(m s<sup>-1</sup>)</th><th>RMS<br>(m s<sup>-1</sup>)</th><th>χ<sup>2</sup><br>–</th></tr></thead><tbody><tr><td>All</td><td>09-09-19</td><td>27-01-22</td><td>201</td><td>13.6 ± 2.5</td><td>2.6</td><td>2.1</td><td>21.3 ± 4.5</td><td>2.4</td><td>1.8</td></tr><tr><td>S<sup>RV</sup><sub>1</sub></td><td>11-09-19</td><td>13-03-20</td><td>53</td><td>10.8 ± 2.9</td><td>2.6</td><td>1.9</td><td>13.7 ± 2.3</td><td>2.4</td><td>1.8</td></tr><tr><td>S<sup>RV</sup><sub>2</sub></td><td>26-08-20</td><td>04-03-21</td><td>72</td><td>9.8 ± 3.0</td><td>2.3</td><td>2.1</td><td>12.0 ± 3.0</td><td>2.0</td><td>1.5</td></tr><tr><td>S<sup>RV</sup><sub>3</sub></td><td>13-08-21</td><td>27-01-22</td><td>76</td><td>8.6 ± 3.3</td><td>2.7</td><td>3.4</td><td>15.2 ± 2.2</td><td>3.0</td><td>2.9</td></tr><tr><td>S<sup>RV</sup><sub>4</sub></td><td>12-10-21</td><td>27-01-22</td><td>58</td><td>8.6 ± 2.8</td><td>2.8</td><td>4.0</td><td>15.2 ± 2.2</td><td>2.9</td><td>2.8</td></tr></tbody></table>

**Notes.** Number of data points *N*. The amplitude of the RV variation is defined as the peak-to-valley difference of the GP model. The goodness of the fit is described with the RMS of the residuals alongside the reduced $\chi^2$.

the expected stellar rotation period of $P_{\text{rot}} = 34.4 \pm 0.5$ days is found in their periodograms.

Within the SOPHIE spectral domain, there is the well-known $H\alpha$ line (Kürster et al. 2003; Bonfils et al. 2007; Boisse et al. 2009) and the CaII H&K lines, which are measured to compute the Mount Wilson S index (Wilson 1968; Baliunas et al. 1995). The so-called S index is a measure of the emission at the core of the CaII H&K lines, located at 3922 Å and 3968 Å, and for the $H\alpha$ is the flux at the core of the line at 6562 Å. To compute this index, we follow Boisse et al. (2009), who define it as the ratio between the flux in a defined range of continuum and the core of the absorption lines.

Figure 9 shows the time series of the $H\alpha$ and S indices of our SOPHIE data along with their GLS periodograms. The time series of $H\alpha$ shows a peak of periodicity at 33.7 d and in the case of the S index, long-term trends dominate the periodogram.

The correlation or anti-correlation between the RVs and the activity indicators usually provides information about the stellar origin of the signal (e.g., Queloz et al. 2001; Boisse et al. 2011). However, the lack of correlation can be due to phase shifts between the RVs and the activity indicators (e.g., Bonfils et al. 2007; Dumusque et al. 2014b). For example, Collier Cameron et al. (2019) measured a temporal lag of 1 and 3 days between the maxima in the RVs and the maxima of the FWHM and BIS, respectively. In the present work, we do not explore the possibility of the phase lags in our data set but we rather warn the reader about it.

The SOPHIE BIS, FWHM, and contrast are not correlated or anti-correlated with the RVS, and the S index and $H\alpha$ are correlated with the RVs with Pearson's coefficients of 0.4. We computed the p-value associated with this correlation in order to prove the significance. The p-value is defined as the probability that the observed correlation is a false positive under a true null hypothesis. In this case, the null hypothesis is that there is no correlation between the variables. The correlations between the RVs and the $H\alpha$ and S indices have p-values of $8 \times 10^{-4}$ and $9 \times 10^{-4}$, respectively. This means that the correlations we find are statistically significant.

We applied a GPR following the same procedure as in Sect. 5.1 to study the periodicity of the signal, as we have seen how the longitudinal magnetic field and the RVs follow a quasi-periodic trend. This behavior could have an effect in the classical GLS periodogram and hide the real period of the signal. The computation of the FWHM and the contrast by the SOPHIE DRS do not provide their error bars. Therefore, we assumed equal

error bars for all the data points of 0.1% for the contrast, and 0.01 km s$^{-1}$ for the FWHM.

We used the same priors as in Sect. 6, which are listed in Table 5. To obtain the posterior distribution from the MCMC we used 50 walkers and 5000 steps after a burn-in phase of 500 steps. The best-fitted GP models for the SOPHIE activity indicators are shown in Fig. E.1. The means of the posterior distributions of the GP hyperparameters are listed in Table E.1 and the corner plots of the posterior distributions are shown in Fig. E.1.

In the case of the SOPHIE activity indicators, the quasi-periodic GP fits a model with a periodicity in agreement with the expected rotation period of $P_{\text{rot}} = 34.4 \pm 0.5$ days. This is particularly interesting for the BIS, FWHM, and CCF contrast, because their GLS periodograms do not exhibit significant peaks of periodicity at the rotation period. The decay time of the GP model of the FWHM, BIS, and $H\alpha$ are consistent with an average timescale decay of at least two rotation cycles, as also seen in the longitudinal magnetic field $B_{\ell}$.

## 7.2. NIR activity indicators

Within the SPIRou LBL framework, we can obtain the differential line width (dLW, Zechmeister et al. 2018) and the chromatic velocity slope. The differential line width corresponds to the second derivative of the spectral profile and can be expressed in units of FWHM. We follow the definition by Artigau et al. (2022) using this quantity as FWHM$_{LBL}$. The chromatic velocity slope is defined as the RV gradient as a function of wavelength (Zechmeister et al. 2018).

The time series and the GLS periodograms of the FWHM$_{LBL}$ and the chromatic velocity slope are shown in Fig. 10 along with their phase-folded and correlation with RVs plots. When phase-folding FWHM$_{LBL}$ by the stellar rotation period of 34.4 days, it seems to be modulated in time. Moreover, FWHM$_{LBL}$ is slightly correlated with the RVs with a Pearson's coefficient of 0.3 and a p-value of $4 \times 10^{-5}$, which means that the found correlation is significant. The highest peak of the FWHM$_{LBL}$ GLS periodogram is at 34.4 days but with FAP slightly above 1%, which corresponds to the stellar rotation period. The chromatic velocity slope indicator is not correlated or anti-correlated with the RVs, but the most significant peak of periodicity in its periodogram is located at 31.3 days, which is close to the rotation period.

As for the SOPHIE activity indicators, we also computed a GP model in the SPIRou activity indicators including the FWHM, chromatic velocity slope, and pEW of the NIR spectral lines. The best-fit GP models for the FWHM and the chromatic velocity slope are shown in Fig. E.2 and the means of the posterior distribution of the hyperparameters are listed in Table E.1. The corner plots of the posterior distributions are shown in Figs. E.4 and E.5. The period measured in the FWHM and the chromatic velocity slope are in agreement – within the error bars – with the expected rotation period of $34.4 \pm 0.5$. Similar to the activity indicators in the optical, the measured decay time of the signal is consistent with two stellar rotation cycles. The results of the NIR spectral lines are discussed in the following subsection.

## 7.3. Diagnostics on NIR spectral lines

Observations in the infrared domain introduce two challenges. On the one hand, we lack well-identified and characterized spectral lines in this domain that can be used as stellar activity tracers. On the other hand, M dwarf spectra do not always show a clear flux continuum within the whole NIR wavelength range, because
<table><caption>**Table 7.** Properties and periodicities of the NIR spectral lines within the SPIRou domain.</caption><thead><tr><th rowspan="2">Line</th><th colspan="2">Line information</th><th colspan="4">Periodicities</th></tr><tr><th>Central wavelength<br>(Å)</th><th>Continuum range<br>(Å)</th><th>All data<br>(days)</th><th>S<sup>RV</sup><sub>1</sub><br>(days)</th><th>S<sup>RV</sup><sub>2</sub><br>(days)</th><th>S<sup>RV</sup><sub>3</sub><br>(days)</th></tr></thead><tbody><tr><td>Ti I</td><td>10 498.99</td><td>10 490 – 10 510</td><td>4.7</td><td>39.6</td><td>8.2</td><td>9.4</td></tr><tr><td>He I</td><td>10 833.25</td><td>10 820 – 10 845</td><td>6.2</td><td>1.6</td><td>27.2</td><td>25.3</td></tr><tr><td>Fe I</td><td>11 693.18</td><td>11 684 – 11 704</td><td><strong>44.7</strong></td><td>30.8</td><td>39.5</td><td>33.4</td></tr><tr><td>K I</td><td>12 435.58</td><td>12 410 – 12 460</td><td>40.5</td><td>12.3</td><td>39.5</td><td>39.5</td></tr><tr><td>Al I</td><td>13 154.32</td><td>13 140 – 13 170</td><td><strong>35.0</strong></td><td>2.0</td><td>35.5</td><td><strong>41.4</strong></td></tr></tbody></table>

**Notes.** Wavelength values are in vacuum. Peaks of the GLS periodograms with FAP < 1% are shown in bold. We searched for periodicities related to the stellar rotation period in a range between 1 *d* and 100 *d*.

their spectra have a large number of absorption lines and hamper the proper measurement of equivalent widths.

We selected NIR absorption lines of chemical species that are proven activity indicators in optical or NIR wavelengths for Sun-like stars or M dwarfs in order to test their potential as activity indicators for SPIRou. These absorption lines are the titanium (Ti I) line at 10 499 Å (Spina et al. 2020), the helium triplet (He I) at 10 833 Å (Schöfer et al. 2019; Fuhrmeister et al. 2019, 2020), the iron (Fe I) line at 11 692 Å (Yana Galarza et al. 2019; Cretignier et al. 2020), the potassium line (K I) at 12 435 Å (Barrado y Navascués et al. 2001; Robertson et al. 2016; Terrien et al. 2022), and the aluminum (Al I) line at 13 154 Å (Spina et al. 2020). The specific central wavelength of these lines is listed in Table 7.

The equivalent width (EW) of spectral lines can be used as a proxy for the strength of the chromospheric activity and is defined as follows Gray (2008):

$$ EW = \int \left(1 - \frac{F_{\lambda}}{F_{0}}\right) d\lambda, \qquad (4) $$

where $F_{\lambda}$ is the flux of the spectral line and $F_{0}$ is the continuum flux surrounding the spectral line.

As the continuum in the M dwarfs spectra is unknown, we used a pseudo-continuum as in Schöfer et al. (2019), and therefore we measured the pEW. To do so, we first stacked the four consecutive exposures of each SPIRou observation to obtain one single high-S/N spectrum per epoch. We then used a Python re-implementation of the continuum-fitting routine from IRAF<fnref n="9" /> to obtain the pseudo-continuum. This pseudo-continuum is defined as the local continuum of a spectral zone of 2–5 Å, where the spectral line is at the center. The details of the wavelength range that we used to compute the pseudo-continuum are listed in Table 7.

To each spectrum per epoch, we fitted a Voigt profile in the absorption lines within an MCMC framework using the emcee (Foreman-Mackey et al. 2013) Python package. The parameters of the Voigt profile are the Gaussian width, the Cauchy-Lorentz scale, the amplitude, and the center of the line. The fitted model of the Voigt profile corresponds to $F_{\lambda}$ in Eq. (4). We used the posterior distributions of the fitted parameters to compute the uncertainties on the pEWs using the bootstrap method. To look for periodicities in the time series of the tested spectral lines, we computed the GLS periodogram. We then considered any peak with a FAP of below 1% as significant. Table 7 shows the highest peaks of the periodicities in the GLS periodograms for each line

within a range of 1 day to 100 days; this excludes periodicities coming from the 1-day alias and long-term variability.

Figure 10 displays the pEW time series of the Al I, Ti I, K I, Fe I, and He I lines along with their GLS periodograms, the phase-folded plot, and the correlation plot with RVs. The highest peak of the periodicity of each line is listed in Table 7 for all the data and for each of the RV seasons. Only Al I exhibits peaks of periodicities with a FAP of below 1% at 35 days, which is close to the expected stellar rotation of 34.4 days. Moreover, Al I shows several peaks at between 30 and 40 days, which may be related to the differential rotation. The lines best correlated with the RVs are Al I and Fe I with Pearson's coefficients of 0.3 and 0.2, respectively. However, the correlation between the RVs and Al I is more significant, with a p-value of $4 \times 10^{-4}$ versus the p-value of 0.01 for the correlation with Fe I.

Excluding Al I, the remaining spectral lines tested, namely Fe I, Ti I, K I, and He I, do not seem to be modulated by the stellar rotation and show small or no correlation with the RVs, at least from the analysis of their periodograms. The Fe I line in particular shows a peak of periodicity at 44.7 days, which is similar to the expected stellar rotation at the poles ($P_{pol}$ = 45.5 days). This line could be tracing active surface features at latitudes close to the stellar pole.

The best-fitting GP models of the tested spectral lines are shown in Fig. E.2 and the final posterior distributions of the hyperparameters are listed in Table E.1. Most of the spectral lines have periodicities of longer than the measured rotation period of 34.4 ± 0.5, except for He I. However, suggested by the previous results within this work, it seems that the He I line is not sensitive to stellar activity. Moreover, it is clear that the GPs do not fully model the variability in the tested spectral lines, and we observe high values of reduced $\chi^2$. Underestimation of the uncertainties on the pEW computation could play a role in this; we will investigate this further in future works.

Thanks to the analysis of the GLS periodograms and the GPR, we can see that the periodicity of the Al I and Fe I lines as computed using these two techniques is consistent; however, the decay time is shorter than expected and does not cover two complete cycles.

## 7.4. Seasonal analysis

As stellar activity shows a quasi-periodic behavior, the activity signals should be consistent over short timescales, which means we should be able to measure their periodicities with the periodograms in a seasonal analysis. We observe that the strongest peaks of the GLS periodograms of the optical activity indicators are not consistent over the three seasons of data previously

<footnote n="9">⁹ Image Reduction and Analysis Facility, https://iraf-community.github.io</footnote>
defined as S<sub>1</sub><sup>RV</sup>, S<sub>2</sub><sup>RV</sup>, and S<sub>3</sub><sup>RV</sup>. The GLS periodograms for each of the three subsets are displayed in Fig. D.1. For the CCF contrast and the BIS, the periodicity close to the rotation period gains significance towards the second subset, but still has a high FAP. On the other hand, the peak at the stellar rotation period decreases in significance towards the 2021/2022 observations. During S<sub>1</sub><sup>RV</sup>, the Hα time series shows an important power excess at a period slightly longer than the expected rotation period. During S<sub>2</sub><sup>RV</sup> and S<sub>3</sub><sup>RV</sup>, Hα keeps some power excess around the rotation period but this is not significant.

The SPIRou activity indicators, the FWHM<sub>LBL</sub> and the chromatic velocity slope, do not show the same seasonal behavior. As shown in Fig. D.2, the FWHM<sub>LBL</sub> has a peak of periodicity that is unrelated to the stellar rotation in S<sub>1</sub><sup>RV</sup>; during S<sub>2</sub><sup>RV</sup> and S<sub>3</sub><sup>RV</sup> we can see some power excess close to the rotation period but with a FAP of above 10%. On the other hand, the chromatic velocity slope has a highly significant peak of periodicity that is slightly shorter than the rotation period in S<sub>3</sub><sup>RV</sup>, but no peaks in S<sub>1</sub><sup>RV</sup> or S<sub>3</sub><sup>RV</sup>. The periodicity seen in S<sub>3</sub><sup>RV</sup> may be affected by the window function of that particular season (see Fig. D.2).

The tested NIR spectral lines also do not keep consistent periodicities across the three subsets of RVs observations S<sub>1</sub><sup>RV</sup>, S<sub>2</sub><sup>RV</sup>, and S<sub>3</sub><sup>RV</sup>. The strongest peaks of the GLS periodograms of Al I, Ti I, K I, and Fe I change from season to season and are close to the rotation period during at least one season (see Table 7). However, this latter periodicity is ~39 d instead of the expected P<sub>rot</sub>. This may be related to active features in the stellar surface at mid latitudes where the rotation should be greater than at the equator. Regarding the Fe I line in particular, although when considering the whole data set its periodicity is larger than P<sub>rot</sub>, the periodicity of each RV season is between 30 and 40 days, which could originate from the stellar rotation modulation. This variability in the periodicity observed for the spectral lines may suggest high temporal variability in the stellar activity manifesting as highly variable active features.

Although we see hints of seasonal behavior in the optical and NIR activity indicators in this work, most of their periodicities are not significant and therefore we cannot describe trends of the stellar activity during each season. For example, Hα is clearly modulated by the stellar rotation during S<sub>1</sub><sup>RV</sup> which could be due to high chromospheric activity levels. However, we do not see this behavior in other activity indicators.

## 8. Stellar activity modeling

Rajpaul et al. (2015) describes a multi-dimensional GP framework within which the stellar activity in the RV time series can be modeled simultaneously with the Keplerian signal using the information from activity indicators. In summary, this framework assumes that the stellar activity signals and the RVs can be modeled by the same latent function, namely G(t), and its time derivative, Ġ(t). Here, G(t) is related to the area of the visible stellar disk covered by active regions and Ġ(t) describes the evolution of these active regions over time. Therefore, the RV variations induced by activity can be expressed as a linear combination of G(t) and Ġ(t). Physically, this linear combination will account for the convective blueshift suppression and the evolution of spots on the stellar surface.

In principle, the ancillary time series, such as the activity indicators, should also be defined as linear combinations of G(t) and Ġ(t) if they are affected by both processes: the area covered by active regions and its time evolution on the stellar surface (e.g., BIS). However, some activity indicators, such as logG'<sub>HK</sub>

and S<sub>HK</sub>, can be described only by G(t) as they do not account for the time evolution of spots (see Rajpaul et al. 2015 and Barragán et al. 2022a for more details).

We used the Python code pyaneti<sup>10</sup> (Barragán et al. 2019b, 2022a) to apply the multi-dimensional GP framework in our SOPHIE and SPIRou data. As the code was built for planet RV and transit modeling and there is no planet signal detected in this system, we set up a planet signal with an amplitude of K = 0 m s<sup>-1</sup>. Then only the stellar activity signal is considered in the modeling. The kernel of the GPR is set to be the same as in Eq. (2) of Sect. 5. However, the amplitude A in Eq. (2) in the multidimensional framework is replaced by the amplitude C<sub>i</sub>, as follows:

$$ 
\begin{aligned}
\Delta RV &= C_0 G(t) + C_1 \dot{G}(t) \\
\Theta_1 &= C_2 G(t) + C_3 \dot{G}(t) \\
\Theta_2 &= C_4 G(t) + C_5 \dot{G}(t) \\
\dots \\
\Theta_N &= C_i G(t) + C_{i+1} \dot{G}(t),
\end{aligned} 
\quad (5) $$

where Θ<sub>1,...,N</sub> are the activity indicators. In this way, it is guaranteed that the hyperparameters of the kernel will be shared, except for the amplitudes C<sub>i</sub>. We follow the procedure as in Sect. 5 to obtain the posterior distributions.

We modeled the stellar activity of the SOPHIE RVs using the Hα index and the BIS as ancillary time series, as follows:

$$ 
\begin{aligned}
\Delta RV_{\text{SOPHIE}} &= C_0 G(t) + C_1 \dot{G}(t) \\
H\alpha &= C_2 G(t) \\
BIS &= C_4 G(t) + C_5 \dot{G}(t).
\end{aligned} 
\quad (6) $$

We chose the Hα instead of the S index because the former shows a peak in periodicity at the stellar rotation period and is clearly correlated with the RVs (see Fig. 9). Details of the computation of the activity indicators are provided in Sect. 7.

The model obtained from the multidimensional GPR is shown in Fig. 11 for the RVs, the Hα, and the BIS time series. The priors and posteriors are listed in Table 8 and Fig. F.1 shows the corner plot of the posterior distributions. The RV jitter is 2.7 ± 0.6 m s<sup>-1</sup> and the rotation period is determined at 35.4 ± 1.0 days. The residuals of the RVs have a scatter of 2.9 m s<sup>-1</sup> and do not show any significant periodicity in the GLS periodogram.

In the case of the SPIRou data set, we used the FWHM<sub>LBL</sub> as ancillary time series as it is the activity indicator that is best correlated to the RVs and shows a significant peak at the stellar rotation period (see Fig. 10). The model is as follows:

$$ 
\begin{aligned}
\Delta RV_{\text{SPIRou}} &= C_0 G(t) + C_1 \dot{G}(t) \\
\text{FWHM}_{\text{LBL}} &= C_2 G(t).
\end{aligned} 
\quad (7) $$

We only used FWHM<sub>LBL</sub> because there is no analog to the BIS or any other activity indicator that could trace the temporal evolution of the active regions in the SPIRou data set. The model obtained for the RVs and the FWHM<sub>LBL</sub> is shown in Fig. 12 and the priors and posteriors are listed in Table 8. Figure F.2 shows a corner plot of the posterior distributions. The rotation period is found at 37.7<sup>+3.0</sup><sub>-2.8</sub> days, the RV jitter is 2.7 ± 0.3 m s<sup>-1</sup>, and the RMS of the RV residuals is 3.1 m s<sup>-1</sup>. Also, the GLS periodogram of the RV residuals does not show any periodicity.

<footnote n="10"><sup>10</sup> https://github.com/oscaribv/pyaneti</footnote>
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Multi-dimensional GP regression of the SOPHIE RVs using the $H\alpha$ index and the BIS as ancillary time series. The black line indicates the GP model and the gray areas are the $1\sigma$ and $2\sigma$ levels. The residuals of the GP model are shown below each time series. The mean value has been removed from all the time series.</figcaption></figure>

We tested the $B_\ell$ as an ancillary time series together with the FWHM<sub>LBL</sub>. The $B_\ell$ is a proven stellar activity tracer and we could therefore use its information to improve the modeling of the stellar activity signal. The stellar rotation period obtained from this test is in better agreement with the value obtained in Sect. 5 than the period obtained when only using the FWHM<sub>LBL</sub> as an ancillary time series, with a value of $P_{\text{rot}} = 34.6 \pm 0.8$ days. However, the scatter of the residuals is higher with a value of $3.7 \text{ m s}^{-1}$.

In Fig. 11, we see the complexity in the structures describing the SOPHIE RVs, which is due to the high harmonic complexity of the signal. We can probe this with the smoothing factor $\beta$ value from the GP model. The smoothing factor $\beta$ is higher for the SOPHIE RVs ($\beta = 0.81_{-0.20}^{+0.29}$) than for SPIRou ($\beta = 0.62_{-0.18}^{+0.31}$).

The high harmonic complexity seen in the optical can be empirically proven with the amplitudes of $G(t)$ and $\dot{G}(t)$ (Barragán et al. 2022b). The function $\dot{G}(t)$ describes the area covered by active regions on the stellar surface as a function of time (Aigrain et al. 2012; Rajpaul et al. 2015). We see that the amplitude of $G(t)$, $C_0 = 3.6_{-1.0}^{+1.6}$ is smaller than the amplitude of $\dot{G}(t)$, $C_1 = 24.2_{-8.9}^{+12.3}$, which explains the observed complexity.

However, we see a different result in the NIR. The amplitudes of $G(t)$ and $\dot{G}(t)$ are similar with values of $C_0 = 3.8_{-0.8}^{+1.5}$ and $C_1 = 3.2_{-2.3}^{+4.7}$. This could be an indication that the function $\dot{G}(t)$ is not fully describing the origin of the RV variations and might be evidence of a different process dominating the stellar activity in the NIR, such as the Zeeman effect.

## 9. Detection of planet RV signals

To date, no planet has been found to orbit Gl 205. From our RV analysis of Sect. 6, we do not observe periodicities that could be attributable to Keplerian signals. Considering the total time span of our SOPHIE and SPIRou observations of almost 900 days, we can discard long-period planets with orbital periods longer than 450 days. In particular, the detection of long-period planets would be highly affected by the two gaps of observations of ~100 days each (see Fig. 1). Further observation campaigns of Gl 205 may yet lead to the discovery of long-period planets around this star.

The well-defined seasons of RV observations in Fig. 1 last 184, 209, and 167 days, respectively, and therefore our highest
<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Multidimensional GP regression of the SPIRou RVs using the FWHM<sub>LBL</sub> as an ancillary time series. The black line indicates the GP model and the gray areas are the $1\sigma$ and $2\sigma$ levels. The residuals of the GP model are shown below each time series. The mean value has been removed from all the time series.</figcaption></figure>

<table><caption>**Table 8.** Priors and best-fit hyperparameters of the multi-dimensional GP model using a quasi-periodic kernel in the SOPHIE and SPIRou data.</caption><thead><tr><th>Parameter</th><th>Prior</th><th>SOPHIE<br>posterior</th><th>SPIRou<br>posterior</th></tr></thead><tbody><tr><td>Mean RV (m s<sup>-1</sup>)</td><td>U(-∞, +∞)</td><td>8729.6 ± 1.5</td><td>9345.9 ± 1.0</td></tr><tr><td>RV jitter (m s<sup>-1</sup>)</td><td>U(0, 100)</td><td>2.7 ± 0.6</td><td>2.7 ± 0.3</td></tr><tr><td>Mean Hα</td><td>U(-∞, +∞)</td><td>0.267 ± 0.009</td><td></td></tr><tr><td>Hα jitter</td><td>U(0, 100)</td><td>0.007 ± 0.001</td><td></td></tr><tr><td>Mean BIS (m s<sup>-1</sup>)</td><td>U(-∞, +∞)</td><td>1.8 ± 1.0</td><td></td></tr><tr><td>BIS jitter (m s<sup>-1</sup>)</td><td>U(0, 100)</td><td>2.0 ± 0.5</td><td></td></tr><tr><td>Mean FWHM<sub>LBL</sub> (m s<sup>-1</sup>)</td><td>U(-∞, +∞)</td><td></td><td>6070 ± 7</td></tr><tr><td>FWHM<sub>LBL</sub> jitter (m s<sup>-1</sup>)</td><td>U(0, 100)</td><td></td><td>20.7 ± 1.7</td></tr><tr><td>C<sub>0</sub></td><td>U(0, 100)</td><td>3.6<sup>+1.6</sup><sub>-1.0</sub></td><td>3.8<sup>+1.5</sup><sub>-0.8</sub></td></tr><tr><td>C<sub>1</sub></td><td>U(0, 100)</td><td>24.2<sup>+12.3</sup><sub>-8.9</sub></td><td>3.2<sup>+4.7</sup><sub>-2.3</sub></td></tr><tr><td>C<sub>2</sub></td><td>U(0, 100)</td><td>17.1<sup>+6.5</sup><sub>-4.0</sub></td><td>22.4<sup>+7.8</sup><sub>-5.3</sub></td></tr><tr><td>C<sub>4</sub></td><td>U(0, 100)</td><td>9.5<sup>+0.8</sup><sub>-0.7</sub></td><td></td></tr><tr><td>C<sub>5</sub></td><td>U(0, 100)</td><td>-17.5<sup>+7.9</sup><sub>-12.3</sub></td><td></td></tr><tr><td>decay time l [d]</td><td>U(30, 150)</td><td>86<sup>+10</sup><sub>-15</sub></td><td>37<sup>+11</sup><sub>-5</sub></td></tr><tr><td>smoothing factor β</td><td>U(0.1, 2.0)</td><td>0.81<sup>+0.29</sup><sub>-0.20</sub></td><td>0.62<sup>+0.31</sup><sub>-0.18</sub></td></tr><tr><td>rotation period P<sub>rot</sub> (d)</td><td>U(20, 50)</td><td>35.4 ± 1.0</td><td>37.7<sup>+3.0</sup><sub>-2.8</sub></td></tr><tr><td>RMS of RV residuals (m/s)</td><td></td><td>2.9</td><td>3.1</td></tr><tr><td>χ<sup>2</sup></td><td></td><td>1.02</td><td>1.01</td></tr></tbody></table>

**Notes.** The symbol *U*(*a*, *b*) defines a uniform prior, with *a* and *b* being the minimum and maximum limits, respectively.

detection sensitivity is given for planets with periods of less than 100 days and with amplitudes greater than 4.8 m s⁻¹, which corresponds to the scatter in the RVs. Moreover, we do not observe periodicity peaks in the GLS periodograms of the SOPHIE and SPIRou RVs that are not related to the stellar rotation period.

The observed peaks of periodicity in the RVs GLS periodogram are not consistent during the three seasons of observations and are always correlated with one or more of the activity indicators studied here. Furthermore, the procedure to filter the stellar activity described in Sect. 8 allows us to discard remnant Keplerian periodic signals in the RVs. In both cases, for the SOPHIE and SPIRou data sets, the multidimensional GP absorbs most of the activity-induced RV signal, leaving an RMS of the residuals of ~3 m s⁻¹. The GLS periodogram of these RV residuals does not exhibit any residual periodicity with a low FAP that could be attributed to a Keplerian signal.

Furthermore, we tested the coherence and consistency of the periodicities using the Stacked Bayesian generalized Lomb-Scargle (Stacked-BGLS) periodogram from Mortier & Collier Cameron (2017). This algorithm is designed to help disentangle short-lived incoherent signals produced by stellar activity from long-lived consistent Keplerian signals.

The Stacked-BGLS periodogram for the RVs from SOPHIE and SPIRou are shown together and independently in Fig. 13. We observed that there are no strong signals in the periodograms at periods of between 0 and 100 days. However, the activity close to 34 days is accompanied by other signals at close periods and shows some level of variable probability, as expected for stellar activity signals.

# 10. Discussion and conclusions

In this work, we present the long-term stellar activity and magnetic field characterization of the early, moderately active M dwarf Gl 205 using quasi-simultaneous optical and NIR high-resolution spectra from the SOPHIE/OHP and SPIRou/CFHT spectrographs. We used the SPIRou circularly polarized spectra to measure the Zeeman signature produced by the presence of the
<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Stacked Bayesian generalized Lomb–Scargle periodogram of the SOPHIE, SPIRou, and SOPHIE+SPIRou RVs. The periodogram of each observation is stacked along the y-axis. The color scale indicates the logarithm of the probability, where dark red indicates higher probability.</figcaption></figure>

magnetic field. The longitudinal magnetic field $B_ℓ$ is modulated by the stellar rotation with a period of $P_{\text{rot}} = 34.4 \pm 0.5$ days, which is in agreement with previous works (Kiraga & Stepien 2007; Bonfils et al. 2013; Hébrard et al. 2016). The comparison of the periodicities found in $B_ℓ$ with those suggested by various activity indicators reinforces the accuracy of $B_ℓ$ in constraining $P_{\text{rot}}$ compared to the standard activity indicators for early M dwarfs.

We applied the Zeeman Doppler imaging technique to reconstruct the large-scale magnetic field map of the star Gl 205 over three seasons of observations. We observe a temporal evolution of the magnetic field topology, while the strength of the large-scale field remains mostly constant. Moreover, we confirmed the differential rotation over the stellar surface of our target, which could explain the disagreement between its rotation period measured by previous studies and that suggested by different activity proxies.

We derived RVs using the CCF technique with the SOPHIE data and the line-by-line method with the SPIRou data. Both RV data sets are comparable in amplitude and scatter with variations due to stellar activity. We applied a quasi-periodic GPR in the SOPHIE and SPIRou RVs. The periodicities found in the RV data sets are in agreement with that of the longitudinal magnetic field $B_ℓ$.

We filtered the activity-induced RV variations using a multi-dimensional GPR framework. For the optical RVs, we used the $H\alpha$ index and the BIS as ancillary time series and, for the NIR, we used the FWHM<sub>LBL</sub>. The obtained models fit the RVs down to the noise level, with a scatter of $2.7 \ \text{m s}^{-1}$ in the residuals for both instruments. As we did not detect residual periodic signals, we can rule out a priori the presence of Keplerian signatures in the given data set. Moreover, we used TESS photometry to discard the presence of transit events in the two available sectors. Nevertheless, we cannot completely rule out the presence of planets around Gl 205 because of the flexibility of the GPs. At certain periods, in particular close to the stellar rotation period, the GP could still absorb the Keplerian signature in the RVs (Rajpaul et al. 2015).

## 10.1. Magnetic field temporal evolution

The long-term follow-up of more than two years of data presented in this work allows us to characterize the evolution of the large-scale magnetic field of Gl 205. The topology evolves significantly from a poloidal dipolar field with a strength of $\langle B \rangle = 12$ G in 2019 to a dominantly toroidal field in 2022 with a similar field strength of $\langle B \rangle = 12.3$ G. This change in the topology occurred quickly, that is, within two months, during our last set of SPIRou observations.

Previous studies of moderately active early M dwarf samples revealed significant diversity in magnetic field topologies (Donati et al. 2008; Hébrard et al. 2016; Martioli et al. 2022), suggesting that there is no unique topology pattern for these stars. This diversity is especially visible in partly convective, slowly rotating M dwarfs. Long-term follow-up campaigns of early M dwarfs, such as Gl 410 and Gl 846 (Donati et al. 2008; Hébrard et al. 2016), revealed a temporal evolution of the topology of the magnetic field. The variability of the magnetic field seen in early M dwarfs is at least partly attributable to their differential rotation over the stellar surface (Donati et al. 2008; Morin et al. 2008; Hébrard et al. 2016).

Although most of the seasons observed in this study and that of Hébrard et al. (2016) exhibit poloidal fields for the case of Gl 205, the fact that the field is dominantly toroidal during at least one of our seasons is unexpected, given that previous studies of partly convective stars (Petit et al. 2008; See et al. 2016) suggested that only fast rotators could generate toroidal dominated fields.

Our results of the large-scale magnetic field during the four seasons of $B_ℓ$ (see Table 4) are mostly consistent with the previous magnetic field study by Hébrard et al. (2016); however, in the present work we have been able to characterize the long-term evolution of the large-scale field. The long-term SPIRou monitoring was key to revealing the true nature of the magnetic field evolution and its intrinsic variability.

## 10.2. Diagnostics on optical and NIR activity indicators

Previous works using photometric, spectroscopic, and spectropolarimetric data for Gl 205 measured similar rotation period values ($P_{\text{rot}} \sim 33$ d). Nevertheless, the variability due to the stellar rotation are not always clear in all the activity indicator time series. Bonfils et al. (2013) showed that the strongest peaks at $P_{\text{rot}}$ of the $H\alpha$ and S index periodograms are seasonal, meaning that they become significant when the data set is restricted to one observation season of ~300 days. Hébrard et al. (2016) measured the periodicity of $B_ℓ$, RV, FWHM, and $H\alpha$, finding values of between 33.46 and 41.9 days, where the best result was given by the $B_ℓ$ with a rotation period of 33.63 ± 0.37 days. This discrepancy between activity indicators was explained by the presence of differential rotation over the stellar surface, but the HARPS-Pol and NARVAL data were not sufficient to confirm this rotation.

We studied stellar activity indicators of Gl 205 in the optical and NIR domains. In the optical domain, the $H\alpha$ index from SOPHIE exhibits a good correlation with the RVs and is modulated by the stellar rotation period as seen by Bonfils et al. (2013) in HARPS data. This is expected, because the $H\alpha$ index
is a well-studied activity indicator for M dwarfs and its correlation with stellar rotation is clear in large samples (e.g., Delfosse et al. 1998; Newton et al. 2017; Jeffers et al. 2018). In particular, the Carmencita sample of Jeffers et al. (2018) shows that $H\alpha$ in inactive (pEW($H\alpha$) > -0.5 Å) early M dwarfs, including Gl 205, tend to exhibit modulations at larger periods from 10 to 100 days.

The other optical indicators included in this work do not show clear signs of tracing activity. Even though most of them tend to have incipient peaks of periodicities close to the stellar rotation in their GLS periodograms, these peaks are not significant or do not persist for longer than one season. For example, the BIS becomes significantly modulated by the stellar rotation during the season $S_2^{RV}$. However, the periodicity is slightly longer than the rotation period, suggesting that the traced active features were closer to the poles. It is not surprising that the BIS is not anti-correlated with the RVs (Desort et al. 2007) in the case of Gl 205, because the SOPHIE resolution is larger than the $v\sin i$ of 0.7 km s⁻¹.

In the NIR domain, we find that the FWHM<sub>LBL</sub> is the best tracer of activity in our data set as it is correlated with the RV with a Pearson's coefficient of 0.4, and its periodogram has the strongest peak at the stellar rotation period. Previous studies found the same trend for other M dwarfs (e.g., Klein et al. 2021; Zicher et al. 2022). Moreover, we tested the performance of NIR spectral lines as possible activity tracers. The He I triplet at 10 833 Å has been extensively studied in the CARMENES sample of M dwarfs (Fuhrmeister et al. 2019, 2020) showing that its variability is correlated with $H\alpha$ only for active cases. In our work, we obtained large error bars from the fitting of the He I line, which is due to the large number of absorption lines around it, making it difficult to obtain a good Voigt profile. The index we derived does not appear to be modulated by the stellar rotation period.

The K line at 12 435 Å is studied in detail by Fuhrmeister et al. (2022). The authors find that the K I is rarely correlated or anti-correlated with $H\alpha$ because this line is less sensitive to chromospheric variability than $H\alpha$. Our results agree with the findings of Fuhrmeister et al. (2022) given that the K I line is not correlated with the RV variations and is not modulated by $P_{rot}$.

Only the Al I line at 13 154 Å shows several peaks of periodicity within the range of the expected differential rotation, suggesting that this line could be tracing several active features across the stellar latitudes. Even though the Ti I at 10 499 Å and Fe I at 11 693 Å lines are not modulated by the stellar rotation, they appear to be affected by the magnetic field, because they follow approximately the same changes in phase. Al I and Ti I have not previously been studied as activity tracers for Sun-like or M dwarf stars, and further analysis in a larger sample would be of great benefit to prove their capacity as activity tracers as they are located in relatively clear parts of the spectrum with few telluric lines.

Interestingly, the analysis of several activity indicators over time reveals changes in the periodicities over different seasons of observations, which may indicate an evolution in the distribution of active features over the stellar surface of Gl 205. All the activity indicators studied in this work reveal different information, suggesting they may be tracing various activity processes in the stellar surface, and demonstrate the complexity of the stellar activity of Gl 205.

As explained above, the lack of correlation or anti-correlation between the RVs and the activity indicators could

be the result of time lags between the time series. For that reason, we cannot completely discard the activity indicators that are not correlated with the RVs as activity tracers. Nevertheless, the existence of the correlation is evidence of stellar nature of the signal.

On the other hand, the GLS periodogram is the standard tool for detecting periodicities in time series; however, its efficiency depends on the number of data points, time span, and sampling, among other characteristics. The absence of a significant peak in periodicity by the periodogram does not directly imply that the indicator is not rotationally modulated. In this work, we compared the periodicities found in the GLS periodograms with the period of the quasi-periodic GP applied in all the activity indicators. In some cases where the GLS periodogram does not find periodicities in the dataset, the GP succeeded in constraining a periodicity close to the expected rotation period. However, in these particular cases, we know the expected stellar rotation and can therefore properly constrain the period of the GP. In blind searches of the stellar rotation period, the results may differ.

The periods found in the GP models are all within the ranges of the differential rotation of the star, namely 32–45.5 days, except for He I. In particular, these periods tend to be longer than the period measured in the longitudinal magnetic field. This may be a geometrical effect due to the inclination of Gl 205, which is 60 deg. From the point of view of the observer, there is one pole hidden, which means that we are more likely to observe active regions with longer rotation periods at greater latitudes.

Finding NIR spectral lines that could function as stellar activity tracers remains an open challenge for stellar characterization and exoplanet searches in this domain. This is particularly important for less active M dwarfs, because chromospheric lines may be better tracers of stellar activity than CCF-based activity indicators (Lafarga et al. 2021). A complementary route to face this challenge consists in identifying subsets of lines that are more or less sensitive to the activity jitter (Bellotti et al. 2022).

## 10.3. Optical and NIR activity-induced RV jitter

In this study, we confirm that the observed RV variations of Gl 205 in the optical and NIR are due to stellar activity, as previously stated. The RV periodicities over the three RV subdatasets and over the whole time series are always linked to one or more of the activity indicators in both domains.

Hébrard et al. (2016) defined the total RV jitter as the sum of the jitter due to the rotational modulation and the jitter originated from random components (e.g., short-lived spots). Although we are close to reaching the noise level of the data set with the multi-dimensional GP modeling, the remaining jitter of 2.7 ± 0.3 m s⁻¹ for SPIRou and 2.7 ± 0.6 m s⁻¹ for SOPHIE is equal to the random component contribution proposed for Gl 205, which is of the order of $J_f = 2.7$ m s⁻¹. This jitter contribution is greater than the average error bar of 1.9 m s⁻¹.

In the GP models applied to the RVs, we observe that the model of the optical RVs is consistent over time and resembles a sinusoid with its first harmonic. The NIR RVs, on the other hand, show a more complex behavior, with a rapidly evolving signal. The amplitudes measured from the GP model show a slightly larger value for the NIR RVs. Suárez Mascareño et al. (2018) measured the optical RV variation due to the rotational modulation and obtained a mean value of 4.4 m s⁻¹ for early M dwarfs and 5.6 m s⁻¹ if only targets from HADES are considered. Our results are in agreement with theirs; the SOPHIE
RV semi-amplitudes through all the seasons are between 4.3 and
$5.4\ \text{m s}^{-1}$<fnref n="11" />.

We find that the optical and NIR RVs of Gl 205 are similar
in amplitude and jitter level for the whole time series, mean-
ing that there is at least a small advantage to observing at NIR
wavelengths, which is expected for early, moderately active or
inactive M dwarfs. Reiners et al. (2010) showed that the gain
in precision toward the NIR is seen especially for M4 and later
spectral types. In the case of early M dwarfs, the best precision
is obtained in the V band, although the highest S/N is reached
in the J band. Moreover, this effect may not be exclusively
dependent on the spectral type only, but on levels of activity
as well.

Robertson et al. (2020) analyzed a small sample of fast-
rotating M dwarfs to compare their optical and NIR RVs. The
late M dwarf Gl 3959 exhibits comparable optical and NIR
RVs from the HIRES and HPF spectrographs, while some other
targets showed comparable RVs only during one season of obser-
vations. Active M dwarfs, such as Gl 388 and AU Mic, show
stronger RV jitter in the optical than in the NIR, while less active
stars, such as Gl 205, do not. The strength of the magnetic field
may play a key role in this phenomenon, because it likely impacts
the brightness contrast between magnetic and nonmagnetic fea-
tures. We show here that Gl 205 has a much weaker field than
those of active M dwarfs. We plan to extend the work of quasi-
simultaneous observations with SOPHIE and SPIRou to more
M dwarfs in order to help understand the multiwavelength stellar
activity contribution in RVs studies.

Acknowledgements. P.C.Z. thanks Oscar Barragán for the helpful discussions about pyaneti and the multi-dimensional GPs. This work is based on observations collected with the SOPHIE spectrograph on the 1.93 m telescope at the Observatoire de Haute-Provence (CNRS), France. We thank the staff of the Observatoire de Haute-Provence for their support at the 1.93 m telescope and on SOPHIE. This work is based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated from the summit of Maunakea by the National Research Council of Canada, the Institut National des Sciences de l'Univers of the Centre National de la Recherche Scientifique of France, and the University of Hawaii. The observations at the Canada-France-Hawaii Telescope were performed with care and respect from the summit of Maunakea which is a significant cultural and historic site. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of MaunaKea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. Based on observations obtained with SPIRou, an international project led by Institut de Recherche en Astrophysique et Planétologie, Toulouse, France. This paper includes data collected by the TESS mission that are publicly available from the Mikulski Archive for Space Telescopes (MAST). We acknowledge funding from the French National Research Agency (ANR) under contract number ANR-18-CE31-0019 (SPlaSH). P.C.Z. thanks the LSSTC Data Science Fellowship Program, which is funded by LSSTC, NSF Cybertraining Grant #1829740, the Brinson Foundation, and the Moore Foundation; her participation in the program has benefited this work. B.K. acknowledges funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (Grant agreement No. 865624, GPRV). A.A.K. acknowledges that this work was supported by FCT – Fundação para a Ciência e a Tecnologia through national funds and by FEDER through COMPETE2020 – Programa Operacional Competitividade e Internacionalização by these grants: UIDB/04434/2020 & UIDP/04434/2020; PTDC/FIS-AST/32113/2017 & POCI-01-0145-FEDER-032113; PTDC/FIS-AST/28953/2017 & POCI-01-0145-FEDER-028953; IF/00849/2015/CP1273/CT0003. J.F.D., C.M. and P.I.C. acknowledge funding from the European Research Council under the H2020 and innovation program (grant #740651 NewWorlds). X.D., X.B., T.F., and A.C. acknowledge a funding from the French National Research Agency in the framework of the Investissements d'Avenir program (ANR-15-IDEX-02), through the funding of the "Origin of Life" project of the Grenoble-Alpes University. S.H. acknowledges CNES funding through the grant 837319.

**References**

Aigrain, S., Pont, F., & Zucker, S. 2012, MNRAS, 419, 3147

Allende Prieto, C., Barklem, P. S., Lambert, D. L., & Cunha, K. 2004, A&A, 420, 183

Andersen, J. M., & Korhonen, H. 2015, MNRAS, 448, 3053

Angus, R., Morton, T., Aigrain, S., Foreman-Mackey, D., & Rajpaul, V. 2018, MNRAS, 474, 2094

Angus, R., Morton, T. D., Foreman-Mackey, D., et al. 2019, AJ, 158, 173

Antoniadis-Karnavas, A., Sousa, S. G., Delgado-Mena, E., et al. 2020, A&A, 636, A9

Artigau, É., Cadieux, C., Cook, N. J., et al. 2022, AJ, 164, 84

Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33

Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123

Astudillo-Defru, N., Bonfils, X., Delfosse, X., et al. 2015, A&A, 575, A119

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

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

Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42

Barrado y Navascués, D., García López, R. J., Severino, G., & Gomez, M. T. 2001, A&A, 371, 652

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

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

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

Bellotti, S., Petit, P., Morin, J., et al. 2022, A&A, 657, A107

Bertaux, J. L., Lallement, R., Ferron, S., Boonne, C., & Bodichon, R. 2014, A&A, 564, A46

Boisse, I., Moutou, C., Vidal-Madjar, A., et al. 2009, A&A, 495, 959

Boisse, I., Bouchy, F., Hébrard, G., et al. 2011, A&A, 528, A4

Bonfils, X., Mayor, M., Delfosse, X., et al. 2007, A&A, 474, 293

Bonfils, X., Delfosse, X., Udry, S., et al. 2013, A&A, 549, A109

Bortle, A., Fausey, H., Ji, J., et al. 2021, AJ, 161, 230

Bouchy, F., Pepe, F., & Queloz, D. 2001, A&A, 374, 733

Bouchy, F., Hébrard, G., Udry, S., et al. 2009a, A&A, 505, 853

Bouchy, F., Isambert, J., Lovis, C., et al. 2009b, in EAS Publications Series, 37, ed. P. Kern, 247

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

Brown, S. F., Donati, J. F., Rees, D. E., & Semel, M. 1991, A&A, 250, 463

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

Carmona, A., Delfosse, X., Bellotti, S., et al. 2023, A&A, submitted [arXiv:2303.16712]

Carolo, E., Desidera, S., Gratton, R., et al. 2014, A&A, 567, A48

Claret, A., & Bloemen, S. 2011, A&A, 529, A75

Collier Cameron, A., Mortier, A., Phillips, D., et al. 2019, MNRAS, 487, 1082

Cook, N. J., Artigau, É., Doyon, R., et al. 2022, PASP, 134, 114509

Courcol, B., Bouchy, F., Pepe, F., et al. 2015, A&A, 581, A38

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

Cretignier, M., Dumusque, X., Allart, R., Pepe, F., & Lovis, C. 2020, A&A, 633, A76

Cretignier, M., Dumusque, X., & Pepe, F. 2022, A&A, 659, A68

Cristofari, P. I., Donati, J. F., Masseron, T., et al. 2022a, MNRAS, 516, 3802

Cristofari, P. I., Donati, J. F., Masseron, T., et al. 2022b, MNRAS, 511, 1893

Davenport, J. R. A., Hawley, S. L., Hebb, L., et al. 2014, ApJ, 797, 122

Davis, A. B., Cisewski, J., Dumusque, X., Fischer, D. A., & Ford, E. B. 2017, ApJ, 846, 59

Delfosse, X., Forveille, T., Perrier, C., & Mayor, M. 1998, A&A, 331, 581

Delisle, J. B., Unger, N., Hara, N. C., & Ségransan, D. 2022, A&A, 659, A182

Delrez, L., Gillon, M., Queloz, D., et al. 2018, SPIE Conf. Ser., 10700, 1070011

Desidera, S., Gratton, R. G., Endl, M., et al. 2004, A&A, 420, L27

Desort, M., Lagrange, A. M., Galland, F., Udry, S., & Mayor, M. 2007, A&A, 473, 983

Díaz, R. F., Delfosse, X., Hobson, M. J., et al. 2019, A&A, 625, A17

Donati, J. F., & Brown, S. F. 1997, A&A, 326, 1135

Donati, J. F., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, MNRAS, 291, 658

Donati, J. F., Mengel, M., Carter, B. D., et al. 2000, MNRAS, 316, 699

Donati, J. F., Howarth, I. D., Jardine, M. M., et al. 2006, MNRAS, 370, 629

Donati, J. F., Morin, J., Petit, P., et al. 2008, MNRAS, 390, 545

Donati, J. F., Kouach, D., Moutou, C., et al. 2020, MNRAS, 498, 5684

Dressing, C. D., & Charbonneau, D. 2013, ApJ, 767, 95

Dressing, C. D., & Charbonneau, D. 2015, ApJ, 807, 45

Ducati, J. R. 2002, VizieR Online Data Catalog

<footnote n="11">¹¹ We approximate the semi-amplitude as half of the amplitudes mea-
sured in Table 6.</footnote>
Dumusque, X., Boisse, I., & Santos, N. C. 2014a, ApJ, 796, 132

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

Dungee, R., van Saders, J., Gaidos, E., et al. 2022, ApJ, 938, 118

Faria, J. P., Suárez Mascareño, A., Figueira, P., et al. 2022, A&A, 658, A115

Feinstein, A., Montet, B., & Ansdell, M. 2020a, J. Open Source Softw., 5, 2347

Feinstein, A. D., Montet, B. T., Ansdell, M., et al. 2020b, AJ, 160, 219

Finociety, B., & Donati, J. F. 2022, MNRAS, 516, 5887

Folsom, C. P., Petit, P., Bouvier, J., et al. 2016, MNRAS, 457, 580

Foreman-Mackey, D. 2016, J. Open Source Softw., 1, 24

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306

Fouqué, P., Martioli, E., Donati, J. F., et al. 2023, ApJL, 927, L11

Fuhrmeister, B., Czesla, S., Hildebrandt, L., et al. 2019, A&A, 632, A24

Fuhrmeister, B., Czesla, S., Hildebrandt, L., et al. 2020, A&A, 640, A52

Fuhrmeister, B., Czesla, S., Nagel, E., et al. 2022, A&A, 657, A125

Gaia Collaboration. 2020, VizieR Online Data Catalog: I/350

García-Mejia, J., Charbonneau, D., Fabricant, D., et al. 2020, SPIE Conf. Ser., 11445, 114457R

Gilbertson, C., Ford, E. B., Jones, D. E., & Stenning, D. C. 2020, ApJ, 905, 155

Gillon, M., Jehin, E., Magain, P., et al. 2011, in European Physical Journal Web of Conferences, 11, 06002

Gomes da Silva, J., Santos, N. C., Bonfils, X., et al. 2011, A&A, 534, A30

Gray, D. F. 2008, The Observation and Analysis of Stellar Photospheres (Cambridge, UK: Cambridge University Press)

Günther, M. N., Zhan, Z., Seager, S., et al. 2020, AJ, 159, 60

Haywood, M. 2008, MNRAS, 388, 1175

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

Hébrard, É. M., Donati, J. F., Delfosse, X., et al. 2014, MNRAS, 443, 2599

Hébrard, É. M., Donati, J. F., Delfosse, X., et al. 2016, MNRAS, 461, 1465

Henry, T. J., Jao, W.-C., Subasavage, J. P., et al. 2006, AJ, 132, 2360

Hippke, M., & Heller, R. 2019, A&A, 623, A39

Hippke, M., David, T. J., Mulders, G. D., & Heller, R. 2019, AJ, 158, 143

Hobson, M. J., Díaz, R. F., Delfosse, X., et al. 2018, A&A, 618, A103

Hobson, M. J., Delfosse, X., Astudillo-Defru, N., et al. 2019, A&A, 625, A18

Hobson, M. J., Bouchy, F., Cook, N. J., et al. 2021, A&A, 648, A48

Holmberg, J., Nordström, B., & Andersen, J. 2009, A&A, 501, 941

Horne, K. 1986, PASP, 98, 609

Huélamo, N., Figueira, P., Bonfils, X., et al. 2008, A&A, 489, L9

Jeffers, S. V., Schöfer, P., Lamert, A., et al. 2018, A&A, 614, A76

Jones, D. E., Stenning, D. C., Ford, E. B., et al. 2017, ArXiv e-prints [arXiv:1711.01318]

Kanodia, S., & Wright, J. 2018, RNAAS, 2, 4

Khata, D., Mondal, S., Das, R., & Baug, T. 2021, MNRAS, 507, 1869

Kiraga, M., & Stepien, K. 2007, Acta Astron., 57, 149

Klein, B., & Donati, J. F. 2020, MNRAS, 493, L92

Klein, B., Donati, J.-F., Moutou, C., et al. 2021, MNRAS, 502, 188

Klein, B., Zicher, N., Kavanagh, R. D., et al. 2016, MNRAS, 512, 5067

Koen, C., Kilkenny, D., van Wyk, F., & Marang, F. 2010, MNRAS, 403, 1949

Kopparapu, R. K. 2013, ApJ, 767, L8

Kotani, T., Tamura, M., Suto, H., et al. 2014, SPIE Conf. Ser., 9147, 91474

Kürster, M., Endl, M., Rouesnel, F., et al. 2003, A&A, 403, 1077

Kurucz, R. L. 1993, SYNTHETHE spectrum synthesis programs and line data (Cambridge, MA: Smithsonian Astrophysical Observatory)

Lafarga, M., Ribas, I., Reiners, A., et al. 2021, A&A, 652, A28

Lomb, N. R. 1976, Ap&SS, 39, 447

Mahadevan, S., Ramsey, L., Bender, C., et al. 2012, SPIE Conf. Ser., 8446, 84461S

Mahmud, N. I., Crockett, C. J., Johns-Krull, C. M., et al. 2011, ApJ, 736, 123

Maldonado, J., Affer, L., Micela, G., et al. 2015, A&A, 577, A132

Martín, E. L., Guenther, E., Zapatero Osorio, M. R., Bouy, H., & Wainscoat, R. 2006, ApJ, 644, L75

Martioli, E., Hébrard, G., Fouqué, P., et al. 2022, A&A, 660, A86

Morin, J., Donati, J. F., Petit, P., et al. 2008, MNRAS, 390, 567

Morin, J., Donati, J. F., Petit, P., et al. 2010, MNRAS, 407, 2269

Morley, C. V., Kreidberg, L., Rustamkulov, Z., Robinson, T., & Fortney, J. J. 2017, ApJ, 850, 121

Mortier, A., & Collier Cameron, A. 2017, A&A, 601, A110

Neves, V., Bonfils, X., Santos, N. C., et al. 2012, A&A, 538, A25

Neves, V., Bonfils, X., Santos, N. C., et al. 2014, A&A, 568, A121

Newton, E. R., Irwin, J., Charbonneau, D., et al. 2016, ApJ, 834, 85

Nicholson, B. A., & Aigrain, S. 2022, MNRAS, 515, 5251

Nutzman, P., & Charbonneau, D. 2008, PASP, 120, 317

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

Perruchot, S., Kohler, D., Bouchy, F., et al. 2008, SPIE Conf. Ser., 7014, 70140J

Petit, P., Donati, J. F., & Collier Cameron, A. 2002, MNRAS, 334, 374

Petit, P., Dintrans, B., Solanki, S. K., et al. 2008, MNRAS, 388, 80

Pinamonti, M., Sozzetti, A., Maldonado, J., et al. 2022, A&A, 664, A65

Queloz, D., Henry, G. W., Sivan, J. P., et al. 2001, A&A, 379, 279

Quirrenbach, A., Amado, P. J., Ribas, I., et al. 2018, SPIE Conf. Ser., 10702, 107020W

Rabus, M., Lachaume, R., Jordán, A., et al. 2019, MNRAS, 484, 2674

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

Reddy, B. E., Lambert, D. L., & Allende Prieto, C. 2006, MNRAS, 367, 1329

Reiners, A. 2012, Living Rev. Solar Phys., 9, 1

Reiners, A., Bean, J. L., Huber, K. F., et al. 2010, ApJ, 710, 432

Reiners, A., Shulyak, D., Anglada-Escudé, G., et al. 2013, A&A, 552, A103

Reylé, C., Jardine, K., Fouqué, P., et al. 2021, A&A, 650, A201

Ricker, G. R. 2014, J. Am. Assoc. Variable Star Observers (JAAVSO), 42, 234

Roberts, S., Osborne, M., Ebden, M., et al. 2012, Philos. Trans. Roy. Soc. Lond. A, 371, 20110550

Robertson, P., Bender, C., Mahadevan, S., Roy, A., & Ramsey, L. W. 2016, ApJ, 832, L12

Robertson, P., Stefansson, G., Mahadevan, S., et al. 2020, ApJ, 897, 125

Sabotta, S., Schlecker, M., Chaturvedi, P., et al. 2021, A&A, 653, A114

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

Scargle, J. D. 1982, ApJ, 263, 835

Schöfer, P., Jeffers, S. V., Reiners, A., et al. 2019, A&A, 623, A44

Schweitzer, A., Passegger, V. M., Cifuentes, C., et al. 2019, A&A, 625, A68

See, V., Jardine, M., Vidotto, A. A., et al. 2016, MNRAS, 462, 4442

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

Semel, M. 1989, A&A, 225, 456

Skillling, J., & Bryan, R. K. 1984, MNRAS, 211, 111

Souto, D., Cunha, K., Smith, V. V., et al. 2022, ApJ, 927, 123

Spina, L., Nordlander, T., Casey, A. R., et al. 2020, ApJ, 895, 52

Suárez Mascareño, A., Rebolo, R., González Hernández, J. I., et al. 2018, A&A, 612, A89

Suárez Mascareño, A., Faria, J. P., Figueira, P., et al. 2020, A&A, 639, A77

Terrien, R. C., Keen, A., Oda, K., et al. 2022, ApJ, 927, L11

Unno, W. 1956, PASJ, 8, 108

Walkowicz, L. M., Basri, G., Batalha, N., et al. 2011, AJ, 141, 50

Wilson, O. C. 1968, ApJ, 153, 221

Yana Galarza, J., Meléndez, J., Lorenzo-Oliveira, D., et al. 2019, MNRAS, 490, L86

Yu, L., Donati, J. F., Grankin, K., et al. 2019, MNRAS, 489, 5556

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

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

Zicher, N., Barragán, O., Klein, B., et al. 2022, MNRAS, 512, 3060

<footnote n="1">¹ Aix Marseille Univ, CNRS, CNES, LAM, Marseille, France e-mail: pia.cortes@lam.fr</footnote>

<footnote n="2">² LSSTC DSFP Fellow</footnote>

<footnote n="3">³ Sub-department of Astrophysics, Department of Physics, University of Oxford, Oxford OX1 3RH, UK</footnote>

<footnote n="4">⁴ Laboratório Nacional de Astrofísica, Rua Estados Unidos 154, 37504-364 Itajubá, MG, Brazil</footnote>

<footnote n="5">⁵ Univ. de Toulouse, CNRS, IRAP, 14 av. Belin, 31400 Toulouse, France</footnote>

<footnote n="6">⁶ Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, 4150-762 Porto, Portugal</footnote>

<footnote n="7">⁷ Departamento de Física e Astronomia, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre, 4169-007 Porto, Portugal</footnote>

<footnote n="8">⁸ Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France</footnote>

<footnote n="9">⁹ Université de Montréal, Département de Physique, IREX, Montréal, QC H3C 3J7, Canada</footnote>

<footnote n="10">¹⁰ International Center for Advanced Studies (ICAS) and ICIFI (CONICET), ECyT-UNSAM, Campus Miguelete, 25 de Mayo y Francia, (1650), Buenos Aires, Argentina</footnote>

<footnote n="11">¹¹ Science Division, Directorate of Science, European Space Research and Technology Centre (ESA/ESTEC), Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands</footnote>
<footnote n="12">12 Institut d'Astrophysique de Paris, CNRS, UMR 7095, Sorbonne Université, 98 bis bd Arago, 75014 Paris, France</footnote>

<footnote n="13">13 Observatoire du Mont-Mégantic, Université de Montréal, Montréal, QC H3C 3J7, Canada</footnote>

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

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

<footnote n="16">16 Laboratoire J.-L. Lagrange, Observatoire de la Côte d'Azur (OCA), Universite de Nice-Sophia Antipolis (UNS), CNRS, Campus Valrose, 06108 Nice Cedex 2, France</footnote>

<footnote n="17">17 Department of Physics & Astronomy, McMaster University, 1280 Main Street West, Hamilton, ON L8S 4K1, Canada</footnote>

<footnote n="18">18 LUPM, Université de Montpellier, CNRS, Place Eugène Bataillon, 34095 Montpellier, France</footnote>

<footnote n="19">19 Canada-France-Hawaii Telescope, CNRS, Kamuela, HI 96743, USA</footnote>

<footnote n="20">20 Observatoire Astronomique de l'Université de Genève, Chemin Pegasi 51b, 1290 Versoix, Switzerland</footnote>

<footnote n="21">21 LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, 5 place Jules Janssen, 92195 Meudon, France</footnote>

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

<footnote n="23">23 Astrophysics Group, University of Exeter, Exeter EX4 2QL, UK</footnote>
Appendix A: RVs and activity indicators

The log of the SOPHIE and SPIRou observations with their activity indicators are listed in Tables A.1 and A.2, respectively.

<table><caption>Table A.1. SOPHIE RVs and activity indicators. The values of the RVs are given as the variation from the mean RV = 8727.9 m s⁻¹.</caption><thead><tr><th>BJD-2400000 d</th><th>ΔRV m s<sup>-1</sup></th><th>σ<sub>RV</sub> m s<sup>-1</sup></th><th>FWHM km s<sup>-1</sup></th><th>contrast %CCF</th><th>BIS m s<sup>-1</sup></th><th>σ<sub>BIS</sub> m s<sup>-1</sup></th><th>S<sub>MW</sub> -</th><th>σ<sub>S</sub> -</th><th>H<sub>α</sub> -</th><th>σ<sub>Hα</sub> -</th></tr></thead><tbody><tr><td>58794.5891</td><td>-1.5</td><td>2.1</td><td>4.77</td><td>21.87</td><td>13.07</td><td>2.87</td><td>1.61</td><td>0.02</td><td>0.2519</td><td>0.0007</td></tr><tr><td>58796.5420</td><td>4.7</td><td>2.3</td><td>4.80</td><td>21.16</td><td>8.53</td><td>3.14</td><td>1.63</td><td>0.01</td><td>0.2499</td><td>0.0009</td></tr><tr><td>58805.5787</td><td>-6.9</td><td>2.5</td><td>4.80</td><td>21.53</td><td>9.48</td><td>3.40</td><td>1.62</td><td>0.03</td><td>0.2659</td><td>0.0010</td></tr><tr><td>58820.5323</td><td>-7.8</td><td>1.7</td><td>4.75</td><td>21.14</td><td>10.76</td><td>2.29</td><td>1.54</td><td>0.02</td><td>0.2457</td><td>0.0006</td></tr><tr><td>58821.5421</td><td>-9.2</td><td>2.4</td><td>4.82</td><td>20.87</td><td>9.65</td><td>3.30</td><td>1.66</td><td>0.03</td><td>0.2445</td><td>0.0009</td></tr><tr><td>58825.5023</td><td>-6.0</td><td>1.8</td><td>4.79</td><td>21.38</td><td>10.43</td><td>2.33</td><td>1.59</td><td>0.02</td><td>0.2522</td><td>0.0007</td></tr><tr><td>58832.5041</td><td>0.7</td><td>1.6</td><td>4.70</td><td>20.93</td><td>14.07</td><td>2.10</td><td>1.52</td><td>0.02</td><td>0.2514</td><td>0.0005</td></tr><tr><td>58851.4754</td><td>-4.0</td><td>2.8</td><td>4.83</td><td>20.24</td><td>24.34</td><td>3.81</td><td>1.49</td><td>0.03</td><td>0.2539</td><td>0.0011</td></tr><tr><td>58854.4361</td><td>-3.0</td><td>1.7</td><td>4.72</td><td>20.96</td><td>10.73</td><td>2.27</td><td>1.51</td><td>0.02</td><td>0.2461</td><td>0.0006</td></tr><tr><td>58859.4193</td><td>-9.0</td><td>2.1</td><td>4.81</td><td>21.25</td><td>6.20</td><td>2.84</td><td>1.43</td><td>0.02</td><td>0.2394</td><td>0.0008</td></tr><tr><td>58860.4112</td><td>-4.8</td><td>1.8</td><td>4.76</td><td>21.25</td><td>5.56</td><td>2.37</td><td>1.48</td><td>0.02</td><td>0.2401</td><td>0.0007</td></tr><tr><td>58878.3905</td><td>3.1</td><td>2.5</td><td>4.83</td><td>20.56</td><td>14.43</td><td>3.42</td><td>1.71</td><td>0.03</td><td>0.2745</td><td>0.0010</td></tr><tr><td>58892.3602</td><td>-3.0</td><td>2.3</td><td>4.81</td><td>20.95</td><td>8.26</td><td>3.12</td><td>1.66</td><td>0.03</td><td>0.2559</td><td>0.0009</td></tr><tr><td>58894.3182</td><td>-1.4</td><td>1.6</td><td>4.74</td><td>21.36</td><td>6.25</td><td>2.10</td><td>1.59</td><td>0.02</td><td>0.2434</td><td>0.0006</td></tr><tr><td>58895.2799</td><td>3.7</td><td>1.7</td><td>4.71</td><td>20.99</td><td>11.49</td><td>2.23</td><td>1.63</td><td>0.03</td><td>0.2456</td><td>0.0006</td></tr><tr><td>58903.2928</td><td>0.6</td><td>1.8</td><td>4.79</td><td>21.38</td><td>9.51</td><td>2.38</td><td>1.68</td><td>0.02</td><td>0.2521</td><td>0.0007</td></tr><tr><td>58904.3100</td><td>-0.7</td><td>1.7</td><td>4.76</td><td>21.27</td><td>9.18</td><td>2.23</td><td>2.02</td><td>0.02</td><td>0.2593</td><td>0.0006</td></tr><tr><td>58905.3351</td><td>-1.0</td><td>2.7</td><td>4.86</td><td>20.87</td><td>9.21</td><td>3.71</td><td>1.67</td><td>0.03</td><td>0.2571</td><td>0.0012</td></tr><tr><td>58911.2803</td><td>9.1</td><td>2.8</td><td>4.87</td><td>20.62</td><td>12.12</td><td>3.89</td><td>1.71</td><td>0.02</td><td>0.2629</td><td>0.0012</td></tr><tr><td>58920.2862</td><td>-6.6</td><td>1.7</td><td>4.77</td><td>21.35</td><td>19.21</td><td>2.24</td><td>1.60</td><td>0.02</td><td>0.2732</td><td>0.0006</td></tr><tr><td>58921.3072</td><td>-10.1</td><td>2.6</td><td>4.79</td><td>20.38</td><td>9.59</td><td>3.56</td><td>1.64</td><td>0.03</td><td>0.2707</td><td>0.0010</td></tr><tr><td>58922.2891</td><td>-4.1</td><td>2.3</td><td>4.84</td><td>21.12</td><td>17.30</td><td>3.11</td><td>1.71</td><td>0.02</td><td>0.2680</td><td>0.0009</td></tr><tr><td>59147.5761</td><td>3.2</td><td>1.6</td><td>4.78</td><td>21.84</td><td>16.70</td><td>2.15</td><td>1.67</td><td>0.02</td><td>0.2645</td><td>0.0005</td></tr><tr><td>59149.5806</td><td>-0.7</td><td>2.0</td><td>4.87</td><td>22.49</td><td>11.77</td><td>2.70</td><td>1.85</td><td>0.02</td><td>0.2664</td><td>0.0008</td></tr><tr><td>59152.5861</td><td>-4.6</td><td>1.7</td><td>4.83</td><td>22.35</td><td>12.68</td><td>2.27</td><td>1.70</td><td>0.01</td><td>0.2649</td><td>0.0006</td></tr><tr><td>59154.5745</td><td>-4.4</td><td>1.8</td><td>4.82</td><td>22.11</td><td>13.54</td><td>2.37</td><td>1.77</td><td>0.02</td><td>0.2658</td><td>0.0007</td></tr><tr><td>59159.6351</td><td>-3.3</td><td>1.6</td><td>4.76</td><td>22.02</td><td>9.15</td><td>2.18</td><td>1.63</td><td>0.02</td><td>0.2481</td><td>0.0005</td></tr><tr><td>59171.5482</td><td>6.6</td><td>1.6</td><td>4.80</td><td>22.13</td><td>11.01</td><td>2.15</td><td>1.74</td><td>0.02</td><td>0.2688</td><td>0.0006</td></tr><tr><td>59186.5302</td><td>-1.5</td><td>1.8</td><td>4.84</td><td>22.13</td><td>16.25</td><td>2.36</td><td>1.72</td><td>0.02</td><td>0.2665</td><td>0.0007</td></tr><tr><td>59194.4799</td><td>-4.1</td><td>1.4</td><td>4.81</td><td>21.93</td><td>4.33</td><td>1.80</td><td>1.82</td><td>0.01</td><td>0.2718</td><td>0.0005</td></tr><tr><td>59196.4731</td><td>1.4</td><td>2.5</td><td>4.89</td><td>20.80</td><td>9.19</td><td>3.43</td><td>1.72</td><td>0.02</td><td>0.2668</td><td>0.0010</td></tr><tr><td>59206.4607</td><td>-2.8</td><td>3.4</td><td>4.86</td><td>19.43</td><td>5.55</td><td>4.69</td><td>1.60</td><td>0.04</td><td>0.2788</td><td>0.0013</td></tr><tr><td>59240.3519</td><td>-3.3</td><td>2.3</td><td>4.94</td><td>21.16</td><td>5.49</td><td>3.14</td><td>1.80</td><td>0.01</td><td>0.2750</td><td>0.0010</td></tr><tr><td>59245.4175</td><td>1.8</td><td>2.8</td><td>4.87</td><td>20.12</td><td>13.67</td><td>3.90</td><td>1.99</td><td>0.03</td><td>0.2948</td><td>0.0011</td></tr><tr><td>59249.3216</td><td>8.7</td><td>1.6</td><td>4.84</td><td>21.50</td><td>10.24</td><td>2.15</td><td>1.96</td><td>0.01</td><td>0.2875</td><td>0.0007</td></tr><tr><td>59263.3251</td><td>-2.7</td><td>1.4</td><td>4.82</td><td>21.99</td><td>15.39</td><td>1.72</td><td>1.72</td><td>0.01</td><td>0.2498</td><td>0.0005</td></tr><tr><td>59271.3509</td><td>5.4</td><td>1.4</td><td>4.83</td><td>21.90</td><td>10.31</td><td>1.77</td><td>1.62</td><td>0.01</td><td>0.2656</td><td>0.0005</td></tr><tr><td>59272.3046</td><td>5.8</td><td>1.7</td><td>4.83</td><td>21.27</td><td>10.45</td><td>2.25</td><td>1.72</td><td>0.02</td><td>0.2718</td><td>0.0007</td></tr><tr><td>59279.3161</td><td>4.0</td><td>1.6</td><td>4.82</td><td>21.31</td><td>6.78</td><td>2.14</td><td>1.89</td><td>0.02</td><td>0.2859</td><td>0.0006</td></tr><tr><td>59475.6741</td><td>0.8</td><td>1.6</td><td>4.89</td><td>22.60</td><td>15.41</td><td>2.09</td><td>1.90</td><td>0.01</td><td>0.2731</td><td>0.0006</td></tr><tr><td>59486.6681</td><td>8.8</td><td>1.6</td><td>4.91</td><td>22.31</td><td>11.11</td><td>2.09</td><td>2.07</td><td>0.01</td><td>0.2984</td><td>0.0006</td></tr><tr><td>59488.6528</td><td>3.8</td><td>1.4</td><td>4.90</td><td>22.42</td><td>15.67</td><td>1.87</td><td>2.11</td><td>0.01</td><td>0.3043</td><td>0.0005</td></tr><tr><td>59500.6615</td><td>0.4</td><td>1.9</td><td>4.87</td><td>22.00</td><td>10.30</td><td>2.54</td><td>1.70</td><td>0.01</td><td>0.2534</td><td>0.0008</td></tr><tr><td>59502.6818</td><td>8.9</td><td>1.6</td><td>4.85</td><td>22.28</td><td>12.14</td><td>2.12</td><td>1.72</td><td>0.01</td><td>0.2512</td><td>0.0006</td></tr><tr><td>59505.5964</td><td>3.1</td><td>1.6</td><td>4.84</td><td>22.15</td><td>13.09</td><td>2.12</td><td>1.85</td><td>0.01</td><td>0.2820</td><td>0.0006</td></tr><tr><td>59506.6446</td><td>3.8</td><td>1.5</td><td>4.85</td><td>22.45</td><td>18.09</td><td>1.98</td><td>1.71</td><td>0.01</td><td>0.2632</td><td>0.0006</td></tr><tr><td>59511.6283</td><td>-0.6</td><td>1.8</td><td>4.86</td><td>22.02</td><td>13.26</td><td>2.35</td><td>1.76</td><td>0.01</td><td>0.2600</td><td>0.0007</td></tr><tr><td>59513.6577</td><td>2.5</td><td>1.9</td><td>4.87</td><td>21.83</td><td>7.76</td><td>2.61</td><td>1.70</td><td>0.01</td><td>0.2606</td><td>0.0008</td></tr><tr><td>59525.5697</td><td>1.3</td><td>1.8</td><td>4.89</td><td>21.76</td><td>10.40</td><td>2.45</td><td>1.85</td><td>0.01</td><td>0.2800</td><td>0.0008</td></tr><tr><td>59532.5356</td><td>-0.5</td><td>2.0</td><td>4.87</td><td>21.86</td><td>17.64</td><td>2.65</td><td>1.88</td><td>0.01</td><td>0.2808</td><td>0.0008</td></tr><tr><td>59536.5829</td><td>-2.0</td><td>1.5</td><td>4.81</td><td>22.22</td><td>13.27</td><td>1.98</td><td>1.71</td><td>0.01</td><td>0.2598</td><td>0.0008</td></tr><tr><td>59537.4909</td><td>-0.2</td><td>1.6</td><td>4.80</td><td>21.84</td><td>12.61</td><td>2.19</td><td>1.62</td><td>0.01</td><td>0.2499</td><td>0.0005</td></tr><tr><td>59538.5451</td><td>3.0</td><td>1.5</td><td>4.80</td><td>22.21</td><td>12.87</td><td>1.92</td><td>1.64</td><td>0.01</td><td>0.2668</td><td>0.0005</td></tr><tr><td>59547.5166</td><td>-6.2</td><td>2.7</td><td>4.88</td><td>20.98</td><td>12.45</td><td>3.68</td><td>1.77</td><td>0.02</td><td>0.2550</td><td>0.0011</td></tr><tr><td>59554.5152</td><td>9.0</td><td>2.0</td><td>4.87</td><td>21.38</td><td>12.03</td><td>2.65</td><td>1.85</td><td>0.02</td><td>0.2845</td><td>0.0008</td></tr><tr><td>59562.5315</td><td>-6.0</td><td>1.5</td><td>4.75</td><td>21.09</td><td>14.40</td><td>1.96</td><td>1.64</td><td>0.02</td><td>0.2650</td><td>0.0005</td></tr><tr><td>59587.4143</td><td>4.6</td><td>2.3</td><td>4.87</td><td>20.92</td><td>8.97</td><td>3.09</td><td>1.55</td><td>0.02</td><td>0.2683</td><td>0.0010</td></tr><tr><td>59594.3647</td><td>4.1</td><td>1.4</td><td>4.77</td><td>21.41</td><td>9.97</td><td>1.83</td><td>1.59</td><td>0.02</td><td>0.2686</td><td>0.0005</td></tr><tr><td>59595.3839</td><td>-1.8</td><td>1.6</td><td>4.81</td><td>21.35</td><td>6.96</td><td>2.18</td><td>1.68</td><td>0.02</td><td>0.2671</td><td>0.0006</td></tr><tr><td>59597.4072</td><td>3.0</td><td>1.6</td><td>4.80</td><td>21.54</td><td>15.26</td><td>2.04</td><td>1.66</td><td>0.01</td><td>0.2644</td><td>0.0006</td></tr><tr><td>59601.3991</td><td>3.9</td><td>1.9</td><td>4.83</td><td>21.42</td><td>12.38</td><td>2.59</td><td>1.75</td><td>0.02</td><td>0.2634</td><td>0.0008</td></tr><tr><td>59607.3530</td><td>3.7</td><td>1.6</td><td>4.79</td><td>21.67</td><td>9.40</td><td>2.04</td><td>1.69</td><td>0.02</td><td>0.2641</td><td>0.0006</td></tr></tbody></table>
<table><caption>**Table A.2.** SPIRou RVs and activity indicators. The values of the RVs are given as the variation from the mean RV = 9229.7 m s⁻¹.</caption><thead><tr><th>BJD-2400000 d</th><th>ΔRV m s<sup>-1</sup></th><th>σ<sub>RV</sub> m s<sup>-1</sup></th><th>FWHM<sub>LBL</sub> km s<sup>-1</sup></th><th>σ<sub>FWHM</sub> km s<sup>-1</sup></th><th>chrom. velocity slope m s<sup>-1</sup></th><th>σ<sub>chrom.</sub> m s<sup>-1</sup></th><th>B<sub>ℓ</sub> G</th><th>σ<sub>B</sub> G</th><th>Phase</th><th>Season</th></tr></thead><tbody><tr><td>58738.1365*</td><td>-1.3</td><td>1.6</td><td>6.006</td><td>0.009</td><td>-16.0</td><td>3.0</td><td>5.70</td><td>1.19</td><td>0.00</td><td>1</td></tr><tr><td>58745.1192</td><td>2.2</td><td>2.2</td><td>5.941</td><td>0.009</td><td>-8.1</td><td>4.6</td><td>0.06</td><td>1.16</td><td>0.20</td><td>1</td></tr><tr><td>58746.0224</td><td>7.8</td><td>2.2</td><td>5.946</td><td>0.010</td><td>-9.7</td><td>4.7</td><td>-1.18</td><td>1.33</td><td>0.23</td><td>1</td></tr><tr><td>58751.1257</td><td>-3.6</td><td>1.9</td><td>5.951</td><td>0.008</td><td>-18.5</td><td>4.0</td><td>-4.12</td><td>1.14</td><td>0.38</td><td>1</td></tr><tr><td>58752.1188</td><td>-5.6</td><td>2.1</td><td>5.929</td><td>0.009</td><td>-17.7</td><td>4.5</td><td>1.28</td><td>1.76</td><td>0.41</td><td>1</td></tr><tr><td>58753.1103</td><td>1.3</td><td>2.1</td><td>5.972</td><td>0.009</td><td>-19.9</td><td>4.1</td><td>0.07</td><td>1.22</td><td>0.44</td><td>1</td></tr><tr><td>58760.1240</td><td>-1.1</td><td>1.9</td><td>5.961</td><td>0.008</td><td>-26.3</td><td>3.9</td><td>5.55</td><td>1.04</td><td>0.64</td><td>1</td></tr><tr><td>58762.1541</td><td>5.6</td><td>2.1</td><td>5.928</td><td>0.009</td><td>-19.2</td><td>4.5</td><td>5.61</td><td>1.30</td><td>0.70</td><td>1</td></tr><tr><td>58763.1257</td><td>3.1</td><td>2.1</td><td>5.942</td><td>0.009</td><td>-21.7</td><td>4.4</td><td>4.51</td><td>1.18</td><td>0.73</td><td>1</td></tr><tr><td>58766.1035</td><td>-1.8</td><td>2.3</td><td>5.905</td><td>0.010</td><td>-9.9</td><td>4.8</td><td>2.23</td><td>1.21</td><td>0.81</td><td>1</td></tr><tr><td>58769.0896</td><td>-0.2</td><td>2.1</td><td>5.935</td><td>0.009</td><td>-4.4</td><td>4.5</td><td>5.79</td><td>1.36</td><td>0.90</td><td>1</td></tr><tr><td>58770.1405</td><td>0.3</td><td>2.0</td><td>5.945</td><td>0.009</td><td>-10.7</td><td>4.3</td><td>7.59</td><td>1.14</td><td>0.93</td><td>1</td></tr><tr><td>58771.0739</td><td>-1.1</td><td>2.3</td><td>5.932</td><td>0.010</td><td>-1.7</td><td>4.9</td><td>5.75</td><td>1.49</td><td>0.96</td><td>1</td></tr><tr><td>58772.1394</td><td>2.0</td><td>2.0</td><td>5.972</td><td>0.009</td><td>-7.0</td><td>4.1</td><td>4.25</td><td>1.26</td><td>0.99</td><td>1</td></tr><tr><td>58773.0344</td><td>0.8</td><td>2.1</td><td>5.959</td><td>0.009</td><td>-6.0</td><td>4.4</td><td>5.78</td><td>1.37</td><td>1.01</td><td>1</td></tr><tr><td>58788.1049</td><td>-2.3</td><td>2.0</td><td>5.922</td><td>0.009</td><td>-9.5</td><td>4.2</td><td>-2.68</td><td>1.19</td><td>1.45</td><td>1</td></tr><tr><td>58789.1097</td><td>-6.3</td><td>1.9</td><td>5.938</td><td>0.008</td><td>-19.7</td><td>3.9</td><td>1.73</td><td>1.05</td><td>1.48</td><td>1</td></tr><tr><td>58791.0882</td><td>0.8</td><td>2.1</td><td>5.898</td><td>0.009</td><td>-24.1</td><td>4.6</td><td>4.03</td><td>1.18</td><td>1.54</td><td>1</td></tr><tr><td>58792.1294</td><td>-0.4</td><td>2.1</td><td>5.930</td><td>0.009</td><td>-12.7</td><td>4.5</td><td>4.26</td><td>1.32</td><td>1.57</td><td>1</td></tr><tr><td>58793.1331</td><td>2.0</td><td>2.1</td><td>5.957</td><td>0.009</td><td>-4.0</td><td>4.3</td><td>5.60</td><td>1.12</td><td>1.60</td><td>1</td></tr><tr><td>58794.1219</td><td>-1.7</td><td>2.3</td><td>5.976</td><td>0.010</td><td>-32.0</td><td>5.1</td><td>5.15</td><td>1.41</td><td>1.63</td><td>1</td></tr><tr><td>58795.1056</td><td>6.6</td><td>1.9</td><td>5.980</td><td>0.008</td><td>-10.7</td><td>4.0</td><td>4.69</td><td>1.07</td><td>1.66</td><td>1</td></tr><tr><td>58796.1638*</td><td>10.7</td><td>1.8</td><td>5.994</td><td>0.011</td><td>-9.7</td><td>3.8</td><td>5.66</td><td>2.13</td><td>1.69</td><td>1</td></tr><tr><td>58796.9823</td><td>5.2</td><td>2.0</td><td>6.018</td><td>0.009</td><td>-14.1</td><td>4.5</td><td>5.73</td><td>1.24</td><td>1.71</td><td>1</td></tr><tr><td>58801.0872</td><td>9.7</td><td>2.0</td><td>5.972</td><td>0.009</td><td>-10.9</td><td>4.4</td><td>2.77</td><td>1.18</td><td>1.83</td><td>1</td></tr><tr><td>58802.0572</td><td>14.6</td><td>2.0</td><td>5.965</td><td>0.009</td><td>-13.1</td><td>4.5</td><td>2.16</td><td>1.22</td><td>1.86</td><td>1</td></tr><tr><td>58823.0087</td><td>-3.7</td><td>2.2</td><td>5.943</td><td>0.009</td><td>-3.2</td><td>4.7</td><td>1.32</td><td>1.42</td><td>2.47</td><td>1</td></tr><tr><td>58825.0963</td><td>-1.4</td><td>2.0</td><td>5.975</td><td>0.009</td><td>-0.8</td><td>4.2</td><td>4.52</td><td>1.15</td><td>2.53</td><td>1</td></tr><tr><td>58826.0288</td><td>-1.9</td><td>2.0</td><td>5.976</td><td>0.009</td><td>-4.5</td><td>4.1</td><td>5.20</td><td>1.15</td><td>2.56</td><td>1</td></tr><tr><td>58827.0595</td><td>-5.3</td><td>2.2</td><td>5.930</td><td>0.010</td><td>-18.0</td><td>4.6</td><td>5.81</td><td>1.11</td><td>2.58</td><td>1</td></tr><tr><td>58827.9807</td><td>-2.2</td><td>2.1</td><td>5.943</td><td>0.009</td><td>-18.2</td><td>4.3</td><td>2.76</td><td>1.15</td><td>2.61</td><td>1</td></tr><tr><td>58829.0723</td><td>-4.4</td><td>2.1</td><td>5.959</td><td>0.009</td><td>-2.7</td><td>4.4</td><td>3.00</td><td>1.18</td><td>2.64</td><td>1</td></tr><tr><td>58829.9988</td><td>4.1</td><td>2.0</td><td>5.947</td><td>0.008</td><td>-3.6</td><td>3.9</td><td>4.29</td><td>1.02</td><td>2.67</td><td>1</td></tr><tr><td>59088.1459</td><td>-0.0</td><td>1.6</td><td>5.986</td><td>0.007</td><td>3.8</td><td>2.9</td><td>1.11</td><td>0.91</td><td>10.17</td><td>2</td></tr><tr><td>59089.0973</td><td>0.3</td><td>1.6</td><td>5.959</td><td>0.007</td><td>1.6</td><td>3.0</td><td>2.59</td><td>0.87</td><td>10.20</td><td>2</td></tr><tr><td>59090.0986</td><td>-1.5</td><td>1.7</td><td>5.949</td><td>0.008</td><td>-0.7</td><td>3.3</td><td>2.87</td><td>0.86</td><td>10.23</td><td>2</td></tr><tr><td>59091.1334</td><td>0.6</td><td>1.9</td><td>5.925</td><td>0.008</td><td>20.8</td><td>3.7</td><td>1.87</td><td>0.89</td><td>10.26</td><td>2</td></tr><tr><td>59092.1231</td><td>-0.9</td><td>1.7</td><td>5.952</td><td>0.007</td><td>4.4</td><td>3.1</td><td>3.62</td><td>0.86</td><td>10.29</td><td>2</td></tr><tr><td>59093.1251</td><td>-1.0</td><td>2.0</td><td>5.897</td><td>0.009</td><td>8.2</td><td>4.0</td><td>2.71</td><td>0.82</td><td>10.32</td><td>2</td></tr><tr><td>59094.1339</td><td>-2.5</td><td>1.9</td><td>5.900</td><td>0.008</td><td>14.4</td><td>3.9</td><td>3.80</td><td>1.02</td><td>10.35</td><td>2</td></tr><tr><td>59095.1265</td><td>2.8</td><td>1.7</td><td>5.945</td><td>0.008</td><td>3.2</td><td>3.3</td><td>2.57</td><td>0.89</td><td>10.38</td><td>2</td></tr><tr><td>59096.1231</td><td>0.2</td><td>1.7</td><td>5.988</td><td>0.007</td><td>4.3</td><td>3.2</td><td>1.37</td><td>1.02</td><td>10.41</td><td>2</td></tr><tr><td>59097.1475</td><td>-2.6</td><td>2.1</td><td>5.903</td><td>0.009</td><td>-0.0</td><td>4.1</td><td>1.94</td><td>1.32</td><td>10.44</td><td>2</td></tr><tr><td>59098.1493</td><td>-6.5</td><td>1.6</td><td>5.968</td><td>0.007</td><td>-7.4</td><td>3.0</td><td>0.74</td><td>0.84</td><td>10.47</td><td>2</td></tr><tr><td>59099.1228</td><td>-6.3</td><td>1.9</td><td>5.914</td><td>0.008</td><td>9.6</td><td>3.8</td><td>1.47</td><td>1.02</td><td>10.49</td><td>2</td></tr><tr><td>59101.1095</td><td>-3.0</td><td>1.9</td><td>5.923</td><td>0.008</td><td>0.0</td><td>3.9</td><td>-0.89</td><td>0.98</td><td>10.55</td><td>2</td></tr><tr><td>59102.0900</td><td>-0.0</td><td>1.7</td><td>5.954</td><td>0.007</td><td>-7.3</td><td>3.3</td><td>1.75</td><td>0.99</td><td>10.58</td><td>2</td></tr><tr><td>59103.0880</td><td>-2.5</td><td>1.9</td><td>5.920</td><td>0.008</td><td>2.0</td><td>3.7</td><td>0.68</td><td>0.92</td><td>10.61</td><td>2</td></tr><tr><td>59111.1241</td><td>-1.0</td><td>1.9</td><td>5.960</td><td>0.008</td><td>-4.6</td><td>4.0</td><td>0.85</td><td>1.14</td><td>10.84</td><td>2</td></tr><tr><td>59112.1125</td><td>-0.1</td><td>1.8</td><td>5.974</td><td>0.008</td><td>6.5</td><td>3.3</td><td>2.15</td><td>1.04</td><td>10.87</td><td>2</td></tr><tr><td>59114.0475</td><td>2.5</td><td>2.1</td><td>5.933</td><td>0.009</td><td>2.9</td><td>4.2</td><td>0.27</td><td>0.92</td><td>10.93</td><td>2</td></tr><tr><td>59115.1229</td><td>0.9</td><td>1.8</td><td>5.983</td><td>0.008</td><td>-11.9</td><td>3.2</td><td>-0.84</td><td>0.99</td><td>10.96</td><td>2</td></tr><tr><td>59116.0896</td><td>4.1</td><td>1.8</td><td>5.957</td><td>0.008</td><td>0.1</td><td>3.5</td><td>-4.41</td><td>1.11</td><td>10.99</td><td>2</td></tr><tr><td>59118.0749*</td><td>3.7</td><td>1.8</td><td>5.914</td><td>0.008</td><td>4.3</td><td>3.6</td><td>-4.26</td><td>0.92</td><td>11.04</td><td>2</td></tr><tr><td>59120.1385</td><td>5.5</td><td>1.9</td><td>5.964</td><td>0.012</td><td>-3.6</td><td>4.2</td><td>-2.66</td><td>2.06</td><td>11.10</td><td>2</td></tr><tr><td>59121.0153</td><td>-4.6</td><td>1.8</td><td>5.931</td><td>0.008</td><td>6.0</td><td>3.5</td><td>-2.10</td><td>0.90</td><td>11.13</td><td>2</td></tr><tr><td>59122.1230</td><td>-0.4</td><td>1.8</td><td>5.929</td><td>0.008</td><td>6.3</td><td>3.4</td><td>1.95</td><td>1.17</td><td>11.16</td><td>2</td></tr><tr><td>59123.0174</td><td>-5.4</td><td>2.0</td><td>5.903</td><td>0.009</td><td>1.1</td><td>4.0</td><td>3.13</td><td>1.01</td><td>11.19</td><td>2</td></tr><tr><td>59127.0178</td><td>-4.9</td><td>1.8</td><td>5.921</td><td>0.008</td><td>3.6</td><td>3.4</td><td>2.37</td><td>0.90</td><td>11.30</td><td>2</td></tr></tbody></table>
<table><caption>Table A.2. SPIRou RVs and activity indicators continued.</caption><thead><tr><th>BJD-2400000<br>d</th><th>ΔRV<br>m s<sup>-1</sup></th><th>σ<sub>RV</sub><br>m s<sup>-1</sup></th><th>FWHM<sub>LBL</sub><br>km s<sup>-1</sup></th><th>σ<sub>FWHM</sub><br>km s<sup>-1</sup></th><th>chrom. velocity slope<br>m s<sup>-1</sup></th><th>σ<sub>chrom.</sub><br>m s<sup>-1</sup></th><th>B<sub>ℓ</sub><br>G</th><th>σ<sub>B</sub><br>G</th><th>Phase</th><th>Season</th></tr></thead><tbody><tr><td>59128.0016</td><td>-4.3</td><td>1.8</td><td>5.918</td><td>0.008</td><td>1.7</td><td>3.3</td><td>2.55</td><td>0.84</td><td>11.33</td><td>2</td></tr><tr><td>59129.0533</td><td>-3.7</td><td>1.8</td><td>5.945</td><td>0.008</td><td>7.0</td><td>3.4</td><td>2.72</td><td>1.05</td><td>11.36</td><td>2</td></tr><tr><td>59129.9930</td><td>-1.9</td><td>1.9</td><td>5.897</td><td>0.008</td><td>7.1</td><td>3.7</td><td>2.36</td><td>0.86</td><td>11.39</td><td>2</td></tr><tr><td>59131.1038</td><td>-1.1</td><td>1.8</td><td>5.917</td><td>0.008</td><td>5.7</td><td>3.4</td><td>1.31</td><td>0.92</td><td>11.42</td><td>2</td></tr><tr><td>59154.1087</td><td>-9.4</td><td>1.7</td><td>5.933</td><td>0.007</td><td>0.4</td><td>3.0</td><td>1.26</td><td>0.75</td><td>12.09</td><td></td></tr><tr><td>59154.9730</td><td>-8.1</td><td>1.7</td><td>5.925</td><td>0.007</td><td>-9.5</td><td>3.0</td><td>3.66</td><td>0.80</td><td>12.12</td><td></td></tr><tr><td>59156.1429</td><td>-6.5</td><td>1.8</td><td>6.000</td><td>0.013</td><td>-10.2</td><td>3.6</td><td>4.49</td><td>1.15</td><td>12.15</td><td></td></tr><tr><td>59157.0273</td><td>-8.2</td><td>1.8</td><td>5.924</td><td>0.008</td><td>-1.0</td><td>3.6</td><td>3.42</td><td>1.24</td><td>12.18</td><td></td></tr><tr><td>59158.0274</td><td>-9.3</td><td>1.8</td><td>5.921</td><td>0.008</td><td>-2.9</td><td>3.5</td><td>6.36</td><td>1.18</td><td>12.21</td><td></td></tr><tr><td>59159.0282</td><td>-3.0</td><td>1.7</td><td>5.937</td><td>0.007</td><td>4.5</td><td>3.2</td><td>7.64</td><td>0.84</td><td>12.24</td><td></td></tr><tr><td>59207.9651</td><td>6.7</td><td>2.0</td><td>5.900</td><td>0.009</td><td>1.5</td><td>4.0</td><td>1.63</td><td>1.15</td><td>13.66</td><td></td></tr><tr><td>59208.9774</td><td>4.2</td><td>1.8</td><td>5.942</td><td>0.008</td><td>0.6</td><td>3.4</td><td>0.71</td><td>0.95</td><td>13.69</td><td></td></tr><tr><td>59209.9278</td><td>1.3</td><td>1.6</td><td>5.961</td><td>0.007</td><td>7.3</td><td>3.0</td><td>2.87</td><td>0.80</td><td>13.71</td><td></td></tr><tr><td>59211.8520</td><td>1.4</td><td>1.8</td><td>5.962</td><td>0.007</td><td>14.1</td><td>3.1</td><td>4.12</td><td>0.92</td><td>13.77</td><td></td></tr><tr><td>59212.9753</td><td>4.1</td><td>1.7</td><td>5.961</td><td>0.007</td><td>2.7</td><td>3.2</td><td>1.49</td><td>1.01</td><td>13.80</td><td></td></tr><tr><td>59214.0266</td><td>2.1</td><td>1.7</td><td>5.979</td><td>0.007</td><td>5.8</td><td>3.2</td><td>2.87</td><td>0.93</td><td>13.83</td><td></td></tr><tr><td>59214.9368</td><td>-1.7</td><td>1.7</td><td>5.953</td><td>0.007</td><td>-0.1</td><td>3.1</td><td>3.91</td><td>0.91</td><td>13.86</td><td></td></tr><tr><td>59215.8867</td><td>1.1</td><td>2.0</td><td>5.985</td><td>0.009</td><td>0.4</td><td>4.2</td><td>-0.06</td><td>1.39</td><td>13.89</td><td></td></tr><tr><td>59216.8987</td><td>5.4</td><td>1.6</td><td>5.938</td><td>0.007</td><td>-0.4</td><td>3.0</td><td>1.14</td><td>0.77</td><td>13.92</td><td></td></tr><tr><td>59217.9291</td><td>4.9</td><td>1.6</td><td>5.963</td><td>0.007</td><td>10.5</td><td>2.9</td><td>-0.82</td><td>0.83</td><td>13.95</td><td></td></tr><tr><td>59218.9165</td><td>1.9</td><td>1.7</td><td>5.945</td><td>0.007</td><td>5.4</td><td>3.2</td><td>-0.64</td><td>0.84</td><td>13.98</td><td></td></tr><tr><td>59219.9436</td><td>5.0</td><td>1.7</td><td>5.943</td><td>0.007</td><td>11.0</td><td>3.1</td><td>-0.10</td><td>0.84</td><td>14.01</td><td></td></tr><tr><td>59220.9034</td><td>0.8</td><td>1.6</td><td>5.947</td><td>0.007</td><td>6.6</td><td>3.0</td><td>2.01</td><td>0.85</td><td>14.03</td><td></td></tr><tr><td>59221.9382</td><td>1.8</td><td>1.6</td><td>5.945</td><td>0.007</td><td>22.3</td><td>2.9</td><td>0.75</td><td>0.83</td><td>14.06</td><td></td></tr><tr><td>59222.9549</td><td>-1.7</td><td>1.6</td><td>5.953</td><td>0.007</td><td>20.6</td><td>2.9</td><td>3.87</td><td>0.84</td><td>14.09</td><td></td></tr><tr><td>59267.7684</td><td>-3.6</td><td>1.7</td><td>5.942</td><td>0.007</td><td>7.3</td><td>3.1</td><td>0.02</td><td>0.97</td><td>15.40</td><td></td></tr><tr><td>59271.7394</td><td>-5.7</td><td>1.6</td><td>5.957</td><td>0.007</td><td>9.2</td><td>2.8</td><td>-0.86</td><td>0.84</td><td>15.51</td><td></td></tr><tr><td>59273.8932</td><td>-2.6</td><td>1.8</td><td>5.956</td><td>0.008</td><td>-0.9</td><td>3.4</td><td>-0.40</td><td>1.06</td><td>15.57</td><td></td></tr><tr><td>59275.8955</td><td>-5.6</td><td>1.9</td><td>5.994</td><td>0.008</td><td>-5.5</td><td>3.8</td><td>0.19</td><td>1.25</td><td>15.63</td><td></td></tr><tr><td>59276.8024</td><td>-5.4</td><td>1.5</td><td>5.964</td><td>0.007</td><td>16.5</td><td>2.9</td><td>-2.27</td><td>0.77</td><td>15.66</td><td></td></tr><tr><td>59277.7671</td><td>-5.3</td><td>1.6</td><td>5.951</td><td>0.007</td><td>-1.4</td><td>3.0</td><td>-1.68</td><td>0.74</td><td>15.69</td><td></td></tr><tr><td>59297.7410</td><td>-6.3</td><td>1.7</td><td>5.896</td><td>0.008</td><td>4.7</td><td>3.4</td><td>2.83</td><td>0.81</td><td>16.27</td><td></td></tr><tr><td>59440.1423</td><td>0.5</td><td>1.8</td><td>5.930</td><td>0.008</td><td>0.6</td><td>3.4</td><td>0.33</td><td>1.07</td><td>20.41</td><td>3</td></tr><tr><td>59442.1181</td><td>-0.4</td><td>1.8</td><td>5.934</td><td>0.008</td><td>3.8</td><td>3.5</td><td>1.99</td><td>0.94</td><td>20.46</td><td>3</td></tr><tr><td>59443.1094</td><td>-1.9</td><td>1.6</td><td>5.975</td><td>0.007</td><td>4.2</td><td>3.0</td><td>2.30</td><td>0.87</td><td>20.49</td><td>3</td></tr><tr><td>59444.1392</td><td>-0.4</td><td>1.7</td><td>5.949</td><td>0.007</td><td>-4.1</td><td>3.2</td><td>3.95</td><td>0.93</td><td>20.52</td><td>3</td></tr><tr><td>59445.1268</td><td>2.1</td><td>1.6</td><td>5.959</td><td>0.007</td><td>-6.8</td><td>2.9</td><td>6.81</td><td>0.87</td><td>20.55</td><td>3</td></tr><tr><td>59447.1009</td><td>1.5</td><td>1.7</td><td>5.972</td><td>0.007</td><td>-6.5</td><td>3.1</td><td>7.99</td><td>0.87</td><td>20.61</td><td>3</td></tr><tr><td>59448.1311</td><td>6.5</td><td>1.8</td><td>5.942</td><td>0.008</td><td>-9.3</td><td>3.4</td><td>5.45</td><td>1.00</td><td>20.64</td><td>3</td></tr><tr><td>59449.1436</td><td>0.6</td><td>1.8</td><td>5.954</td><td>0.008</td><td>-11.1</td><td>3.4</td><td>8.44</td><td>0.90</td><td>20.67</td><td>3</td></tr><tr><td>59450.1248</td><td>3.1</td><td>1.7</td><td>5.973</td><td>0.007</td><td>-2.1</td><td>3.1</td><td>4.68</td><td>0.98</td><td>20.70</td><td>3</td></tr><tr><td>59451.1236</td><td>2.6</td><td>1.9</td><td>5.951</td><td>0.008</td><td>-8.3</td><td>3.4</td><td>4.54</td><td>0.95</td><td>20.73</td><td>3</td></tr><tr><td>59452.1097</td><td>1.7</td><td>1.7</td><td>5.971</td><td>0.007</td><td>0.7</td><td>3.0</td><td>3.13</td><td>0.86</td><td>20.76</td><td>3</td></tr><tr><td>59453.1165</td><td>9.1</td><td>1.7</td><td>5.981</td><td>0.007</td><td>1.5</td><td>3.1</td><td>6.72</td><td>0.92</td><td>20.78</td><td>3</td></tr><tr><td>59473.0971</td><td>3.0</td><td>1.6</td><td>5.964</td><td>0.007</td><td>0.7</td><td>3.0</td><td>-2.72</td><td>0.89</td><td>21.37</td><td>3</td></tr><tr><td>59475.0965</td><td>6.6</td><td>1.8</td><td>5.932</td><td>0.008</td><td>-0.5</td><td>3.4</td><td>-2.60</td><td>1.09</td><td>21.42</td><td>3</td></tr><tr><td>59477.1203</td><td>6.5</td><td>1.9</td><td>5.922</td><td>0.008</td><td>2.4</td><td>3.7</td><td>1.47</td><td>0.86</td><td>21.48</td><td>3</td></tr><tr><td>59478.0777</td><td>-0.2</td><td>1.9</td><td>5.919</td><td>0.008</td><td>7.0</td><td>3.8</td><td>3.05</td><td>0.93</td><td>21.51</td><td>3</td></tr><tr><td>59480.0106</td><td>7.1</td><td>1.7</td><td>5.955</td><td>0.007</td><td>0.0</td><td>3.1</td><td>4.14</td><td>0.85</td><td>21.57</td><td>3</td></tr><tr><td>59481.0283</td><td>8.0</td><td>1.7</td><td>5.954</td><td>0.008</td><td>-4.3</td><td>3.2</td><td>4.29</td><td>0.97</td><td>21.60</td><td>3</td></tr><tr><td>59482.0279</td><td>9.6</td><td>1.8</td><td>5.945</td><td>0.008</td><td>-3.3</td><td>3.6</td><td>6.76</td><td>1.34</td><td>21.62</td><td>3</td></tr><tr><td>59502.1252</td><td>-1.7</td><td>1.9</td><td>5.892</td><td>0.008</td><td>-9.2</td><td>3.5</td><td>-3.36</td><td>0.99</td><td>22.21</td><td></td></tr><tr><td>59504.1295</td><td>-3.5</td><td>1.9</td><td>5.971</td><td>0.008</td><td>2.1</td><td>3.9</td><td>-1.66</td><td>1.19</td><td>22.27</td><td></td></tr><tr><td>59506.1212</td><td>-16.5</td><td>2.5</td><td>5.904</td><td>0.011</td><td>8.9</td><td>5.4</td><td>-2.18</td><td>1.46</td><td>22.33</td><td></td></tr><tr><td>59507.1298</td><td>-2.2</td><td>1.8</td><td>5.931</td><td>0.008</td><td>-4.6</td><td>3.5</td><td>-2.24</td><td>0.91</td><td>22.35</td><td></td></tr><tr><td>59509.1337</td><td>-13.4</td><td>2.0</td><td>5.962</td><td>0.009</td><td>3.2</td><td>4.3</td><td>-0.33</td><td>1.19</td><td>22.41</td><td></td></tr><tr><td>59510.1171</td><td>-5.5</td><td>2.1</td><td>5.950</td><td>0.009</td><td>15.0</td><td>4.4</td><td>-1.46</td><td>1.22</td><td>22.44</td><td></td></tr><tr><td>59511.1047</td><td>-5.4</td><td>1.7</td><td>5.918</td><td>0.007</td><td>0.7</td><td>3.2</td><td>0.15</td><td>0.93</td><td>22.47</td><td></td></tr><tr><td>59513.0782</td><td>-6.6</td><td>1.5</td><td>5.961</td><td>0.007</td><td>-5.2</td><td>2.7</td><td>0.50</td><td>0.78</td><td>22.53</td><td></td></tr><tr><td>59514.0972</td><td>-1.4</td><td>1.6</td><td>5.946</td><td>0.007</td><td>1.1</td><td>2.9</td><td>0.82</td><td>0.88</td><td>22.56</td><td></td></tr></tbody></table>
<table><caption>**Table A.2.** SPIRou RVs and activity indicators continued.</caption><thead><tr><th>BJD-2400000<br>d</th><th>ΔRV<br>m s<sup>-1</sup></th><th>σ<sub>RV</sub><br>m s<sup>-1</sup></th><th>FWHM<sub>LBL</sub><br>km s<sup>-1</sup></th><th>σ<sub>FWHM</sub><br>km s<sup>-1</sup></th><th>chrom. velocity slope<br>m s<sup>-1</sup></th><th>σ<sub>chrom.</sub><br>m s<sup>-1</sup></th><th>B<sub>ℓ</sub><br>G</th><th>σ<sub>B</sub><br>G</th><th>Phase</th><th>Season</th></tr></thead><tbody><tr><td>59515.1549</td><td>-2.8</td><td>1.8</td><td>5.924</td><td>0.008</td><td>-3.0</td><td>3.2</td><td>0.33</td><td>0.82</td><td>22.59</td><td></td></tr><tr><td>59516.0071</td><td>2.5</td><td>1.7</td><td>5.970</td><td>0.007</td><td>0.6</td><td>2.8</td><td>0.83</td><td>0.78</td><td>22.61</td><td></td></tr><tr><td>59535.0426</td><td>-7.0</td><td>2.0</td><td>5.926</td><td>0.009</td><td>0.4</td><td>4.1</td><td>-6.65</td><td>1.31</td><td>23.17</td><td>4</td></tr><tr><td>59537.9797</td><td>-1.3</td><td>1.7</td><td>5.906</td><td>0.007</td><td>-15.9</td><td>3.2</td><td>-3.38</td><td>1.05</td><td>23.25</td><td>4</td></tr><tr><td>59539.0540</td><td>-4.0</td><td>1.7</td><td>5.925</td><td>0.007</td><td>-3.6</td><td>3.3</td><td>-1.83</td><td>0.93</td><td>23.28</td><td>4</td></tr><tr><td>59540.0405</td><td>-0.1</td><td>1.7</td><td>5.925</td><td>0.007</td><td>-9.8</td><td>3.1</td><td>-1.41</td><td>0.91</td><td>23.31</td><td>4</td></tr><tr><td>59541.0569</td><td>0.1</td><td>1.6</td><td>5.923</td><td>0.007</td><td>-3.0</td><td>3.1</td><td>-0.24</td><td>0.76</td><td>23.34</td><td>4</td></tr><tr><td>59558.0071</td><td>7.6</td><td>1.7</td><td>5.985</td><td>0.007</td><td>3.7</td><td>2.9</td><td>0.96</td><td>0.80</td><td>23.83</td><td>4</td></tr><tr><td>59560.0476</td><td>5.2</td><td>1.6</td><td>5.996</td><td>0.007</td><td>14.7</td><td>3.0</td><td>1.01</td><td>0.93</td><td>23.89</td><td>4</td></tr><tr><td>59561.0577</td><td>9.7</td><td>1.6</td><td>5.977</td><td>0.007</td><td>14.2</td><td>3.1</td><td>0.82</td><td>0.94</td><td>23.92</td><td>4</td></tr><tr><td>59562.0182</td><td>3.4</td><td>1.7</td><td>5.974</td><td>0.008</td><td>14.3</td><td>3.3</td><td>0.48</td><td>1.06</td><td>23.95</td><td>4</td></tr><tr><td>59563.8546</td><td>2.7</td><td>1.6</td><td>5.952</td><td>0.007</td><td>2.3</td><td>2.8</td><td>0.55</td><td>0.76</td><td>24.00</td><td>4</td></tr><tr><td>59564.9371</td><td>3.8</td><td>1.7</td><td>5.950</td><td>0.007</td><td>7.4</td><td>3.1</td><td>0.19</td><td>0.89</td><td>24.03</td><td>4</td></tr><tr><td>59566.9795</td><td>-0.3</td><td>1.8</td><td>5.896</td><td>0.008</td><td>7.8</td><td>3.4</td><td>-1.53</td><td>0.92</td><td>24.09</td><td>4</td></tr><tr><td>59585.9682</td><td>1.3</td><td>1.9</td><td>5.916</td><td>0.008</td><td>6.9</td><td>3.2</td><td>-3.13</td><td>0.77</td><td>24.65</td><td>4</td></tr><tr><td>59587.9584</td><td>2.1</td><td>1.8</td><td>5.919</td><td>0.008</td><td>4.8</td><td>3.5</td><td>-0.20</td><td>0.89</td><td>24.70</td><td>4</td></tr><tr><td>59588.9581</td><td>2.1</td><td>1.8</td><td>5.925</td><td>0.008</td><td>1.2</td><td>3.4</td><td>-1.16</td><td>1.07</td><td>24.73</td><td>4</td></tr><tr><td>59589.8836</td><td>-0.7</td><td>1.9</td><td>5.912</td><td>0.008</td><td>5.8</td><td>3.8</td><td>-3.13</td><td>0.99</td><td>24.76</td><td>4</td></tr><tr><td>59590.7756</td><td>2.5</td><td>1.6</td><td>5.954</td><td>0.007</td><td>13.7</td><td>2.8</td><td>-0.73</td><td>0.81</td><td>24.79</td><td>4</td></tr><tr><td>59591.8296</td><td>2.5</td><td>1.6</td><td>5.960</td><td>0.007</td><td>12.8</td><td>2.9</td><td>1.69</td><td>0.86</td><td>24.82</td><td>4</td></tr><tr><td>59592.7883</td><td>3.3</td><td>1.7</td><td>5.930</td><td>0.007</td><td>14.1</td><td>3.2</td><td>1.57</td><td>0.86</td><td>24.84</td><td>4</td></tr><tr><td>59593.9250</td><td>7.8</td><td>1.5</td><td>5.961</td><td>0.007</td><td>10.1</td><td>2.7</td><td>2.23</td><td>0.74</td><td>24.88</td><td>4</td></tr><tr><td>59594.9754</td><td>5.8</td><td>1.6</td><td>5.952</td><td>0.007</td><td>21.2</td><td>2.8</td><td>1.81</td><td>0.80</td><td>24.91</td><td>4</td></tr><tr><td>59596.9323</td><td>-1.5</td><td>1.6</td><td>5.939</td><td>0.007</td><td>11.8</td><td>2.9</td><td>2.60</td><td>0.81</td><td>24.96</td><td>4</td></tr><tr><td>59597.9986</td><td>-0.9</td><td>1.6</td><td>5.923</td><td>0.007</td><td>16.1</td><td>3.1</td><td>2.66</td><td>0.93</td><td>25.00</td><td>4</td></tr><tr><td>59598.9765</td><td>0.2</td><td>1.6</td><td>5.929</td><td>0.007</td><td>6.8</td><td>3.0</td><td>1.15</td><td>0.84</td><td>25.02</td><td>4</td></tr><tr><td>59599.9542</td><td>-3.6</td><td>1.7</td><td>5.913</td><td>0.007</td><td>-8.2</td><td>3.2</td><td>0.40</td><td>0.82</td><td>25.05</td><td>4</td></tr><tr><td>59601.9558</td><td>0.4</td><td>1.7</td><td>5.973</td><td>0.007</td><td>4.3</td><td>3.3</td><td>-1.90</td><td>1.13</td><td>25.11</td><td>4</td></tr><tr><td>59603.9449</td><td>-3.7</td><td>1.7</td><td>5.909</td><td>0.007</td><td>1.1</td><td>3.2</td><td>-3.85</td><td>0.84</td><td>25.17</td><td>4</td></tr><tr><td>59604.9396</td><td>-1.2</td><td>1.6</td><td>5.953</td><td>0.007</td><td>-3.3</td><td>2.9</td><td>-1.71</td><td>0.87</td><td>25.20</td><td>4</td></tr><tr><td>59606.8864</td><td>-2.4</td><td>1.7</td><td>5.938</td><td>0.007</td><td>0.5</td><td>3.0</td><td>-1.20</td><td>0.92</td><td>25.25</td><td>4</td></tr><tr><td>59607.9038</td><td>-4.6</td><td>1.7</td><td>5.978</td><td>0.007</td><td>-1.7</td><td>2.9</td><td>1.27</td><td>1.60</td><td>25.28</td><td>4</td></tr><tr><td>59609.8721</td><td>-4.1</td><td>1.5</td><td>5.969</td><td>0.008</td><td>-5.7</td><td>2.8</td><td>6.21</td><td>1.46</td><td>25.34</td><td>4</td></tr></tbody></table>

**Notes.** *Night with more than one spectropolarimetric sequence. The RV and activity indicators are computed with the nightly average spectra.
**Appendix B: Longitudinal magnetic field *B*<sub>ℓ</sub> GP posterior distributions**

The posterior distributions of the GP model applied in the longitudinal magnetic field are shown in Fig. B.1.

<figure><img src="image_15.png" /><figcaption>**Fig. B.1.** Corner plot of the GP hyperparameters on the *B*<sub>*ℓ*</sub> time series sampled from our MCMC framework generated with the corner Python package (Foreman-Mackey 2016). The blue line over the posterior distributions indicates the mean and the dashed lines show the 0.16 and 0.84 quantile uncertainties.</figcaption></figure>
**Appendix C: ZDI inversion of the Stokes V LSD profiles**

The best ZDI fit to the time series of Stokes V LSD profiles observed in September – December 2019, August – October 2020, August – September 2021, and November 2021 – January 2022 are shown in Figure C.1 and Figure C.2.

<figure><img src="image_16.png" /><figcaption>**Fig. C.1.** Stokes V LSD profiles (thin black lines) and best-fitting ZDI profiles (thick red lines). Profiles observed between September and December 2019 (Season S₁ in Table 4) are shown in the left-hand panel whereas the profiles observed between August and October 2020 (Season S₂ in Table 4) are shown in the right-hand panel. The stellar rotation cycle computed using the first observation (BJD = 2 458 738.128) as a reference time and assuming a rotation period of 34.4 d (see Table 3) is given on the right-hand side of each observation. The ±1σ error bars are indicated on the left-hand side of each profile.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Fig. C.2.** Same as Figure C.1, but for the observations collected in Season 3 (August to September 2021, left-hand panel), and Season 4 (November 2021 and January 2022, right-hand panel). We note that we use the same rotation period and reference time as in Figure C.1 to compute the stellar rotation phases.</figcaption></figure>
**Appendix D: Periodogram of the activity indicators**

The GLS periodograms of all the SOPHIE and SPIRou activity indicators, for each of the three subsets of observations $S_1^{RV}$, $S_2^{RV}$, and $S_3^{RV}$, are shown in Figures D.1 and D.2.

<figure><img src="image_18.png" /><figcaption>**Fig. D.1.** Generalized Lomb-Scargle periodograms of the SOPHIE RVs, activity indicators, and window function of the three seasons of RV observations: $S_1^{RV}$, $S_2^{RV}$, and $S_3^{RV}$. The vertical color lines mark the $P_{\text{rot}}$ at 34.4 d in yellow, the first harmonic at $P_{\text{rot}}/2$ in light blue, and the time range of differential rotation in light green between $P_{\text{eq}} = 32.0$ d and $P_{\text{pol}} = 45.5$ days. The black horizontal lines indicate the FAP levels of 10% (solid), 1% (dashed), and 0.1% (dotted).</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>Fig. D.2. Same as Figure D.1 for SPIRou data.</figcaption></figure>
**Appendix E: GP model of the activity indicators time series**

The best-fit hyperparameters of the GP model applied in the SOPHIE and SPIRou activity indicators are listed in Table E.1. Figures E.1 and E.2 show the best GP model for each activity indicator, and Figures E.3, E.4, and E.5 show the corner plots of the posterior distributions.

<figure><img src="image_20.png" /><figcaption>**Fig. E.1.** Time series of the SOPHIE RVs, FWHM, CCF contrast, BIS, S index, and $H\alpha$. The GP model of each time series is shown in blue with the $3\sigma$ uncertainties in light blue.</figcaption></figure>
<table><caption>Table E.1. Best-fit hyperparameters of the GP model using a quasi-periodic kernel in the SOPHIE and SPIRou activity indicators.</caption><thead><tr><th rowspan="3">Parameter</th><th colspan="6">SOPHIE Posterkerns</th><th colspan="6">SPIRou Posterkerns</th></tr><tr><th>FWHM</th><th>contrast</th><th>BIS</th><th>S index</th><th>Hα</th><th>FWHMLBL</th><th>chromatic<br>velocity slope</th><th>pEW K</th><th>pEW Ti</th><th>pEW Al</th><th>pEW Fe</th><th>pEW He</th></tr><tr><th>km s<sup>-1</sup></th><th>% CCF</th><th>ms<sup>-1</sup></th><th></th><th></th><th></th><th>km s<sup>-1</sup></th><th>ms<sup>-1</sup></th><th>Å</th><th>Å</th><th>Å</th><th>Å</th><th>Å</th></tr></thead><tbody><tr><td>mean value μ</td><td>4.82 ± 0.02</td><td>21.6 ± 0.3</td><td>11.7 ± 0.9</td><td>1.73 ± 0.6</td><td>0.265 ± 0.005</td><td>5.945 ± 0.005</td><td>-1.1 ± 3.6</td><td>0.403 ± 0.005</td><td>0.289 ± 0.002</td><td>1.48 ± 0.02</td><td>0.71 ± 0.01</td><td>0.83 ± 0.01</td></tr><tr><td>jitter σ</td><td>0.04 ± 0.05</td><td>0.37 ± 0.05</td><td>1.8 ± 0.8</td><td>0.09 ± 0.01</td><td>0.007 ± 0.001</td><td>0.021 ± 0.002</td><td>4.9 ± 0.5</td><td>0.0111 ± 0.0008</td><td>0.008 ± 0.001</td><td>0.043 ± 0.004</td><td>0.011 ± 0.001</td><td>0.039 ± 0.003</td></tr><tr><td>amplitude α</td><td>0.04 ± 0.01</td><td>0.63 ± 0.16</td><td>2.3 ± 1.1</td><td>0.14 ± 0.4</td><td>0.014 ± 0.003</td><td>0.015 ± 0.005</td><td>8.7 ± 2.0</td><td>0.013 ± 0.004</td><td>0.005 ± 0.002</td><td>0.08 ± 0.02</td><td>0.015 ± 0.003</td><td>0.03 ± 0.01</td></tr><tr><td>decay time l [d]</td><td>82 ± 36</td><td>48<sup>+24</sup><sub>-10</sub></td><td>96 ± 38</td><td>67<sup>+34</sup><sub>-21</sub></td><td>88<sup>+38</sup><sub>-27</sub></td><td>71<sup>+55</sup><sub>-29</sub></td><td>71<sup>+40</sup><sub>-22</sub></td><td>62 ± 24</td><td>81 ± 35</td><td>51<sup>+29</sup><sub>-12</sub></td><td>41 ± 8</td><td>77 ± 35</td></tr><tr><td>smoothing factor β</td><td>1.5 ± 0.4</td><td>1.6 ± 0.3</td><td>0.7 ± 0.5</td><td>1.4 ± 0.4</td><td>0.7 ± 0.3</td><td>0.5 ± 0.4</td><td>1.3 ± 0.4</td><td>1.5 ± 0.4</td><td>1.3 ± 0.5</td><td>0.6 ± 0.2</td><td>1.7 ± 0.3</td><td>1.2 ± 0.6</td></tr><tr><td>rotation period P<sub>rot</sub> [d]</td><td>32.1 ± 3.9</td><td>37.4 ± 3.5</td><td>36.7 ± 3.5</td><td>33.8 ± 4.4</td><td>37 ± 2.3</td><td>36.7 ± 4.1</td><td>37.1 ± 4.5</td><td>42.2 ± 4.9</td><td>37.1 ± 6.3</td><td>37.9 ± 2.1</td><td>36.9 ± 7.0</td><td>27.0 ± 4.9</td></tr><tr><td>RMS of residuals</td><td>0.03</td><td>0.31</td><td>2.8</td><td>0.07</td><td>0.005</td><td>0.02</td><td>5.6</td><td>0.01</td><td>0.01</td><td>0.04</td><td>0.01</td><td>0.04</td></tr><tr><td>χ<sup>2</sup></td><td>12.1</td><td>10.7</td><td>1.3</td><td>1.8</td><td>23.4</td><td>7.1</td><td>2.7</td><td>63.6</td><td>23.1</td><td>282</td><td>28.6</td><td>158.5</td></tr></tbody></table>
<figure><img src="image_21.png" /><figcaption>**Fig. E.2.** Time series of the SPIRou longitudinal magnetic field $B_\ell$, RVs, FWHM<sub>LBL</sub>, chromatic velocity slope, and the pEW of K I, Ti I, Al I, Fe I, and He I. The GP model of each time series is shown in red with the $3\sigma$ uncertainties in light red.</figcaption></figure>
<figure><img src="image_22.png" /></figure>

SOPHIE RV

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

BIS

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

CCF contrast

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

Hα

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

S index

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

FWHM

**Fig. E.3.** Corner plot of the posterior distributions of the GP regression applied in the SOPHIE RVs, Hα, BIS, S index, CCF contrast, and FWHM, independently.
<figure><img src="image_28.png" /><figcaption>**Fig. E.4.** Corner plot of the posterior distributions of the GP regression applied in the SPIRou RVs, FWHM, chromatic velocity slope, Al I, Ti I, independently.</figcaption></figure>
<figure><img src="image_29.png" /><figcaption>Fig. E.5. Same as Figure E.4, but for Fe I and He I.</figcaption></figure>
# Appendix F: Multi-dimensional GP corner plots

The posterior distributions of the multidimensional GP model applied in the SOPHIE and SPIRou data are shown in Figures F.1 and F.2.

<figure><img src="image_30.png" /><figcaption>**Fig. F.1.** Corner plot with the posterior distributions of the multidimensional GP applied in the SOPHIE RVs, $H\alpha$, and BIS time series.</figcaption></figure>
<figure><img src="image_31.png" /><figcaption>Fig. F.2. Corner plot with the posterior distributions of the multidimensional GP applied in the SPIRou RVs and FWHM<sub>LBL</sub> time series.</figcaption></figure>