Magnetic Fields in M Dwarf Members of the Pleiades Open Cluster Using APOGEE spectra

Fábio Wanderley Observatório Nacional/MCTIC, R. Gen. José Cristino, 77, 20921-400, Rio de Janeiro, Brazil Katia Cunha Observatório Nacional/MCTIC, R. Gen. José Cristino, 77, 20921-400, Rio de Janeiro, Brazil Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721-0065, USA Oleg Kochukhov Department of Physics and Astronomy, Uppsala University, Box 516, S-75120 Uppsala, Sweden Verne V. Smith NSF’s NOIRLab, 950 N. Cherry Ave. Tucson, AZ 85719 USA Diogo Souto Departamento de Física, Universidade Federal de Sergipe, Av. Marcelo Deda Chagas, S/N Cep 49.107-230, São Cristóvão, SE, Brazil Lyra Cao Department of Astronomy, The Ohio State University, Columbus, OH 43210, USA Kevin Covey Department of Physics &\&& Astronomy, Western Washington University, Bellingham, WA, 98225, USA Steven R. Majewski Department of Astronomy, University of Virginia, Charlottesville, VA 22904-4325, USA Cintia Martinez Instituto de Astronomía y Física del Espacio (CONICET-UBA), C.C. 67 Sucursal 28, C1428EHA, Buenos Aires, Argentina Philip S. Muirhead Institute for Astrophysical Research and Department of Astronomy, Boston University, 725 Commonwealth Ave., Boston, MA 02215, USA Center for Interdisciplinary Exploration and Research in Astrophysics (CIERA) and Department of Physics and Astronomy, Northwestern University, 1800 Sherman Ave., Evanston, IL 60201, USA Marc Pinsonneault Department of Astronomy, The Ohio State University, Columbus, OH 43210, USA C. Allende Prieto Instituto de Astrofísica de Canarias, E-38205 La Laguna, Tenerife, Spain Departamento de Astrofísica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, Spain Keivan G. Stassun Department of Physics and Astronomy, Vanderbilt University, 6301 Stevenson Center Ln., Nashville, TN 37235, USA
Abstract

Average magnetic field measurements are presented for 62 M-dwarf members of the Pleiades open cluster, derived from Zeeman-enhanced Fe I lines in the H-band. An MCMC methodology was employed to model magnetic filling factors using SDSS-IV APOGEE high-resolution spectra, along with the radiative transfer code SYNMAST, MARCS stellar atmosphere models, and the APOGEE DR17 spectral line list. There is a positive correlation between mean magnetic fields and stellar rotation, with slow-rotator stars (Rossby number, Ro>>>0.13) exhibiting a steeper slope than rapid-rotators (Ro<<<0.13). However, the latter sample still shows a positive trend between Ro and magnetic fields, which is given by <<<B>>> = 1604 ×\times× Ro-0.20. The derived stellar radii, when compared with physical isochrones, show that on average, our sample shows radius inflation, with median enhanced radii ranging from +3.0%percent\%% to +7.0%percent\%%, depending on the model. There is a positive correlation between magnetic field strength and radius inflation, as well as with stellar spot coverage, correlations that together indicate that stellar spot-filling factors generated by strong magnetic fields might be the mechanism that drives radius inflation in these stars. We also compare our derived magnetic fields with chromospheric emission lines (Hα𝛼\alphaitalic_α, Hβ𝛽\betaitalic_β and Ca II K), as well as with X-ray and Hα𝛼\alphaitalic_α to bolometric luminosity ratios, and find that stars with higher chromospheric and coronal activity tend to be more magnetic.

Near Infrared astronomy(1093) — Open star clusters(1160) — M dwarf stars(982) — Stellar activity(1580) — Stellar magnetic fields(1610)
journal: ApJsoftware: Matplotlib (Hunter 2007), Numpy (Harris et al. 2020), Synmast (Kochukhov et al. 2010), Turbospectrum (Plez 2012)

1 Introduction

Quantitative characterization of magnetic fields provides a deeper understanding of stellar physics. As a star evolves, its stellar wind interacts with the magnetosphere, generating a torque that converts kinetic energy into magnetic energy, reducing the stellar angular momentum and resulting in a slowing of its rotational velocity (Kawaler, 1988). This process of magnetic braking over time enables gyrochronology to estimate the age of a star based on its stellar rotation (Skumanich, 1972; Barnes, 2003), with younger stars tending to have higher magnetic fields and activity than older stars of similar Teff. Another effect caused by magnetic fields is the heating of the stellar chromosphere and coronae, causing the emission of, respectively, strong UV and X-ray non-thermal radiation (Hawley et al., 2014; Astudillo-Defru et al., 2017; Newton et al., 2017).

Magnetic fields are especially important in M dwarf stars, as these stars have longer spin-down timescales, maintaining magnetic fields for longer periods than hotter stars (Newton et al., 2016). This means that any surrounding exoplanets will experience high-energy fluxes for longer periods, which can impact the habitability of these systems. In addition, since M dwarfs are cool, their habitable zones are located closer when compared to hotter stars, which increases the incident flux, as well as the probability of orbit locking, which, in turn, also impacts habitability. Nonetheless, M dwarf stars represent around 70%percent7070\%70 % of the stars of our galaxy (Salpeter, 1955; Reid & Gizis, 1997), and some of them can live trillions of years on the main sequence, which gives life many opportunities and time to form and evolve around these stars. A deep understanding of their magnetic fields and implications for their environment is fundamental for the understanding of these stars and the habitability of their exo-planetary systems.

One method to study stellar magnetic fields is through the Zeeman effect. Magnetically sensitive lines (e.g., having high effective Landé g-factors) when subject to magnetic fields are split into components, which we can observe as an additional line broadening. This broadening scales with the square of the line central wavelength, therefore spectral lines located at longer wavelengths are more sensitive to Zeeman splitting than lines with the same effective Landé-g factors located in the bluer part of the spectrum (see Kochukhov 2021 and references therein)

Due to other broadening mechanisms, such as Doppler broadening, stellar rotation, and instrumental broadening, we may not be able to resolve the Zeeman splitting in the spectrum, and what we measure is the Zeeman intensification of an affected line (Stift & Leone, 2003; Basri et al., 1992; Basri & Marcy, 1994).

There are two main ways to characterize stellar magnetic fields, one considering large-scale and the other small-scale magnetic fields. Large-scale analyses provide a topological view of the stellar magnetic field, separating it into components based on different orientations. The technique used for this characterization, Zeeman-Doppler Imaging (ZDI, Kochukhov 2016) is based on the analysis of the circular polarization of the line, described by the Stokes V𝑉Vitalic_V parameter, and therefore spectro-polarimetric data are needed for this technique. The small-scale magnetic approach models the total intensity of the field and is based on the Stokes I𝐼Iitalic_I parameter, and no polarimetric stellar data are needed. The review by Kochukhov (2021) discusses the current state of M dwarf magnetic field studies and presents a compilation of large and small-scale magnetic field measurements from the literature. Below we mention the results of a few of these studies.

The first work to model the magnetic field for an M dwarf star in the literature was Saar & Linsky (1985), who used Ti I lines from a high-resolution (Rsimilar-to\sim45,000) K-band spectrum of the flare star AD Leo and found a mean magnetic field of 3.8 kG. Many works followed, modeling the Zeeman effect in M dwarfs spectral lines in the optical and infrared and finding magnetic fields ranging from zero up to 8 kG (Johns-Krull & Valenti, 1996; Shulyak et al., 2011, 2014, 2017, 2019; Reiners et al., 2022; Cristofari et al., 2023a, b; Han et al., 2023). Reiners et al. (2022) (hereafter, R22) determined magnetic fields for a large sample of M dwarfs using CARMENES spectra, and analyzed how magnetic fields are related to parameters such as magnetic flux, activity, and Rossby numbers, and how these distributions change when comparing the saturated and non-saturated regimes. Some works studied magnetic fields deriving large-scale magnetic fields from Stokes V, as well as small-scale fields, and found that the latter is considerably greater than the one obtained from circular polarization, which indicates that most of the star’s magnetic field is probably stored in small structures at its surface (Phan-Bao et al., 2009; Kochukhov & Lavail, 2017; Kochukhov & Shulyak, 2019; Kochukhov & Reiners, 2020).

All of the studies mentioned above explored magnetic fields in M dwarfs from the Galactic field star population and these can in principle have different ages and metallicities. Stars from an open cluster, on the contrary, originate from the same molecular cloud and are expected to have approximately the same age and chemistry, making them great benchmarks with which to study stellar evolution and atomic diffusion, but these are also excellent benchmarks for studying stellar magnetic fields. Because cluster stars form at the same time and share the same chemical composition, metallicity and age dependencies are removed by their inter-comparison, and this allows for an investigation of magnetic fields primarily as a function of other stellar properties, such as effective temperatures or rotational periods (Souto et al., 2021). The recent work of Wanderley et al. (2023) used APOGEE near-infrared spectra (Majewski et al., 2017; Abdurro’uf et al., 2022) to derive atmospheric parameters and metallicities, and study radius inflation in a sample of M dwarf stars members of the young Hyades open cluster. (Wanderley et al., 2023) found that these stars are on average inflated, and this may be caused by stellar magnetic fields.

In this work, we use spectral lines affected by Zeeman broadening present in the SDSS APOGEE spectra to derive average magnetic fields for a sample of 62 M dwarfs from the young (age = 112 ±plus-or-minus\pm± 5 Myr; Dahm 2015), near-solar metallicity (Soderblom et al., 2009) Pleiades open cluster. This is the first study to derive magnetic fields for a sample of M dwarf members of an open cluster, and also the first to derive magnetic fields for M dwarfs based on APOGEE H-band spectra.

This paper is organized as follows: Section 2 presents the APOGEE data and sample selection. In Section 3 we present the methodology employed to derive average stellar magnetic fields for the Pleiades M dwarf star sample. In Section 4, we discuss the results which include the relation between magnetic fields and stellar rotation, comparisons with the literature, radius inflation, and analysis of activity indicators. Finally, Section 5 summarizes the conclusions.

