# Evolving magnetic lives of Sun-like stars - I.
## Characterisation of the large-scale magnetic field with Zeeman-Doppler imaging

S Bellotti, T Lüftinger, S.Boro Saikia, C.P Folsom, P Petit, J Morin, M

Güdel, J.F Donati, E Alecian

► To cite this version:

S Bellotti, T Lüftinger, S.Boro Saikia, C.P Folsom, P Petit, et al.. Evolving magnetic lives of Sun-like stars - I. Characterisation of the large-scale magnetic field with Zeeman-Doppler imaging. Astronomy & Astrophysics - A&A, 2025, 700, pp.A282. 10.1051/0004-6361/202555502. hal-05271734

HAL Id: hal-05271734

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

Submitted on 23 Sep 2025

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# Evolving magnetic lives of Sun-like stars

## I. Characterisation of the large-scale magnetic field with Zeeman-Doppler imaging

S. Bellotti<sup>1,2,*</sup>①, T. Lüftinger<sup>3</sup>, S. Boro Saikia<sup>4,†</sup>, C. P. Folsom<sup>5</sup>, P. Petit<sup>2</sup>, J. Morin<sup>6</sup>, M. Güdel<sup>4,‡</sup>, J.-F. Donati<sup>2</sup>, and E. Alecian<sup>7</sup>⑦

Received 13 May 2025 / Accepted 14 July 2025

### ABSTRACT

*Context.* Planets orbiting young, solar-type stars are embedded in a more energetic environment than that of the solar neighbourhood. They experience harsher conditions due to enhanced stellar magnetic activity and wind shaping the secular evolution of a planetary atmosphere.

*Aims.* This study is dedicated to the characterisation of the magnetic activity of eleven Sun-like stars, with ages between 0.2 and 6.1 Gyr and rotation periods between 4.6 and 28.7 d. Based on a sub-sample of six stars, we aim to study the large-scale magnetic field, which we then use to simulate the associated stellar wind and environment. Finally, we want to determine the conditions during the early evolution of planetary habitability.

*Methods.* We analysed high-resolution spectropolarimetric data collected in 2018 and 2019 with Narval. We computed activity diagnostics from chromospheric lines such as Ca II H&K, Hα, and the Ca II infrared triplet, as well as the longitudinal magnetic field from circularly polarised least-squares deconvolution profiles. For six stars exhibiting detectable circular polarisation signals, we reconstructed the large-scale magnetic field at the photospheric level by means of Zeeman-Doppler imaging (ZDI).

*Results.* In agreement with previous studies, we found a global decrease in the activity indices and longitudinal field with increasing age and rotation period. The large-scale magnetic field of the six sub-sample stars displays a strength between 1 and 25 G and reveals substantial contributions from different components such as poloidal (40–90%), toroidal (10–60%), dipolar (30–80%), and quadrupolar (10–40%), with distinct levels of axisymmetry (6–84%) and short-term variability of the order of months. Ultimately, this implies that exoplanets tend to experience a broad variety of stellar magnetic environments after their formation.

**Key words.** techniques: polarimetric – stars: activity – stars: magnetic field

## 1. Introduction

Stellar magnetic activity powers outflows and eruptive processes such as flares, coronal mass ejections, and energetic particle events. The collective action of these phenomena determines the space weather planets orbit within and plays a major role in the processing of planetary magnetospheres and upper atmospheres (e.g. Lammer et al. 2003, 2012; Rugheimer et al. 2015). Magnetic interactions in the outer stellar layers (chromosphere, transition region, and corona) produce energetic X-ray and extreme ultraviolet (EUV) radiation that can heat and ionise upper planetary atmospheres, triggering their evaporation and photochemistry (Cecchi-Pestellini et al. 2006; Murray-Clay et al. 2009; Owen & Jackson 2012; Tsai et al. 2023; Van Looveren et al. 2025). A similar effect is also induced by stellar cosmic rays

(Rodgers-Lee et al. 2021a,b; Raeside et al. 2025). Ultimately, the evolution of stellar activity during the history of a star dictates any observed planetary atmosphere (Penz & Micela 2008; Sanz-Forcada et al. 2011; Locci et al. 2019; Allan & Vidotto 2019). This is fundamental, for example, to explain the atmospheric evolution of Solar System planets (Lammer et al. 2003; Wood 2004; Kulikov et al. 2007; Airapetian & Usmanov 2016).

In this context, studying Sun-like stars at different ages is relevant to understanding planetary atmospheric evolution at early stages of formation, as well as the consequences of such extreme conditions for the habitability of Earth-like planets (See et al. 2014; Airapetian et al. 2020). Observations of young stars indicate an enhanced magnetic activity and a harsher space environment compared to the current Sun (Wood 2006), implying more frequent and more energetic phenomena (Wood 2004, 2006; Güdel 2007). For instance, X-ray and extreme ultraviolet (EUV) emission can be two or three orders of magnitude larger than

<footnote n="1">¹ Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands</footnote>

<footnote n="2">² Institut de Recherche en Astrophysique et Planétologie, Université de Toulouse, CNRS, IRAP/UMR 5277, 14 avenue Edouard Belin, F-31400 Toulouse, France</footnote>

<footnote n="3">³ Science Division, Directorate of Science, European Space Research and Technology Centre (ESA/ESTEC), Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands</footnote>

<footnote n="4">⁴ University of Vienna, Department of Astrophysics, Türkenschanzstrasse 17, A-1180 Vienna, Austria</footnote>

<footnote n="5">⁵ Tartu Observatory, University of Tartu, Observatoorium 1, Tõravere 61602, Estonia</footnote>

<footnote n="6">⁶ Laboratoire Univers et Particules de Montpellier, Université de Montpellier, CNRS, F-34095 Montpellier, France</footnote>

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

<footnote>* Corresponding author: bellotti@strw.leidenuniv.nl</footnote>
in the case of our Sun (Guinan et al. 2003; Ribas et al. 2005; Tu et al. 2015; Johnstone et al. 2021; Ketzer & Poppenhaeger 2023).

The nature of stellar activity is dictated by the intensity and structure of the magnetic field, which, in turn, depends on the rotation and internal structure of the star. Knowing the intensity of the magnetic fields on young planetary host stars, which can be several tens to hundreds of times stronger than that of the Sun (e.g. Folsom et al. 2018), is a pivotal precondition to understanding whether planets can become habitable in the early evolution of a planetary system. Moreover, the habitability conditions are not solely determined by the star's temperature, but also by its magnetic field (Vidotto et al. 2013; Cockell et al. 2016), which can influence the survival of the planetary atmosphere (see Fig. 5 in Van Looveren et al. 2025, for a recent example).

The level of stellar magnetic activity is evaluated via proxies such as photometric variability, X-ray luminosity, and emission in the cores of chromospheric lines. The latter are typically Ca II H&K, Hα, and Ca II infrared triplet lines (Wilson 1968; Gizis et al. 2002; Busà et al. 2007). A detailed characterisation of the stellar magnetic field can be carried out with spectropolarimetry, an observational technique that detects magnetically induced polarisation in spectral lines owing to the Zeeman effect (Zeeman 1897; Donati et al. 1997; Landstreet 2009). By collecting a time series of polarised spectra, it is possible to reconstruct the geometry of the large-scale magnetic field at the surface of the star by means of tomographic inversion, that is, Zeeman-Doppler imaging (ZDI; Semel 1989; Donati & Brown 1997). The resulting maps of the magnetic field can then be used as boundary conditions to extrapolate the field lines outward and model the stellar wind via magnetohydrodynamical simulations (e.g. Vidotto et al. 2014a, 2023; Nicholson et al. 2016; Boro Saikia et al. 2020; Folsom et al. 2020; Alvarado-Gómez et al. 2022; Evensberget et al. 2021, 2022, 2023).

The ZDI technique has been extensively applied to study the strength and configuration of the large-scale magnetic field of Sun-like stars of different ages and rotation periods (Donati & Collier Cameron 1997; Petit et al. 2008; Rosén et al. 2016; Folsom et al. 2016, 2018; Willamo et al. 2022). Distinct trends have emerged: for instance, slowly rotating Sun-like stars tend to exhibit poloidal-dominated field geometries, as shown by Petit et al. (2008) and Brown et al. (2022), together with a decrease in chromospheric activity. Furthermore, the average magnetic field strength decreases with increasing age and rotation period (e.g. Vidotto et al. 2014a; Folsom et al. 2018; Willamo et al. 2022), which is in agreement with trends of the total, unsigned magnetic field measurements performed by theoretical modelling of Zeeman broadening (Reiners 2012; Kochukhov et al. 2020). Repeated reconstructions of the large-scale magnetic field with ZDI have also revealed the presence of long-term trends and magnetic cycles for these stars (e.g. Fares et al. 2009; Rosén et al. 2016; Boro Saikia et al. 2016; Jeffers et al. 2017; Alvarado-Gómez et al. 2018; Bellotti et al. 2025).

In this work, we characterise the magnetic field and activity level for a sample of Sun-like stars whose ages span between 0.2 and 6.1 Gyr and with their rotation periods between 4.6 and 28.7 d. Thus, we are set to continue the exploration of magnetic activity properties across a wide stellar parameter space. For a sub-sample of six stars having detectable Zeeman signatures in circular polarisation and suitable temporal sampling of their rotation, we reconstructed the large-scale magnetic field geometry with ZDI. These six stars have ages between 0.2 and 1.6 Gyr

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Rotational evolutionary tracks computed by Johnstone et al. (2015a,b) and used in Tu et al. (2015), with our sample stars overplotted. Magenta data points indicate stars with a ZDI map both in the literature and presented in this work, cyan data points indicate stars with a first ZDI reconstruction in this work, and yellow data points indicate stars without a ZDI map. Some of the error bars are smaller than the data point size. The curves represent the 10th (red), 50th (green), and 90th (blue) percentiles of the stellar rotational distribution, of the envelope (solid lines) and core (dotted lines). The small horizontal lines are the observational constraints for the corresponding percentiles. The dashed black line is a time-dependent saturation relation for the stellar dipolar field, the wind mass loss rate and the X-ray luminosity. For more information, see the work of Tu et al. (2015).</figcaption></figure>

and rotation periods between 4.6 and 12.4 d. In this age range, the relation between rotation period and age is not unique, thus, monitoring many stars in this age range is important. The magnetic field maps resulting from this study will be used in a subsequent work for stellar wind and environment modelling.

The paper is structured as follows. In Sect. 2, we describe the sample of Sun-like stars examined in this study and their spectropolarimetric observations performed with the high-resolution spectropolarimeter Narval. In Sect. 3, we outline the computation of canonical activity indices as well as the longitudinal magnetic field. In Sect. 4, we recall the principles and assumptions of ZDI. Finally, we discuss our results star-by-star and general trends in Sects. 5 and 6, presenting our conclusions in Sect. 7.

## 2. Observations

The targets of our study are 11 G-type stars. According to the literature, they exhibit considerably distinct activity levels, along with likely differences in the magnetic field structure, wind properties, and high-energy output. We placed our stars on the rotation evolutionary tracks of Tu et al. (2015) to highlight this point, as illustrated in Fig. 1. The cited authors showed how the temporal evolution of stellar rotation, along with the associated X-ray and EUV luminosity, has a significant impact on the atmospheric evolution of exoplanets. Our stellar sample covers the three rotational evolution tracks (slow, medium, and fast), providing good coverage of the non-unique rotational history of Sun-like stars during the first few hundred million years. The Sun could have taken any of these three paths, so this coverage of the rotational evolution is necessary to remark on the magnetic history of the Sun.
<table><caption>**Table 1.** Properties of our stars.</caption><thead><tr><th>Name</th><th>B<br>[mag]</th><th>V<br>[mag]</th><th>Dist<br>[pc]</th><th>T<sub>eff</sub><br>[K]</th><th>log g<br>[dex]</th><th>M<br>[M<sub>☉</sub>]</th><th>R<br>[R<sub>☉</sub>]</th><th>Age<br>[Gyr]</th><th>P<sub>rot</sub><br>[d]</th><th>v<sub>eq</sub> sin i<br>[km s<sup>-1</sup>]</th><th>Inc<br>[°]</th></tr></thead><tbody><tr><td>HD 1835</td><td>7.05<sup>(1)</sup></td><td>6.39<sup>(1)</sup></td><td>21.33<sup>(8)</sup></td><td>5837<sup>(9)</sup></td><td>4.47<sup>(9)</sup></td><td>0.980<sup>(9)</sup></td><td>0.958<sup>(9)</sup></td><td>0.40±0.10<sup>(11)</sup></td><td>7.68<sup>(13)</sup></td><td>7.0<sup>(9)</sup></td><td>65<sup>(24)</sup></td></tr><tr><td>HD 28205</td><td>7.95<sup>(2)</sup></td><td>7.40<sup>(2)</sup></td><td>46.98<sup>(8)</sup></td><td>6234<sup>(10)</sup></td><td>4.34<sup>(10)</sup></td><td>1.210<sup>(10)</sup></td><td>1.228<sup>(10)</sup></td><td>0.63±0.05<sup>(12)</sup></td><td>5.87<sup>(15)</sup></td><td>10.6*</td><td>...</td></tr><tr><td>HD 30495</td><td>6.14<sup>(3)</sup></td><td>5.50<sup>(3)</sup></td><td>13.23<sup>(8)</sup></td><td>5759<sup>(9)</sup></td><td>4.39<sup>(9)</sup></td><td>0.860<sup>(9)</sup></td><td>0.983<sup>(9)</sup></td><td>0.40±0.10<sup>(11)</sup></td><td>11.36±0.17<sup>(16)</sup></td><td>4.1<sup>(16)</sup></td><td>...</td></tr><tr><td>HD 43162</td><td>7.07<sup>(1)</sup></td><td>6.37<sup>(1)</sup></td><td>16.69<sup>(8)</sup></td><td>5617<sup>(10)</sup></td><td>4.54<sup>(10)</sup></td><td>0.990<sup>(10)</sup></td><td>0.889<sup>(10)</sup></td><td>0.32±0.04<sup>(12)</sup></td><td>8.50<sup>+1.3(17)</sup><sub>-0.8</sub></td><td>7.0<sup>†</sup></td><td>70<sup>*</sup></td></tr><tr><td>HD 82443</td><td>7.78<sup>(3)</sup></td><td>7.05<sup>(3)</sup></td><td>18.07<sup>(8)</sup></td><td>5287<sup>(10)</sup></td><td>4.57<sup>(10)</sup></td><td>0.910<sup>(10)</sup></td><td>0.819<sup>(10)</sup></td><td>0.25±0.05<sup>(12)</sup></td><td>5.34±0.07<sup>(14)</sup></td><td>6.5<sup>(14)</sup></td><td>60<sup>(14)</sup></td></tr><tr><td>HD 114710</td><td>4.84<sup>(4)</sup></td><td>4.25<sup>(5)</sup></td><td>9.19<sup>(8)</sup></td><td>6075<sup>(9)</sup></td><td>4.57<sup>(9)</sup></td><td>1.540<sup>(9)</sup></td><td>1.063<sup>(9)</sup></td><td>1.60±0.40<sup>(11)</sup></td><td>12.35<sup>(18)</sup></td><td>4.1<sup>(19)</sup></td><td>70<sup>*</sup></td></tr><tr><td>HD 149026*</td><td>8.75<sup>(1)</sup></td><td>8.14<sup>(1)</sup></td><td>76.22<sup>(8)</sup></td><td>6084<sup>(10)</sup></td><td>4.17<sup>(10)</sup></td><td>1.140<sup>(10)</sup></td><td>1.460<sup>(10)</sup></td><td>2.78±0.50<sup>(11)</sup></td><td>&lt; 12.30 ± 1.20</td><td>6.0<sup>(20)</sup></td><td>...</td></tr><tr><td>HD 189733*</td><td>8.58<sup>(6)</sup></td><td>7.65<sup>(6)</sup></td><td>19.78<sup>(8)</sup></td><td>5023<sup>(10)</sup></td><td>4.58<sup>(10)</sup></td><td>0.776<sup>(10)</sup></td><td>0.840<sup>(10)</sup></td><td>0.40±0.10<sup>(11)</sup></td><td>11.94±0.16<sup>(21)</sup></td><td>2.97<sup>(21)</sup></td><td>75<sup>(21)</sup></td></tr><tr><td>HD 190406</td><td>6.39<sup>(1)</sup></td><td>5.79<sup>(1)</sup></td><td>17.77<sup>(8)</sup></td><td>5932<sup>(9)</sup></td><td>4.45<sup>(9)</sup></td><td>1.131<sup>(9)</sup></td><td>1.052<sup>(9)</sup></td><td>1.70±0.40<sup>(11)</sup></td><td>13.94<sup>(18)</sup></td><td>4.0*</td><td>...</td></tr><tr><td>HD 206860*</td><td>5.29<sup>(7)</sup></td><td>4.70<sup>(7)</sup></td><td>18.13<sup>(8)</sup></td><td>5974<sup>(9)</sup></td><td>4.47<sup>(9)</sup></td><td>1.085<sup>(9)</sup></td><td>1.002<sup>(9)</sup></td><td>0.20±0.10<sup>(11)</sup></td><td>4.55±0.03<sup>(22)</sup></td><td>10.1<sup>(22)</sup></td><td>75<sup>(22)</sup></td></tr><tr><td>HD 219828*</td><td>8.78<sup>(1)</sup></td><td>8.01<sup>(1)</sup></td><td>72.53<sup>(8)</sup></td><td>5807<sup>(10)</sup></td><td>4.02<sup>(10)</sup></td><td>1.040<sup>(10)</sup></td><td>1.651<sup>(10)</sup></td><td>6.08<sup>+0.9(11)</sup><sub>-0.6</sub></td><td>28.70±0.60<sup>(23)</sup></td><td>3.7<sup>(23)</sup></td><td>...</td></tr></tbody></table>