2 APOGEE Data and Sample Selection

We determined average magnetic fields by analyzing near-infrared (λ𝜆\lambdaitalic_λ1.51μ𝜇\muitalic_μm to λ𝜆\lambdaitalic_λ1.69μ𝜇\muitalic_μm), high-resolution (average resolution of R22,500similar-to𝑅22500R\sim 22,500italic_R ∼ 22 , 500) spectra of Pleiades M dwarf stars observed by the SDSS IV APOGEE survey (Majewski et al., 2017; Blanton et al., 2017). As part of SDSS IV, the APOGEE spectra analyzed here were obtained at the 2.5-m telescope located at APO in the northern hemisphere (Bowen & Vaughan, 1973; Gunn et al., 2006; Wilson et al., 2019).

An initial sample of Pleiades members was obtained from Heyl et al. (2022) and confirmed using the membership analysis in Cantat-Gaudin et al. (2020). We adopted a threshold of 80%percent8080\%80 % of minimum membership probability from Cantat-Gaudin et al. (2020) for a star to be considered a member of the Pleiades. We cross matched this sample with APOGEE DR17 (Abdurro’uf et al., 2022), and selected for M dwarfs with 4.7<MKs<6.24.7subscript𝑀subscript𝐾𝑠6.24.7<M_{K_{s}}<6.24.7 < italic_M start_POSTSUBSCRIPT italic_K start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUBSCRIPT < 6.2 (Mann et al., 2015, 2016), using 2MASS KSsubscript𝐾𝑆K_{S}italic_K start_POSTSUBSCRIPT italic_S end_POSTSUBSCRIPT magnitudes (Skrutskie et al., 2006) and distances from Bailer-Jones et al. (2021). All magnitudes were corrected for extinction using the mean extinction to the Pleiades of AV = 0.12 (Stauffer et al., 2007), and the relations from Wang & Chen (2019). To remove binary stars from the sample, we considered only stars with a scatter in APOGEE radial velocity smaller than 1 km s-1. We also removed stars that presented large Gaia DR3 (Gaia Collaboration, 2022) RUWE numbers (RUWE >>> 1.4), as this can indicate the presence of an unresolved companion (Belokurov et al., 2020). In addition, stars without a vsini𝑖\sin{i}roman_sin italic_i measurement in DR17 were also removed from the sample.

We also analyzed the distribution of distances, proper motions (Gaia DR3, Gaia Collaboration 2022), and radial velocities (from APOGEE DR17) of the selected targets to check for outliers, but we found none. Finally, we also removed from the sample stars that had noisy and problematic APOGEE spectra. Our final sample of Pleiades members analyzed in this study is composed of 62 M dwarfs.

Figure 1 presents the Gaia and 2MASS CMDs (color-magnitude diagram) of the selected targets (top panels), their distribution in space, and proper motions (middle panels), as well as histograms for their star distances (d𝑑ditalic_d), and radial velocities (RV𝑅𝑉RVitalic_R italic_V) (bottom panels). The average distance, radial velocity, and proper motions along α𝛼\alphaitalic_α (right ascension) and δ𝛿\deltaitalic_δ (declination) for our Pleiades M dwarf sample are, respectively: <<<d>>>=135.12±plus-or-minus\pm±2.19 pc, <<<RV>>>=5.86±plus-or-minus\pm±2.19 km s-1, <<< μαcosδsubscript𝜇𝛼𝛿\mu_{\alpha}\cos{\delta}italic_μ start_POSTSUBSCRIPT italic_α end_POSTSUBSCRIPT roman_cos italic_δ >>>=18.20±plus-or-minus\pm±0.97 mas yr-1 and <<< μδsubscript𝜇𝛿\mu_{\delta}italic_μ start_POSTSUBSCRIPT italic_δ end_POSTSUBSCRIPT >>>=45.35±limit-from45.35plus-or-minus-45.35\pm- 45.35 ±1.10 mas yr-1. These results are in good agreement with measurements from Lodieu et al. (2019) of respectively: 135.15±plus-or-minus\pm±0.43 pc, 5.67 ±plus-or-minus\pm± 2.93 kms1𝑘𝑚superscript𝑠1km\,s^{-1}italic_k italic_m italic_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT, 19.5 mas yr-1 and 45.545.5-45.5- 45.5 mas yr-1.

In the top panels of Figure 1 we show several isochrones from the literature: a MIST isochrone (Choi et al., 2016), a DARTMOUTH isochrone (Dotter et al., 2008), a PARSEC isochrone (Bressan et al., 2012; Nguyen et al., 2022), a BHAC15 isochrone (Baraffe et al., 2015), and two SPOTS isochrones from Somers et al. (2020), one with spots covering 20%percent\%% of the stellar photosphere and another with a spot coverage of 80%percent\%%. All isochrones shown are for solar metallicity and an age of 100 Myr, which is roughly the estimated age for the Pleiades open cluster.

The Pleiades study by Covey et al. (2016) found significant scatter in the K vs J--Ks diagram of Pleiades stars (see Figure 3 in their paper). That study also compared the observed colors of Pleiades stars with physical models and found that the V--K colors in rapidly rotating stars present a positive offset for the same V magnitude if compared to slow rotators, which was interpreted as being due to binarity or a dependency of photospheric/spot properties on rotation rate. The photometric data in the Gaia G vs GBP{}_{\rm BP}-start_FLOATSUBSCRIPT roman_BP end_FLOATSUBSCRIPT -GRP CMD, shown in the top left panel of Figure 1, shows a clear offset, with most isochrones that do not consider stellar spots presenting bluer colors for the same magnitudes, while the SPOTS isochrone associated with an 80%percent\%% photospheric spot coverage presents an excellent match to the photometric data of the selected stars. We note that this pattern is not seen in the 2MASS CMD, where all isochrones present very small variations, even for different spot fractions, which might be related to the lower photometric spot contrast in the infrared when compared to the visible spectrum.

Refer to caption
Figure 1: From left to right, the top panels present respectively Gaia and 2MASS CMDs, the blue points are the selected M dwarf members of the Pleiades open cluster. Several 100 Myr solar metallicity isochrones are shown: MIST, DARTMOUTH, PARSEC, BHAC15, and SPOTS. Two SPOTS isochrones are shown, one for a photosphere spot coverage of 20%percent\%% and another for 80%percent\%%. The middle panels present respectively the distribution in space (right ascension and declination) and proper motions (from Gaia DR3) of the target stars. The bottom panels present respectively the distance (Bailer-Jones et al., 2021) and radial velocity histograms. The middle and bottom panels also present the mean and standard deviations for the parameters.

3 Methodology and Results

To derive average magnetic fields from Zeeman intensified lines we used the APOGEE line list (Smith et al., 2021) and model atmospheres from the MARCS grid (Gustafsson et al., 2008). The average magnetic field modeling in this study was done using the SYNMAST spectral synthesis code (Kochukhov et al., 2010), which computes the effects of magnetic fields on stellar spectra. This code uses polarised radiative transfer calculations to derive IQUV local Stokes parameters for a given magnetic field vector (in radial, meridional, and azimuthal orientations). In this study, we assumed a radial magnetic field and used SYNMAST to calculate intensity at seven limb angles. Then another code was used to perform disk integration, converting these intensity fluxes into density fluxes, which can be compared to APOGEE spectra.

To search for the best iron lines in the APOGEE region that can be used as magnetic field indicators for M dwarf stars, we compared two SYNMAST syntheses, one computed for 0 kG (no magnetic field) and another for 3 kG. We selected four Fe I lines as best indicators: λ𝜆\lambdaitalic_λ15207.526 Å, λ𝜆\lambdaitalic_λ15294.56 Å, λ𝜆\lambdaitalic_λ15621.654 Å, and λ𝜆\lambdaitalic_λ15631.948 Å. Table 3 presents the selected spectral lines, their central wavelengths in vacuum, the excitation potentials, loggfgf\log{\rm gf}roman_log roman_gf values from the APOGEE line list (Smith et al., 2021), effective Landé-g factors collected from VALD database (Piskunov et al., 1995; Kupka et al., 1999), as well as term designations associated with the upper and lower energy levels. The four selected Fe I lines are considerably sensitive to magnetic fields presenting effective Landé-g factors that range between 1.5similar-toabsent1.5\sim 1.5∼ 1.5 and 1.7similar-toabsent1.7\sim 1.7∼ 1.7. All the selected Fe I lines, except λ𝜆\lambdaitalic_λ15621.654 Å are from the same multiple. These lines were used to compute the magnetic fields for all stars in our sample and gave overall consistent results.