**Notes.** The columns indicate: identifier of the star, B band magnitude, V band magnitude, distance, effective temperature, surface gravity, mass, radius, age, rotation period, projected equatorial velocity ($v_{eq} \sin i$), and stellar inclination. References: (1) Høg et al. (2000), (2) Joner et al. (2006), (3) SIMBAD (Wenger et al. 2000), (4) Ducati (2002), (5) van Belle & von Braun (2009), (6) Koen et al. (2010), (7) Boro Saikia et al. (2015), (8) Gaia Collaboration (2020), (9) Valenti & Fischer (2005), (10) Paegert et al. (2021), (11) Ramírez et al. (2012), (12) Linsky et al. (2020), (13) Güdel (2007), (14) Folsom et al. (2016), (15) Pizzolato et al. (2003), (16) Egeland et al. (2015), (17) Marsden et al. (2014), (18) Wright et al. (2011), (19) Gray & Baliunas (1997), (20) Bonomo et al. (2017), (21) Fares et al. (2010), (22) Boro Saikia et al. (2015), (23) Santos et al. (2016), (24) Rosén et al. (2016). The symbol ⋆ in the name of a star represents a known exoplanet host. The symbol * indicates values of $v_{eq} \sin i$ or inclination that we derived in this work from geometrical considerations. We estimated $v_{eq} \sin i$ as the ratio between the radius and rotation period of the star. For stars like HD 28205 and HD 190406, these values are only references and were not used for the reconstructions because the stars have either no detected Zeeman signatures or too few observations. The symbol † means that we derived the value by means of ZDI optimisation.

Four stars in our sample (HD 1835, HD 82443, HD 189733, and HD 206860) have already been observed with spectropolarimetry, but these observations date back to different epochs. Therefore, our analyses also address the evolution of the magnetic field structure on timescales of many years, as reported previously for other Sun-like stars (Donati et al. 2008; Boro Saikia et al. 2015; Rosén et al. 2016; Willamo et al. 2022; Bellotti et al. 2025).

The properties of our stars are summarised in Table 1. The ages of the stars were mostly extracted from the work of Ramírez et al. (2012) and Linsky et al. (2020). Furthermore, Ramírez et al. (2012) used either isochrones or gyrochronology to estimate the stellar age. The latter method is more accurate for young active stars. For our stars, the ages of HD 149026 and HD 219828 were obtained via isochrones, while for the remaining stars, ages were obtained from gyrochronology.

We analysed optical spectropolarimetric observations collected with Narval in 2018 and 2019¹. Narval is the spectropolarimeter on the 2 m Télescope Bernard Lyot (TBL) at the Pic du Midi Observatory in France (Donati 2003), which operates between 370 and 1050 nm at high resolution ($R = 65,000$). The observations were carried out in circular polarisation mode, providing both unpolarised (Stokes I) and circularly polarised (Stokes V) high-resolution spectra. The data were reduced with the LIBRE-ESPRIT pipeline (Donati et al. 1997) and the continuum-normalised spectra were retrieved from PolarBase (Petit et al. 2014). We provide the full list of observations in Table A.1 and the total number of observations for each star in Table B.1.

The detection and characterisation of Zeeman signatures in circularly polarised light are performed by means of least-squares deconvolution (LSD; Donati et al. 1997; Kochukhov et al. 2010) using the LSDPY code which is part of

the Specpolflow software (Folsom et al. 2025)<fnref n="2" />. This numerical technique produces high signal-to-noise ratio (S/N) line profiles (unpolarised and circularly polarised) from the combination of thousands of photospheric spectral lines included in a synthetic line list. The line lists were produced using the Vienna Atomic Line Database<fnref n="3" /> (VALD, Ryabchikova et al. 2015). They contain information of atomic lines with known Landé factor (indicated by $g_{eff}$ and describing the magnetic sensitivity of a spectral line) and with a depth greater than 40% the level of the unpolarised continuum (following Donati et al. 1997). A summary of the line lists used and their properties is given in Table B.1.

We recorded substantially lower S/N in Stokes V LSD profiles for two observations of HD 206860 on July 26th 2018 and June 22nd 2019, as well as for HD 43162 on December 19th 2018; hence, these observations were excluded from the analysis. For HD 206860, we also noticed that the Stokes V profile on July 28, 2018 featured several oscillations around the zero flux level and was also excluded. For three stars, namely, HD 149026, HD 190406, and HD 219828, we did not report any clear detections of Zeeman signatures in Stokes V profiles. Indeed, for each star, the false-alarm probability (FAP; see Donati et al. 1997, for more details) of the putative Zeeman signature is greater than $10^{-3}$. It was only in a couple of observations that we were able to see a marginal detection (FAP ~ $10^{-3} - 10^{-4}$). For this reason, these stars were excluded from the spectropolarimetric characterisation analyses, however, they were retained for the measurement of unpolarised magnetic activity indicators, as outlined in the next sections. In Appendix D, we describe the polarisation signatures in Stokes V as well as spurious polarisation signatures in Stokes N. In the next following, the observations have been phased with the following ephemeris,

$$ \text{HJD} = \text{HJD}_0 + \text{P}_{\text{rot}} \cdot n_{\text{cyc}}, \qquad (1) $$

<footnote n="1">¹ The data is available at https://www.polarbase.ovgso.fr/</footnote>

<footnote n="2">² Available at https://github.com/folsomcp/LSDpy</footnote>

<footnote n="3">³ http://vald.astro.uu.se/</footnote>
where HJD₀ is the heliocentric Julian Date reference (the first one of the time series for each star, see Table A.1), P<sub>rot</sub> is the rotation period of the star (see Table 1), and n<sub>cyc</sub> represents the rotation cycle.

### 3. Activity indices

To gauge the magnetic activity level of our stars, we computed time series of activity indices. We used chromospheric spectral lines falling in the optical domain covered by Narval, namely the Ca II H&K lines, Hα, and Ca II infrared triplet lines. In Table C.1, we list the main statistical features of the time series for all three activity indicators.

The log R'<sub>HK</sub> quantifies the amount of flux in the Ca II H&K lines relative to the nearby continuum, without colour dependence and photospheric contribution (Middelkoop 1982; Noyes et al. 1984; Rutten 1984). The recipe to compute the index consists of three steps: (i) measuring the S index, (ii) calibrating it to the Mount Wilson scale, and (iii) converting it to the log R'<sub>HK</sub> index. Following the definition of Vaughan et al. (1978), the S index is measured via

$$S = \frac{aF_H + bF_K}{cF_R + dF_V} + e, \qquad (2)$$

where F<sub>H</sub> and F<sub>K</sub> are the fluxes in two triangular band passes with FWHM = 1.09 Å centred on the cores of the H line (3968.470 Å) and K line (3933.661 Å), whereas F<sub>R</sub> and F<sub>V</sub> are the fluxes within two 20-Å rectangular band passes centred at 3901 and 4001 Å, respectively. The set of coefficients {a, b, c, d, e} is used to calibrate the S index from a specific instrument scale to the Mount Wilson scale, and were estimated by Marsden et al. (2014) for Narval. We then used the formula by Rutten (1984) to convert the S index to the log R'<sub>HK</sub> index.

The Hα index is measured as

$$H\alpha = \frac{F_{Ha}}{H_R + H_V}, \qquad (3)$$

where F<sub>Ha</sub> is the flux within a rectangular band pass of 3.60 Å centred on the Hα line at 6562.85 Å, and H<sub>V</sub> and H<sub>R</sub> are the fluxes within two rectangular bandpasses of 2.2 Å centred on 6558.85 Å and 6567.30 Å (Gizis et al. 2002). Finally, we followed Petit et al. (2013) and Marsden et al. (2014) to measure the Ca II infrared triplet index as

$$CaIIRT = \frac{IR1 + IR2 + IR3}{IR_R + IR_V}, \qquad (4)$$

where IR1, IR2, and IR3 are the fluxes within a rectangular band pass of 2 Å centred on the Ca II lines at 8498.023 Å, 8542.091 Å, and 8662.410 Å, respectively, and IR<sub>R</sub> and IR<sub>V</sub> are the fluxes within two rectangular bandpasses of 5 Å centred on 8704.9 Å and 8475.8 Å.

#### 3.1. Longitudinal magnetic field

We computed the disc-integrated, line-of-sight projected component of the large-scale magnetic field following Rees & Semel (1979). Formally,

$$B_l[G] = \frac{-2.14 \cdot 10^{11} \int vV(v)dv}{\lambda_0 g_{\text{eff}} c \int (I_c - I)dv}, \qquad (5)$$

where λ₀ (in nm) and g<sub>eff</sub> are the normalisation wavelength and Landé factor of the LSD profiles, I<sub>c</sub> is the unpolarised continuum level, v is the radial velocity associated to a point in the spectral line profile in the star's rest frame (in km s⁻¹), and c the speed of light in vacuum (in km s⁻¹). For all our stars, we set the normalisation parameters to λ₀ = 700 nm and g<sub>eff</sub> = 1.24. The velocity range over which the integration is carried out encompass the width of both Stokes I and V LSD profiles. In Table B.1, we list the velocity range for each star.

### 4. Magnetic imaging

We applied ZDI to reconstruct the large-scale magnetic field topology for the stars in our study. The magnetic field vector is expressed as the sum of a poloidal and toroidal component, each described via spherical harmonics formalism (Lehmann & Donati 2022). The algorithm inverts a time series of Stokes V LSD profiles into a magnetic field map while applying a maximum-entropy regularisation scheme (Skilling & Bryan 1984). In practice, ZDI fits the spherical harmonics coefficients α<sub>ℓ,m</sub>, β<sub>ℓ,m</sub>, and γ<sub>ℓ,m</sub> (with ℓ and m being the degree and order of the mode, respectively) in an iterative fashion, until a target reduced χ² is reached (for more information see Skilling & Bryan 1984; Semel 1989; Donati et al. 1997). The coefficients α<sub>ℓ,m</sub>, β<sub>ℓ,m</sub>, and γ<sub>ℓ,m</sub> are complex numbers used to describe the radial poloidal, tangential poloidal, and toroidal field components, respectively. The algorithm searches for the maximum-entropy solution, that is, the magnetic field configuration compatible with the data and with the lowest information content.

We employed the zdipy code described in Folsom et al. (2018) and adopted the weak-field approximation, for which Stokes V is proportional to the derivative of Stokes I with respect to wavelength (e.g. Landi Degl'Innocenti 1992). As outlined in Folsom et al. (2018), the local unpolarised line profiles are modelled with a Voigt kernel. For each star, we performed a χ² minimisation between the median of the observed Stokes I LSD profiles and its model over a grid of line depth, Gaussian width and Lorentzian width values. The optimal values we used in the ZDI reconstruction are listed in Table D.1.

We set the limb darkening coefficient to 0.6964 (Claret & Bloemen 2011) and the maximum degree of spherical harmonic coefficients to l<sub>max</sub> = 10 for every ZDI reconstruction. The latter is consistent with the projected equatorial velocity (v<sub>eq</sub> sin i) for our stars, which is at most 10.6 km s⁻¹. v<sub>eq</sub> sin(i) determines l<sub>max</sub> proportionally and the amount of magnetic field complexity one can image (Hussain et al. 2009). In fact, faster rotation implies more separated Zeeman signatures in radial velocity space, thus limiting polarity cancellation. In our case, we homogeneously set l<sub>max</sub> equal to 10, but a lower value could have been used without changing the results, as most of our stars are slow rotators and the magnetic energy is stored in the low-l degrees.

In cases where the stellar inclination was not available in the literature, it was estimated from geometrical considerations. More precisely, it was computed by comparing the stellar radius provided in the literature with the projected radius, R sin i = P<sub>rot</sub>v<sub>eq</sub> sin i/50.59, where R sin i is measured in solar radii, P<sub>rot</sub> is the stellar rotation period measured in days, and the v<sub>eq</sub> sin i unit is km s⁻¹. When the estimated inclination was larger than 80°, we adopted a value of 70° to conservatively prevent mirroring effects between the northern and southern hemispheres.
The ZDIPY code includes differential rotation as a function of colatitude (θ) expressed as

$$ \Omega(\theta) = \Omega_{eq} - d\Omega \sin^2(\theta), \qquad (6) $$

with $\Omega_{eq} = 2\pi/P_{rot}$ is the rotational frequency at the equator and $d\Omega$ is the differential rotation rate in rad d⁻¹. Following the parameter optimisation outlined in Donati et al. (2000) and Petit et al. (2002), we generated a grid of ($P_{rot}$, $d\Omega$) pairs and searched for the values that minimised the $\chi^2$ distribution between observations and synthetic LSD profiles, at a fixed entropy level. The best-fit parameters were measured by fitting a 2D paraboloid to the $\chi^2$ distribution, while the error bars were obtained from a variation of $\Delta\chi^2 = 1$ away from the minimum (Press et al. 1992; Petit et al. 2002).

Once the stellar input parameters were fixed, it was necessary to determine a target $\chi^2$ that the ZDI algorithm is meant to converge to. This value represents the optimal fit quality of the Stokes V LSD profiles and should be the lowest possible. At the same time, it should be large enough to prevent overfitting the noise in the observations, which would lead to spurious artefacts in the final ZDI map. As outlined in Alvarado-Gómez et al. (2015), finding the target $\chi^2$ is achieved by running ZDI over a grid of $\chi^2$ values, each time recording the entropy at convergence. In this way, we can observe the change in the information content of the ZDI map with $\chi^2$ and identify the target $\chi^2$ as the value corresponding to the maximum of the change rate in the entropy.

## 5. Results

In this section, we describe the computation of activity indices and longitudinal magnetic field for each star individually. The average measurements are reported in Table C.1, while all the values are listed in Table A. We also describe the magnetic field reconstructions with ZDI and (when possible) we compare them with previous maps available in the literature. The ZDI maps are given in Figs. 3 and 4 for the radial, azimuthal, and meridional components of the magnetic field. The properties of the ZDI maps and the results of the differential rotation search are summarised in Table 2. The table reports only the stars for which a Zeeman signature in Stokes V was detectable and with a number of observations producing sufficient rotational phase coverage to perform the tomographic inversion reliably.

We find average magnetic field strengths ranging from 1 to 25 G. The magnetic field topologies are predominantly poloidal (>60%), but there are cases where the toroidal component is dominant (63%) or significant (>20%). Most stars exhibit complex magnetic fields, as the dipolar component accounts for at most 60% of the poloidal energy, and the quadrupolar and octupolar modes store between 10–40% and 3–15% of the energy. The axisymmetry of the reconstructed topologies also varies significantly within our sample, since the fraction of total energy in the corresponding modes spans between 6% and 84%, with two out of six stars featuring relatively axisymmetric magnetic orientations (65–73% and 64–84%). These results are in agreement with the expected magnetic topology for fast-rotating stars with Sun-like properties and interior, that is, with a convective envelope surrounding a radiative core (Petit et al. 2005; Donati & Landstreet 2009; Folsom et al. 2016; Rosén et al. 2016; Folsom et al. 2018; Willamo et al. 2022; Bellotti et al. 2025).

<figure><img src="image_3.png" /><figcaption>Fig. 2. Joint search of differential rotation and equatorial rotation period for HD 1835, HD 82443, and HD 206860. The panels illustrate the $\chi^2$ landscape over a grid of ($P_{rot,eq}$, $d\Omega$) pairs, with the $1\sigma$ and $3\sigma$ contours. The best values are obtained by fitting a 2D paraboloid around the minimum, while their error bars are estimated from the projection of the $1\sigma$ contour on the respective axis (Press et al. 1992).</figcaption></figure>

### 5.1. HD 1835 (BE Cet)

HD 1835 is a G2.5 dwarf (Keenan & McNeil 1989) at a distance of 21.3 pc (Gaia Collaboration 2020). It is located in the Hyades cluster (Montes et al. 2001) and its age was estimated to be 400 Myr (Güdel 2007; Ramírez et al. 2012), and its rotation period to 7.68 d (Güdel 2007).
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Reconstructed large-scale magnetic field map of HD 1835, HD 43162, and HD 82443 in flattened polar view. From the left, the radial, azimuthal, and meridional components of the magnetic field vector are illustrated. The radial ticks are located at the rotational phases when the observations were collected (see Eq. (1)), while the concentric circles represent different stellar latitudes: -30°, +30°, and +60° (dashed lines), as well as the equator (solid line).</figcaption></figure>

We measured $\log R'_{HK}$ values between -4.414 and -4.338, with a median of -4.384 and a standard deviation of 0.025. These values are consistent with those reported in Boro Saikia et al. (2018) and Brown et al. (2022). The $H\alpha$ index varies between 0.321 and 0.328, with a median of 0.326 and standard deviation of 0.002, and the Ca II IRT ranges 0.853–0.890 with a median value of 0.866 ± 0.010.

The longitudinal field is between -5.7 and 2.0 G, with a median value of -2.2 ± 2.5 G (the median error bar is 1.0 G). In comparison, the range reported by Rosén et al. (2016) and Willamo et al. (2022) from 2013 and 2017 HARPS-Pol observations is approximately -5 to 10 G. While consistent, our 2018 range of $B_l$ values is narrower, possibly indicating a weakening of the magnetic field over a timescale of a year.

We constrained differential rotation, as outlined in Sect. 4, finding that the optimised parameters are $P_{\text{rot}} = 7.57 \pm 0.02$ d and $d\Omega = 0.018 \pm 0.008 \text{ rad d}^{-1}$ (see Fig. 2). Using Eq. (6), we estimated a rotation period at the pole of 7.74 ± 0.08 d; thus, the rotation period of 7.68 d reported in Güdel (2007) is consistent with the expected equator-pole range of rotation rates. The initial $\chi_r^2$ was 6.5, which is associated to zero magnetic field map. The target $\chi_r^2$ of the reconstruction improved to 2.0 when solid body rotation was assumed, and down to 1.4 when differential rotation was included. The deviation of $\chi_r^2$ from 1.0 is symbolic of intrinsic variability that occurred during the time span of our observations and which is unaccounted for by the ZDI model.

The time series for the modelled Stokes V profiles are given in Appendix D.

Our ZDI reconstruction of HD 1835 is shown in Fig. 3, exhibiting low-latitude, almost equatorial magnetic structures. The field has a predominantly poloidal component (91% of the total magnetic energy), of which 32%, 28%, and 37% is in the dipolar, quadrupolar and octupolar modes. The field is largely non-axisymmetric, since the axisymmetric component (that is, $\ell \ge 1$ and $m = 0$) accounts for 6% of the total energy. The reconstructed field shows similarities with the 2013 map of Rosén et al. (2016) and the 2017 map of Willamo et al. (2022), both in terms of dominant component and axisymmetry. Compared to Rosén et al. (2016) reconstruction, our map manifests an increased quadrupolar component at the expense of the dipolar component, indicating a higher complexity. Compared to Willamo et al. (2022), our ZDI reconstruction of the radial field features a more extended area covered by negative polarity.

The average magnetic field strength we recovered is 7 G, which is approximately half the value from previous reconstructions. This goes in the same direction as the decreased range of variability of the longitudinal field. For HD 1835, Egeland (2017) reported two superimposed activity cycles with two main periodicities of 7.8 yr and 20.8 yr. In a similar manner, Boro Saikia et al. (2018) estimated periodicities of 9.06 yr and 22 yr. Although there are three large-scale magnetic field maps reconstructed for this star, they are more recent than the temporal baseline of S-index values used by Egeland (2017) and Boro Saikia et al. (2018), which complicates the comparison of the field topology and the phase of the S-index cycle.

## 5.2. HD 28205 (V993 Tau)

HD 28205 is a G0.0 dwarf (Linsky et al. 2020) at a distance of 46.98 pc (Gaia Collaboration 2020). It is located in the Hyades cluster (Kraft 1965) and its age was estimated to be 630 Myr (Linsky et al. 2020), and its rotation period to 5.87 d (Pizzolato et al. 2003).

We measured $\log R'_{HK}$ values between -4.453 and -4.424, with a median of -4.438 ± 0.009 (the median error is 0.036 in comparison). These values are about 0.06 dex larger than the average index given in Brown et al. (2022). The $H\alpha$ index is stable around 0.313 ± 0.001, and the Ca II IRT varies between 0.843 and 0.849 with a median value of 0.848 ± 0.002.

The longitudinal field was measured to be between -2.9 and 0.8 G, with a median value of -1.0 ± 1.3 G (the median error bar is 1.7 G). The value is on the same order of magnitude as the unsigned, average longitudinal field reported by Brown et al. (2022). Since we only had five observations for this star, mostly clustered around rotational phase 0.0, we did not perform temporal analyses or ZDI.

## 5.3. HD 30495 (IX Eri)

HD 30495 is a G2.5 dwarf (Gray et al. 2006) at a distance of 13.24 pc (Gaia Collaboration 2020). It has an age of 400 Myr (Ramírez et al. 2012) and a rotation period of 11.36 d (Egeland et al. 2015).

We measured $\log R'_{HK}$ values between -4.481 and -4.495, with a median of -4.490 ± 0.006 (the median error is 0.030 in comparison). These values are compatible with Boro Saikia et al. (2018) and Brown et al. (2022), but on the higher end of the range listed. This could be explained by the chromospheric activity cycle of the star, composed of two superimposed periodicities of 2 and 12 yr (Egeland et al. 2015;
<table><caption>**Table 2.** Properties of the magnetic maps.</caption><thead><tr><th>Star</th><th>Epoch</th><th>dΩ<br>[rad s<sup>-1</sup>]</th><th>⟨|B<sub>v</sub>|⟩<br>[G]</th><th>|B<sub>max</sub>|<br>[G]</th><th>⟨B<sup>2</sup>⟩<br>[×10<sup>2</sup> G<sup>2</sup>]</th><th>f<sub>pol</sub><br>[%]</th><th>f<sub>tor</sub><br>[%]</th><th>f<sub>dip</sub><br>[%]</th><th>f<sub>quad</sub><br>[%]</th><th>f<sub>oct</sub><br>[%]</th><th>f<sub>axi</sub><br>[%]</th><th>f<sub>axi,pol</sub><br>[%]</th><th>f<sub>axi,tor</sub><br>[%]</th></tr></thead><tbody><tr><td>HD 1835</td><td>2018.77</td><td>0.018 ± 0.008</td><td>7.4</td><td>29.3</td><td>0.78</td><td>91.1</td><td>8.9</td><td>32.1</td><td>28.4</td><td>27.4</td><td>6.4</td><td>6.9</td><td>1.2</td></tr><tr><td>HD 43162</td><td>2019.01</td><td>...</td><td>11.0</td><td>23.3</td><td>1.50</td><td>63.8</td><td>36.2</td><td>84.0</td><td>11.1</td><td>3.5</td><td>63.9</td><td>45.7</td><td>96.0</td></tr><tr><td>HD 43162</td><td>2019.22</td><td>...</td><td>16.4</td><td>36.1</td><td>3.33</td><td>36.9</td><td>63.1</td><td>62.5</td><td>25.8</td><td>7.9</td><td>83.9</td><td>61.8</td><td>96.8</td></tr><tr><td>HD 82443</td><td>2019.22</td><td>0.114 ± 0.022</td><td>23.6</td><td>65.1</td><td>7.05</td><td>63.9</td><td>36.1</td><td>34.3</td><td>46.8</td><td>9.2</td><td>36.1</td><td>4.6</td><td>92.0</td></tr><tr><td>HD 114710</td><td>2019.23</td><td>...</td><td>1.3</td><td>3.4</td><td>0.02</td><td>88.6</td><td>11.4</td><td>59.9</td><td>25.4</td><td>9.8</td><td>31.4</td><td>25.0</td><td>80.9</td></tr><tr><td>HD 189733</td><td>2018.89</td><td>0.110<sup>†</sup></td><td>18.8</td><td>56.5</td><td>5.50</td><td>52.1</td><td>47.9</td><td>67.7</td><td>16.1</td><td>8.2</td><td>43.6</td><td>1.5</td><td>89.5</td></tr><tr><td>HD 206860</td><td>2018.57</td><td>0.109 ± 0.004</td><td>20.8</td><td>47.7</td><td>5.10</td><td>85.5</td><td>14.5</td><td>46.8</td><td>16.8</td><td>13.5</td><td>64.6</td><td>61.0</td><td>85.7</td></tr><tr><td>HD 206860</td><td>2019.55</td><td>0.051 ± 0.042</td><td>15.8</td><td>35.3</td><td>2.80</td><td>69.2</td><td>30.8</td><td>59.4</td><td>21.0</td><td>10.5</td><td>73.4</td><td>65.2</td><td>91.9</td></tr></tbody></table>

**Notes.** The following quantities are listed: star's name, epoch of observations, differential rotation rate, mean unsigned magnetic strength, maximum unsigned magnetic strength, total reconstructed magnetic energy, poloidal and toroidal magnetic energy as a fraction of the total energy, dipolar, quadrupolar, and octupolar magnetic energy as a fraction of the poloidal energy, axisymmetric magnetic energy as a fraction of the total energy, poloidal axisymmetric energy as a fraction of the poloidal energy, and toroidal axisymmetric energy as a fraction of the toroidal energy. The symbol † indicates that the value from Fares et al. (2017) was used.

Brandenburg et al. (2017). More precisely, extrapolating the trend shown in Fig. 1 of Egeland et al. (2015), we note that the 2015 observations analysed by Brown et al. (2022) lie at a lower activity state than our 2018 observations. The $H\alpha$ index we computed is constant around 0.323 ± 0.001, and the Ca II IRT is stable around 0.855 ± 0.001.

The longitudinal field was measured to be between 0.4 and 3.1 G, with a median value of 1.6 ± 1.1 G (the median error bar is 1.1 G). The range is compatible with the unsigned, average longitudinal field reported by Brown et al. (2022). Since we only have three observations for this star, we did not perform any temporal analyses or ZDI.

## 5.4. HD 43162 (V352 CMa)

HD 43162 A is a G6.5 dwarf (Linsky et al. 2020), and it is the primary component of a triple system at a distance of 16.69 pc (Chini et al. 2014; Gaia Collaboration 2020). It has an age of 320 Myr (Linsky et al. 2020) and a rotation period of 8.50 d (Marsden et al. 2014).

We measured $\log R'_{HK}$ values between -4.398 and -4.315, with a median of -4.364 ± 0.023. While these values indicate a higher activity level than reported in Marsden et al. (2014) and Lehtinen et al. (2016) by about 0.06 dex, they are compatible with the measurements by Boro Saikia et al. (2018) and Brown et al. (2022). The $H\alpha$ index varies between 0.336 and 0.342, with a median of 0.339 ± 0.002, and the Ca II IRT ranges between 0.891 and 0.919 with a median value of 0.903 ± 0.007. The longitudinal field was measured to be between -4.9 and 7.2 G, with a median value of 2.4 ± 3.3 G (the median error bar is 1.5 G). This range encompasses the value reported in Marsden et al. (2014) as well as the unsigned average in Brown et al. (2022).

We carried out a ZDI reconstructions for both epochs of HD 43162. The maps are shown in Fig. 3 and the model Stokes V profiles in Fig. D.1. In this case, reliable Stokes V models of the LSD profiles could only be obtained with an optimisation of $v_{eq} \sin i$ based on $\chi^2$ minimisation. The idea is similar to the grid search used to find differential rotation (see Sect. 4); however, in this case, only a grid over $v_{eq} \sin i$ was used, while all the other parameters were kept fixed. We obtained a value of 7 km s$^{-1}$, which is lower than the value of 9.6 km s$^{-1}$ estimated by Valenti & Fischer (2005), but consistent with previous measurements by Zejda et al. (2012). We then searched for the optimised values of rotation period and differential rotation jointly, but the results were inconclusive; hence, we performed ZDI reconstructions assuming solid body rotation. When searching for the optimised value of rotation period alone, we found a $\chi^2$ minimum at the expected rotation period of 8.5 d (Lehtinen et al. 2020) for both epochs; hence, we used this value for the ZDI maps.

For the 2019.01 epoch, we set a target $\chi^2$ to 1.7 from an initial value of 6.2. We found an average field strength of 11 G, with the topology being 64% poloidal and 36% toroidal. The dipolar mode dominates the poloidal component by storing 84% of the magnetic energy, while the quadrupolar and octupolar account for 11% and 4%, respectively. The axisymmetric component is 64%. For 2019.22, the target $\chi^2$ is 1.5, from an initial value of 11.2. The average field strength is 16 G, and the field configuration is predominantly toroidal (63%). The poloidal component is still largely dipolar (63%), but the quadrupolar and octupolar components increased to 26% and 8%.

No polarity reversal is seen between the two epochs. While we cannot exclude an evolution of the field on short timescales, the difference in dominant topology between the two epochs could also be due to the scarce phase coverage affecting the 2019.01 epoch. Furthermore, we observe that the Stokes V signatures are asymmetric: a larger positive lobe than the negative one (see Fig. D.1). The spurious polarisation signature in Stokes N is negligible compared to the Stokes V signature; hence, it does not affect its shape and symmetry. The asymmetry in Stokes V has been encountered before for solar-like stars such as $\xi$ Boo A (Petit et al. 2005; Morgenthaler et al. 2012) as well as more massive stars such as $\theta$ Leo and $\epsilon$ UMa (Blazère et al. 2016). We note that this might be due to vertical gradients of photospheric velocity and magnetic field strength (López Ariste 2002). Interestingly, we only found this feature in the case of HD 43162.

## 5.5. HD 82443 (DX Leo)

HD 82443 is a G9.0 dwarf (Gray et al. 2003), and it is the primary of a binary system with an M6 dwarf (Lépine & Bongiorno 2007) at a distance of 18.07 pc (Gaia Collaboration 2020). It is a member of the Her-Lyr moving group (Gaidos 1998; Eisenbeiss et al. 2013), with an age of 250 Myr (Folsom et al. 2016; Linsky et al. 2020) and a rotation period of $P_{rot} = 5.377$ d (Folsom et al. 2016).
We measured $\log R'_{\text{HK}}$ values between -4.170 and -4.115, with a median of -4.147 ± 0.016 (the median error bar is 0.011). These values are 0.05 dex larger than the average values reported in Folsom et al. (2016), Boro Saikia et al. (2018), and Brown et al. (2022). The $H\alpha$ index varies between 0.375 and 0.383, with a median of 0.380 ± 0.003, and the Ca II IRT ranges 0.974–0.995 with a median value of 0.986 ± 0.006. The values of both indices are ~0.015 units higher than the values measured by (Folsom et al. 2016).

The longitudinal field was measured to be between -10.9 and 1.3 G, with a median value of -4.1 ± 8.4 G (the median error bar is 1.3 G). In comparison, Folsom et al. (2016) measured the longitudinal field to be between ±20 G, and Brown et al. (2022) obtained an average $B_l$ of 11 G in absolute value.

We applied ZDI and fit the Stokes V time series down to $\chi_r^2 = 2.2$ from an initial value of 35.1, when assuming solid body rotation. We then searched for differential rotation and constrained $P_{\text{rot}} = 5.460 \pm 0.026$ d and $d\Omega = 0.114 \pm 0.022$ rad d$^{-1}$ (see Fig. 2). We caution that measurements of the differential rotation with less than 15 observations across the stellar rotation may be biased due to rotational phase gaps (Petit et al. 2002). The rotation period is consistent with the value constrained with ZDI by Folsom et al. (2016) within $1\sigma$. Our differential rotation rate is ~2$\sigma$ larger than the lower limit estimated from photometric light curves by Messina et al. (1999). Using Eq. (6), we estimated a rotation period at the pole of 6.070 ± 0.130 d. By assuming our $P_{\text{rot}}$, $d\Omega$ pair, the target $\chi_r^2$ of the ZDI reconstruction was improved to 1.7.

The magnetic field topology is shown in Fig. 3 and the properties listed in Table 2. We obtained an average magnetic field strength of 23 G, which is consistent with the reconstruction of Folsom et al. (2016). The magnetic field topology is predominantly poloidal (64%), with 34%, 47%, and 9% of the magnetic energy stored in the dipolar, quadrupolar, and octupolar modes, respectively. The field is also non-axisymmetric (36%). In comparison, the reconstruction of Folsom et al. (2016) using 2013 data exhibits a predominantly dipolar (71%) and non-axisymmetric (8%) topology, indicating an increase in complexity with our recent observations. Considering the toroidal component specifically, the axisymmetric energy fraction increases from 38% in Folsom et al. (2016) to 92% in our reconstruction. This evolution occurred over a timescale of approximately 5 yr, which is similar on the same order of magnitude as the photometric activity cycle reported by Baliunas et al. (1995), Messina et al. (1999), Lehtinen et al. (2016). If we assume that the dynamo processes operating in the stellar interior are Sun-like, we could attribute the evolution of the axisymmetric-toroidal fraction to a correlated, equatorward distribution of starspots.

## 5.6. HD 114710 (β Com)

HD 114710 is a G0.0 dwarf (Keenan & McNeil 1989) at a distance of 9.19 pc (Gaia Collaboration 2020) and part of a visual binary system (Mason et al. 2001). It has an age of 1.60 Gyr (Ramírez et al. 2012) and a rotation period of $P_{\text{rot}} = 12.35$ d (Wright et al. 2011).

We measured $\log R'_{\text{HK}}$ values between -4.633 and -4.591, with a median of -4.609 ± 0.011 (the median error bar is comparable to the standard deviation). These values are compatible with the measurements reported by Boro Saikia et al. (2018) and Brown et al. (2022). The $H\alpha$ index varies between 0.312 and 0.315, with a median of 0.313 ± 0.001, and the Ca II IRT ranges 0.793–0.813 with a median value of 0.798 ± 0.006. The longi-

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Reconstructed large-scale magnetic field map of HD 114710, HD 189733, and HD 206860 for the two epochs. The format is the same as Fig. 3.</figcaption></figure>

tudinal field was measured to be between -2.7 and 0.7 G, with a median value of -0.3 ± 1.0 G (where the median error bar is 0.5 G).

The Stokes V time series was fit down to $\chi_r^2 = 1.6$ from an initial value of 5.5 assuming solid body rotation, since the differential rotation search was inconclusive. The maps and the Stokes V fits are shown in Fig. 4 and Fig. D.1, respectively. We caution that the amplitude of the spurious polarisation signature in Stokes N is of the same order of magnitude as the Zeeman signature in Stokes V for this star, but the shape and amplitude of Stokes V seems to be unaffected by the presence (or absence) of the Stokes N signature (see also Sect. 2). The properties of the field topology are listed in Table 2. The magnetic field is predominantly poloidal (89%), with 60%, 25%, and 10% of the magnetic energy stored in the dipolar, quadrupolar, and octupolar modes, respectively. The field is also non-axisymmetric (31% in the axisymmetric components) and exhibits an average field strength of 1.3 G.

## 5.7. HD 149026 (HIP 808308)

HD 149026 is a G0.0 dwarf (Sato et al. 2005) at a distance of 76.21 pc (Gaia Collaboration 2020). Its age was estimated to be 2.78 Gyr (Ramírez et al. 2012). This star is known to host a Saturn-like planet with a mass of 0.38 M$_{\text{Jup}}$ on a 2.87 d orbit (Sato et al. 2005; Bonomo et al. 2017). While there is no measurement of the rotation period for this star in the literature, we estimated a value of 12.3 ± 1.2 d using the radius and the projected equatorial velocity (see Table 1).
From the three Narval observations, we measured $log R'_{HK}$ values between -4.903 and -4.886, with a median of -4.900 ± 0.010. The value is consistent with Boro Saikia et al. (2018) and Brown et al. (2022). The $H\alpha$ index is stable at 0.308, and the Ca II IRT ranges 0.763-0.774 with a median value of 0.768 ± 0.004.

Overall, this star is the least active of our sample. We did not detect any evident Zeeman signatures from the circularly polarised spectra of this star, but only one observation led to a marginal detection with a longitudinal magnetic field of 3.4 ± 2.0 G. From observations with non-detections, we can still use Eq. (5) to compute a $B_l$ value with error bar to determine the upper limits. The two non-detection observations led to 0.35 ± 2.1 G and 2.2 ± 2.2 G; thus, we placed a $3\sigma$ upper limit of 6 G to the longitudinal field detection (see Table C.1).

## 5.8. HD 189733 (V452 Vul)

HD 189733 is a K2.0 dwarf (Gray et al. 2003) at a distance of 19.78 pc (Gaia Collaboration 2020), and it is part of a visual binary system (Mason et al. 2001). The age is 400 Myr (Ramírez et al. 2012), while the rotation period is 11.94 d (Fares et al. 2010). The star hosts a hot Jupiter on a 2.22 d orbit (Bouchy et al. 2005).

We measured $log R'_{HK}$ values between -4.485 and -4.449, with a median of -4.473 ± 0.011 (the median error bar is 0.026). Our range is compatible with the value computed by Brown et al. (2022) and 0.05 dex larger than that reported by Boro Saikia et al. (2018). The $H\alpha$ index varies between 0.362 and 0.364, with a median of 0.363 ± 0.001, and the Ca II IRT ranges 0.862-0.880 with a median value of 0.867 ± 0.005.

The longitudinal field was measured to be between -4.3 and 1.8 G, with a median value of -0.3 ± 2.0 G (the median error bar is 1.1 G). In comparison, Brown et al. (2022) measured an average $B_l$ of 4.5 G in absolute value.

We applied ZDI and fit the Stokes V time series down to $\chi_r^2 = 1.8$ from an initial value of 15.1, assuming solid body rotation. The maps are shown in Fig. 4 and the Stokes V models in Fig. D.1. The search for differential rotation was inconclusive, owing to both a sparse time series and the poor phase coverage (mostly between 0.0 and 0.5). Thus, we attempted to carry out ZDI reconstructions with values of differential rotation rates obtained from previous ZDI reconstructions: $d\Omega = 0.146 \pm 0.049 \text{ rad d}^{-1}$ (Fares et al. 2010) and $d\Omega = 0.110 \pm 0.050 \text{ rad d}^{-1}$ from Fares et al. (2017). In both cases, the value of the target $\chi_r^2$ was improved to ~1.7. Compared to the ZDI map obtained assuming solid body rotation, the maps including $d\Omega$ differ by less than 0.5 G in average field strength, less than 3% in poloidal energy fraction, and less than 0.5% for the axisymmetric energy fraction. The same was found for the dipolar and octupolar components, while the quadrupolar energy fraction decreased from 30 to 20%, implying a simpler magnetic field configuration when differential rotation was accounted for. We decided to fix $d\Omega = 0.110 \text{ rad d}^{-1}$ from Fares et al. (2017) since it was derived using the largest number of observations and for a data set in 2013, which is closest in time to our observations.

The properties of the magnetic topology are listed in Table 2. The field has an average field strength of 18 G, and it has 52% of the magnetic energy in the poloidal component, which is predominantly dipolar (68%). The quadrupolar and octupolar components account for 16% and 8% of the energy, while the axisymmetric component is 44%. Even though our map is reconstructed without an ideal phase coverage of the stellar rotation, our results are consistent with the previous reconstructions of

Moutou et al. (2007), Fares et al. (2010, 2017), with an equatorial magnetic spot of positive polarity in the radial component and an azimuthal component dominated by a negative polarity (see Fig. 3 in Fares et al. 2017). These authors studied the evolution of the magnetic topology of HD 189733 from a nine-year spectropolarimetric monitoring. In addition, although a magnetic cycle could not be constrained, they corroborated the fast evolution of the topology over a few stellar rotations. This was manifested mainly in the variations in axisymmetry and poloidal-to-toroidal energy fraction, but the toroidal component was always the dominant one.

## 5.9. HD 190406 (GJ 779)

HD 190406 is a G0.0 dwarf (Gray et al. 2006) at a distance of 17.77 pc (Gaia Collaboration 2020), and it is part of a visual binary system (Mason et al. 2001). The estimated age is 1.70 Gyr (Ramírez et al. 2012) and the rotation period is $P_{rot} = 13.94 \text{ d}$ (Wright et al. 2011).

We measured $log R'_{HK}$ values between -4.719 and -4.663, with a median of -4.697 ± 0.019, which is consistent with the values reported by Boro Saikia et al. (2018) and Brown et al. (2022). The $H\alpha$ index varies between 0.312 and 0.314, and the Ca II IRT ranges 0.789-0.800 with a median value of 0.795 ± 0.003. We did not detect any Stokes V signatures in the eight Narval observations of this star. The average longitudinal field is 0.9 G and the average uncertainty is 1.0 G; thus, we set a $3\sigma$ upper limit of 3 G (see Table C.1).

## 5.10. HD 206860 (HN Peg)

HD 206860 is a G0.0 dwarf (Gray et al. 2001) at a distance of 18.13 pc (Gaia Collaboration 2020), and it is part of a visual binary system (Mason et al. 2001). It is a member of the Her-Lyr moving group (Eisenbeiss et al. 2013), with an age of 200 Myr (Güdel 2007; Ramírez et al. 2012) and a rotation period of $P_{rot} = 4.55 \text{ d}$ (Boro Saikia et al. 2015). The star hosts a 16-M$_{Jup}$ planet at a distance of 773 au (Luhman et al. 2007).

We measured $log R'_{HK}$ values between -4.356 and -4.294, with a median of -4.334 ± 0.017 (the median error bar is similar to the standard deviation). These values are at most 0.07 dex larger than the values reported in Marsden et al. (2014), Boro Saikia et al. (2015, 2018), and Brown et al. (2022). The $H\alpha$ index varies between 0.326 and 0.331, with a median of 0.327 ± 0.001, which is 0.012 and 0.005 units larger than the values of Marsden et al. (2014) and Boro Saikia et al. (2015), respectively. The Ca II IRT ranges 0.895-0.929 with a median value of 0.910 ± 0.008, which is comparable to the range of values from both Marsden et al. (2014) and Boro Saikia et al. (2015).

The longitudinal field was measured to be between 2.5 and 10.5 G, with a median value of 5.2 ± 3.4 G (the median error bar is 1.5 G). In comparison, Boro Saikia et al. (2015) measured longitudinal field values between ±14 G, and Brown et al. (2022) obtained an average $B_l$ of 6.4 G in absolute value terms.

The ZDI reconstructions were performed for 2018 and 2019 epochs of HD 206860. The maps are shown in Fig. 4 and the Stokes V models in Fig. D.1. Assuming solid body rotation, the target $\chi_r^2$ is 2.0 and 1.6 starting from 11.7 and 13.0 for the 2018 and 2019 epochs, respectively. We then constrained differential rotation via $\chi^2$ minimisation (see Fig. 2). In 2018, we obtained $P_{rot} = 4.723 \pm 0.003 \text{ d}$ and $d\Omega = 0.109 \pm 0.004 \text{ rad d}^{-1}$ (implying a rotation period at the pole of 5.145±0.017 d), while in 2019 we obtained $P_{rot} = 4.650 \pm 0.017 \text{ d}$ and $d\Omega = 0.051 \pm 0.042 \text{ rad d}^{-1}$
(corresponding to $4.832 \pm 0.157$ d at the pole). Using these values, the target $\chi_r^2$ was improved to 1.3 and 1.1 for the two epochs. The values of $d\Omega$ differ by roughly a factor of two, but they are reasonably consistent within uncertainties. The error bar on the differential rotation rate extracted from the 2019 epoch is larger owing to the incomplete phase coverage of the observations used in the analysis.

We note that the differential rotation rates we constrained for 2018 and 2019 are a factor of 2 and 4 smaller, respectively, than those measured by Boro Saikia et al. (2015). If we were to perform the ZDI reconstruction using $d\Omega = 0.22 \text{ rad d}^{-1}$ from Boro Saikia et al. (2015), the $\chi_r^2$ improvement would not be as substantial (of about 0.1–0.2) as when using our derived values. This may suggest an evolution of the differential shear at the surface of the star, which was also recorded for other stars (AB Dor, LQ Hya, and HR 1099; Donati et al. 2003), and for the Sun in correlation with the magnetic cycle (Poljančić Beljan et al. 2022).

The properties of the field topology are listed in Table 2. In 2018, the magnetic field is predominantly poloidal (86%), with 47%, 17%, and 14% of the magnetic energy stored in the dipolar, quadrupolar, and octupolar modes, respectively, while the axisymmetric component accounts for 65% of the total magnetic energy. The situation is similar for 2019, the main difference being the decrease in the poloidal component to 69% and the corresponding increase in the toroidal component. The poloidal component also becomes more dipolar (60%). In terms of average magnetic field strength, the 2018 epoch featured a value of 21 G, while in 2019 it decreased to 16 G.

In general, the maps of HD 206860 are consistent with both the ZDI reconstructions done by Rosén et al. (2016), with Narval data between 2007–2011 and HARPS-Pol in 2013, and by Boro Saikia et al. (2015), with the Narval 2007–2011 time series. Indeed, they also feature a primarily poloidal, dipolar, and axisymmetric configuration, with the positive polarity predominantly located in the northern hemisphere. The maps also share similarities in terms of complexity, with magnetic patches extending down to low latitudes. More specifically, our reconstructions from 2018 and 2019 resemble the previous ZDI maps based on the 2009 data.

### 5.11. HD 219828 (HIP 115100)

HD 219828 is a G0.0 dwarf (Moore & Paddock 1950) at a distance of 72.56 pc (Gaia Collaboration 2020). The estimated age is 6.08 Gyr (Ramírez et al. 2012) and the rotation period is $P_{\text{rot}} = 28.70$ d (Santos et al. 2016). This is the most distant, oldest, and slowest rotator in our sample. The star hosts a hot Neptune (mass of 0.06 M$_{\text{Jup}}$) and a long-period planet (mass of 15 M$_{\text{Jup}}$) at a distance of 0.05 and 5.8 au, respectively (Santos et al. 2016).

We measured $\log R'_{\text{HK}}$ values between –5.000 and –4.966, with a median of –4.986 ± 0.011, which is 0.1 dex lower than what was reported by Boro Saikia et al. (2018). The $H\alpha$ index varies between 0.316 and 0.318, and the Ca II IRT ranges 0.747–0.783 with a median value of 0.758 ± 0.010. We did not clearly detect Stokes V signatures in the ten Narval observations of this star. For one of them, there is a marginal detection and its associated longitudinal magnetic field is 3.6 ± 2.0 G. The average longitudinal field is 1 G and the average uncertainty is 1.8 G, so we set a $3\sigma$ upper limit of 5 G.

## 6. Activity trends

As summarised in Table 1, the stars in our sample have ages between 0.2 and 6.1 Gyr, while the rotation period ranges

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

**Fig. 5.** Trends of activity indices and reconstructed magnetic field geometry with age. From the top: $\log R'_{\text{HK}}$, average magnetic field strength, poloidal fraction of the total energy, dipolar fraction of the poloidal energy, and axisymmetric fraction of the total energy. Properties measured for the same star but at different epochs are connected by a solid line. Our measurements are represented by ‘X’ and measurements by Folsom et al. (2016, 2018, and references therein) and Bellotti et al. (2025) are indicated with circles. The solar values from the work of Egeland (2017) and Vidotto et al. (2018) are also included.

between 4.6 and 28.7 d. Although the main purpose of this work is to provide magnetic field maps for modelling the stellar wind and environment in a subsequent paper, it is interesting to contextualise our measurements of activity indicators, longitudinal field, and magnetic field properties with previous studies. In particular, we discuss the trends with respect to stellar age and rotation period below.

Studies of activity indices around cool stars on the main sequence by Pace (2013), Marsden et al. (2014), and Brown et al. (2022) showed a decrease in stellar activity with increasing stellar age and rotation period. We noted a consistent behaviour (see Fig. 5 and Table C.1), where larger average values of activity indicators and longitudinal field towards younger and faster rotating stars. Quantitatively, we observe for example that our median longitudinal field decreases from 5 G to less than 1 G over from 0.1 to 4 Gyr (or from 5 d to 15 d rotation period) and our $\log R'_{\text{HK}}$ decreases from –4.2 to –5.0.

Turning to the large-scale magnetic field geometry, Vidotto et al. (2014b), Folsom et al. (2016, 2018), Rosén et al. (2016) reported trends across different ages and rotation periods for Sun-like stars similar to those examined in this work. They
reported a clear trend of decreasing average magnetic field with increasing age and rotation period, which we consistently observed for our stars as well. The average field strength of our six stars exhibits a decrease from 25 to 10 G between 0.2 and 0.4 Gyr. The oldest star in this sub-sample (at 1.6 Gyr) has an average magnetic field strength of 1 G. We show this trend in Fig. 5, where we complement our measurements with those of Petit et al. (2008), Folsom et al. (2016, 2018), and Bellotti et al. (2025).

Figure 5 also shows the variation with age and rotation period of the magnetic energy fraction in the poloidal, dipolar, and axisymmetric components. The addition of our measurements did not reveal clear trends in dipolar nor axisymmetric energy fraction, as already reported in Folsom et al. (2018). For the poloidal energy fraction, we expect the magnetic field of slow rotators and old stars to have a predominantly poloidal geometry, while faster-rotating and younger stars to have mixed geometries (Petit et al. 2008; Folsom et al. 2016). Our measurements are consistent with these expectations.

Overall, although the age estimate is not used in a ZDI reconstruction, inhomogeneous or inaccurate age measurements can increase scatter. Moreover, the search of these trends is complicated by temporal variations in the magnetic field in the form of intrinsic variability (short-term) or magnetic cycles (long-term), as revealed also for other Sun-like stars by spectropolarimetry (e.g. Boro Saikia et al. 2016; Mengel et al. 2016; Fares et al. 2017; Jeffers et al. 2022; Bellotti et al. 2025). For two of the stars examined here, HD 43162 and HD 206860, we reconstructed two maps corresponding to two distinct observational epochs. The magnetic properties reconstructed with ZDI are different at these epochs, which increases the scatter as illustrated in Fig. 5. A more thorough search for trends of the magnetic geometry over age and characterisation of this scatter would require multiple ZDI maps for the same star over different years to include such a temporal variability (Jeffers et al. 2023).

## 7. Conclusions

In this study, we analysed the spectropolarimetric time series of eleven Sun-like stars collected with Narval between 2018 and 2019. Our aim has been to characterise the magnetic properties of these stars for subsequent modelling of their stellar environments, which is topical in light of space-based missions dedicated to planetary atmospheric characterisation (e.g. *JWST* currently and *Ariel* in the future), as well as for direct imaging (e.g. *LIFE*; Quanz et al. 2018, 2022 and the Habitable World Observatory; National Academies of Sciences and Engineering, Medicine 2021). Indeed, knowledge of the stellar magnetic environment that exoplanets are embedded in feeds back to, for instance, the interpretation of atmospheric signatures due to star-planet interactions (e.g. Carolan et al. 2021; Gupta et al. 2023), as well as habitability assessments (e.g. Vidotto et al. 2013; Lingam & Loeb 2019; Van Looveren et al. 2024).

The stellar ages and rotation periods of our stars span different regions of the parameter space, from 0.2 to 6.1 Gyr and from 4.6 to 28.7 d, respectively, allowing us to sample stellar activity and magnetic field properties at different stages on the main sequence. As a first-order characterisation of the stars' magnetic activity, we computed chromospheric activity indices such as $log R'_{HK}$, $H\alpha$, and the Ca II infrared triplet from unpolarised spectra as well as the longitudinal magnetic field $B_l$ from circularly polarised spectra. Overall, we noted trends expected from previous work: in particular, stellar activity decreasing with

increasing age and rotation period (e.g. Marsden et al. 2014; Brown et al. 2022).

For the six stars with a detectable circular polarisation signature and large number of observations, we reconstructed the large-scale magnetic field configuration via ZDI. This was the first magnetic field reconstruction for two stars: HD 43162 and HD 114710. The properties of the large-scale magnetic field are in good agreement in terms of field strength and complexity with stars of similar properties. We found average magnetic field strengths ranging from 1 to 25 G, where the lowest values belong to old or slow rotating stars, consistently with what is expected for this stellar parameter space (see e.g. Rosén et al. 2015; Folsom et al. 2016, 2018). Most stars exhibit complex magnetic fields that are predominantly poloidal and with a significant (>10%) toroidal component and with distinct degrees of axisymmetry.

In agreement with See et al. (2015) and Folsom et al. (2018), we found that the poloidal-dominated topologies have various levels of axisymmetry, whereas the toroidal-dominated topologies are mostly axisymmetric. We also observed a decreasing trend of the average field strength with increasing age and rotation period. We did not record evident trends among the magnetic field properties, such as the poloidal fraction or axisymmetric fraction with respect to stellar age and rotation period, which could be due to temporal variations in the field topology, namely, in the magnetic cycles.

Our results stress the importance of spectropolarimetric studies to investigate the magnetic field properties of exoplanet-hosting stars, given the variety of properties of the large-scale magnetic field for Sun-like stars over time. This likely translates into a variety of stellar wind properties and offers a range of diverse magnetic environments in which exoplanets could be embedded. Ultimately, this study comprises a systematic attempt to characterise the evolution of environmental conditions and habitability of exoplanets from the perspective of stellar magnetism.

## Data availability

All analysed spectropolarimetric observations collected with Narval are available in Polarbase. They can be found at https://www.polarbase.ovgso.fr/

**Acknowledgements.** S.B. acknowledges funding by the Dutch Research Council (NWO) under the project "Exo-space weather and contemporaneous signatures of star-planet interactions" (with project number OCENW.M.22.215 of the research programme "Open Competition Domain Science- M"). S.B. also acknowledges funding from the SCI-S department of the European Space Agency (ESA), under the Science Faculty Research fund E/0429-03. S.B.S. acknowledges funding from the Austrian Science Fund (FWF) Lise Meitner project M2829-N. C.P.F. acknowledges funding from the European Union's Horizon Europe research and innovation program under grant agreement No. 101079231 (EXOHOST), and from the United Kingdom Research and Innovation Horizon Europe Guarantee Scheme (grant No. 10051045). This work used the Dutch national e-infrastructure with the support of the SURF Cooperative using grant nos. EINF-2218 and EINF-5173. Based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated 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. This work has made use of the VALD database, operated at Uppsala University, the Institute of Astronomy RAS in Moscow, and the University of Vienna; Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2013, 2018); NumPy (van der Walt et al. 2011); Matplotlib: Visualization with Python (Hunter 2007); SciPy (Virtanen et al. 2020) and PyAstronomy (Czesla et al. 2019).

## References

Airapetian, V. S., & Usmanov, A. V. 2016, ApJ, 817, L24
Airapetian, V. S., Barnes, R., Cohen, O., et al. 2020, Int. J. Astrobiol., 19, 136

Allan, A., & Vidotto, A. A. 2019, MNRAS, 490, 3760

Alvarado-Gómez, J. D., Hussain, G. A. J., Grunhut, J., et al. 2015, A&A, 582, A38

Alvarado-Gómez, J. D., Hussain, G. A. J., Drake, J. J., et al. 2018, MNRAS, 473, 4326

Alvarado-Gómez, J. D., Cohen, O., Drake, J. J., et al. 2022, ApJ, 928, 147

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

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

Bagnulo, S., Landolfi, M., Landstreet, J. D., et al. 2009, PASP, 121, 993

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

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

Bellotti, S., Fares, R., Vidotto, A. A., et al. 2023, A&A, 676, A139

Bellotti, S., Petit, P., Jeffers, S. V., et al. 2025, A&A, 693, A269

Blazère, A., Petit, P., Lignières, F., et al. 2016, A&A, 586, A97

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

Boro Saikia, S., Jeffers, S. V., Petit, P., et al. 2015, A&A, 573, A17

Boro Saikia, S., Jeffers, S. V., Morin, J., et al. 2016, A&A, 594, A29

Boro Saikia, S., Lüftinger, T., Jeffers, S. V., et al. 2018, A&A, 620, L11

Boro Saikia, S., Jin, M., Johnstone, C. P., et al. 2020, A&A, 635, A178

Bouchy, F., Udry, S., Mayor, M., et al. 2005, A&A, 444, L15

Brandenburg, A., Mathur, S., & Metcalfe, T. S. 2017, ApJ, 845, 79

Brown, E. L., Jeffers, S. V., Marsden, S. C., et al. 2022, MNRAS, 514, 4300

Busà, I., Aznar Cuadrado, R., Terranegra, L., Andretta, V., & Gomez, M. T. 2007, A&A, 466, 1089

Carolan, S., Vidotto, A. A., Villarreal D'Angelo, C., & Hazra, G. 2021, MNRAS, 500, 3382

Cecchi-Pestellini, C., Ciaravella, A., & Micela, G. 2006, A&A, 458, L13

Chini, R., Fuhrmann, K., Barr, A., et al. 2014, MNRAS, 437, 879

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

Cockell, C. S., Bush, T., Bryce, C., et al. 2016, Astrophysics, 16, 89

Czesla, S., Schröter, S., Schneider, C. P., et al. 2019, Astrophysics Source Code Library [record ascl:1906.010]

Donati, J. F. 2003, in Solar Polarization, eds. J. Trujillo-Bueno, & J. Sanchez Almeida, Astronomical Society of the Pacific Conference Series, 307, 41

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

Donati, J. F., & Collier Cameron, A. 1997, MNRAS, 291, 1

Donati, J. F., & Landstreet, J. D. 2009, Annu. Rev. Astron. Astrophys., 47, 333

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., Collier Cameron, A., & Petit, P. 2003, MNRAS, 345, 1187

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

Ducati, J. R. 2002, CDS/ADC Collection of Electronic Catalogues, 2237

Egeland, R. 2017, Ph.D. Thesis, Montana State University, Bozeman

Egeland, R., Metcalfe, T. S., Hall, J. C., & Henry, G. W. 2015, ApJ, 812, 12

Eisenbeiss, T., Ammler-von Eiff, M., Roell, T., et al. 2013, A&A, 556, A53

Evensberget, D., Carter, B. D., Marsden, S. C., Brookshaw, L., & Folsom, C. P. 2021, MNRAS, 506, 2309

Evensberget, D., Carter, B. D., Marsden, S. C., et al. 2022, MNRAS, 510, 5226

Evensberget, D., Marsden, S. C., Carter, B. D., et al. 2023, MNRAS, 398, 1383

Fares, R., Donati, J. F., Moutou, C., et al. 2009, MNRAS, 398, 1383

Fares, R., Donati, J. F., Moutou, C., et al. 2010, MNRAS, 406, 409

Fares, R., Bourrier, V., Vidotto, A. A., et al. 2017, MNRAS, 471, 1246

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

Folsom, C. P., Bouvier, J., Petit, P., et al. 2018, MNRAS, 474, 4956

Folsom, C. P., Ó Fionnagáin, D., Fossati, L., et al. 2020, A&A, 633, A48

Folsom, C. P., Erba, C., Petit, V., et al. 2025, arXiv e-prints [arXiv:2505.18476]

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

Gaidos, E. J. 1998, PASP, 110, 1259

Gizis, J. E., Reid, I. N., & Hawley, S. L. 2002, AJ, 123, 3356

Gray, D. F., & Baliunas, S. L. 1997, ApJ, 475, 303

Gray, R. O., Napier, M. G., & Winkler, L. I. 2001, AJ, 121, 2148

Gray, R. O., Corbally, C. J., Garrison, R. F., McFadden, M. T., & Robinson, P. E. 2003, AJ, 126, 2048

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

Güdel, M. 2007, Liv. Rev. Sol. Phys., 4, 3

Guinan, E. F., Ribas, I., & Harper, G. M. 2003, ApJ, 594, 561

Gupta, S., Basak, A., & Nandy, D. 2023, ApJ, 953, 70

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

Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90

Hussain, G. A. J., Collier Cameron, A., Jardine, M. M., et al. 2009, MNRAS, 398, 189

Jeffers, S. V., Boro Saikia, S., Barnes, J. R., et al. 2017, MNRAS, 471, L96

Jeffers, S. V., Cameron, R. H., Marsden, S. C., et al. 2022, A&A, 661, A152

Jeffers, S. V., Kiefer, R., & Metcalfe, T. S. 2023, Space Sci. Rev., 219, 54

Johnstone, C. P., Güdel, M., Lüftinger, T., Toth, G., & Brott, I. 2015a, A&A, 577, A27

Johnstone, C. P., Güdel, M., Brott, I., & Lüftinger, T. 2015b, A&A, 577, A28

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

Joner, M. D., Taylor, B. J., Laney, C. D., & van Wyk, F. 2006, AJ, 132, 111

Keenan, P. C., & McNeil, R. C. 1989, ApJS, 71, 245

Ketzet, L., & Poppenhaeger, K. 2023, MNRAS, 518, 1683

Kochukhov, O., Makaganiuk, V., & Piskunov, N. 2010, A&A, 524, A5

Kochukhov, O., Hackman, T., Lehtinen, J. J., & Wehrhahn, A. 2020, A&A, 635, A142

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

Kraft, R. P. 1965, ApJ, 142, 681

Kulikov, Y. N., Lammer, H., Lichtenegger, H. I. M., et al. 2007, Space Sci. Rev., 129, 207

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

Lammer, H., Güdel, M., Kulikov, Y., et al. 2012, Earth Planets Space, 64, 179

Landi Degl'Innocenti, E. 1992, in Magnetic field measurements, eds. F. Sanchez, M. Collados, & M. Vazquez, 71

Landstreet, J. D. 2009, in Observing and Modelling Stellar Magnetic Fields 1. Basic Physics and Simple Models, eds. C. Neiner, & J. P. Zahn, EAS Publications Series, 39, 1

Lehmann, L. T., & Donati, J. F. 2022, MNRAS, 514, 2333

Lehtinen, J., Jetsu, L., Hackman, T., Kajatkari, P., & Henry, G. W. 2016, A&A, 588, A38

Lehtinen, J. J., Spada, F., Käpylä, M. J., Olspert, N., & Käpylä, P. J. 2020, Nat. Astron., 4, 658

Lépine, S., & Bongiorno, B. 2007, AJ, 133, 889

Lingam, M., & Loeb, A. 2019, Int. J. Astrobiol., 18, 527

Linsky, J. L., Wood, B. E., Youngblood, A., et al. 2020, ApJ, 902, 3

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

López Ariste, A. 2002, ApJ, 564, 379

Luhman, K. L., Patten, B. M., Marengo, M., et al. 2007, ApJ, 654, 570

Marsden, S. C., Petit, P., Jeffers, S. V., et al. 2014, MNRAS, 444, 3517

Mason, B. D., Wycoff, G. L., Hartkopf, W. I., Douglass, G. G., & Worley, C. E. 2001, AJ, 122, 3466

Mathias, P., Aurière, M., López Ariste, A., et al. 2018, A&A, 615, A116

Mengel, M. W., Fares, R., Marsden, S. C., et al. 2016, MNRAS, 459, 4325

Messina, S., Guinan, E. F., Lanza, A. F., & Ambruster, C. 1999, A&A, 347, 249

Middelkoop, F. 1982, A&A, 107, 31

Montes, D., López-Santiago, J., Gálvez, M. C., et al. 2001, MNRAS, 328, 45

Moore, J. H., & Paddock, G. F. 1950, ApJ, 112, 48

Morgenthaler, A., Petit, P., Saar, S., et al. 2012, A&A, 540, A138

Moutou, C., Donati, J. F., Savalle, R., et al. 2007, A&A, 473, 651

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

National Academies of Sciences and Engineering, Medicine. 2021, Pathways to Discovery in Astronomy and Astrophysics for the 2020s (Washington, DC: The National Academies Press)

Nicholson, B. A., Vidotto, A. A., Mengel, M., et al. 2016, MNRAS, 459, 1907

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

Owen, J. E., & Jackson, A. P. 2012, MNRAS, 425, 2931

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

Paegert, M., Stassun, K. G., Collins, K. A., et al. 2021, arXiv e-prints [arXiv:2108.04778]

Penz, T., & Micela, G. 2008, A&A, 479, 579

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

Petit, P., Donati, J. F., Aurière, M., et al. 2005, MNRAS, 361, 837

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

Petit, P., Aurière, M., Konstantinova-Antova, R., et al. 2013, in The Environments of the Sun and the Stars, eds. J. P. Rozelot, & C. E. Neiner (Berlin Springer Verlag), Lecture Notes in Physics, 857, 231

Petit, P., Louge, T., Théado, S., et al. 2014, PASP, 126, 469

Pizzolato, N., Maggio, A., Micela, G., Sciortino, S., & Ventura, P. 2003, A&A, 397, 147

Poljančić Beljan, I., Jurdana-Šepić, R., Jurkić, T., et al. 2022, A&A, 663, A24

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

Quanz, S. P., Kammerer, J., Defrère, D., et al. 2018, in Optical and Infrared Interferometry and Imaging VI, eds. M. J. Creech-Eakman, P. G. Tuthill, & A. Mérand, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 10701, 1070111

Quanz, S. P., Ottiger, M., Fontanet, E., et al. 2022, A&A, 664, A21

Raeside, S. R., Rodgers-Lee, D., & Rimmer, P. B. 2025, A&A, 697, A26

Ramírez, I., Fish, J. R., Lambert, D. L., & Allende Prieto, C. 2012, ApJ, 756, 46

Rees, D. E., & Semel, M. D. 1979, A&A, 74, 1

Reiners, A. 2012, Liv. Rev. Sol. Phys., 9, 1

Ribas, I., Guinan, E. F., Güdel, M., & Audard, M. 2005, ApJ, 622, 680

Rodgers-Lee, D., Taylor, A. M., Vidotto, A. A., & Downes, T. P. 2021a, MNRAS, 504, 1519
Rodgers-Lee, D., Vidotto, A. A., & Mesquita, A. L. 2021b, MNRAS, 508, 4696

Rosén, L., Kochukhov, O., & Wade, G. A. 2015, ApJ, 805, 169

Rosén, L., Kochukhov, O., Hackman, T., & Lehtinen, J. 2016, A&A, 593, A35

Rugheimer, S., Kaltenegger, L., Segura, A., Linsky, J., & Mohanty, S. 2015, ApJ, 809, 57

Rutten, R. G. M. 1984, A&A, 130, 353

Ryabchikova, T., Piskunov, N., Kurucz, R. L., et al. 2015, Phys. Scr., 90, 054005

Santos, N. C., Santerne, A., Faria, J. P., et al. 2016, A&A, 592, A13

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

Sato, B., Fischer, D. A., Henry, G. W., et al. 2005, ApJ, 633, 465

See, V., Jardine, M., Vidotto, A. A., et al. 2014, A&A, 570, A99

See, V., Jardine, M., Vidotto, A. A., et al. 2015, MNRAS, 453, 4301

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

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

Tsai, S.-M., Lee, E. K. H., Powell, D., et al. 2023, Nature, 617, 483

Tu, L., Johnstone, C. P., Güdel, M., & Lammer, H. 2015, A&A, 577, L3

Valenti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141

van Belle, G. T., & von Braun, K. 2009, ApJ, 694, 1085

van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Comput. Sci. Eng., 13, 22

Van Looveren, G., Güdel, M., Boro Saikia, S., & Kislyakova, K. 2024, A&A, 683, A153

Van Looveren, G., Boro Saikia, S., Herbort, O., et al. 2025, A&A, 694, A310

Vaughan, A. H., Preston, G. W., & Wilson, O. C. 1978, PASP, 90, 267

Vidotto, A. A., Jardine, M., Morin, J., et al. 2013, A&A, 557, A67

Vidotto, A. A., Jardine, M., Morin, J., et al. 2014a, MNRAS, 438, 1162

Vidotto, A. A., Gregory, S. G., Jardine, M., et al. 2014b, MNRAS, 441, 2361

Vidotto, A. A., Lehmann, L. T., Jardine, M., & Pevtsov, A. A. 2018, MNRAS, 480, 477

Vidotto, A. A., Bourrier, V., Fares, R., et al. 2023, A&A, 678, A152

Virtanen, P., Gommers, R., Burovski, E., et al. 2020, https://doi.org/10.5281/zenodo.595738

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

Willamo, T., Lehtinen, J. J., Hackman, T., et al. 2022, A&A, 659, A71

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

Wood, B. E. 2004, Liv. Rev. Sol. Phys., 1, 2

Wood, B. E. 2006, Space Sci. Rev., 126, 3

Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 2011, ApJ, 743, 48

Zeeman, P. 1897, Nature, 55, 347

Zejda, M., Paunzen, E., Baumann, B., Mikulášek, Z., & Liška, J. 2012, A&A, 548, A97
# Appendix A: Journal of observations

<table><caption>**Table A.1.** List of median activity indices and standard deviations for our stars.</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2018</td></tr><tr><td>Aug 19</td><td>01:58:56</td><td>0.00</td><td>4x600.0</td><td>450</td><td>4.3</td><td>-4.406 ± 0.022</td><td>0.324 ± 0.001</td><td>0.867 ± 0.003</td><td>-4.6 ± 1.0</td></tr><tr><td>Aug 20</td><td>01:06:47</td><td>0.13</td><td>4x600.0</td><td>534</td><td>3.2</td><td>-4.413 ± 0.018</td><td>0.323 ± 0.001</td><td>0.860 ± 0.002</td><td>-3.7 ± 0.8</td></tr><tr><td>Aug 21</td><td>01:35:49</td><td>0.26</td><td>4x600.0</td><td>453</td><td>4.2</td><td>-4.410 ± 0.024</td><td>0.323 ± 0.001</td><td>0.867 ± 0.002</td><td>1.0 ± 1.0</td></tr><tr><td>Aug 22</td><td>01:27:16</td><td>0.39</td><td>4x600.0</td><td>597</td><td>2.9</td><td>-4.370 ± 0.015</td><td>0.324 ± 0.001</td><td>0.871 ± 0.002</td><td>0.8 ± 0.7</td></tr><tr><td>Aug 23</td><td>01:52:38</td><td>0.52</td><td>4x600.0</td><td>461</td><td>4.2</td><td>-4.350 ± 0.020</td><td>0.326 ± 0.001</td><td>0.883 ± 0.002</td><td>1.7 ± 0.9</td></tr><tr><td>Aug 24</td><td>00:26:20</td><td>0.64</td><td>4x600.0</td><td>439</td><td>4.6</td><td>-4.354 ± 0.024</td><td>0.328 ± 0.001</td><td>0.890 ± 0.003</td><td>1.0 ± 1.0</td></tr><tr><td>Aug 27</td><td>01:53:18</td><td>1.04</td><td>4x600.0</td><td>509</td><td>3.4</td><td>-4.413 ± 0.022</td><td>0.323 ± 0.001</td><td>0.864 ± 0.002</td><td>-2.1 ± 0.9</td></tr><tr><td>Aug 28</td><td>02:46:09</td><td>1.18</td><td>4x600.0</td><td>534</td><td>3.2</td><td>-4.414 ± 0.017</td><td>0.321 ± 0.001</td><td>0.855 ± 0.002</td><td>-3.2 ± 0.8</td></tr><tr><td>Sep 22</td><td>23:41:26</td><td>4.55</td><td>4x600.0</td><td>440</td><td>4.1</td><td>-4.352 ± 0.022</td><td>0.326 ± 0.001</td><td>0.876 ± 0.003</td><td>2.0 ± 1.0</td></tr><tr><td>Sep 24</td><td>23:16:43</td><td>4.81</td><td>4x600.0</td><td>367</td><td>5.2</td><td>-4.364 ± 0.027</td><td>0.327 ± 0.002</td><td>0.870 ± 0.003</td><td>-3.2 ± 1.2</td></tr><tr><td>Sep 25</td><td>24:04:11</td><td>4.94</td><td>4x600.0</td><td>473</td><td>3.5</td><td>-4.379 ± 0.019</td><td>0.325 ± 0.001</td><td>0.866 ± 0.002</td><td>-5.7 ± 0.9</td></tr><tr><td>Sep 27</td><td>00:30:02</td><td>5.07</td><td>4x600.0</td><td>418</td><td>3.9</td><td>-4.392 ± 0.024</td><td>0.327 ± 0.001</td><td>0.864 ± 0.003</td><td>-3.3 ± 1.1</td></tr><tr><td>Sep 28</td><td>00:25:54</td><td>5.20</td><td>4x600.0</td><td>605</td><td>2.9</td><td>-4.386 ± 0.013</td><td>0.327 ± 0.001</td><td>0.857 ± 0.002</td><td>0.3 ± 0.7</td></tr><tr><td>Oct 03</td><td>23:21:19</td><td>5.98</td><td>4x600.0</td><td>367</td><td>5.0</td><td>-4.383 ± 0.029</td><td>0.328 ± 0.002</td><td>0.863 ± 0.003</td><td>-4.2 ± 1.2</td></tr><tr><td>Oct 05</td><td>23:34:60</td><td>6.24</td><td>4x600.0</td><td>616</td><td>3.0</td><td>-4.401 ± 0.014</td><td>0.326 ± 0.001</td><td>0.853 ± 0.002</td><td>-2.3 ± 0.7</td></tr><tr><td>Oct 08</td><td>21:31:13</td><td>6.62</td><td>4x600.0</td><td>416</td><td>4.5</td><td>-4.338 ± 0.022</td><td>0.328 ± 0.001</td><td>0.885 ± 0.003</td><td>0.6 ± 1.1</td></tr></tbody></table>

**Notes.** Columns (1 and 2) Date and universal time of the observations; (3) Rotational cycle of the observations found using Eq. 1; (4) Exposure time of a polarimetric sequence, (5) S/N at 650 nm per polarimetric sequence; and (6) RMS noise level of Stokes V signal in units of an unpolarised continuum.

<table><caption>**Table A.2.** Same as Table A.1 for HD 28205</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2018</td></tr><tr><td>Nov 17</td><td>01:50:24</td><td>0.00</td><td>4x800.0</td><td>435</td><td>3.8</td><td>-4.443 ± 0.023</td><td>0.313 ± 0.001</td><td>0.847 ± 0.003</td><td>-0.1 ± 1.2</td></tr><tr><td>Nov 28</td><td>23:50:36</td><td>2.03</td><td>4x800.0</td><td>282</td><td>6.5</td><td>-4.438 ± 0.041</td><td>0.313 ± 0.002</td><td>0.848 ± 0.004</td><td>-1.0 ± 1.9</td></tr><tr><td>Dec 04</td><td>23:23:05</td><td>3.05</td><td>4x800.0</td><td>288</td><td>5.9</td><td>-4.437 ± 0.040</td><td>0.313 ± 0.002</td><td>0.849 ± 0.004</td><td>-2.9 ± 1.9</td></tr><tr><td>Dec 06</td><td>24:28:10</td><td>3.40</td><td>4x800.0</td><td>323</td><td>6.1</td><td>-4.453 ± 0.036</td><td>0.313 ± 0.002</td><td>0.843 ± 0.003</td><td>0.8 ± 1.7</td></tr><tr><td>Dec 10</td><td>22:51:31</td><td>4.07</td><td>4x800.0</td><td>406</td><td>4.4</td><td>-4.424 ± 0.023</td><td>0.314 ± 0.002</td><td>0.848 ± 0.003</td><td>-2.0 ± 1.3</td></tr></tbody></table>

<table><caption>**Table A.3.** Same as Table A.1 for HD 30495</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2018</td></tr><tr><td>Dec 05</td><td>00:46:09</td><td>0.00</td><td>4x400.0</td><td>445</td><td>3.7</td><td>-4.495 ± 0.030</td><td>0.323 ± 0.001</td><td>0.856 ± 0.002</td><td>0.4 ± 1.1</td></tr><tr><td>Dec 06</td><td>23:38:31</td><td>0.17</td><td>4x400.0</td><td>418</td><td>4.2</td><td>-4.490 ± 0.030</td><td>0.323 ± 0.001</td><td>0.854 ± 0.002</td><td>1.6 ± 1.2</td></tr><tr><td>Dec 10</td><td>24:10:40</td><td>0.53</td><td>4x400.0</td><td>654</td><td>2.6</td><td>-4.481 ± 0.020</td><td>0.323 ± 0.001</td><td>0.854 ± 0.002</td><td>3.0 ± 0.7</td></tr></tbody></table>

<table><caption>**Table A.4.** Same as Table A.1 for HD 43162</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2018</td></tr><tr><td>Dec 18</td><td>00:25:56</td><td>0.00</td><td>4x400.0</td><td>275</td><td>7.6</td><td>-4.348 ± 0.047</td><td>0.339 ± 0.002</td><td>0.902 ± 0.004</td><td>1.5 ± 1.8</td></tr><tr><td>Dec 18</td><td>01:51:37</td><td>0.01</td><td>4x400.0</td><td>315</td><td>7.2</td><td>-4.353 ± 0.039</td><td>0.338 ± 0.002</td><td>0.902 ± 0.003</td><td>-0.8 ± 1.5</td></tr><tr><td>Dec 19</td><td>23:06:57</td><td>0.23</td><td>4x400.0</td><td>136</td><td>17.4</td><td>-4.371 ± 0.187</td><td>0.336 ± 0.004</td><td>0.891 ± 0.009</td><td>4.4 ± 3.7</td></tr><tr><td colspan="10" style="text-align: center;">2019</td></tr><tr><td>Jan 03</td><td>23:02:49</td><td>1.99</td><td>4x400.0</td><td>377</td><td>6.0</td><td>-4.397 ± 0.059</td><td>0.338 ± 0.002</td><td>0.898 ± 0.003</td><td>1.6 ± 1.2</td></tr><tr><td>Jan 04</td><td>00:48:02</td><td>2.00</td><td>4x400.0</td><td>321</td><td>7.3</td><td>-4.398 ± 0.095</td><td>0.338 ± 0.002</td><td>0.898 ± 0.004</td><td>2.4 ± 1.5</td></tr><tr><td>Jan 04</td><td>22:32:52</td><td>2.11</td><td>4x400.0</td><td>351</td><td>6.3</td><td>-4.364 ± 0.046</td><td>0.340 ± 0.002</td><td>0.905 ± 0.004</td><td>4.4 ± 1.3</td></tr><tr><td>Jan 05</td><td>10:08:20</td><td>2.12</td><td>4x400.0</td><td>239</td><td>9.5</td><td>-4.368 ± 0.104</td><td>0.340 ± 0.003</td><td>0.904 ± 0.005</td><td>7.2 ± 2.0</td></tr><tr><td>Jan 05</td><td>22:22:25</td><td>2.23</td><td>4x400.0</td><td>329</td><td>6.7</td><td>-4.343 ± 0.034</td><td>0.341 ± 0.002</td><td>0.909 ± 0.003</td><td>4.3 ± 1.4</td></tr></tbody></table>
<table><caption>Table A.5. Same as Table A.1 for HD 82443</caption><tr><td>Jan 06</td><td>00:41:12</td><td>2.24</td><td>4x400.0</td><td>296</td><td>7.6</td><td>-4.332 ± 0.041</td><td>0.341 ± 0.002</td><td>0.911 ± 0.004</td><td>1.6 ± 1.6</td></tr><tr><td>Jan 06</td><td>22:44:04</td><td>2.34</td><td>4x400.0</td><td>320</td><td>7.1</td><td>-4.323 ± 0.034</td><td>0.342 ± 0.002</td><td>0.918 ± 0.003</td><td>5.8 ± 1.5</td></tr><tr><td>Jan 07</td><td>10:00:37</td><td>2.36</td><td>4x400.0</td><td>314</td><td>7.1</td><td>-4.315 ± 0.037</td><td>0.342 ± 0.002</td><td>0.919 ± 0.003</td><td>5.8 ± 1.5</td></tr><tr><td>Mar 20</td><td>19:24:42</td><td>10.92</td><td>4x600.0</td><td>342</td><td>6.3</td><td>-4.377 ± 0.044</td><td>0.339 ± 0.002</td><td>0.903 ± 0.003</td><td>2.4 ± 1.4</td></tr><tr><td>Mar 21</td><td>19:24:26</td><td>11.03</td><td>4x600.0</td><td>378</td><td>5.4</td><td>-4.379 ± 0.030</td><td>0.337 ± 0.002</td><td>0.896 ± 0.003</td><td>4.2 ± 1.2</td></tr><tr><td>Mar 22</td><td>19:29:24</td><td>11.15</td><td>4x600.0</td><td>462</td><td>4.5</td><td>-4.370 ± 0.025</td><td>0.338 ± 0.001</td><td>0.898 ± 0.002</td><td>2.4 ± 1.0</td></tr><tr><td>Mar 23</td><td>19:34:30</td><td>11.27</td><td>4x600.0</td><td>324</td><td>6.9</td><td>-4.391 ± 0.049</td><td>0.338 ± 0.002</td><td>0.896 ± 0.003</td><td>-2.7 ± 1.4</td></tr><tr><td>Mar 25</td><td>19:31:48</td><td>11.50</td><td>4x600.0</td><td>232</td><td>10.4</td><td>-4.347 ± 0.046</td><td>0.339 ± 0.002</td><td>0.908 ± 0.004</td><td>-4.9 ± 2.2</td></tr><tr><td>Mar 26</td><td>19:23:28</td><td>11.62</td><td>4x600.0</td><td>323</td><td>7.1</td><td>-4.342 ± 0.035</td><td>0.341 ± 0.002</td><td>0.911 ± 0.003</td><td>-4.3 ± 1.4</td></tr><tr><td>Mar 27</td><td>19:21:34</td><td>11.74</td><td>4x600.0</td><td>325</td><td>6.9</td><td>-4.353 ± 0.033</td><td>0.341 ± 0.002</td><td>0.908 ± 0.003</td><td>4.1 ± 1.5</td></tr><tr><td>Mar 28</td><td>19:43:50</td><td>11.86</td><td>4x600.0</td><td>483</td><td>4.4</td><td>-4.371 ± 0.026</td><td>0.339 ± 0.001</td><td>0.896 ± 0.002</td><td>4.3 ± 1.0</td></tr></table>

<table><caption>Table A.6. Same as Table A.1 for HD 114710</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2019</td></tr><tr><td>Mar 08</td><td>21:40:16</td><td>0.00</td><td>4x900.0</td><td>374</td><td>5.4</td><td>-4.170 ± 0.017</td><td>0.377 ± 0.002</td><td>0.987 ± 0.003</td><td>6.6 ± 2.6</td></tr><tr><td>Mar 11</td><td>21:59:23</td><td>0.56</td><td>4x900.0</td><td>450</td><td>3.6</td><td>-4.130 ± 0.011</td><td>0.381 ± 0.001</td><td>0.995 ± 0.002</td><td>-10.9 ± 1.3</td></tr><tr><td>Mar 20</td><td>23:10:32</td><td>2.24</td><td>4x900.0</td><td>425</td><td>4.6</td><td>-4.143 ± 0.012</td><td>0.380 ± 0.001</td><td>0.984 ± 0.002</td><td>15.6 ± 1.3</td></tr><tr><td>Mar 21</td><td>22:36:06</td><td>2.42</td><td>4x900.0</td><td>484</td><td>3.6</td><td>-4.133 ± 0.011</td><td>0.380 ± 0.001</td><td>0.986 ± 0.002</td><td>-6.7 ± 1.2</td></tr><tr><td>Mar 22</td><td>22:38:54</td><td>2.61</td><td>4x900.0</td><td>518</td><td>3.8</td><td>-4.159 ± 0.010</td><td>0.375 ± 0.001</td><td>0.974 ± 0.002</td><td>-5.5 ± 1.1</td></tr><tr><td>Mar 23</td><td>22:45:25</td><td>2.80</td><td>4x900.0</td><td>488</td><td>3.8</td><td>-4.161 ± 0.011</td><td>0.375 ± 0.001</td><td>0.982 ± 0.002</td><td>-4.1 ± 1.2</td></tr><tr><td>Mar 24</td><td>22:36:58</td><td>2.98</td><td>4x900.0</td><td>556</td><td>3.1</td><td>-4.147 ± 0.009</td><td>0.377 ± 0.001</td><td>0.983 ± 0.002</td><td>-7.1 ± 1.0</td></tr><tr><td>Mar 25</td><td>23:23:19</td><td>3.17</td><td>4x900.0</td><td>297</td><td>7.7</td><td>-4.147 ± 0.022</td><td>0.380 ± 0.002</td><td>0.987 ± 0.004</td><td>8.7 ± 2.2</td></tr><tr><td>Mar 26</td><td>22:47:25</td><td>3.36</td><td>4x900.0</td><td>382</td><td>5.4</td><td>-4.115 ± 0.015</td><td>0.383 ± 0.002</td><td>0.993 ± 0.003</td><td>3.1 ± 1.6</td></tr></tbody></table>

<table><caption>Table A.7. Same as Table A.1 for HD 149026</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2019</td></tr><tr><td>Mar 12</td><td>01:33:25</td><td>0.00</td><td>4x400.0</td><td>1132</td><td>1.4</td><td>-4.624 ± 0.010</td><td>0.313 ± 0.001</td><td>0.795 ± 0.001</td><td>0.0 ± 0.5</td></tr><tr><td>Mar 21</td><td>00:34:41</td><td>0.73</td><td>4x400.0</td><td>1113</td><td>1.5</td><td>-4.609 ± 0.009</td><td>0.314 ± 0.001</td><td>0.797 ± 0.001</td><td>0.6 ± 0.5</td></tr><tr><td>Mar 21</td><td>24:02:23</td><td>0.80</td><td>4x400.0</td><td>854</td><td>1.8</td><td>-4.608 ± 0.013</td><td>0.314 ± 0.001</td><td>0.797 ± 0.001</td><td>0.6 ± 0.7</td></tr><tr><td>Mar 22</td><td>24:04:39</td><td>0.89</td><td>4x400.0</td><td>1117</td><td>1.4</td><td>-4.610 ± 0.010</td><td>0.312 ± 0.001</td><td>0.793 ± 0.001</td><td>-0.6 ± 0.5</td></tr><tr><td>Mar 23</td><td>24:11:45</td><td>0.97</td><td>4x400.0</td><td>971</td><td>1.8</td><td>-4.614 ± 0.012</td><td>0.314 ± 0.001</td><td>0.799 ± 0.001</td><td>-0.3 ± 0.6</td></tr><tr><td>Mar 25</td><td>01:23:38</td><td>1.05</td><td>4x400.0</td><td>1220</td><td>1.4</td><td>-4.610 ± 0.009</td><td>0.313 ± 0.000</td><td>0.796 ± 0.001</td><td>-1.3 ± 0.5</td></tr><tr><td>Mar 26</td><td>00:47:24</td><td>1.13</td><td>4x400.0</td><td>859</td><td>1.9</td><td>-4.594 ± 0.013</td><td>0.315 ± 0.001</td><td>0.800 ± 0.001</td><td>-2.1 ± 0.7</td></tr><tr><td>Mar 27</td><td>01:22:59</td><td>1.21</td><td>4x400.0</td><td>784</td><td>2.3</td><td>-4.591 ± 0.015</td><td>0.314 ± 0.001</td><td>0.804 ± 0.001</td><td>0.2 ± 0.8</td></tr><tr><td>Mar 28</td><td>01:11:13</td><td>1.29</td><td>4x400.0</td><td>1196</td><td>1.3</td><td>-4.600 ± 0.009</td><td>0.313 ± 0.000</td><td>0.799 ± 0.001</td><td>0.7 ± 0.5</td></tr><tr><td>Mar 29</td><td>01:29:45</td><td>1.38</td><td>4x400.0</td><td>1050</td><td>1.6</td><td>-4.609 ± 0.010</td><td>0.313 ± 0.001</td><td>0.797 ± 0.001</td><td>-0.4 ± 0.5</td></tr><tr><td>Apr 19</td><td>22:56:18</td><td>3.15</td><td>4x400.0</td><td>887</td><td>1.9</td><td>-4.606 ± 0.013</td><td>0.313 ± 0.001</td><td>0.813 ± 0.001</td><td>-2.7 ± 0.7</td></tr><tr><td>Apr 27</td><td>23:15:53</td><td>3.80</td><td>4x400.0</td><td>1154</td><td>1.4</td><td>-4.633 ± 0.010</td><td>0.312 ± 0.001</td><td>0.810 ± 0.001</td><td>-0.8 ± 0.5</td></tr></tbody></table>

<table><caption>Table A.8. Same as Table A.1 for HD 189733</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2018</td></tr><tr><td>Aug 20</td><td>20:51:13</td><td>0.00</td><td>4x900.0</td><td>327</td><td>6.9</td><td>-4.903 ± 0.091</td><td>0.308 ± 0.002</td><td>0.768 ± 0.003</td><td>0.3 ± 2.1</td></tr><tr><td>Aug 21</td><td>20:39:25</td><td>0.05</td><td>4x900.0</td><td>342</td><td>6.7</td><td>-4.886 ± 0.080</td><td>0.308 ± 0.002</td><td>0.763 ± 0.003</td><td>3.4 ± 2.0</td></tr><tr><td>Aug 22</td><td>20:52:59</td><td>0.11</td><td>4x900.0</td><td>306</td><td>7.0</td><td>-4.902 ± 0.099</td><td>0.309 ± 0.002</td><td>0.774 ± 0.004</td><td>2.2 ± 2.3</td></tr></tbody></table>

<table><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2018</td></tr><tr><td>Nov 14</td><td>18:06:32</td><td>0.00</td><td>4x800.0</td><td>367</td><td>5.0</td><td>-4.471 ± 0.023</td><td>0.364 ± 0.002</td><td>0.867 ± 0.003</td><td>1.8 ± 1.1</td></tr><tr><td>Nov 15</td><td>18:06:59</td><td>0.08</td><td>4x800.0</td><td>409</td><td>4.3</td><td>-4.479 ± 0.021</td><td>0.362 ± 0.002</td><td>0.867 ± 0.002</td><td>1.5 ± 0.9</td></tr><tr><td>Nov 16</td><td>20:22:51</td><td>0.18</td><td>4x800.0</td><td>366</td><td>5.4</td><td>-4.485 ± 0.036</td><td>0.363 ± 0.002</td><td>0.866 ± 0.003</td><td>-1.7 ± 1.0</td></tr><tr><td>Nov 17</td><td>18:00:16</td><td>0.25</td><td>4x800.0</td><td>410</td><td>4.4</td><td>-4.484 ± 0.020</td><td>0.363 ± 0.001</td><td>0.862 ± 0.002</td><td>0.9 ± 0.9</td></tr></tbody></table>
<table><caption>Table A.9. Same as Table A.1 for HD 190406</caption><tr><td>Nov 28</td><td>18:00:01</td><td>1.17</td><td>4x800.0</td><td>206</td><td>9.5</td><td>-4.475 ± 0.059</td><td>0.363 ± 0.003</td><td>0.868 ± 0.005</td><td>-0.7 ± 1.9</td></tr><tr><td>Dec 04</td><td>18:11:16</td><td>1.68</td><td>4x800.0</td><td>310</td><td>6.7</td><td>-4.449 ± 0.029</td><td>0.364 ± 0.002</td><td>0.880 ± 0.003</td><td>-2.5 ± 1.3</td></tr><tr><td>Dec 10</td><td>17:55:26</td><td>2.18</td><td>4x800.0</td><td>401</td><td>5.0</td><td>-4.469 ± 0.019</td><td>0.363 ± 0.002</td><td>0.867 ± 0.003</td><td>0.0 ± 0.9</td></tr><tr><td>Dec 11</td><td>18:00:29</td><td>2.26</td><td>4x800.0</td><td>332</td><td>5.5</td><td>-4.472 ± 0.035</td><td>0.364 ± 0.002</td><td>0.870 ± 0.003</td><td>-4.3 ± 1.2</td></tr></table>

<table><caption>Table A.10. Same as Table A.1 for HD 206860</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2018</td></tr><tr><td>Jul 09</td><td>24:07:36</td><td>0.00</td><td>4x400.0</td><td>562</td><td>3.0</td><td>-4.704 ± 0.031</td><td>0.313 ± 0.001</td><td>0.795 ± 0.002</td><td>-1.2 ± 0.8</td></tr><tr><td>Jul 18</td><td>00:47:09</td><td>0.58</td><td>4x400.0</td><td>621</td><td>2.7</td><td>-4.714 ± 0.028</td><td>0.313 ± 0.001</td><td>0.789 ± 0.002</td><td>2.3 ± 0.8</td></tr><tr><td>Jul 19</td><td>01:14:59</td><td>0.65</td><td>4x400.0</td><td>536</td><td>3.4</td><td>-4.719 ± 0.036</td><td>0.312 ± 0.001</td><td>0.790 ± 0.002</td><td>0.9 ± 0.9</td></tr><tr><td>Jul 22</td><td>24:13:41</td><td>0.93</td><td>4x400.0</td><td>599</td><td>3.0</td><td>-4.698 ± 0.029</td><td>0.313 ± 0.001</td><td>0.795 ± 0.002</td><td>0.1 ± 0.8</td></tr><tr><td>Jul 23</td><td>23:04:34</td><td>1.00</td><td>4x400.0</td><td>423</td><td>4.6</td><td>-4.695 ± 0.044</td><td>0.314 ± 0.001</td><td>0.796 ± 0.002</td><td>2.1 ± 1.2</td></tr><tr><td>Jul 25</td><td>23:55:40</td><td>1.15</td><td>4x400.0</td><td>382</td><td>4.7</td><td>-4.694 ± 0.049</td><td>0.313 ± 0.002</td><td>0.794 ± 0.003</td><td>-0.1 ± 1.3</td></tr><tr><td>Jul 28</td><td>00:46:38</td><td>1.29</td><td>4x400.0</td><td>339</td><td>5.3</td><td>-4.666 ± 0.054</td><td>0.312 ± 0.002</td><td>0.800 ± 0.003</td><td>1.7 ± 1.5</td></tr><tr><td>Jul 29</td><td>00:38:55</td><td>1.36</td><td>4x400.0</td><td>516</td><td>3.4</td><td>-4.663 ± 0.032</td><td>0.313 ± 0.001</td><td>0.799 ± 0.002</td><td>1.2 ± 0.9</td></tr></tbody></table>

<table><caption>Table A.11. Same as Table A.1 for HD 219828</caption><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2018</td></tr><tr><td>Jul 18</td><td>03:22:53</td><td>0.00</td><td>4x300.0</td><td>487</td><td>3.4</td><td>-4.294 ± 0.015</td><td>0.330 ± 0.001</td><td>0.926 ± 0.003</td><td>8.7 ± 1.3</td></tr><tr><td>Jul 19</td><td>02:20:40</td><td>0.21</td><td>4x300.0</td><td>492</td><td>3.9</td><td>-4.320 ± 0.015</td><td>0.329 ± 0.001</td><td>0.919 ± 0.002</td><td>3.2 ± 1.3</td></tr><tr><td>Jul 23</td><td>03:30:05</td><td>1.10</td><td>4x300.0</td><td>449</td><td>3.6</td><td>-4.295 ± 0.016</td><td>0.331 ± 0.001</td><td>0.929 ± 0.003</td><td>7.6 ± 1.4</td></tr><tr><td>Jul 24</td><td>01:25:49</td><td>1.30</td><td>4x300.0</td><td>348</td><td>5.2</td><td>-4.332 ± 0.026</td><td>0.328 ± 0.002</td><td>0.915 ± 0.003</td><td>5.9 ± 1.9</td></tr><tr><td>Jul 24</td><td>03:14:51</td><td>1.32</td><td>4x500.0</td><td>457</td><td>3.8</td><td>-4.331 ± 0.018</td><td>0.328 ± 0.001</td><td>0.912 ± 0.003</td><td>7.6 ± 1.4</td></tr><tr><td>Jul 26</td><td>02:20:24</td><td>1.75</td><td>4x300.0</td><td>356</td><td>4.8</td><td>-4.338 ± 0.026</td><td>0.327 ± 0.002</td><td>0.908 ± 0.003</td><td>6.0 ± 1.9</td></tr><tr><td>Jul 26</td><td>02:56:22</td><td>1.75</td><td>4x500.0</td><td>442</td><td>4.0</td><td>-4.335 ± 0.019</td><td>0.327 ± 0.001</td><td>0.910 ± 0.003</td><td>4.0 ± 1.5</td></tr><tr><td>Jul 27</td><td>21:55:06</td><td>2.15</td><td>4x500.0</td><td>357</td><td>5.0</td><td>-4.307 ± 0.027</td><td>0.330 ± 0.002</td><td>0.920 ± 0.003</td><td>4.8 ± 1.9</td></tr><tr><td>Jul 29</td><td>01:53:50</td><td>2.40</td><td>4x300.0</td><td>437</td><td>4.2</td><td>-4.336 ± 0.018</td><td>0.327 ± 0.001</td><td>0.909 ± 0.003</td><td>10.4 ± 1.5</td></tr><tr><td>Jul 29</td><td>02:29:02</td><td>2.41</td><td>4x500.0</td><td>539</td><td>3.1</td><td>-4.339 ± 0.014</td><td>0.327 ± 0.001</td><td>0.910 ± 0.002</td><td>9.9 ± 1.2</td></tr><tr><td>Sep 11</td><td>00:12:31</td><td>12.06</td><td>4x300.0</td><td>442</td><td>3.8</td><td>-4.323 ± 0.017</td><td>0.327 ± 0.001</td><td>0.910 ± 0.003</td><td>2.6 ± 1.4</td></tr><tr><td>Sep 19</td><td>23:50:46</td><td>14.03</td><td>4x300.0</td><td>373</td><td>4.5</td><td>-4.338 ± 0.024</td><td>0.327 ± 0.002</td><td>0.912 ± 0.003</td><td>7.5 ± 1.7</td></tr><tr><td>Sep 20</td><td>21:28:42</td><td>14.23</td><td>4x300.0</td><td>373</td><td>4.4</td><td>-4.339 ± 0.022</td><td>0.326 ± 0.002</td><td>0.910 ± 0.003</td><td>4.2 ± 1.7</td></tr><tr><td>Sep 22</td><td>21:56:32</td><td>14.68</td><td>4x300.0</td><td>413</td><td>3.8</td><td>-4.309 ± 0.018</td><td>0.328 ± 0.001</td><td>0.921 ± 0.003</td><td>-2.1 ± 1.5</td></tr><tr><td>Sep 26</td><td>23:48:02</td><td>15.57</td><td>4x300.0</td><td>325</td><td>5.2</td><td>-4.310 ± 0.026</td><td>0.330 ± 0.002</td><td>0.914 ± 0.004</td><td>1.2 ± 2.0</td></tr><tr><td colspan="10" style="text-align: center;">2019</td></tr><tr><td>Jun 22</td><td>01:50:36</td><td>74.49</td><td>4x600.0</td><td>555</td><td>3.0</td><td>-4.348 ± 0.013</td><td>0.328 ± 0.001</td><td>0.908 ± 0.002</td><td>9.2 ± 1.1</td></tr><tr><td>Jul 12</td><td>01:43:13</td><td>78.89</td><td>4x600.0</td><td>667</td><td>2.4</td><td>-4.348 ± 0.010</td><td>0.327 ± 0.001</td><td>0.900 ± 0.002</td><td>-2.5 ± 0.9</td></tr><tr><td>Jul 16</td><td>02:41:07</td><td>79.77</td><td>4x600.0</td><td>681</td><td>2.5</td><td>-4.356 ± 0.010</td><td>0.327 ± 0.001</td><td>0.901 ± 0.002</td><td>4.9 ± 0.9</td></tr><tr><td>Jul 21</td><td>02:46:06</td><td>80.87</td><td>4x600.0</td><td>690</td><td>2.5</td><td>-4.349 ± 0.011</td><td>0.327 ± 0.001</td><td>0.896 ± 0.002</td><td>1.3 ± 0.9</td></tr><tr><td>Jul 22</td><td>00:54:03</td><td>81.08</td><td>4x600.0</td><td>554</td><td>3.0</td><td>-4.331 ± 0.013</td><td>0.328 ± 0.001</td><td>0.909 ± 0.002</td><td>4.9 ± 1.1</td></tr><tr><td>Jul 23</td><td>00:48:36</td><td>81.30</td><td>4x600.0</td><td>622</td><td>2.5</td><td>-4.326 ± 0.011</td><td>0.327 ± 0.001</td><td>0.909 ± 0.002</td><td>0.6 ± 0.9</td></tr><tr><td>Jul 25</td><td>02:04:19</td><td>81.75</td><td>4x600.0</td><td>619</td><td>2.7</td><td>-4.356 ± 0.011</td><td>0.326 ± 0.001</td><td>0.895 ± 0.002</td><td>6.0 ± 1.0</td></tr><tr><td>Jul 29</td><td>02:09:57</td><td>82.63</td><td>4x600.0</td><td>557</td><td>2.8</td><td>-4.338 ± 0.012</td><td>0.328 ± 0.001</td><td>0.906 ± 0.002</td><td>6.3 ± 1.0</td></tr><tr><td>Aug 02</td><td>02:20:10</td><td>83.51</td><td>4x600.0</td><td>516</td><td>3.6</td><td>-4.322 ± 0.014</td><td>0.329 ± 0.001</td><td>0.914 ± 0.002</td><td>5.5 ± 1.2</td></tr></tbody></table>

<table><thead><tr><th>Date</th><th>UT<br>[hh:mm:ss]</th><th>n<sub>cyc</sub></th><th>t<sub>exp</sub><br>[s]</th><th>S/N</th><th>σ<sub>LSD</sub><br>[10<sup>-5</sup>I<sub>c</sub>]</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>I</sub><br>[G]</th></tr></thead><tbody><tr><td colspan="10" style="text-align: center;">2018</td></tr><tr><td>Sep 19</td><td>23:04:02</td><td>0.00</td><td>4x800.0</td><td>257</td><td>9.1</td><td>-5.000 ± 0.147</td><td>0.317 ± 0.002</td><td>0.783 ± 0.004</td><td>0.8 ± 1.9</td></tr><tr><td>Sep 22</td><td>22:43:05</td><td>0.10</td><td>4x800.0</td><td>260</td><td>7.6</td><td>-4.994 ± 0.135</td><td>0.317 ± 0.002</td><td>0.765 ± 0.004</td><td>1.5 ± 1.8</td></tr><tr><td>Sep 24</td><td>22:20:31</td><td>0.17</td><td>4x800.0</td><td>184</td><td>13.8</td><td>-4.984 ± 0.257</td><td>0.316 ± 0.003</td><td>0.764 ± 0.006</td><td>1.2 ± 2.7</td></tr><tr><td>Sep 26</td><td>01:00:51</td><td>0.21</td><td>4x800.0</td><td>283</td><td>7.8</td><td>-4.998 ± 0.116</td><td>0.316 ± 0.002</td><td>0.761 ± 0.004</td><td>0.5 ± 1.6</td></tr><tr><td>Sep 28</td><td>01:22:59</td><td>0.28</td><td>4x800.0</td><td>344</td><td>6.6</td><td>-4.987 ± 0.087</td><td>0.318 ± 0.002</td><td>0.760 ± 0.003</td><td>-0.1 ± 1.3</td></tr><tr><td>Oct 03</td><td>24:17:46</td><td>0.49</td><td>4x800.0</td><td>242</td><td>9.9</td><td>-4.966 ± 0.153</td><td>0.317 ± 0.003</td><td>0.757 ± 0.004</td><td>3.6 ± 2.0</td></tr><tr><td>Oct 06</td><td>00:31:03</td><td>0.56</td><td>4x800.0</td><td>353</td><td>6.3</td><td>-4.985 ± 0.084</td><td>0.318 ± 0.002</td><td>0.752 ± 0.003</td><td>0.5 ± 1.3</td></tr><tr><td>Oct 08</td><td>22:34:57</td><td>0.66</td><td>4x800.0</td><td>257</td><td>9.1</td><td>-4.987 ± 0.136</td><td>0.317 ± 0.002</td><td>0.747 ± 0.004</td><td>1.0 ± 1.8</td></tr></tbody></table>
<table><tr><td>Oct 22</td><td>22:29:16</td><td>1.15</td><td>4x800.0</td><td>282</td><td>8.0</td><td>-4.986 ± 0.116</td><td>0.317 ± 0.002</td><td>0.755 ± 0.004</td><td>0.8 ± 1.7</td></tr><tr><td>Oct 24</td><td>22:07:20</td><td>1.22</td><td>4x800.0</td><td>268</td><td>9.3</td><td>-4.966 ± 0.137</td><td>0.317 ± 0.002</td><td>0.749 ± 0.004</td><td>0.7 ± 1.8</td></tr></table>

## Appendix B: Stokes LSD profiles of our stars

This appendix contains examples of LSD profiles for the stars in our sample (see Fig. B.1). The vertical dotted line in the plots indicates the radial velocity of the star, computed by fitting the Stokes *I* profile with a Voigt kernel and a linear component to account for residuals of continuum normalisation. The radial velocity is obtained as the centroid of the fit. The properties of the line lists adopted to perform LSD are listed in Table B.1. The number of lines reported in Table B.1 accounts for the removal of the wavelength intervals affected by telluric lines or containing the $H\alpha$ line, following Bellotti et al. (2022). These intervals are: [627,632], [655.5,657], [686,697], [716,734], [759,770], [813,835], and [895,986] nm.

The LIBRE-ESPRIT pipeline provides the null spectrum (Stokes *N*) as well, obtained by dividing polarimetric sub-exposures with the same polarisation state (see Donati et al. 1997, for more details). Least-squares deconvolution was used on the null spectrum to retrieve the corresponding average profile, which is a useful check for the presence of spurious polarisation signatures and the overall noise level in the LSD output (Donati et al. 1997; Bagnulo et al. 2009). As shown in Fig. B.1, except for HD 43162, HD 189733, and HD 206860, the Stokes *N* LSD profile for all other stars in our study exhibits a systematically positive signature at the radial velocity of the star. Surprisingly, the Stokes *N* signal is persistent even when subsets of the line lists are used, more precisely when magnetically sensitive ($g_{eff} > 1.2$) or red ($\lambda > 500$ nm) lines are used. The latter subset was indeed used by Folsom et al. (2016) and Bellotti et al. (2023) to mitigate the Stokes *N* signal attributed to an imperfect background subtraction during data reduction. Moreover, we note that the Stokes *N* signal in the data sets of a certain star may be present for the same night when observations of HD 206860 were collected and for which this type of signal was not reported. Following Mathias et al. (2018), we searched for (nearly) simultaneous ESPaDOnS observations of our stars to check whether the Stokes *V* signal was present for both instruments; whereas the Stokes *N* signal appeared only for Narval, possibly pointing at some instrumental effect. However, no ESPaDOnS observations are present for our stars or close enough in time to our Narval observations. These combined aspects make the origin of the signal more difficult to decipher.

Although the source of the Stokes *N* signal is not constrained, we noticed that the amplitude of Stokes *V* is not significantly affected by the presence of the Stokes *N* signal (see Fig. B.1). In practice, we observed that i) for some stars, the amplitude of Stokes *V* did not vary appreciably throughout the time series, while the Stokes *N* amplitude oscillated visibly in strength, ii) vice versa, for some other stars the Stokes *V* LSD profile varied while Stokes *N* did not, and iii) the Stokes *N* signal was observed even for stars whose Stokes *V* signature was not detected. Given the lack of correlation between the Stokes *V* and *N* signals, we proceeded to perform the magnetic characterisation of the stars.