After the selection of diagnostic lines for measuring magnetic fields, the next step in our analysis was to generate a grid of synthetic spectra which was used in the analysis of each star. We adopted the DR17 ASPCAP (APOGEE Stellar Parameter and Chemical Abundances Pipeline, García Pérez et al. 2016; Abdurro’uf et al. 2022) Teff values for each star, along with an approximate logg𝑔\log{g}roman_log italic_g (similar-to\sim 4.7 – 4.8), depending on the stellar Teff. The ASPCAP results used in this work were computed with the Turbospectrum (Plez, 2012) code instead of Synspec (Hubeny & Lanz, 2011), however, we note that there are no significant differences between both sets of results for M dwarfs. The spectra of M dwarfs are relatively insensitive to the microturbulent velocity parameter, as noted by Souto et al. (2017, 2020), who found that a value of 1 km s-1 provides good fits to the observations, and we adopt this value in our analysis. All measurable Fe I lines in the APOGEE spectra of M dwarfs have high effective Landé-g factors, which makes them not suitable to be used as rotational broadening indicators. Therefore, we used a sample of OH lines, which are insensitive to magnetic fields. This approach is similar to the one adopted by Johns-Krull et al. (2004); Johns-Krull (2007); Yang et al. (2008); Yang & Johns-Krull (2011); Lavail et al. (2019), who derived mean magnetic fields for T-Tauri stars from K-band near-infrared spectra, using magnetically insensitive CO lines to measure non-magnetic broadening. To derive projected rotational velocities (vsini𝑖\sin{i}roman_sin italic_i) for the stars we used the radiative transfer code Turbospectrum (Plez, 2012); we adopted vsini𝑖{\sin{i}}roman_sin italic_i threshold of 3 km s-1 given the spectral resolution of the APOGEE spectra. We computed a grid of synthetic spectra, for the adopted stellar parameters for each star and Fe I line, with metallicities ranging from -0.75\leq[Fe/H]\leq+0.5 in steps of 0.25, and magnetic field values from 0 to 12 kG in steps of 2 kG, considering only the radial component. We then convolved the synthetic spectra with a rotational profile for the adopted vsini𝑖\sin{i}roman_sin italic_i as well as a Gaussian profile corresponding to the spectrum LSF (see Wilson et al. 2019; Nidever et al. 2015). Each synthesis was fitted to the DR17 normalized APOGEE spectrum (García Pérez et al., 2016) and was subject to small wavelength shifts when needed.

We employed Monte Carlo and Markov Chain (MCMC) to model the observed spectra and derive magnetic fields for the Pleiades M dwarfs. MCMC is a powerful tool, not only because it provides best fits to observations but also because it gives realistic and well-defined uncertainties based on the posterior distribution. Our methodology considers that the surface of the star can be divided into different components, each associated with a different magnetic field value. The filling factor describes the fraction of the stellar surface associated with a specific <<<B>>>. Many works in the literature used MCMC to derive average magnetic fields from filling factor determinations (Lavail et al., 2019; Kochukhov & Reiners, 2020; Hahlin et al., 2021; Hahlin & Kochukhov, 2022; Reiners et al., 2022; Cristofari et al., 2023a, b; Hahlin et al., 2023; Pouilly et al., 2023).

We developed a methodology to derive magnetic fields, that employs the python-code emcee (Foreman-Mackey et al., 2013) and finds the combination of metallicity, and 6 filling factors (2–12 kG in a 2 kG step) that best fits the four selected Fe I lines at the same time. We note that the non-magnetic filling factor f0 is given by 1fnsubscript𝑓𝑛-\sum f_{n}- ∑ italic_f start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT. For each entry in the posterior distribution, we calculate an average magnetic field in Gauss units given by:

<B>=nfn×1000n,n=[2,4,6,8,10,12]formulae-sequenceexpectation𝐵subscript𝑛subscript𝑓𝑛1000𝑛𝑛24681012<B>=\sum_{n}f_{n}\times 1000n\,\,,\,\,n=[2,4,6,8,10,12]< italic_B > = ∑ start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT italic_f start_POSTSUBSCRIPT italic_n end_POSTSUBSCRIPT × 1000 italic_n , italic_n = [ 2 , 4 , 6 , 8 , 10 , 12 ] (1)

The adopted average magnetic field, along with the lower and upper uncertainties, are given respectively by the median, the 16thsuperscript16𝑡16^{th}16 start_POSTSUPERSCRIPT italic_t italic_h end_POSTSUPERSCRIPT, and 84thsuperscript84𝑡84^{th}84 start_POSTSUPERSCRIPT italic_t italic_h end_POSTSUPERSCRIPT percentiles of the posterior distribution. In Table 1 we provide the filling factors and mean magnetic fields for two stars, the ones having the highest and lowest <<<B>>> in our sample.

In Figure 2, we illustrate the methodology by presenting the best fits, as well as the corner plot for the star 2M03511207+2355575. The left panels show fits for the four lines, where black dots are the observed spectrum, and the blue and red lines represent respectively the synthesis with the derived magnetic field, and the result with the same metallicity but no magnetic field. The right panel shows the corner plot describing the results. It presents the median and uncertainties (from 16thsuperscript16𝑡16^{th}16 start_POSTSUPERSCRIPT italic_t italic_h end_POSTSUPERSCRIPT and 84thsuperscript84𝑡84^{th}84 start_POSTSUPERSCRIPT italic_t italic_h end_POSTSUPERSCRIPT percentiles) of the modeled filling factors and the metallicity. We also show the posterior distribution for the derived average magnetic field.

Refer to caption
Refer to caption
Figure 2: Mean magnetic field results for star 2M03511207+2355575. The left panels show the four Fe I lines used in the modeling: black dots are the observed APOGEE spectrum, red dashed lines are synthetic profiles computed without magnetic field, and dark blue lines are our best-fits obtained from the MCMC modeling. The right panel is a corner plot that presents the median and uncertainties (from 16thsuperscript16𝑡16^{th}16 start_POSTSUPERSCRIPT italic_t italic_h end_POSTSUPERSCRIPT and 84thsuperscript84𝑡84^{th}84 start_POSTSUPERSCRIPT italic_t italic_h end_POSTSUPERSCRIPT percentiles) of the derived parameters, which includes filling factors from 2 to 12 kG in steps of 2 kG. We also show the obtained average magnetic field (in kG) and its uncertainties.

The derived mean magnetic fields for the studied Pleiades M dwarf stars range from similar-to\sim1.0 to similar-to\sim4.2 kG, with a median±plus-or-minus\pm±MAD of 3.0±plus-or-minus\pm±0.3 kG. Table 2 presents the derived mean magnetic fields for the sample stars, along with the adopted Teff and vsini𝑖\sin{i}roman_sin italic_i, and the SNR of the APOGEE spectra analyzed. This table also contains some quantities that will be discussed in Section 4: rotational periods, and Rossby numbers the ratio between the derived magnetic fields and the magnetic field limit based on kinetic-to-magnetic energy conversion, magnetic fluxes, and activity indicators: X-ray to bolometric ratios and Hα𝛼\alphaitalic_α, Hβ𝛽\betaitalic_β and Ca II K equivalent widths.

4 Discussion

4.1 Magnetic Fields, Rossby Numbers and Stellar Rotation

We collected rotational periods for 53 stars in our sample from the works by Hartman et al. (2010), Rebull et al. (2016) and Covey et al. (2016), who analyzed data respectively from the HATNet (Bakos et al., 2004), K2 (Howell et al., 2014), and PTF (Law et al., 2009; Rau et al., 2009) surveys. The rotational periods of the sample stars include both slow and rapid rotators, encompassing rotational periods similar-to\sim0.7 to 17.3 days.