<table><caption>**Table B.1.** Synthetic line lists associated to our stars and used with LSD.</caption><thead><tr><th>Name</th><th>Mask</th><th>N<sub>lines</sub></th><th>N<sub>obs</sub></th><th>⟨S/N⟩<sub>V</sub></th><th>v<sub>range</sub><br>[km s<sup>-1</sup>]</th></tr></thead><tbody><tr><td>HD 1835</td><td>6000_4.5</td><td>6960</td><td>16</td><td>479</td><td>±20</td></tr><tr><td>HD 28205</td><td>6000_4.5</td><td>6980</td><td>5</td><td>346</td><td>...</td></tr><tr><td>HD 30495</td><td>6000_4.5</td><td>7270</td><td>3</td><td>505</td><td>±20</td></tr><tr><td rowspan="2">HD 43162</td><td>5500_4.5</td><td>5340</td><td>11</td><td>298</td><td>±16</td></tr><tr><td>5500_4.5</td><td>5340</td><td>8</td><td>359</td><td>±16</td></tr><tr><td>HD 82443</td><td>5250_4.5</td><td>6180</td><td>9</td><td>441</td><td>±20</td></tr><tr><td>HD 114710</td><td>6000_4.5</td><td>7250</td><td>12</td><td>1028</td><td>±22</td></tr><tr><td>HD 149026</td><td>6000_4.0</td><td>4430</td><td>3</td><td>325</td><td>...</td></tr><tr><td>HD 189733</td><td>5000_4.5</td><td>6615</td><td>8</td><td>350</td><td>±17</td></tr><tr><td>HD 190406</td><td>6000_4.5</td><td>6960</td><td>8</td><td>497</td><td>...</td></tr><tr><td rowspan="2">HD 206860</td><td>6000_4.5</td><td>7260</td><td>17</td><td>396</td><td>±22</td></tr><tr><td>6000_4.5</td><td>7260</td><td>10</td><td>560</td><td>±22</td></tr><tr><td>HD 219828</td><td>6000_4.0</td><td>4430</td><td>10</td><td>273</td><td>...</td></tr></tbody></table>

**Notes.** The columns indicate: identifier of the star, properties of the line list in the form temperature[K]_logarithm surface gravity[cm s<sup>-1</sup>], number of lines used in LSD, number of observations, mean signal-to-noise ratio of the Stokes *V* LSD profile, velocity range centred on the radial velocity of the star encompassing both the absorption of the Stokes *I* LSD profiles and the lobes of Stokes *V* LSD profiles. When the Stokes *V* LSD profiles does not exhibit a clear Zeeman signature, the velocity range is not shown. For two stars, HD 43162 and HD 206860, we have two epochs so we list the values for each.

## Appendix C: Activity indices measurements

In this Appendix, we provide the list of measurements for the following activity indices: log $R'_{HK}$, $H\alpha$, Ca II IRT, and $B_I$. In Table C.1, we provide the average value for each index, together as the uncertainty computed as the standard deviation of the measurements.
<table><caption>**Table C.1.** List of median activity indices and standard deviations for our stars.</caption><thead><tr><th>Name</th><th>log R'<sub>HK</sub></th><th>Hα</th><th>Ca II IRT</th><th>B<sub>l</sub> [G]</th></tr></thead><tbody><tr><td>HD 1835</td><td>-4.384 ± 0.025</td><td>0.326 ± 0.002</td><td>0.866 ± 0.010</td><td>-2.2 ± 2.5</td></tr><tr><td>HD 28205</td><td>-4.438 ± 0.009</td><td>0.313 ± 0.001</td><td>0.848 ± 0.002</td><td>-1.0 ± 1.3</td></tr><tr><td>HD 30495</td><td>-4.490 ± 0.006</td><td>0.323 ± 0.001</td><td>0.855 ± 0.001</td><td>1.6 ± 1.1</td></tr><tr><td>HD 43162</td><td>-4.364 ± 0.023</td><td>0.339 ± 0.002</td><td>0.903 ± 0.007</td><td>2.4 ± 3.3</td></tr><tr><td>HD 82443</td><td>-4.147 ± 0.016</td><td>0.380 ± 0.003</td><td>0.986 ± 0.006</td><td>-4.1 ± 8.4</td></tr><tr><td>HD 114710</td><td>-4.609 ± 0.016</td><td>0.313 ± 0.001</td><td>0.798 ± 0.006</td><td>-0.3 ± 1.0</td></tr><tr><td>HD 149026</td><td>-4.900 ± 0.010</td><td>0.308 ± 0.002</td><td>0.768 ± 0.004</td><td>&lt; 6</td></tr><tr><td>HD 189733</td><td>-4.473 ± 0.011</td><td>0.363 ± 0.001</td><td>0.867 ± 0.005</td><td>-0.3 ± 2.0</td></tr><tr><td>HD 190406</td><td>-4.697 ± 0.019</td><td>0.317 ± 0.002</td><td>0.795 ± 0.003</td><td>&lt; 3</td></tr><tr><td>HD 206860</td><td>-4.334 ± 0.017</td><td>0.327 ± 0.001</td><td>0.910 ± 0.008</td><td>5.2 ± 3.4</td></tr><tr><td>HD 219828</td><td>-4.986 ± 0.011</td><td>0.317 ± 0.002</td><td>0.758 ± 0.010</td><td>&lt; 5</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Fig. B.1.** Least-squares deconvolution profiles for the young, Sun-like stars examined in this work. Each panel corresponds to a different star and contains the Stokes V (solid blue line) and Stokes N (dashed blue line) profiles. The vertical dotted line indicates the radial velocity of the star. We observe that some stars have a clearly detected Zeeman signature (e.g. HD 1835, HD 82443), while others present a flat Stokes V profile (HD 149026, HD 219828), hence they were excluded from tomographic inversion.</figcaption></figure>
## Appendix D: ZDI models of Stokes V profiles