The Rossby number (Ro) is an important indicator of stellar activity, it is given by the ratio between the rotational period and the convective turnover time. The convective turnover time, τ𝜏\tauitalic_τ, is defined as the time that it takes a convective cell to traverse the convective envelope of a star. Due to both their deeper convective envelopes and slower convective velocities, cooler stars are expected to have longer convective turnover times, and also to be able to sustain dynamos longer than hotter stars. For τ𝜏\tauitalic_τ measured in units of days, we used the relation τ=12.3×(L/L)0.5𝜏12.3superscript𝐿subscript𝐿direct-product0.5\tau=12.3\times(L/L_{\odot})^{-0.5}italic_τ = 12.3 × ( italic_L / italic_L start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT - 0.5 end_POSTSUPERSCRIPT, derived by R22. Bolometric luminosities were determined using distances from Bailer-Jones et al. (2021), 2MASS Kssubscript𝐾𝑠K_{s}italic_K start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT magnitudes, and bolometric corrections from photometric calibrations from Mann et al. (2015, 2016). This calibration is based on 2MASS J𝐽Jitalic_J and V𝑉Vitalic_V band magnitudes collected from Stauffer et al. (2007); Lasker et al. (2008); Zacharias et al. (2012); Muirhead et al. (2018). We also employed the same method to account for extinction described in Section 2.

The mean magnetic fields derived for the sample stars versus their projected rotational velocities, vsini𝑖\sin{i}roman_sin italic_i (both of which are quantities derived from the APOGEE spectra), are shown in Figure 3 as blue symbols. The blue circles represent stars for which we can estimate the vsini𝑖\sin{i}roman_sin italic_i, while the blue triangles are stars having vsini𝑖\sin{i}roman_sin italic_i values up to 3 km s-1, the adopted threshold in vsini𝑖\sin{i}roman_sin italic_i that can be estimated in this study. We also show in this figure literature results from R22 as black x’s. The mean magnetic fields for the Pleiades stars generally overlap with the M dwarfs in R22 within the overlapping vsini𝑖\sin{i}roman_sin italic_i range, showing just a small systematic difference in <<<B>>> when compared to R22.

Refer to caption
Figure 3: Distribution of the derived average magnetic fields versus adopted projected rotational velocities, vsini𝑖\sin{i}roman_sin italic_i, for the studied M dwarfs in the Pleiades are represented by blue circles and triangles for respectively stars with vsini>3𝑖3\sin{i}>3roman_sin italic_i > 3 km s-1 and with vsini3𝑖3\sin{i}\leq 3roman_sin italic_i ≤ 3 km s-1. Results from Reiners et al. (2022) are also shown (black x’s).

Figure 4 shows our results of magnetic fields versus rotational periods for the Pleiades M dwarfs (red and blue open circles), along with results from the literature for comparison. We compiled a total of 281 average magnetic field measurements from the literature. The results from Shulyak et al. (2017) and Shulyak et al. (2019) are represented by green squares in Figure 4. Results from Cristofari et al. (2023a) and Cristofari et al. (2023b), (results from radial-tangential macroturbulent velocity profiles and fitted logg𝑔\log{g}roman_log italic_g) are represented by cyan diamonds. Results from R22 are for stars with Teff<4000subscript𝑇eff4000T_{\rm eff}<4000italic_T start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT < 4000 K and are represented by black x’s. Finally, other results from the literature (orange triangles), are from the following works: Shulyak et al. (2011); Kochukhov & Lavail (2017); Kochukhov et al. (2001, 2009); Kochukhov & Shulyak (2019); Johns-Krull & Valenti (2000); Afram et al. (2009); Shulyak et al. (2014); Saar & Linsky (1985); Saar (1994); Phan-Bao et al. (2009); Saar (1996); Johns-Krull & Valenti (1996).

We can divide the <<<B>>>–Prot plane in Figure 4 into two regions that correspond to the saturated and unsaturated regimes. Red and blue open circles in this figure represent respectively stars in our sample that are slow rotators (Ro>>>0.13, non-saturated regime) and rapid rotators (Ro<<<0.13, saturated regime). We added a vertical dashed line at Protsubscript𝑃rotP_{\rm rot}italic_P start_POSTSUBSCRIPT roman_rot end_POSTSUBSCRIPT=7 days to Figure 4, which is an estimated threshold analogous to the separation based on Rossby numbers. Stars in the non-saturated regime are expected to have lower activity levels and also to present a greater dependency between Rossby numbers and magnetic fields. Stars in the saturated regime are in the limit of kinetic to magnetic energy conversion, and show a much flatter relation, although there is still some dependency between Rossby numbers and magnetic fields.

Most of the stars in our sample have rotational periods roughly between 1 and 10 days. Overall, stars in the saturated regime (blue open circles) follow an approximately constant <<<B>>> as a function of Protsubscript𝑃rotP_{\rm rot}italic_P start_POSTSUBSCRIPT roman_rot end_POSTSUBSCRIPT, showing just a modest inclination in the trend. It is clear that our results for this population overlap with results from the literature in this regime. For sample stars in the unsaturated regime (open red circles), <<<B>>> values begin to show a de-saturation signal (at Prot similar-to\sim7 days), with the slowest rotating stars of our sample presenting a steeper negative relation between rotational periods and magnetic fields than for faster rotators. Although our sample does not include very slowly rotating stars, our results generally overlap with literature values and are in line with the steep relation between magnetic fields and rotational periods in the literature, which extends to Prot similar-to\sim 100 days. It is interesting to note that despite the fact that our sample and that from R22 both have the same upper Teffsubscript𝑇effT_{\rm eff}italic_T start_POSTSUBSCRIPT roman_eff end_POSTSUBSCRIPT limit of 4000 K, the latter sample reaches much greater rotational periods and lower magnetic fields. The possible main reason for this is that R22 results are for field M dwarfs, while our work considers only M dwarf star members of the very young Pleiades open cluster. M dwarfs from the field may have had much more time to lose their magnetic fields and angular momentum when compared to the young stars from our sample.

Refer to caption
Figure 4: Mean magnetic fields versus rotational periods for Pleiades M dwarfs rapid (open blue circles, Ro<<<0.13) and slow rotators (open red circles, Ro>>>0.13), along with data compiled from the literature. Green squares, black x’s, and cyan diamonds are data from respectively Shulyak et al. (2017, 2019), Reiners et al. (2022), and Cristofari et al. (2023a, b). Orange triangles show results from other works, see text for references.

The dependency between rotation and magnetic fields can be explained by the dynamo mechanism that converts kinetic into magnetic energy. The total magnetic field of the star is limited by its available kinetic energy, and this maximum magnetic field is described as Bkin. Here we adopt the relation from Reiners et al. (2009) for Bkin (in units of Gauss):

Bkin=4800×(ML2R7)1/6subscript𝐵kin4800superscript𝑀superscript𝐿2superscript𝑅716B_{\rm kin}=4800\times\bigg{(}\frac{ML^{2}}{R^{7}}\bigg{)}^{1/6}italic_B start_POSTSUBSCRIPT roman_kin end_POSTSUBSCRIPT = 4800 × ( divide start_ARG italic_M italic_L start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_R start_POSTSUPERSCRIPT 7 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT 1 / 6 end_POSTSUPERSCRIPT (2)

where L𝐿Litalic_L is the derived luminosity (in L units), M𝑀Mitalic_M is the stellar mass (in M units), and R𝑅Ritalic_R is the radius (in R units) derived using the Stefan-Boltzmann equation, with the adopted Teff. We derived magnetic fluxes (ϕBsubscriptitalic-ϕB\phi_{\rm B}italic_ϕ start_POSTSUBSCRIPT roman_B end_POSTSUBSCRIPT) by multiplying the obtained average magnetic field by 4πR24𝜋superscript𝑅24\pi R^{2}4 italic_π italic_R start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Masses were determined from isochrones, following the methodology discussed in Wanderley et al. (2023). In summary, we selected SPOTS isochrones from Somers et al. (2020) for 100 Myr and solar metallicity, interpolated isochrones associated with different photospheric spots fractions (fspots), and adopted as the mass of the star the point in the interpolated isochrone plane with closest Teff and luminosity of the star.

We also estimated masses using DARTMOUTH (Dotter et al., 2008), MIST (Choi et al., 2016) and BHAC15 (Baraffe et al., 2015) isochrones, and found that, as expected, the choice of isochrone does not change significantly the Bkin results.

Figure 5 presents our results for <<<B>>> as a function of Rossby numbers (left panel), <<<B>>>/Bkin (middle panel), and ϕBsubscriptitalic-ϕB\phi_{\rm B}italic_ϕ start_POSTSUBSCRIPT roman_B end_POSTSUBSCRIPT (right panel) as a function of rotational periods. Stars with Ro>>>0.13 are represented by open red circles, and stars with Ro<<<0.13 are represented by open blue circles. The grey dashed vertical line in the left panel at Ro=0.13 is the threshold that corresponds to the separation between the saturated and the non-saturated regimes.

For stars in the saturated regime in our sample (open blue circles), we derived relations between their mean magnetic fields with Rossby numbers, and between <<<B>>>/Bkin and ϕBsubscriptitalic-ϕB\phi_{\rm B}italic_ϕ start_POSTSUBSCRIPT roman_B end_POSTSUBSCRIPT with rotational periods by using non-linear least squares to fit power-law functions. The obtained relations are found below and these are shown as solid blue lines in Figure 5.

<B>=1604G×Ro0.20expectation𝐵1604𝐺𝑅superscript𝑜0.20<B>=1604G\times Ro^{-0.20}< italic_B > = 1604 italic_G × italic_R italic_o start_POSTSUPERSCRIPT - 0.20 end_POSTSUPERSCRIPT (3)
<B>/Bkin=1.13×Prot0.21expectation𝐵subscript𝐵kin1.13superscriptsubscript𝑃rot0.21<B>/B_{\rm kin}=1.13\times P_{\rm rot}^{-0.21}< italic_B > / italic_B start_POSTSUBSCRIPT roman_kin end_POSTSUBSCRIPT = 1.13 × italic_P start_POSTSUBSCRIPT roman_rot end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 0.21 end_POSTSUPERSCRIPT (4)
ϕB=6.01×1025Mx×Prot0.23subscriptitalic-ϕB6.01superscript1025𝑀𝑥superscriptsubscript𝑃rot0.23\phi_{\rm B}=6.01\times 10^{25}Mx\times P_{\rm rot}^{-0.23}italic_ϕ start_POSTSUBSCRIPT roman_B end_POSTSUBSCRIPT = 6.01 × 10 start_POSTSUPERSCRIPT 25 end_POSTSUPERSCRIPT italic_M italic_x × italic_P start_POSTSUBSCRIPT roman_rot end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 0.23 end_POSTSUPERSCRIPT (5)

For comparison, we also show in Figure 5 similar relations derived by R22 for their sample of rapid and slow rotators (represented by cyan and red dashed lines, respectively).

As previously discussed in Figure 4, the results in the left panel of Figure 5 show a clear difference in the trends between the saturated and unsaturated regimes. The relation obtained here for the saturated stars (Equation 3) is less steep than for the unsaturated ones and similar to the one derived in R22. Given the small number of slow-rotating stars and small range in Rossby number of our sample, we did not compute a best-fit relation for the unsaturated stars but the comparison of our results (open red circles) with the relation in R22 (red dashed line) shows good agreement. The results in the middle panel of Figure 5 show that overall most of our stars present <<<B>>>/Bkin ratios near 1 (represented by the grey horizontal dashed line), indicating that they are near the peak of magnetic energy production based on their available kinetic energy. There is a small trend with rotational periods (or Rossby numbers) showing an inverse correlation between <<<B>>>/Bkin and rotational periods. This is illustrated by our derived relation for rapid rotators (Equation 4), which is very similar to the one derived by R22 (dashed blue line), presenting a similar slope, and being almost the same as our relation. The results in the right panel of Figure 5 show an overall similar behavior as found for <<<B>>> versus Rossby number, with saturated stars presenting a smaller dependency with Prot than non-saturated ones. This is illustrated by the blue solid line showing the relation between magnetic fluxes and rotational periods for the rapid rotator sample (Equation 5). The well-defined relation obtained here for the saturated regime for the Pleiades M dwarfs is different from the results in R22, who found a large scatter in magnetic fluxes for fast-rotating stars and without a well-defined relation. Such difference in the results may come from the fact that R22 studied field stars with different ages, while our sample has a well-constrained age of similar-to\sim 100 million years and, the well-defined trend for the saturated regime may be an indication of a strong connection between magnetic fluxes and stellar ages. The rotational periods and computed Rossby numbers, magnetic fluxes and <<<B>>>/Bkin for the stars are presented in Table 2.

Refer to caption
Figure 5: Comparison of different magnetic field indicators as a function of Rossby numbers (left panel) and rotational periods (middle and right panels). Stars with Ro<<<0.13 are represented by open blue circles, and stars with Ro>>>0.13 are represented by open red circles. The left panel presents the distribution of the obtained magnetic fields, the middle panel presents the ratio between the obtained average magnetic fields and the magnetic field limit based on the stellar kinetic energy and the right panel presents the magnetic flux of the stars. The left panel presents a vertical grey line to represent Ro=13. The horizontal grey line in the middle panel represents the point where the magnetic field is at its maximum physical limit, based on kinetic to magnetic energy conversion, where <<<B>>>= Bkin. We also show relations between rotational and magnetic parameters derived by this work (solid lines) and Reiners et al. (2022) (dashed lines), for slow (red lines) and/or rapid rotators (cyan lines).

4.2 Radius Inflation in the Pleiades M dwarfs

Several works in the literature have compared stellar radii measurements with predicted radii from stellar isochrones that do not include magnetic fields and have found that stars, M dwarf stars in particular, have larger radii than predicted by the models (Reiners et al., 2012; Jackson et al., 2016, 2018, 2019; Jeffers et al., 2018; Kesseli et al., 2018; Jaehnig et al., 2019; Wanderley et al., 2023). Stellar magnetic fields can have the effect of decreasing convective efficiency and/or generating stellar spots, which reduces the stellar photospheric temperature and cause the star to inflate (Chabrier et al., 2007). The term radius inflation (Rfrac, fractional radius inflation) refers to the fractional difference between the radius obtained from measurements and isochrone models (Riso).

In this study, we follow the same methodology as discussed in Wanderley et al. (2023) to measure radius inflation for our sample of Pleiades M dwarfs, using as baseline MIST, DARTMOUTH, and BHAC15 isochrones with solar-metallicities and ages of 100 Myr. As discussed in Section 4.1, stellar radii were derived from the Stefan-Boltzmann equation, using adopted Teff along with luminosities derived from photometric relations. The median±plus-or-minus\pm±MAD radius inflation obtained in this work is 7.0±plus-or-minus\pm±1.5%percent\%%, 6.5±plus-or-minus\pm±1.4%percent\%% and 5.4 ±plus-or-minus\pm±1.3%percent\%%, respectively, based upon the MIST, DARTMOUTH, and BHAC15 isochrones. A median radius inflation of similar-to\sim5–7%percent\%% found here for the Pleiades M dwarfs is considerably larger than the median radius inflation of similar-to\sim2%percent\%% that was obtained by Wanderley et al. (2023) for M-dwarf members of the older (age=625 ±plus-or-minus\pm± 50 Myr, Perryman et al. 1998) Hyades open cluster. Greater radius inflation for younger M dwarfs, is generally in line with expectations from gyrochronology that stars from younger clusters should present higher average magnetic fields if compared to older stars with the same masses. It should be kept in mind, however, that while the sample of Wanderley et al. (2023) included M dwarfs beyond the fully convective threshold, our Pleiades sample is composed only of partially convective M dwarfs, and models including magnetic fields predict that radius inflation can be inhibited in fully convective stars (Feiden et al., 2015).

Refer to caption
Figure 6: Derived average magnetic fields and radius inflation for the sample 62 stars, considering different sets of isochrones. From left to right radius inflation is based on respectively MIST, DARTMOUTH, BHAC15, and SPOTS isochrones.

In addition, in this study, we also computed radius inflation using as a baseline SPOTS isochrone models from Somers et al. (2020). These models consider that stellar spots change the internal structure of the star by suppressing convection, reducing the photospheric temperature, and as a consequence inflating the star. As previously, to derive the fractional radius inflation (Rfrac) and photospheric spot fractions (fspots), we adopted SPOTS isochrones for solar metallicity, and an age of 100 Myr and the same methodology to measure radius inflation presented in Wanderley et al. (2023) (We refer to this previous study for details).

Using SPOTS models, we find a median radius inflation for our sample of 3.0±plus-or-minus\pm±1.2%percent\%%, which is smaller when compared to that obtained with either the MIST, DARTMOUTH or BHAC15 isochrones. Less radius inflation is expected when using Somers et al. (2020) isochrones as a baseline due to the decrease in luminosity for more spotted models. The median radius inflation found here for the Pleiades M dwarfs is larger than the median radius inflation of 1.0±plus-or-minus\pm±0.5%percent\%% reported by Wanderley et al. (2023) for the Hyades.

Figure 6 presents our results for the mean magnetic fields and derived radius inflation (Rfrac) for the sample M dwarfs. Panels from left to right show the radius inflation for the MIST, DARTMOUTH, BHAC15, and SPOTS isochrones. The main result from all panels is that overall Pleiades M dwarfs having higher mean magnetic fields have more inflated radii. The work of Feiden et al. (2015) applied modifications in the DARTMOUTH models to include the effects of magnetic fields and found a positive correlation between magnetic field strength and radius inflation when considering either rotational or turbulent dynamos. This finding agrees with our results, as for all studied isochrones, stars with higher magnetic fields tend to be more inflated. There is, however, a tendency for <<<B>>> to have a weaker dependency with radius inflation, in particular for the SPOTS models with Rfrac between 0.02 and 0.10. Figure 7 presents the derived photospheric stellar spot fractions as a function of average magnetic fields. For the majority of the Pleiades M dwarfs, there is a positive correlation between stellar spot fractions and magnetic field, although with a small number of outliers. Finally, we note that there is one clear outlier star in Figures 6 and 7, having high radius inflation and low <<<B>>>. This star has the lowest Ks mag from our sample (Ks mag = 10.42; Skrutskie et al. (2006)). Its position the J-Ks vs Ks diagram (see Figure 1) might hint that it is not a member of the Pleiades cluster. However, this star was found in Cantat-Gaudin et al. (2020) to have a 100%percent\%% probability of being a member of the Pleiades. Note that our fspots may progressively increase for <<<B>>> larger than roughly 2500 G. Although, Cao & Pinsonneault (2022) previously found constant fspots at low Rossby number (Ro <<< 0.21) in their analysis of Pleiades dwarfs with their values of fspots modeled from a two-component surface defined by a star-spot filling factor and a star-spot temperature contrast.

Refer to caption
Figure 7: Fractional stellar photospheric spots coverage from SPOTS models (Somers et al., 2020), as a function of our derived average magnetic fields.

4.3 Magnetic Fields and Activity

Stellar activity is a term that represents the stellar variability phenomenon which is mainly a consequence of strong magnetic fields. Stellar variability can be caused by stellar spots, which are temporary cooler regions in the stellar photosphere where the strong magnetic fields suppress convection. As the star rotates these cold spots can reduce the observed stellar flux, and result in periodic changes in the magnitude of the star. Another important mechanism that can generate stellar variability is high-energy non-thermal emission. Strong magnetic fields are responsible for heating the coronae and chromosphere of active stars, which results in considerably larger amounts of high-energy radiation emission than expected from their black-body profiles. Heating of the coronae results in excess in the X-ray stellar output that can be measured as the ratio between the X-ray and bolometric luminosity, while heating of the chromosphere is often studied from magnetic sensitive emission lines, such as Hα𝛼\alphaitalic_α, Hβ𝛽\betaitalic_β, and Ca II H and K lines. Since coronal and chromospheric non-thermal emissions are produced by the effect of stellar magnetic fields, they are also variable, and change according to the stellar magnetic cycle.

To study the relation between activity and magnetic fields in the Pleiades M dwarfs, we cross-matched our sample with the targets in the X-ray studies by Núñez & Agüeros (2016) and Wright et al. (2011), and found 31 stars in common with the first study (not considering 17 stars having only upper limit measurements), and 19 stars in common with the latter study; the X-ray to bolometric luminosity ratios in those works are X-ray flux measurements from Einstein and ROSAT observations (Micela et al., 1990; Stauffer et al., 1994; Micela et al., 1996, 1999; Stelzer et al., 2000). The results are summarized in the top panel of Figure 8 where we show the logarithm of the X-ray to bolometric luminosity ratios versus our derived magnetic fields. There is a clear correlation between the two variables, with a Pearson correlation coefficient of 0.46 for the sample in Núñez & Agüeros (2016) (shown as maroon circles) and 0.35 for the sample in Wright et al. (2011) (shown as orange circles), indicating that stars with stronger magnetic fields tend to present greater X-ray to bolometric luminosity ratios.

Refer to caption
Figure 8: Derived average magnetic fields as a function of the ratios between X-ray to bolometric luminosities for stars of our sample in common with the study of Wright et al. (2011) (orange circles), and Núñez & Agüeros (2016) (maroon circles). The same stars are connected by dashed lines.

Concerning chromospheric activity the work of Fang et al. (2018) measured equivalent widths for Hα𝛼\alphaitalic_α, Hβ𝛽\betaitalic_β, and Ca II K emission lines for stars in open clusters using LAMOST DR3 spectra (Cui et al., 2012). We cross-matched our Pleiades sample with their database and found 33 stars in common with measured Hα𝛼\alphaitalic_α, 31 with measured Hβ𝛽\betaitalic_β, and 19 with measured Ca II K. The three panels in Figure 9 present our derived average magnetic field measurements versus the total equivalent width for the Ca II K line (top panel), Hβ𝛽\betaitalic_β (middle panel), and Hα𝛼\alphaitalic_α (bottom panel), all equivalent widths are in units of Å. Similarly to what was found for the X-ray to bolometric ratios, there is a correlation between the mean magnetic fields for the Pleiades M dwarfs and the equivalent widths of the emission lines for all three lines, although the results for the Ca II K line are less clear.

Refer to caption
Figure 9: Derived average magnetic fields as a function of equivalent widths of magnetic sensitive chromospheric emission lines for stars of our sample that are in common with the work of Fang et al. (2018). The panels present from top to bottom: Ca II K, Hβ𝛽\betaitalic_β and Hα𝛼\alphaitalic_α emission lines.

Figure 10 shows the Hα𝛼\alphaitalic_α to bolometric luminosity ratios as a function of mean magnetic field for stars in our sample having Hα𝛼\alphaitalic_α emission strength measurements (EW) reported in Table 2. To convert Hα𝛼\alphaitalic_α equivalent widths into L/Lbol we employed the methodology described in Stassun et al. (2012). Similarly to what is seen for the relation between EW and <<<B>>> (bottom panel of Figure 9), there is a clear correlation between Hα𝛼\alphaitalic_α to bolometric luminosity ratios and the derived magnetic fields, with a Pearson correlation coefficient of 0.78.

Refer to caption
Figure 10: Derived average magnetic fields as a function of Hα𝛼\alphaitalic_α to bolometric luminosities. We employed the methodology described in Stassun et al. (2012) to convert Hα𝛼\alphaitalic_α equivalent widths (taken from Fang et al. (2018)) into L/Lbol.

Overall, it is reassuring that in this study we find that the more magnetic stars in our sample tend also to be more active, and this is based on our magnetic field measurements and independent results from activity indicators from other works in the literature. The X-ray and Hα𝛼\alphaitalic_α to bolometric luminosity ratios as well as the equivalent widths for Hα𝛼\alphaitalic_α, Hβ𝛽\betaitalic_β and Ca II K of the stars are presented in Table 2.

5 Conclusions

We used the SYNMAST code Kochukhov et al. (2010), along with a Markov Chain Monte Carlo (MCMC) methodology, to compute synthetic spectra and analyze magnetically sensitive Fe I lines to derive average magnetic fields for 62 M dwarf members of the young (age = 112 ±plus-or-minus\pm± 5 Myr; Dahm 2015) and nearby Pleiades open cluster. This analysis is based on the SDSS IV APOGEE spectra (Majewski et al., 2017; Abdurro’uf et al., 2022), the APOGEE line list (Smith et al., 2021), effective Landé-g factors from VALD (Piskunov et al., 1995; Kupka et al., 1999), and model atmospheres from the MARCS grid (Gustafsson et al., 2008).

A search was carried out to find the best Fe I lines in the APOGEE spectral region that were sensitive to Zeeman broadening and could be used as diagnostic lines for measuring mean magnetic fields in M dwarfs, with four Fe I lines identified and selected as the best indicators: λ𝜆\lambdaitalic_λ15207.526 Å, λ𝜆\lambdaitalic_λ15294.56 Å, λ𝜆\lambdaitalic_λ15621.654 Å, and λ𝜆\lambdaitalic_λ15631.948 Å. The derived mean magnetic fields for the studied Pleiades M dwarf stars range from similar-to\sim 1.0 to similar-to\sim 4.2 kG, with a median±plus-or-minus\pm±MAD of 3.0±plus-or-minus\pm±0.3 kG, not reaching the lowest magnetic field and rotation levels reported in other studies that explored field stars of similar masses, which is probably explained by the young age of the cluster.

The derived mean magnetic field measurements in the Pleiades M dwarfs were used to study correlations with Rossby number (Ro) and stellar rotation. The Rossby number is given by the ratio between the rotational period and the convective turnover time, being an important indicator of stellar activity. We find a clear trend that more magnetic stars have, on average, higher projected rotational velocities, vsini𝑖\sin{i}roman_sin italic_i. The Rossby number was used to separate our sample into rapid (Ro<<<0.13) and slow rotators (Ro>>>0.13). Overall, our results for <<<B>>> versus Prot and B versus Rossby number overlap with results from the literature for field stars, and indicate that the population of stars with Ro<<<0.13 exhibit a steeper relation between magnetic field and rotational period, or Rossby number, when compared to stars with Ro>>>0.13. However, even for Ro>>>0.13, there remains a shallow trend between Rossby number and magnetic field, which is given by: <<<B>>> = 1604×1604\times1604 ×Ro-0.20.

For this sample of Pleiades M dwarfs, we also investigated the ratio between mean magnetic fields and the maximum magnetic field limit (<<<B>>>/Bkin) that is reached based on the kinetic-to-magnetic energy conversion. It is found that most of the studied Pleiades M dwarfs are at the limit of kinetic to magnetic energy conversion, or, are in the saturated regime, having <<<B>>>/Bkin{}_{\rm kin}\approxstart_FLOATSUBSCRIPT roman_kin end_FLOATSUBSCRIPT ≈1. Unlike previous results in the literature for field stars in the saturated regime, the computed stellar magnetic fluxes ϕBsubscriptitalic-ϕB\phi_{\rm B}italic_ϕ start_POSTSUBSCRIPT roman_B end_POSTSUBSCRIPT as a function of Protrot{\rm rot}roman_rot for the Pleiades M dwarfs show similar trends obtained for <<<B>>> versus Protrot{\rm rot}roman_rot.

Another important effect of magnetic fields is to inflate the stellar radii of cool dwarfs. Radius inflation corresponds to the fractional difference (Rfrac) between the radius obtained from measurements and predictions from isochrone models. In this study, we derived the radii for the studied M dwarfs and used MIST, DARTMOUTH, BHAC15, and SPOTS isochrones as baselines to infer their radius inflation. We obtain a median±plus-or-minus\pm±MAD radius inflation for our sample of respectively 7.0±plus-or-minus\pm±1.5%percent\%%, 6.5±plus-or-minus\pm±1.4%percent\%%, 5.4 ±plus-or-minus\pm±1.3%percent\%% and 3.0 ±plus-or-minus\pm±1.2%percent\%%, being more inflated than M dwarfs in the older Hyades open cluster Wanderley et al. (2023). For the Pleiades, it is noticeable that for SPOTS isochrones there is less radius inflation when compared to non-spotted isochrones, as expected. In addition, our results indicate that more magnetic stars have more inflated radii, showing a correlation between Rfrac and <<<B>>>. For SPOTS models in particular, there is, however, a tendency for <<<B>>> to exhibit a weaker dependency with radius inflation for Rfrac between 0.02 and 0.10. For the majority of the Pleiades M dwarfs in our sample, there is a positive correlation between stellar spot fraction and magnetic field, although with a small number of outliers.

To study the relation between chromospheric stellar activity and magnetic fields in the Pleiades M dwarfs, we compared our derived mean magnetic fields with high-energy non-thermal emission indicators, such as equivalent width measurements of the lines of Hα𝛼\alphaitalic_α, Hβ𝛽\betaitalic_β, and Ca II K, as well as Hα𝛼\alphaitalic_α to bolometric luminosity ratios. For all of these, we find a positive correlation between chromospheric activity and magnetic fields. A positive correlation was also obtained between the mean magnetic field and the ratio between X-ray to bolometric luminosity, which is an important indicator of coronal activity. Overall, it is reassuring that the more magnetic stars in this study also tend to be more active, and this is based on our magnetic field measurements from APOGEE spectra, which is independent of the results for activity indicators obtained from the literature. This study opens a new window into using the APOGEE survey to investigate magnetic fields in cool stars.

We thank the referee for the comments that improved the paper. F.W. acknowledges support from fellowship by Coordenação de Ensino Superior - CAPES. K.C. and V.S. acknowledge that their work is supported, in part, by the National Science Foundation through NSF grant No. AST-2009507. O.K. acknowledges support by the Swedish Research Council (grant agreements no. 2019-03548 and 2023-03667), the Swedish National Space Agency, and the Royal Swedish Academy of Sciences. D.S. thanks the National Council for Scientific and Technological Development – CNPq. Funding for the Sloan Digital Sky Survey IV has been provided by the Alfred P. Sloan Foundation, the U.S. Department of Energy Office of Science, and the Participating Institutions. SDSS-IV acknowledges support and resources from the Center for High-Performance Computing at the University of Utah. The SDSS web site is www.sdss.org. SDSS-IV is managed by the Astrophysical Research consortium for the Participating Institutions of the SDSS Collaboration including the Brazilian Participation Group, the Carnegie Institution for Science, Carnegie Mellon University, the Chilean Participation Group, the French Participation Group, Harvard-Smithsonian Center for Astrophysics, Instituto de Astrofísica de Canarias, The Johns Hopkins University, Kavli Institute for the Physics and Mathematics of the Universe (IPMU) / University of Tokyo, Lawrence Berkeley National Laboratory, Leibniz Institut für Astrophysik Potsdam (AIP), Max-Planck-Institut für Astronomie (MPIA Heidelberg), Max-Planck-Institut für Astrophysik (MPA Garching), Max-Planck-Institut für Extraterrestrische Physik (MPE), National Astronomical Observatory of China, New Mexico State University, New York University, University of Notre Dame, Observatório Nacional / MCTI, The Ohio State University, Pennsylvania State University, Shanghai Astronomical Observatory, United Kingdom Participation Group, Universidad Nacional Autónoma de México, University of Arizona, University of Colorado Boulder, University of Oxford, University of Portsmouth, University of Utah, University of Virginia, University of Washington, University of Wisconsin, Vanderbilt University, and Yale University. Sloan \startlongtable
Table 1: Filling Factors
APOGEE ID f0 f2 f4 f6 f8 f10 f12 <<<B>>>
G
2M03420291+2355538 0.218 ±plus-or-minus\pm±0.087 0.274 ±plus-or-minus\pm±0.152 0.147 ±plus-or-minus\pm±0.134 0.034 ±plus-or-minus\pm±0.153 0.246 ±plus-or-minus\pm±0.22 0.037 ±plus-or-minus\pm±0.049 0.044 ±plus-or-minus\pm±0.098 4201 529511subscriptsuperscriptabsent511529{}^{511}_{529}start_FLOATSUPERSCRIPT 511 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT 529 end_POSTSUBSCRIPT
2M03481801+2353294 0.784 ±plus-or-minus\pm±0.058 0.113 ±plus-or-minus\pm±0.064 0.037 ±plus-or-minus\pm±0.03 0.011 ±plus-or-minus\pm±0.011 0.013 ±plus-or-minus\pm±0.01 0.005 ±plus-or-minus\pm±0.025 0.036 ±plus-or-minus\pm±0.016 1032 215236subscriptsuperscriptabsent236215{}^{236}_{215}start_FLOATSUPERSCRIPT 236 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT 215 end_POSTSUBSCRIPT
\startlongtable
Table 2: Stellar Data
APOGEE ID SNR vsini𝑖\sin{i}roman_sin italic_i Teff Prot Ro <<<B>>> <<<B>>>/Bkin ϕBsubscriptitalic-ϕB\phi_{\rm B}italic_ϕ start_POSTSUBSCRIPT roman_B end_POSTSUBSCRIPT log(LX/Lbol)subscript𝐿Xsubscript𝐿bol\log{(L_{\rm X}/L_{\rm bol})}roman_log ( italic_L start_POSTSUBSCRIPT roman_X end_POSTSUBSCRIPT / italic_L start_POSTSUBSCRIPT roman_bol end_POSTSUBSCRIPT ) log(LHα/Lbol)subscript𝐿𝐻𝛼subscript𝐿bol\log{(L_{H\alpha}/L_{\rm bol})}roman_log ( italic_L start_POSTSUBSCRIPT italic_H italic_α end_POSTSUBSCRIPT / italic_L start_POSTSUBSCRIPT roman_bol end_POSTSUBSCRIPT ) EW EW EWCaIIK
km s-1 K day G 1025superscript102510^{25}10 start_POSTSUPERSCRIPT 25 end_POSTSUPERSCRIPT Mx Å Å Å
2M03420291+2355538 138 17.7 3424 0.89 0.02 4202530+511subscriptsuperscriptabsent511530{}^{+511}_{-530}start_FLOATSUPERSCRIPT + 511 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 530 end_POSTSUBSCRIPT 1.52 9.6 a/-2.92b 5.43 -3.74 6 6.55
2M03424379+2532064 103 7.2 3443 3362615+669subscriptsuperscriptabsent669615{}^{+669}_{-615}start_FLOATSUPERSCRIPT + 669 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 615 end_POSTSUBSCRIPT 1.06 3.81 a/…b
2M03515114+2317414 112 13.4 3452 0.87 0.01 3463567+581subscriptsuperscriptabsent581567{}^{+581}_{-567}start_FLOATSUPERSCRIPT + 581 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 567 end_POSTSUBSCRIPT 1.12 4.64 a/…b 6.1 -3.68 7.2
2M03422864+2501004 131 15.2 3462 1.39 0.02 3850478+449subscriptsuperscriptabsent449478{}^{+449}_{-478}start_FLOATSUPERSCRIPT + 449 end_FLOATSUPERSCRIPT start_POSTSUBSCRIPT - 478 end_POSTSUBSCRIPT 1.21 4.61 a/-3.13b 5.58 -3.72 9.46 24.76
  • The complete Table is available in electronic format.

  • aWright et al. (2011)

  • bNúñez & Agüeros (2016)

\startlongtable
Table 3: Diagnostic lines
λ𝜆\lambdaitalic_λ χexcsubscript𝜒exc\chi_{\rm exc}italic_χ start_POSTSUBSCRIPT roman_exc end_POSTSUBSCRIPT loggfgf\log{\rm gf}roman_log roman_gf Effective Landé-g Lower level term designation Upper level term designation
Å̊𝐴\mathring{A}over̊ start_ARG italic_A end_ARG eV
15207.526 5.3852 0.067 1.532 LS 3d6.(5D).4s (6D).5s e7D LS 3d6.(5D).4s (6D).5p n7D*
15294.560 5.3085 0.523 1.590 LS 3d6.(5D).4s (6D).5s e7D LS 3d6.(5D).4s (6D).5p n7D*
15621.654 5.5392 0.280 1.494 LS 3d6.(5D).4s (6D).5s e5D LS 3d6.(5D).4s (6D).5p t5D*
15631.948 5.3516 -0.032 1.655 LS 3d6.(5D).4s (6D).5s e7D LS 3d6.(5D).4s (6D).5p n7D*

References

  • Abdurro’uf et al. (2022) Abdurro’uf, Accetta, K., Aerts, C., et al. 2022, ApJS, 259, 35
  • Afram et al. (2009) Afram, N., Reiners, A., & Berdyugina, S. V. 2009, in Astronomical Society of the Pacific Conference Series, Vol. 405, Solar Polarization 5: In Honor of Jan Stenflo, ed. S. V. Berdyugina, K. N. Nagendra, & R. Ramelli, 527
  • Astudillo-Defru et al. (2017) Astudillo-Defru, N., Delfosse, X., Bonfils, X., et al. 2017, A&A, 600, A13
  • Bailer-Jones et al. (2021) Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, AJ, 161, 147
  • Bakos et al. (2004) Bakos, G., Noyes, R. W., Kovács, G., et al. 2004, PASP, 116, 266
  • Baraffe et al. (2015) Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, A&A, 577, A42
  • Barnes (2003) Barnes, S. A. 2003, ApJ, 586, 464
  • Basri & Marcy (1994) Basri, G., & Marcy, G. W. 1994, ApJ, 431, 844
  • Basri et al. (1992) Basri, G., Marcy, G. W., & Valenti, J. A. 1992, ApJ, 390, 622
  • Belokurov et al. (2020) Belokurov, V., Penoyre, Z., Oh, S., et al. 2020, MNRAS, 496, 1922
  • Blanton et al. (2017) Blanton, M. R., Bershady, M. A., Abolfathi, B., et al. 2017, AJ, 154, 28
  • Bowen & Vaughan (1973) Bowen, I. S., & Vaughan, A. H., J. 1973, Appl. Opt., 12, 1430
  • Bressan et al. (2012) Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127
  • Cantat-Gaudin et al. (2020) Cantat-Gaudin, T., Anders, F., Castro-Ginard, A., et al. 2020, A&A, 640, A1
  • Cao & Pinsonneault (2022) Cao, L., & Pinsonneault, M. H. 2022, MNRAS, 517, 2165
  • Chabrier et al. (2007) Chabrier, G., Gallardo, J., & Baraffe, I. 2007, A&A, 472, L17
  • Choi et al. (2016) Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102
  • Covey et al. (2016) Covey, K. R., Agüeros, M. A., Law, N. M., et al. 2016, ApJ, 822, 81
  • Cristofari et al. (2023a) Cristofari, P. I., Donati, J. F., Folsom, C. P., et al. 2023a, MNRAS, 522, 1342
  • Cristofari et al. (2023b) Cristofari, P. I., Donati, J. F., Moutou, C., et al. 2023b, MNRAS, 526, 5648
  • Cui et al. (2012) Cui, X.-Q., Zhao, Y.-H., Chu, Y.-Q., et al. 2012, Research in Astronomy and Astrophysics, 12, 1197
  • Dahm (2015) Dahm, S. E. 2015, ApJ, 813, 108
  • Dotter et al. (2008) Dotter, A., Chaboyer, B., Jevremović, D., et al. 2008, ApJS, 178, 89
  • Fang et al. (2018) Fang, X.-S., Zhao, G., Zhao, J.-K., & Bharat Kumar, Y. 2018, MNRAS, 476, 908
  • Feiden et al. (2015) Feiden, G. A., Jones, J., & Chaboyer, B. 2015, in Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, Vol. 18, 18th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, 171
  • Foreman-Mackey et al. (2013) Foreman-Mackey, D., Conley, A., Meierjurgen Farr, W., et al. 2013, emcee: The MCMC Hammer, Astrophysics Source Code Library, record ascl:1303.002
  • Gaia Collaboration (2022) Gaia Collaboration. 2022, VizieR Online Data Catalog: Gaia DR3 Part 1. Main source (Gaia Collaboration, 2022), VizieR On-line Data Catalog: I/355. Originally published in: Astron. Astrophys., in prep. (2022)
  • García Pérez et al. (2016) García Pérez, A. E., Allende Prieto, C., Holtzman, J. A., et al. 2016, AJ, 151, 144
  • Gunn et al. (2006) Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, AJ, 131, 2332
  • Gustafsson et al. (2008) Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951
  • Hahlin & Kochukhov (2022) Hahlin, A., & Kochukhov, O. 2022, A&A, 659, A151
  • Hahlin et al. (2021) Hahlin, A., Kochukhov, O., Alecian, E., Morin, J., & BinaMIcS Collaboration. 2021, A&A, 650, A197
  • Hahlin et al. (2023) Hahlin, A., Kochukhov, O., Rains, A. D., et al. 2023, A&A, 675, A91
  • Han et al. (2023) Han, E., López-Valdivia, R., Mace, G. N., & Jaffe, D. T. 2023, AJ, 166, 4
  • Harris et al. (2020) Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357
  • Hartman et al. (2010) Hartman, J. D., Bakos, G. Á., Kovács, G., & Noyes, R. W. 2010, MNRAS, 408, 475
  • Hawley et al. (2014) Hawley, S. L., Davenport, J. R. A., Kowalski, A. F., et al. 2014, ApJ, 797, 121
  • Heyl et al. (2022) Heyl, J., Caiazzo, I., & Richer, H. B. 2022, ApJ, 926, 132
  • Howell et al. (2014) Howell, S. B., Sobeck, C., Haas, M., et al. 2014, PASP, 126, 398
  • Hubeny & Lanz (2011) Hubeny, I., & Lanz, T. 2011, Synspec: General Spectrum Synthesis Program, Astrophysics Source Code Library, record ascl:1109.022
  • Hunter (2007) Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90
  • Jackson et al. (2018) Jackson, R. J., Deliyannis, C. P., & Jeffries, R. D. 2018, MNRAS, 476, 3245
  • Jackson et al. (2019) Jackson, R. J., Jeffries, R. D., Deliyannis, C. P., Sun, Q., & Douglas, S. T. 2019, MNRAS, 483, 1125
  • Jackson et al. (2016) Jackson, R. J., Jeffries, R. D., Randich, S., et al. 2016, A&A, 586, A52
  • Jaehnig et al. (2019) Jaehnig, K., Somers, G., & Stassun, K. G. 2019, ApJ, 879, 39
  • Jeffers et al. (2018) Jeffers, S. V., Schöfer, P., Lamert, A., et al. 2018, A&A, 614, A76
  • Johns-Krull (2007) Johns-Krull, C. M. 2007, ApJ, 664, 975
  • Johns-Krull & Valenti (1996) Johns-Krull, C. M., & Valenti, J. A. 1996, ApJ, 459, L95
  • Johns-Krull & Valenti (2000) Johns-Krull, C. M., & Valenti, J. A. 2000, in Astronomical Society of the Pacific Conference Series, Vol. 198, Stellar Clusters and Associations: Convection, Rotation, and Dynamos, ed. R. Pallavicini, G. Micela, & S. Sciortino, 371
  • Johns-Krull et al. (2004) Johns-Krull, C. M., Valenti, J. A., & Saar, S. H. 2004, ApJ, 617, 1204
  • Kawaler (1988) Kawaler, S. D. 1988, ApJ, 333, 236
  • Kesseli et al. (2018) Kesseli, A. Y., Muirhead, P. S., Mann, A. W., & Mace, G. 2018, AJ, 155, 225
  • Kochukhov (2016) Kochukhov, O. 2016, in Lecture Notes in Physics, Berlin Springer Verlag, ed. J.-P. Rozelot & C. Neiner, Vol. 914, 177
  • Kochukhov (2021) —. 2021, A&A Rev., 29, 1
  • Kochukhov et al. (2009) Kochukhov, O., Heiter, U., Piskunov, N., et al. 2009, in American Institute of Physics Conference Series, Vol. 1094, 15th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun, ed. E. Stempels, 124
  • Kochukhov & Lavail (2017) Kochukhov, O., & Lavail, A. 2017, ApJ, 835, L4
  • Kochukhov et al. (2010) Kochukhov, O., Makaganiuk, V., & Piskunov, N. 2010, A&A, 524, A5
  • Kochukhov et al. (2001) Kochukhov, O., Piskunov, N., Valenti, J., & Johns-Krull, C. 2001, in Astronomical Society of the Pacific Conference Series, Vol. 248, Magnetic Fields Across the Hertzsprung-Russell Diagram, ed. G. Mathys, S. K. Solanki, & D. T. Wickramasinghe, 219
  • Kochukhov & Reiners (2020) Kochukhov, O., & Reiners, A. 2020, ApJ, 902, 43
  • Kochukhov & Shulyak (2019) Kochukhov, O., & Shulyak, D. 2019, ApJ, 873, 69
  • Kupka et al. (1999) Kupka, F., Piskunov, N., Ryabchikova, T. A., Stempels, H. C., & Weiss, W. W. 1999, A&AS, 138, 119
  • Lasker et al. (2008) Lasker, B. M., Lattanzi, M. G., McLean, B. J., et al. 2008, AJ, 136, 735
  • Lavail et al. (2019) Lavail, A., Kochukhov, O., & Hussain, G. A. J. 2019, A&A, 630, A99
  • Law et al. (2009) Law, N. M., Kulkarni, S. R., Dekany, R. G., et al. 2009, PASP, 121, 1395
  • Lodieu et al. (2019) Lodieu, N., Pérez-Garrido, A., Smart, R. L., & Silvotti, R. 2019, A&A, 628, A66
  • Majewski et al. (2017) Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2017, AJ, 154, 94
  • Mann et al. (2015) Mann, A. W., Feiden, G. A., Gaidos, E., Boyajian, T., & von Braun, K. 2015, ApJ, 804, 64
  • Mann et al. (2016) —. 2016, ApJ, 819, 87
  • Micela et al. (1996) Micela, G., Sciortino, S., Kashyap, V., Harnden, F. R., J., & Rosner, R. 1996, ApJS, 102, 75
  • Micela et al. (1990) Micela, G., Sciortino, S., Vaiana, G. S., et al. 1990, ApJ, 348, 557
  • Micela et al. (1999) Micela, G., Sciortino, S., Harnden, F. R., J., et al. 1999, A&A, 341, 751
  • Muirhead et al. (2018) Muirhead, P. S., Dressing, C. D., Mann, A. W., et al. 2018, AJ, 155, 180
  • Newton et al. (2017) Newton, E. R., Irwin, J., Charbonneau, D., et al. 2017, ApJ, 834, 85
  • Newton et al. (2016) Newton, E. R., Irwin, J., Charbonneau, D., Berta-Thompson, Z. K., & Dittmann, J. A. 2016, ApJ, 821, L19
  • Nguyen et al. (2022) Nguyen, C. T., Costa, G., Girardi, L., et al. 2022, A&A, 665, A126
  • Nidever et al. (2015) Nidever, D. L., Holtzman, J. A., Allende Prieto, C., et al. 2015, AJ, 150, 173
  • Núñez & Agüeros (2016) Núñez, A., & Agüeros, M. A. 2016, ApJ, 830, 44
  • Perryman et al. (1998) Perryman, M. A. C., Brown, A. G. A., Lebreton, Y., et al. 1998, A&A, 331, 81
  • Phan-Bao et al. (2009) Phan-Bao, N., Lim, J., Donati, J.-F., Johns-Krull, C. M., & Martín, E. L. 2009, ApJ, 704, 1721
  • Piskunov et al. (1995) Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 112, 525
  • Plez (2012) Plez, B. 2012, Turbospectrum: Code for spectral synthesis
  • Pouilly et al. (2023) Pouilly, K., Kochukhov, O., Kóspál, Á., et al. 2023, MNRAS, 518, 5072
  • Rau et al. (2009) Rau, A., Kulkarni, S. R., Law, N. M., et al. 2009, PASP, 121, 1334
  • Rebull et al. (2016) Rebull, L. M., Stauffer, J. R., Bouvier, J., et al. 2016, AJ, 152, 114
  • Reid & Gizis (1997) Reid, I. N., & Gizis, J. E. 1997, AJ, 113, 2246
  • Reiners et al. (2009) Reiners, A., Basri, G., & Christensen, U. R. 2009, ApJ, 697, 373
  • Reiners et al. (2012) Reiners, A., Joshi, N., & Goldman, B. 2012, AJ, 143, 93
  • Reiners et al. (2022) Reiners, A., Shulyak, D., Käpylä, P. J., et al. 2022, A&A, 662, A41
  • Saar (1994) Saar, S. H. 1994, in Infrared Solar Physics, ed. D. M. Rabin, J. T. Jefferies, & C. Lindsey, Vol. 154, 493
  • Saar (1996) Saar, S. H. 1996, in Stellar Surface Structure, ed. K. G. Strassmeier & J. L. Linsky, Vol. 176, 237
  • Saar & Linsky (1985) Saar, S. H., & Linsky, J. L. 1985, ApJ, 299, L47
  • Salpeter (1955) Salpeter, E. E. 1955, ApJ, 121, 161
  • Shulyak et al. (2017) Shulyak, D., Reiners, A., Engeln, A., et al. 2017, Nature Astronomy, 1, 0184
  • Shulyak et al. (2014) Shulyak, D., Reiners, A., Seemann, U., Kochukhov, O., & Piskunov, N. 2014, A&A, 563, A35
  • Shulyak et al. (2011) Shulyak, D., Seifahrt, A., Reiners, A., Kochukhov, O., & Piskunov, N. 2011, MNRAS, 418, 2548
  • Shulyak et al. (2019) Shulyak, D., Reiners, A., Nagel, E., et al. 2019, A&A, 626, A86
  • Skrutskie et al. (2006) Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163
  • Skumanich (1972) Skumanich, A. 1972, ApJ, 171, 565
  • Smith et al. (2021) Smith, V. V., Bizyaev, D., Cunha, K., et al. 2021, AJ, 161, 254
  • Soderblom et al. (2009) Soderblom, D. R., Laskar, T., Valenti, J. A., Stauffer, J. R., & Rebull, L. M. 2009, AJ, 138, 1292
  • Somers et al. (2020) Somers, G., Cao, L., & Pinsonneault, M. H. 2020, ApJ, 891, 29
  • Souto et al. (2021) Souto, D., Cunha, K., & Smith, V. V. 2021, ApJ, 917, 11
  • Souto et al. (2017) Souto, D., Cunha, K., García-Hernández, D. A., et al. 2017, ApJ, 835, 239
  • Souto et al. (2020) Souto, D., Cunha, K., Smith, V. V., et al. 2020, ApJ, 890, 133
  • Stassun et al. (2012) Stassun, K. G., Kratter, K. M., Scholz, A., & Dupuy, T. J. 2012, ApJ, 756, 47
  • Stauffer et al. (1994) Stauffer, J. R., Caillault, J. P., Gagne, M., Prosser, C. F., & Hartmann, L. W. 1994, ApJS, 91, 625
  • Stauffer et al. (2007) Stauffer, J. R., Hartmann, L. W., Fazio, G. G., et al. 2007, ApJS, 172, 663
  • Stelzer et al. (2000) Stelzer, B., Neuhäuser, R., & Hambaryan, V. 2000, A&A, 356, 949
  • Stift & Leone (2003) Stift, M. J., & Leone, F. 2003, A&A, 398, 411
  • Wanderley et al. (2023) Wanderley, F., Cunha, K., Souto, D., et al. 2023, ApJ, 951, 90
  • Wang & Chen (2019) Wang, S., & Chen, X. 2019, ApJ, 877, 116
  • Wilson et al. (2019) Wilson, J. C., Hearty, F. R., Skrutskie, M. F., et al. 2019, PASP, 131, 055001
  • Wright et al. (2011) Wright, N. J., Drake, J. J., Mamajek, E. E., & Henry, G. W. 2011, ApJ, 743, 48
  • Yang & Johns-Krull (2011) Yang, H., & Johns-Krull, C. M. 2011, ApJ, 729, 83
  • Yang et al. (2008) Yang, H., Johns-Krull, C. M., & Valenti, J. A. 2008, AJ, 136, 2286
  • Zacharias et al. (2012) Zacharias, N., Finch, C. T., Girard, T. M., et al. 2012, VizieR Online Data Catalog, I/322A