This appendix contains the modelled time series of Stokes V LSD profiles via ZDI. For each star, we show the time series of observed Stokes V profiles, their ZDI models and the rotational cycle computed with Eq. 1. Each profile is shifted vertically for visualisation. In Table D.1, we also include the values of depth, Gaussian width, and Lorentzian width that describe the local Stokes I LSD line profiles. These parameters were optimised via $\chi^2$ minimisation between the observed and modelled Stokes I LSD profiles (see Sect. 4).

<table><caption>**Table D.1.** Optimised Stokes I parameters for ZDI reconstruction.</caption><thead><tr><th>Name</th><th>Depth</th><th>Gaussian width [km s<sup>-1</sup>]</th><th>Lorentzian width [km s<sup>-1</sup>]</th></tr></thead><tbody><tr><td>HD 1835</td><td>1.43</td><td>0.90</td><td>0.36</td></tr><tr><td>HD 43162</td><td>1.18</td><td>0.88</td><td>0.44</td></tr><tr><td>HD 82443</td><td>1.95</td><td>0.55</td><td>0.17</td></tr><tr><td>HD 114710</td><td>0.87</td><td>1.30</td><td>0.76</td></tr><tr><td>HD 189733</td><td>0.85</td><td>1.30</td><td>0.68</td></tr><tr><td>HD 206860</td><td>1.06</td><td>1.02</td><td>0.64</td></tr></tbody></table>

**Notes.** The depth values are relative to the local continuum.
<figure><img src="image_8.png" /><figcaption>**Fig. D.1.** Stokes V profiles for all stars examined in this work. The cycle number is computed using Eq. 1. For HD 43162 and HD 206860, two columns are shown and correspond to two different epochs.</figcaption></figure>