Astronomy &amp; Astrophysics manuscript no. afram˙final

August 2, 2019

## Complexity of magnetic fields on red dwarfs /star

N. Afram 1 and S. V. Berdyugina 1 2 ,

1 Leibniz-Institut f¨r Sonnenphysik (KIS), Freiburg, Germany u

2 Kavli Institute for Theoretical Physics, UC Santa Barbara, CA, USA

Received date; accepted date

## ABSTRACT

Context. Magnetic fields in cool stars can be investigated by measuring Zeeman line broadening and polarization in atomic and molecular lines. Similar to the Sun, these fields are complex and height-dependent. Many molecular lines dominating M-dwarf spectra (e.g., FeH, CaH, MgH, and TiO) are temperature- and Zeeman- sensitive and form at di ff erent atmospheric heights, which makes them excellent probes of magnetic fields on M dwarfs.

Aims. Our goal is to analyze the complexity of magnetic fields in M dwarfs. We investigate how magnetic fields vary with the stellar temperature and how 'surface' inhomogeneities are distributed in height - the dimension that is usually neglected in stellar magnetic studies.

Methods. We have determined e ff ective temperatures of the photosphere and of magnetic features, magnetic field strengths and filling factors for nine M dwarfs (M1-M7). Our χ 2 analysis is based on a comparison of observed and synthetic intensity and circular polarization profiles. Stokes profiles were calculated by solving polarized radiative transfer equations.

Results. Properties of magnetic structures depend on the analyzed atomic or molecular species and their formation heights. Two types of magnetic features similar to those on the Sun have been found: a cooler (starspots) and a hotter (network) one. The magnetic field strength in both starspots and network is within 3 kG to 6 kG, on average it is 5 kG. These fields occupy a large fraction of M dwarf atmospheres at all heights, up to 100%. The plasma β is less than one, implying highly magnetized stars.

Conclusions. A combination of molecular and atomic species and a simultaneous analysis of intensity and circular polarization spectra have allowed us to better decipher the complexity of magnetic fields on M dwarfs, including their dependence on the atmospheric height. This work provides an opportunity to investigate a larger sample of M dwarfs and L-type brown dwarfs.

Key words. Molecular processes - Polarization - Radiative transfer - Line: formation - Stars: Magnetic fields

## 1. Introduction

Magnetic fields are pervasive in a variety of astrophysical objects, ranging from planets to stars, and to galaxies. They strongly influence the object's structure, dynamics, and evolution. This fact is now being incorporated into evolutionary models. Furthermore, increasing attention has been put on the impact of stellar activity on planets, where the host star's magnetic fields perturb the planetary magnetosphere, which is thought to be essential in providing a protective shield for the evolution of life.

the bottom of the convection zone like in the Sun, where at least the cyclic part of the solar dynamo is believed to operate. Thus, the existence of activity on those stars requires an alternative mechanism to generate magnetic fields. It is still unknown how the dynamo mechanism on late-type dwarfs changes as the mass reduces from solar type to fully convective objects. For early-type M dwarfs, it is possible that both the solar type and turbulent dynamos coexist, whereas the fully convective latetype M dwarfs and brown dwarfs can certainly maintain only a turbulent (distributed) dynamo.

One of the fundamental astrophysical challenges is the understanding of the generation, amplification, and complexity of stellar magnetic fields. Stellar surface activity arises in the presence of a convection zone, which appears in F stars and takes over the whole star roughly at a spectral type M3 to M4 (Chabrier &amp; Bara ff e 1997; West et al. 2008). These fully convective low-mass main-sequence stars lack an interface layer at

/star The spectra of the stars analyzed in this work are also available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http: // cdsweb.u-strasbg.fr cgi-bin qcat?J A / / / + A . /

Mdwarfs are the most common type of stars in the galaxy. They are in the center of attention in the exoplanetary research field as they are numerous and can host Earth-size planets in the star's habitable zone. These can be detected even for low-luminosity M dwarfs (e.g., Anglada-Escud´ e et al. 2016; Gillon et al. 2017). Their strong activity also a ff ects the planetary atmospheres and magnetospheres. In addition, the knowledge of the star's magnetic field complexity and topology, and possible changes thereof, can help to reduce the uncertainties in searching for exoplanets around M dwarfs.

c © ESO 2019

The direct detection of stellar magnetism is still a di ffi cult task. Even for our closest neighbor, the Sun, mysteries about the magnetic fields remain unsolved. Considering the amount of information on solar magnetic fields, it is evident that the gaps in our knowledge about magnetic fields on M dwarfs are enormous. An overview of magnetic fields on cool stars can be found in Berdyugina (2005, 2009), Strassmeier (2009, 2011), and Reiners (2012).

Stellar magnetic fields are estimated by using the Zeeman e ff ect in atomic and molecular spectral lines. To fully describe the magnetic field vector, one needs to measure all states of line polarization described by the Stokes parameters I (total intensity or flux), Q &amp; U (linear polarization), and V (circular polarization). In solar observations, all Stokes parameters are measured quasi-simultaneously, but such a sequence of measurements is significantly longer in stellar observations due to flux limitations. Also, since linear polarization caused by magnetic fields is typically smaller than circular polarization, measuring Stokes Q &amp; U is often neglected.

Due to these challenges and a lack of high spectral resolution stellar spectropolarimeters, many previous magnetic field studies were based on an analysis of Stokes I only until recently. The basis of such an analysis is a comparison of Zeeman broadened lines with large and small Land´ e factors, that is, magnetically sensitive and insensitive line profiles in observed spectra (Robinson 1980). A similar technique was applied by Basri et al. (1992) who analyzed the change in line equivalent widths to measure magnetic fields on cool stars. Reiners &amp; Basri (2006) applied this approach to M dwarf spectra of atoms and molecules.

Detecting Stokes V in stellar spectra is considered to be an unambiguous detection of a magnetic field, in the absence of possible systematic errors. To increase the detection sensitivity, various techniques averaging Stokes profiles of thousands of atomic lines were developed (e.g., Donati et al. 1997; Sennhauser et al. 2009). The disadvantage of such global filtering is that the spectral information is lost and only an average atmosphere can be diagnosed. Semel (1989) and Semel et al. (1993) studied Stokes V atomic spectra and proposed to deconvolve their time-series into surface magnetic field maps (Zeeman-Doppler Imaging, ZDI). This technique has been extensively employed in recent years (e.g., Donati &amp; Landstreet 2009; Donati 2003b; Morin et al. 2010).

Using only Stokes I or only Stokes V as well as average atomic spectral lines for measuring and mapping stellar magnetic fields have large limitations (see, e.g., Berdyugina 2009). In particular, these approaches do not fully describe magnetic fields in starspots, whose internal physical parameters remain unaccessible.

While the atomic Zeeman e ff ect can probe warmer regions on A to mid-M type stars, the use of molecular lines is of advantage for studying cool objects starspots in the op-/ tical and infrared wavelength regions, where atomic lines diminish their utility among the forest of molecular features. The power of molecular spectroscopy and spectropolarimetry for studying stellar magnetic fields has been demonstrated during the last two decades. Thorough theoretical and observational studies of magnetic properties of diatomic molecules radicals /

observed in spectra of the Sun and cool stars they have been established as unique probes of physical parameters prevalent in the quiet Sun, sunspots and starspots, including temperature and magnetic field strength (Berdyugina et al. 2000, 2003; Berdyugina 2011). The sensitivity of di ff erent molecules (such as MgH, TiO, FeH, and CaH) for probing magnetic fields in starspots on F, G, K, and M stars was recently demonstrated by Afram &amp; Berdyugina (2015). These developments have led to the first direct detections of magnetic fields in starspots on M dwarfs (Berdyugina et al. 2006b) and on a brown dwarf (Berdyugina et al. 2017; Kuzmychov et al. 2017). Theoretical developments and astrophysical applications have encouraged studies of magnetic properties of diatomic molecules in physics laboratory experiments, such as TiO (Virgo et al. 2005), CrH (Chen et al. 2007), MgH (Zhang &amp; Steimle 2014), etc.

In this paper, we analyze Stokes I and Stokes V together, and a number of atomic and molecular lines simultaneously. Our goal is to shed light on the complexity of the magnetic fields observed on M dwarfs. Here, when we speak of the complexity, we particularly mean how magnetic fields (their strengths and filling factors) are distributed on the surface (local horizontal plane) and within di ff erent layers of the stellar atmosphere (local vertical direction). This 3D structure can be unraveled by using molecular and atomic species with di ff erent formation depths and temperature sensitivities (Berdyugina 2011).

To achieve this goal, we determined the e ff ective (photosphere and spot) temperature, the magnetic field strength, and its filling factors for a sample of M dwarfs with spectral types ranging from M1 to M7. Following Afram et al. (2008), we carried out a χ 2 minimization of discrepancies between observed and modeled intensity and circular polarization spectra featuring seven di ff erent molecular and atomic bands or lines that are magnetically sensitive. The sensitivity of the considered molecules for probing magnetic fields in starspots was studied by Afram &amp; Berdyugina (2015). The right choice of magnetically sensitive lines is very important for the magnetic field determination because the inclusion of many marginally magnetically sensitive lines in any analysis only dilutes the result.

The paper is organized as follows: In Sect. 2, a description of the observational data is given, and in Sect. 3, the employed molecular and atomic data for FeH, CaH, MgH, TiO, Fe I, and Ti I lines are described. Section 4 presents our method: we have modeled the observed line profiles employing our numerical codes based on the full Stokes parameter radiative transfer in atomic and molecular lines in the presence of magnetic fields (Berdyugina et al. 2003). In Sect.4.2, we explain each step of our procedure. An analysis of sensitivities of our results is given in Sect. 4.3. We discuss our results for each parameter separately in Sect. 5, and summarize our conclusions in Sect. 6.

## 2. Observational data

The data used throughout this paper are described in Table 1 and contain a sample of relatively slowly rotating earlyto late-M dwarfs. The calibrated data are taken from the online archive of observations with the spectropolarimeter ESPaDOnS(Donati 2003a) at the 3.6 m Canada-France-Hawaii

Table 1. Red dwarf data from the online archive of the Canadian Astronomy Data Centre. Some of the data used in this paper were previously published by Berdyugina et al. (2006b); Morin et al. (2010); Fouqu´ e et al. (2018).

| Name            | Spectral Type   |   vsin i [ km / s ] | S / N (at 800nm)   | Telescope / Instrument   | Obs. Date UT   | P.I. Name   |
|-----------------|-----------------|---------------------|--------------------|--------------------------|----------------|-------------|
| AU Mic          | M1              |                   7 | ∼ 920              | CFHT / ESPaDOnS          | 2005-07-16     | Berdyugina  |
| GJ 360          | M2              |                   5 | ∼ 380              | CFHT / ESPaDOnS          | 2014-11-11     | Malo        |
| AD Leo          | M3.5            |                   3 | ∼ 770              | CFHT / ESPaDOnS          | 2007-06-25     | Harrington  |
| EV Lac (Gl 873) | M3.5            |                   5 | ∼ 800              | CFHT / ESPaDOnS          | 2005-07-17     | Berdyugina  |
| YZ Cmi          | M4.5            |                   7 | ∼ 330              | CFHT / ESPaDOnS          | 2006-02-07     | Forveille   |
| GJ 1224         | M4.5            |                   0 | ∼ 200              | CFHT / ESPaDOnS          | 2008-06-25     | Morin       |
| GJ 1245B        | M5.5            |                   7 | ∼ 180              | CFHT / ESPaDOnS          | 2007-10-03     | Wade        |
| CN Leo (Gl 406) | M5.5            |                   3 | ∼ 210              | CFHT / ESPaDOnS          | 2008-03-22     | Morin       |
| VB 8            | M7.0            |                   5 | ∼ 100              | CFHT / ESPaDOnS          | 2009-05-04     | Morin       |

Telescope (CFHT). The resolving power for CFHT ESPaDOnS / is approximately 68,000, and the observed wavelength region covers the range from 3700 to 10480 Å, including many molecular (FeH, CaH, MgH, and TiO) and atomic (Fe and Ti) lines.

Stars were chosen to evenly cover the spectral type range from M1 to M7. Properties for the selected targets, such as the ages and activity levels, can be found in the literature (e.g., Gizis et al. 2002; Schmitt &amp; Liefke 2004; Silvestri et al. 2005; Moutou et al. 2017). If several measurements were taken for a given star, we have chosen the one with the strongest polarimetric signal. We analyze here only one phase at which the signal is maximal, that is, a snapshot of the star's rotational phase. The exposure times for the observed spectra are short compared to the rotational periods of the stars with values of about 0.7%-6.0% of the rotational periods (except for VB 8 where the exposure time is ∼ 38% of the rotational period). A next step will be an analysis of the rotational dependence of parameters. The basic information on the selected data is provided in Table 1. These data were previously analyzed and presented in several papers (see references provided in the table caption).

Mand L dwarfs, but they are less blended in contrast to other molecular features. In addition, many lines in these bands become noticeably broadened due to Zeeman splitting already at field strengths of 2 -3kG. Zeeman broadened FeH lines in an active M dwarf were detected by Valenti et al. (2001). They employed the sunspot measurements by Wallace et al. (1998) to model the stellar spectrum, since a theoretical description of the FeH molecule was not available at that time. Berdyugina et al. (2001, 2003) modeled synthetic Stokes profiles of FeH lines and showed the usefulness of the FeH F 4 ∆ -X 4 ∆ system around 1 µ m for diagnosing solar and stellar magnetic fields, once the spin-coupling constants are available.

## 3. Molecular and atomic data

The considered molecular species were previously described in Afram &amp; Berdyugina (2015), but for the sake of completeness, we briefly present them in this section again, with a more detailed overview of FeH, as it is a very sensitive diagnostic for M dwarfs. We also include magnetically sensitive atomic lines, to determine surface temperatures and to measure magnetic fields.

## 3.1. FeH at 8700 and 9900-10000 Å

The molecule FeH is one of the most sensitive indicators of magnetic fields in cool stellar atmospheres. Here, we consider two di ff erent FeH bands of the Wing-Ford system F 4 ∆ -X 4 ∆ : the (1,0) band near 8700 Å and (0,0) band around 1 µ m. The magnetic sensitivity of the 1 µ mband was clearly demonstrated in the sunspot atlas of Wallace et al. (1998). A great advantage of using the Wing-Ford bands of FeH as a Zeeman diagnostic is the fact that numerous lines of the same species are visible in a relatively narrow wavelength region throughout spectra of

The FeH F 4 ∆ -X 4 ∆ system is produced by transitions between two electronic quartet states with Ω = 7 2 , 5 2 , 3 2 , and 1 2 , Ω being the quantum number for the component of the total electronic angular momentum along the internuclear axis of a diatomic molecule. The coupling of the angular momenta is intermediate between limiting Hund's cases (a) and (b) (see Herzberg 1950). The necessary molecular constants, which describe how strongly the di ff erent internal angular momenta are coupled to each other, were unknown until provided by Dulick et al. (2003), allowing further analyses of this particular system. Starting from the available Hamiltonian, Afram et al. (2007, 2008) calculated a perturbed molecular Zeeman e ff ect, computed Land´ e factors of the energy levels and transitions, and calibrated molecular constants to account for an unknown perturbation by comparing observed and modeled polarized sunspot spectra. Polarized spectra are a ff ected more strongly by internal perturbations than intensity spectra. Hence, di ff erences in the Land´ e factors can be better seen in polarization which makes this method more sensitive. In addition, polarization is important for determining the sign of the e ff ective Land´ e factor, which can easily be negative in molecular transitions, as shown by Berdyugina &amp; Solanki (2002). The progress in the theoretical description of this FeH system finally allowed Afram et al. (2009) to start applying this theory for studying magnetic fields on M dwarfs. In this work, we fully exploit its diagnostic capabilities.

The magnetic sensitivity varies strongly among the FeH lines of di ff erent rotational branches and total angular momentum J -numbers within the observed wavelength range. Transitions within P-, Q-, and R-branches have ∆ J = -1 0, ,

and + 1, respectively. Our primary diagnostic was the highly sensitive FeH Q 3 (2.5) line at 9975.476 Å with the e ff ective Land´ e factor g e ff = 0.96 (Afram et al. 2008). Its Stokes I profile changes quickly as the magnetic field strength increases from 1 to 5 kG at T e ff , 2 = 3500K (see Afram et al. 2009). The synthetic profile develops a characteristic, almost rectangular shape until it completely splits into Zeeman components. This particular line is therefore excellent for determining magnetic fields from intensity spectra. The R 1 (2.5) and R 1 (6.5) lines at 9976.397 Å and 9978.721 Å (with g e ff of 0.18 and 0.19, respectively) react only moderately to incremental changes of 1 kG in the magnetic field strength and are almost insensitive to small variations of about 0.1 kG (Afram et al. 2009). A joint analysis of such lines with high and low magnetic sensitivity helps to disentangle the magnetic field strength and its filling factor. However, we emphasize that including too many marginally magnetically sensitive lines into an analysis may significantly dilute the result and increase uncertainties.

Previously, Reiners &amp; Basri (2006) investigated Zeeman broadening in FeH lines of early M dwarfs. They developed a technique to estimate the mean magnetic field from a comparison to reference spectra of early M dwarfs with magnetic field measurements calculated in atomic lines. Reiners &amp; Basri (2007) employed this method for a sample of 22 M stars providing mean magnetic fields (averaged over the visible surface) for M dwarfs of spectral types M2-M9 with the accuracy of ± 1 kG. Wende et al. (2011) identified magnetically sensitive FeH lines in high resolution near-infrared CRIRES spectra through comparison of an active and inactive M dwarf. Laboratory measurements of the Zeeman response of WingFord FeH lines were carried out by Crozet et al. (2012, 2014). Then, they used the newly determined Land´ e factors to deduce the magnetic field in sunspots from Stokes V profiles. Shulyak et al. (2010, 2011b,a, 2014) used FeH to carry out modeling and analysis of magnetic fields in selected M dwarfs with a method similar to Reiners &amp; Basri (2007).

## 3.2. TiO at 7054 Å

Spectra of sunspots and cool stars show many TiO lines. The first astronomical detection of TiO was recorded in spectra of M dwarfs (Fowler 1904) and first observations of TiO molecular bands from starspots was reported by Vogt (1979) for a K2 star. Using TiO bands for analyzing starspot properties was suggested by Ramsey &amp; Nations (1980), which allowed the spot area and temperature to be measured (e.g., Ne ff et al. 1995; O'Neal et al. 1996). Valenti et al. (1998) tested the spectral synthesis of TiO lines for reproducing the optical spectrum of an inactive M dwarf. The first spectropolarimetric models and observations in sunspots were obtained by Berdyugina et al. (2000). Berdyugina et al. (2003) concluded that the TiO γ (0 , 0) R 3 system band head at 7054Å is one of the best molecular diagnostics of the magnetic field in sunspot umbra and starspots. First circular polarization signals in TiO lines were observed and modeled in three M dwarfs (Berdyugina et al. 2006b).

The triplet states of this system are in an intermediate Hund's case (a-b) and require accounting for perturbations for precise measurements (Berdyugina &amp; Solanki 2002). This is the approach we follow in this paper. Inferences of starspot temperature and filling factor on five highly active stars were shown by O'Neal et al. (2004), by fitting TiO bands using spectra of inactive G and K stars to represent the unspotted photospheres of the active stars and spectra of M stars to represent the spots.

## 3.3. MgH at 5200 Å

The molecular Zeeman splitting for MgH was theoretically predicted by Kronig (1928) and Hill (1929), while laboratory measurements of the Zeeman e ff ect were analyzed in Crawford (1934). Laborde (1961) described MgH spectral lines and compared laboratory wavelengths with solar ones. MgH lines were employed for determining the surface gravity (log g ) of Arcturus Bell et al. (1985) and other cool stars (Berdyugina &amp; Savanov 1992). Berdyugina et al. (2000) reported the first spectropolarimetric measurements and model of MgH lines in sunspots.

The A 2 Π state of the MgH A 2 Π -X 2 Σ + (0,0) transition is an intermediate Hund case (a-b) with spin-orbit constants of Y = 5 7, i.e., approaching a pure Hund case (b). Th X . 2 Σ + state is pure case (b). Both states require proper perturbation calculations of the Paschen-Back e ff ect. With this, the MgH A 2 Π -X 2 Σ + system at 5200 Å is a sensitive tool for studying stellar magnetic (Berdyugina et al. 2000, 2005).

## 3.4. CaH at 6938 Å

Calcium monohydride (CaH) is an important astrophysical molecule, and its bands have been used as indicators of cool stars luminosities (Oehman 1934; Mould 1976; Mould &amp; Wallis 1977; Barbuy et al. 1993). Because CaH absorption is an important opacity source in brown dwarfs (e.g., Burrows et al. 2001; Kirkpatrick 2005), CaH bands are valuable for studying brown dwarfs too. The CaH A 2 Π -X 2 Σ + band system is observed in the wavelength region at 6600-7600Å.

The ground state X 2 Σ + is described by a pure Hund case (b), while the excited state A 2 Π is intermediate between Hund cases (a) and (b). In addition, the Paschen-Back e ff ect must be considered for CaH (Berdyugina et al. 2003). Berdyugina et al. (2006a) obtained the first polarimetric observations and synthesis of CaH Stokes profiles and confirmed the diagnostic value of CaH in cool astrophysical sources.

## 3.5. Atomic lines

A few Zeeman sensitive atomic lines are also included in this study: two Fe I lines at 8468.4 Å (e ff ective Land´ e factor g = 2.50) and 8514.1Å ( g = 1.83) and three Ti I lines at 8364.2 Å ( g = 1.43), 8377.9 Å ( g = 0.88), 8382.5 Å ( g = 1.25). The two last Ti I lines are referred throughout the paper as 8380 Å lines. These lines were chosen because they cover a range of Land´ e factors, that helps disentangling the magnetic field strength

9

7

9

7

9

7

.

0

9

7

9

7

.

5

n

9

7

9

7

7

.

0

9

7

9

7

.

5

Fig. 1. Upper plot: Response of the highly magnetically sensitive FeH Q 3 (2.5) line at 9975.476 Å and the marginally sensitive FeH R 1 (2.5) and R 1 (6.5) lines at 9976.397 Å and 9978.721 Å, respectively, to changes of 200 K in the e ff ective temperatures for the range of 2600 K to 3400 K. The line profiles become deeper with cooler temperatures. A magnetic field strength of 3000 G was applied and instrumental (0.3 Å) broadening. Lower plot: The same for changes of log g from 4.0 to 5.0 in steps of 0.5 and an e ff ective temperature of 3000 K.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00076v1-figure-d74d7a33-0456-4de8-9ff9-bb8465e9f8e2.jpg)

from the area filled by that field (filling factor). Also, a number of weak blending molecular lines were included into the line list for computing synthetic spectra.

## 4. Method

Here, we present our method with which we determined e ff ective temperatures of the photosphere and of magnetic features, magnetic field strengths and filling factors for nine M dwarfs (M1-M7). Our χ 2 analysis is based on a comparison of observed and synthetic intensity and circular polarization profiles (Stokes I and V) of many magnetically sensitive atomic and molecular lines.

## 4.1. Modeling Stokes profiles

We calculated synthetic molecular and atomic intensity and polarization profiles in the presence of magnetic fields using the code STOPRO (Solanki 1987; Solanki et al. 1992; Frutiger et al. 2000; Berdyugina et al. 2003, 2005). This code solves the polarized radiative transfer equations for all four Stokes parameters simultaneously under the assumption of the local thermodynamic equilibrium (LTE). The equations are solved for model atmospheres from Allard et al. (2001) assuming solar element abundances. Both perturbed Zeeman and Paschen-Back e ff ects were accounted for molecular transitions. Here, we use 1D models, which in the future can be replaced by 3D MHD models of starspots which are not yet available; even for sunspots, such models are still not fully realistic.

We took into account the following line-broadening e ff ects: intrinsic, thermal (according to the local temperature in the model atmosphere), magnetic, microturbulent (1 km s), rota-/ tional, and instrumental (0.15 Å). For stars with the projected rotational velocity v sin i ≥ 5km s we considered / v sin i values given by Reiners &amp; Basri (2007). For stars with v sin i &lt; 5 km s / we adjusted v sin i values to obtain best fits. Uncertainties in v sin i &lt; 5 result in negligible errors. For example, for the magnetic field strength, the error is smaller than 0.1 kG.

For our analysis, we selected several spectral regions where various atomic and molecular transitions were observed. We divided them into three sets according to their temperature sensitivities:

- 'Atoms' depicts all considered atomic lines, i.e., Fe I at 8468Å and 8514Å, and Ti I at 8364 Å and 8380 Å.
- 'Mol1' depicts the molecular lines of FeH at 9900Å and 9975Å and TiO at 7055Å. In some cases, this set was further split into two subsets: one with the FeH lines only, the other only with the TiO lines, because TiO and FeH lines form at different heights in early M dwarf atmospheres.

'Mol2' depicts the molecular lines of CaH at 6938 Å, MgH at 5200Å, and FeH at 8700Å.

The observed intensity (Stokes I) and circular polarization (Stokes V) profiles were used to determine parameters of M dwarf atmospheres. We employed a model with two components: a warm photosphere and cool spots, both being allowed to be magnetic. Cool magnetic regions are referred to as starspots, while warm magnetic regions are considered to be analogues to the solar photosphere network fields. This model can be described by five parameters as follows:

T phot , photosphere temperature,

B , magnetic field strength,

T mag , magnetic region temperature,

- f I , magnetic field filling factor determined from Stokes I,
- f V, magnetic field filling factor determined from Stokes V.

Here, a filling factor f is a fraction of the visible stellar surface area covered by regions with a magnetic field strength B and temperature T mag . The fraction 1 -f is the stellar surface with the photosphere temperature T phot that is free from cool starspots or other magnetic regions. Because circular polarization of opposite signs for complex, spatially unresolved magnetic fields cancels out, we distinguish between filling factors determined from Stokes I and V. Accordingly, B refers to a module (unsigned) of the magnetic field | B | , when determined from only Stokes I, and to a residual (signed) lineof-sight (LOS) component B LOS, when determined from only Stokes V. When both Stokes I and V are fitted simultaneously (like in our model for most species), B refers to the magnetic field module (unsigned) along the LOS.

We characterize the magnetic field complexity (entanglement) as the di ff erence between the filling factors determined from the intensity and polarization (here, Stokes V):

This di ff erence quantifies the fraction of the magnetic area on the star that is highly entangled and remains hidden if only Stokes V is analyzed.

With this model, we have searched for solutions with and without cool magnetic spots as well as with and without magnetic network fields. As described in the Sections 4.2 and 5, various combinations of these features have been found for different atomic and molecular species. For instance, atomic and mol2 lines have been found to be suitable to probe warmer magnetic regions (network) on M dwarfs, while mol1 lines are sensitive to cooler regions (starspots).

We have also tested a three-component model, assuming the existence of a non-magnetic photosphere, magnetic areas in the photosphere (network) and magnetic cool spots (starspots) simultaneously. Results obtained with such a three-component model and the two-component model described above were very similar.

Since metal hydrides are sensitive to both stellar temperature and luminosity (see Sect. 3), we used FeH lines for evaluating how the model atmosphere e ff ective temperature T e ff and surface gravity log g a ff ect the line profiles (Fig. 1). We found a stronger response to T e ff changes for cooler stars. For log g , the di ff erence between 4.0 and 4.5 is negligible, but it is somewhat larger between 4.5 and 5.0, i.e., also for lower-mass, cooler stars. Considering uncertainties in spectral classification of our targets, we fixed log g = 4 5 and then determined e . ff ective temperatures. In this case, a log g uncertainty of 0.1 dex would lead to errors in T e ff of ∼ 30K at an e ff ective temperature of 3000 K (see also discussion of our results in Sect. 5).

## 4.2. Parameter fitting procedure

We determined the model parameters introduced in the previous section for the observed stars using a χ 2 minimization analysis by fitting the spectral range for each considered set.

Since the model parameters are convolved, we used an iterative procedure. First, we attempted to formally fit all selected sets of atomic and molecular species simultaneously to obtain all the parameters at once. However, we found that the selected parameter set is insu ffi cient for describing all the data, since the parameters refer to a single layer in the atmosphere and do not account for its complexity in height. In contrast, our data carry this height dependent information and, thus, cannot be fitted with parameters describing only one layer within the atmosphere. In addition, magnetic field parameters determined from polarization are influenced by the field entanglement via cancellations. In principle, if the model is complex enough, fitting all parameters at once should be possible, especially when including profile variability due to stellar rotation, i.e., similar to ZDI but in 3D. This is however beyond the scope of the present paper.

Therefore, we adapted our procedure. We have fitted various combinations of the model parameters at the same time using di ff erent sets of data, thus, attempting to probe both the entire atmosphere and its di ff erent layers separately. We arrived at the most optimal procedure for our analysis, a summary of which is given in Table 2.

We first determined T phot assuming B = 0 from the Ti I 8380Å Stokes I two line profiles. These lines were found to indicate systematically higher temperatures than other species. Therefore, it was fixed as a reference temperature for a nonmagnetic photosphere in the subsequent analysis. Since in reality it fits a mix of the photosphere and spot contributions, we actually estimate a mean e ff ective temperature, which could be 100-200K lower than the true photospheric temperature. On the other hand, network magnetic regions could be hotter than the photosphere (like on the Sun) which may help counteract the temperature bias. After fixing T phot , we iteratively determined the remaining model parameters from Stokes I and V profiles of all species as presented in Table 2. We have found that atomic Stokes I and V profiles can be well fitted with a model assuming no cool spots but allowing for network magnetic regions. These model results are listed in Table 3 under the entry 'atoms'.

Examples of best fits to the observed Stokes profiles are shown in Fig. 2 and Fig. 3: the first shows fits for FeH lines only but in all stars, and the second shows fits to spectra for all atomic and molecular species but only for the star AD Leo. For illustration purpose only, the best-fit synthetic spectra were calculated for each atom or molecule individually (since the analysis of all species together does not necessarily yield the best fit for the individual species), and with all the five model parameters obtained at once. The fit in Fig. 2 was obtained including only FeH lines without other blends. Hence, it gives an impression of the dominance of the FeH F 4 ∆ -X 4 ∆ system throughout the extracted wavelength range of about 10 Å. Poorer fits in Stokes V result from lower signal to noise, low or very complex signals, and cancellations in the circular polarization.

## 4.3. Analysis of sensitivities

The sensitivities of the derived parameters were estimated from the average signal-to-noise ratio (S N) level of the spectral re-/ gion as listed in the online data archive files. The sensitivities correspond to the internal error of each parameter, when all other parameters are fixed. The 1 σ errors of all the parameters (photosphere and spot temperatures, filling factors in intensity and circular polarization and field strength) have been estimated from the unnormalized ∆ χ 2 level for the respective number of degrees of freedom. Sensitivities are given in Table 3. They do not include systematic errors due to missing blends and imperfections of the fits. Those systematic errors are di ffi -cult to determine because of the limited number of the parameters considered for describing complex spectra of M dwarfs.

## 5. Results

Main results following the procedure described in Sect.4.2 comprise temperatures of the photosphere and magnetic features, magnetic field strengths, and magnetic filling factors in intensity and circular polarization for all nine stars and the atomic (atoms) and molecular (mol1 and mol2) sets. They are presented in Table 3, visualized in Fig.4, and discussed here.

Table 2. Summary of iterative χ 2 minimization procedure used in this paper. The sets of spectral features (atoms, mol1, mol2) are described in Sect. 4.1. The results are presented in Table 3 and shown in Fig. 4

.

| Species     | Used Stokes   | Assumed                    | Retrieved   | Range, step                    | Probed atmosphere features                                                         |
|-------------|---------------|----------------------------|-------------|--------------------------------|------------------------------------------------------------------------------------|
| Ti I 8380 Å | I             | B = 0                      | T phot (Ti) | 2000K-4000K, 100K              | An average photosphere with the highest temperature.                               |
| mol1        | I             | T phot (Ti) T mag ≤ T phot | T mag f I   | 2000K-4000K, 100K 0.0-1.0, 0.1 | Starspots: TiO and FeH lines were fitted separately (form at di ff erent heights). |
| mol2        | I             | T phot (Ti) T mag ≤ T phot | T mag f I   | 2000K-4000K, 100K 0.0-1.0, 0.1 | Magnetic network.                                                                  |
| atoms       | I, V          | T phot = T mag             | B f V       | 0G-7000G, 100G 0.0-0.2, 0.01   | Magnetic network.                                                                  |
| mol1        | I, V          | T phot (Ti) T mag f I      | B f V       | 0G-7000G, 100G 0.0-0.2, 0.01   | Starspots.                                                                         |
| mol2        | I             | T phot (Ti) T mag f I      | B           | 0G-7000G, 100G                 | Magnetic network. The Stokes V signals are too weak and not used.                  |

## 5.1. The photosphere temperature

In fact, we have not found any lines that could be formed in a non-magnetized atmosphere on these M dwarfs.

As noted above, the photosphere temperature was evaluated from the intensity profiles of the two Ti I 8380 Å lines as the highest inferred value among all species. We conclude that these lines are least contaminated by cool spots (which are approximately 2200 K cooler than the photosphere), while they can still be contaminated by hot areas (which are about 200 K hotter). The evaluated photosphere temperature was fixed for the analysis of other line sets. It changes from 3800 K for an M1 dwarf to 3100K for an M7 dwarf. These are somewhat higher than average e ff ective temperatures (photosphere + cool spots) of M dwarfs determined from observed colors or spectral fitting without accounting for cool spots (3700 K for M1 to 2700K for an M7 (e.g., Rajpurohit et al. 2013). This discrepancy is explained by our finding that spots on M dwarfs are significantly cooler than the photosphere and can cover a large surface area (see Sect. 5.2 and 5.4).

Interestingly, when we fitted T phot using all selected atomic lines, we found even higher values for the three coolest M dwarfs (Fig. 4a and third column in Table 3). Among other possible explanations, this may indicate that network regions are significantly hotter than the photosphere on the coolest dwarfs, while we assumed they are of the same temperature as the photosphere. When we attempted to include such features into our analysis (three-component model), the results were redundant. After T phot was determined as described here, its value was fixed to obtain the other four parameters of our model.

## 5.2. The temperature of magnetic regions

Fitting Stokes I and V of the selected Fe I and Ti I lines using four free model parameters ( T phot was fixed to T phot(Ti) ) resulted in T mag = T phot and a strong magnetic field. Thus, these lines probe a warm, highly magnetized photosphere, i.e., network.

When fitting Stokes I spectra of the CaH, MgH, and FeH 8700Å lines (set mol2; their Stokes V spectra were not included because of low S N) using the three model parameters / (excluding f V), we found that the resulting T mag matched well T phot obtained with the Ti I 8380 Å lines (Fig.4b) for all but two stars. This indicates that the mol2 set lines also form preferentially in the warm photosphere on M dwarfs, in contrast to G and K dwarfs where these lines are excellent diagnostics of cool starspots (Afram &amp; Berdyugina 2015). Hence, the mol2 set CaH, MgH and FeH 8700Å lines can be used to determine T phot of M dwarfs and, therefore, provide an independent justification for our choice of T phot . Also, as is the case for atomic lines, they reveal that M dwarf photospheres are highly magnetized (see Sect. 5.3).

Fitting Stokes I and V of the mol1 set TiO and FeH lines using the four free model parameters revealed for most stars T mag between 2000-2300K with an average of 2100K, which is approximately constant across the M1 to M7 spectral classes (Fig.4b). These cool spots are also highly magnetic (see Sect. 5.3). Here, the di ff erence between T phot and T mag reduces on average from 1700K for an M1 dwarf to about 1000 K for an M7 dwarf. Interestingly, this range of spot temperature di ff erences is typical for sunspots of various sizes (e.g., Berdyugina 2005). For a young, active M8.5 brown dwarf, Berdyugina et al. (2017) and Kuzmychov et al. (2017) determined T phot = 2800K and T mag = 2200K using also atomic and molecular lines. This result agrees with the found trend and extends our finding to ultra-cool dwarfs.

The spot temperature di ff erence with respect to the photosphere determined here for M dwarfs di ff ers from that found from broad-band photometric measurements, as summarized by Berdyugina (2005). This is probably due to a heightdependence of spot parameters as well as the higher sensitivity

Fig. 2. This illustrates the high temperature and magnetic sensitivity of FeH lines which is reproduced by our models. Observed (black) and synthetic (red) intensity and circular polarization profiles for FeH at 9897-9913 Å in the considered sample of M dwarfs. The synthetic profiles were calculated using two components: one non-magnetic ('photospheric') and one magnetic ('spot') atmosphere model (Allard et al. 2001) with best-fit parameters obtained from fits of Stokes I and V spectra for FeH in the considered wavelength region only. All the following parameters were obtained simultaneously in one χ 2 minimization: the photosphere and spot temperatures, the magnetic field strength, the spot filling factor for the intensity f I , and the circular polarization f V.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00076v1-figure-58149e93-94ae-487c-99d1-0d37a53f4f59.jpg)

of molecular lines to the temperature than that of photometry (see Sect. 6).

## 5.3. The magnetic field strength

In summary, thanks to the variety of our spectral diagnostics, we found two kinds of magnetic regions with clearly different temperatures in M dwarf atmospheres as a function of the spectral class: one is cold (starspots), and the other is warm (network).

The magnetic field strength B was determined during the same fitting procedure as described for T mag in Sect. 5.2. It corresponds therefore to an average field strength of a given magnetic structure at a given layer of the atmosphere, as defined by the corresponding species.

The B values found from atomic and molecular lines are confined to the range of 3-6kG (there are only three outliers), with an average of about 5 kG (Fig. 4d). This average is practically constant throughout the considered spec-

Fig. 3. Observed (black) and synthetic (red) intensity and circular polarization profiles for AD Leo in all considered wavelength regions. The synthetic profiles were calculated using two components: one non-magnetic ('photospheric') and one magnetic ('spot') atmosphere model (Allard et al. 2001). Fits to Stokes I and V spectra for each atom and molecule were obtained individually (see Sect. 4.2), but all the parameters were obtained simultaneously: the photosphere and spot temperatures, the magnetic field strength, the spot filling factor for the intensity f I , and the circular polarization f V. In three regions Stokes V was not fitted because of lower signal to noise in this particular wavelength region and observation, while in other regions the signal is stronger, despite the fact that it looks like noise because the region contains many lines.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00076v1-figure-27dc727f-ef99-486f-8686-102f5077bfaf.jpg)

tral class range of M1-M7 dwarfs. Berdyugina et al. (2017) and Kuzmychov et al. (2017) similarly found B ∼ 5 kG in their study of a young, active M8.5 brown dwarf.

## 5.4. The magnetic field filling factors

In summary, a cooler starspot atmosphere probed by the mol1 set and a hotter network atmosphere probed by the atoms and mol2 set are both magnetic, indicating highly magnetized atmospheres on these stars. The constant trend for the magnetic field implies that magnetic structures probably become locally smaller in scale toward lower mass M dwarfs.

The magnetic field filling factors f I and f V were determined during the same fitting procedure as described for T mag in Sect. 5.2. The filling factor f V has been determined for the atom and mol1 sets (Fig. 4e), while for the mol2 set, the circular polarization signals were too weak when compared to the noise.

For the atom and mol2 sets, the photosphere and magnetic region temperatures were found to be the same, except for one case (see Sect. 5.2). Thus, to distinguish the non-magnetic photosphere and network, we analyzed atomic Stokes I and V si-

Table 3. Best fit parameters (including sensitivities as determined in Sect. 4.3) for all considered stars obtained from the χ 2 analyses. For molecules (sets mol1 and mol2), T phot is fixed to T phot (Ti). For atoms (set atoms), T phot = T mag .

| Star sensitivity   | Species   | T phot [K] ± 50 K   | T mag [K] ± 50 K   | B [kG] ± 0.5 kG   | f I ± 0.1   | f V ± 0.02   | /epsilon1 IV   |
|--------------------|-----------|---------------------|--------------------|-------------------|-------------|--------------|----------------|
| AU Mic             | Ti I      | 3700                | -                  | -                 | -           | -            | -              |
| AU Mic             | atoms     | 3800                | -                  | 4.3               | 0.8         | 0.01         | 0.79           |
| AU Mic             | mol1      | -                   | 3700               | 3.7               | 0.1         | 0.01         | 0.09           |
| AU Mic             | mol2      | -                   | 3700               | 3.2               | 1.0         | -            | -              |
| GJ 360             | Ti I      | 3800                | -                  | -                 | -           | -            | -              |
| GJ 360             | atoms     | 3800                | -                  | 5.1               | 0.2         | 0.01         | 0.19           |
| GJ 360             | mol1      | -                   | 2000               | 4.5               | 0.7         | 0.01         | 0.69           |
| GJ 360             | mol2      | -                   | 3800               | 5.0               | 1.0         | -            | -              |
| AD Leo             | Ti I      | 3700                | -                  | -                 | -           | -            | -              |
| AD Leo             | atoms     | 3700                | -                  | 3.5               | 0.9         | 0.05         | 0.85           |
| AD Leo             | mol1      | -                   | 2200               | 5.6               | 0.7         | 0.05         | 0.65           |
| AD Leo             | mol2      | -                   | 3700               | 4.3               | 0.4         | -            | -              |
| EV Lac             | Ti I      | 3600                | -                  | -                 | -           | -            | -              |
| EV Lac             | atoms     | 3600                | -                  | 5.1               | 0.7         | 0.03         | 0.67           |
| EV Lac             | mol1      | -                   | 2000               | 5.6               | 0.7         | 0.04         | 0.66           |
| EV Lac             | mol2      | -                   | 3600               | 6.8               | 0.8         | -            | -              |
| YZ Cmi             | Ti I      | 3400                | -                  | -                 | -           | -            | -              |
| YZ Cmi             | atoms     | 3400                | -                  | 4.3               | 1.0         | 0.01         | 0.99           |
| YZ Cmi             | mol1      | -                   | 2300               | 4.4               | 0.6         | 0.01         | 0.59           |
| YZ Cmi             | mol2      | -                   | 3400               | 5.0               | 1.0         | -            | -              |
| GJ 1224            | Ti I      | 3500                | -                  | -                 | -           | -            | -              |
| GJ 1224            | atoms     | 3500                | -                  | 4.2               | 0.8         | 0.08         | 0.72           |
| GJ 1224            | mol1      | -                   | 2100               | 4.5               | 0.8         | 0.20         | 0.60           |
| GJ 1224            | mol2      | -                   | 2000               | 5.0               | 0.7         | -            | -              |
| GJ 1245B           | Ti I      | 3200                | -                  | -                 | -           | -            | -              |
| GJ 1245B           | atoms     | 3400                | -                  | 4.3               | 1.0         | 0.02         | 0.98           |
| GJ 1245B           | mol1      | -                   | 2300               | 5.7               | 0.7         | 0.02         | 0.68           |
| GJ 1245B           | mol2      | -                   | 3200               | 5.0               | 1.0         | -            | -              |
| CN Leo             | Ti I      | 3100                | -                  | -                 | -           | -            | -              |
| CN Leo             | atoms     | 3400                | -                  | 4.5               | 1.0         | 0.08         | 0.92           |
| CN Leo             | mol1      | -                   | 2000               | 5.0               | 0.7         | 0.05         | 0.65           |
| CN Leo             | mol2      | -                   | 3100               | 5.0               | 1.0         | -            | -              |
| VB 8               | Ti I      | 3100                | -                  | -                 | -           | -            | -              |
| VB 8               | atoms     | 4000                | -                  | 7.0               | 1.0         | 0.20         | 0.80           |
| VB 8               | mol1      | -                   | 3100               | 1.0               | 0.6         | 0.01         | 0.59           |
| VB 8               | mol2      | -                   | 3100               | 5.0               | 1.0         | -            | -              |

multaneously. This resulted in very high f I (between 0.7 and 1.0, except for one star) and very low f V (0.01 to 0.08, with only one star having 0.2). This large di ff erence between f I and f V indicates a strong entanglement of the network magnetic field in M dwarfs. This also explains why the mol2 set lines have very low Stokes V signals.

For the mol1 set, we analyzed results for FeH and TiO separately. These molecules reveal cool magnetic starspots (see Sect. 5.2) with f I up to 0.8. Figure 4c shows that the FeH f I increases from 0.3 to 0.8 for earlier to later M dwarfs, while the TiO f I remains on average constant at about 0.7 (except for one star). We computed linear fits with an outlier-resistant two-variable linear regression for the FeH and TiO values separately for T e ff ≤ 3500 K and for T e ff ≥ 3500 K. The parameters for TiO and for FeH (for T e ff ≤ 3500 K) are: f I TiO , = (0.00008 ± 0.000088) ∗ T e ff + (0.41 ± 0.30) and f I FeH , , = (0.0001 ± 0.00026) ∗ T e ff + (0.36 ± 0.86). The pa- rameters for FeH (for T e ff ≥ 3500 K) are: f I FeH , , = (-0.001 ± 0.0002) ∗ T e ff + (5.5 ± 0.9). The errors of the linear fits show that the di ff erence between the trends is significant for stars with T e ff ≥ 3500 K. At the same time, f V is in the range of 0.01-0.05 for both molecules. This indicates that cool starspots in M dwarfs may be relatively small structures (perhaps similar to magnetic pores on the Sun) with large surface coverage and highly intermittent polarity.

The contrasting behavior between the mol1 FeH and TiO f I can be explained by di ff erences in their formation heights. Figure 5 shows averages of the TiO and FeH line core contribution functions for the atmosphere models with T e ff = 3200K and 3800K. Contribution functions for other temperatures are not plotted, since Figure 4c shows that the di ff erences in filling factors depend on the photosphere temperature. One can see that the formation height of FeH peaks deeper and that of TiO higher in the atmosphere for both temperatures. Contribution

Fig. 4. Temperature and magnetic field best-fit parameters (see Table 3) for di ff erent magnetic regions on M dwarfs which are probed by five groups of spectral species. T mag and f I are obtained from fitting Stokes I profiles only; B is from fits of Stokes I and V together; f V is from fits to Stokes V only; for the atoms, T phot = T mag (see text for details). To avoid overlaps, some symbols are slightly shifted in the horizontal direction. Uncertainties of the parameters are indicated as thin-line crosses in each panel. The linear fits in panel 4c are described in Sect. 5.4.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00076v1-figure-f55a50ab-94db-4427-91e7-74341aca1d58.jpg)

functions of the considered atomic line cores peak at intermediate heights. Hence, the TiO f I probes a higher atmosphere where the magnetic fields fill about 0.7 of it in all M dwarfs, while the FeH f I probes a lower atmosphere where the magnetic field filling factor increases from 0.3 in M1 to 0.7 in M5M7 dwarfs.

## 5.5. The magnetic field complexity

netic field of 5 kG, which we have found in M1-M7 dwarfs, β &lt; 1 throughout the entire atmosphere of all M dwarfs (above the optical depth 1 at 1.2 µ m). This indicates that plasma in M dwarf atmospheres is structured by the magnetic field. On the Sun, this is only observed in the solar chromosphere and corona as well as in sunspot umbra, while in the photosphere outside sunspots plasma motion is controlled by thermodynamic processes.

We have found that in early M dwarfs (M1-M3) the cool spot area is smaller in the lower atmosphere (0.3) and larger in the higher atmosphere (0.8). This is consistent with magnetic spot models expanding from a lower to upper atmosphere (e.g., Rempel et al. 2009). Using the model atmospheres, we have computed the so-called plasma β parameter, which is the ratio of the gas pressure to the magnetic pressure. The magnetic field strength producing the pressure equal to the ambient gas pressure (i.e., when β = 1) is sometimes called the equipartition magnetic field, B eq (not to be confused with the definition arising from the equipartition theorem). For the average mag-

In late M dwarfs (M5-M7), the cool spot area is large (on average 0.7) in both the lower and higher atmosphere. This may indicate that the magnetic flux emergence is larger in these stars, and strong magnetic fields fill gaps in the lower atmosphere. Thus, the entire atmosphere of a late M dwarf is more complex and active. More frequent and more powerful flares observed in late M dwarfs (Yang et al. 2017; Mondrik et al. 2019) support this interpretation of our results.

A sketch of di ff erent magnetic regions in M dwarf atmospheres including spots expanding toward higher layers as well as small-scale, network magnetic fields in the lower atmospheres is given in Figure 6. It can be seen that atomic and

Fig. 5. Contribution functions for selected atomic and molecular species for two atmosphere models with e ff ective temperatures 3200K and 3800K versus the optical depth at 1.2 µ m. Maxima indicate e ff ective formation heights of the lines.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00076v1-figure-63029f05-a473-4c7f-8da8-a8586724a051.jpg)

molecular lines are formed at di ff erent optical depths which allows us to probe changing magnetic strengths and structure size at di ff erent heights.

The complexity /epsilon1 IV of the magnetic field quantity was found to vary with the height in the atmosphere. Also, the height-dependence itself depends on the stellar mass. The largest complexity is seen in network fields, with atomic f I = 0 7 . -1 0 (except for one star) and . f V = 0 01 . -0 2, leading to . about 70-99% of hidden network fields. Cool spots probed by molecular lines are less entangled than the network, implying that, locally, they are larger-scale structures than network fields, similar to the Sun. The complexity range in spots is about 1070%. Furthermore, we have found that the complexity tends to increase toward lower mass M dwarfs.

In contrast, ZDI maps based on only atomic Stokes V (e.g., Donati 2003b; Morin et al. 2010) reveal a very small part of the stellar magnetic flux. Thus, interpretation of this residual flux as representative of the global magnetic field in terms of di ff erent components (toroidal and poloidal fields, and the degree of axisymmetry of the poloidal field) is biased due to the underestimated true complexity of the field. A discussion of tomographic field topology reconstructions and their limitations is given by Berdyugina (2009) and Kochukhov et al. (2017).

Another definition of the complexity of magnetic fields was used by Shulyak et al. (2014): as the minimum number of magnetic field components required to fit the observed Stokes I line profiles. These approaches, based only on one Stokes parameter, cannot fully describe the complexity of magnetic fields on M dwarfs (and other stars). As we have found, the complexity varies in three dimensions, depending on the height in the atmosphere.

## 6. Summary

Low-mass stars with fully convective interiors exhibit magnetic phenomena similar to those observed on the Sun, although they lack an interface layer where dynamo processes are thought to take place. Thus, exploring magnetic fields and their topology and complexity provides a key tool for understanding stellar dynamos. In this work, we have determined temperatures of the photosphere and magnetic regions as well as mean magnetic field strengths and their filling factors for nine M dwarfs (M1M7). We employed observed intensity and circular polarizaion spectra of selected atomic and molecular lines and a χ 2 analysis based on model atmospheres and polarized radiative transfer.

We have found that all studied M dwarfs are highly magnetic with the average field of ∼ 5kG covering up to 100% of the star, similar to the conclusions by Johns-Krull &amp; Valenti (2000). In particular, two kinds of magnetic regions have been distinguished: one is significantly cooler than the photosphere (starspots) and has a filling factor of about 60-80% in the upper atmosphere, while the other one is as hot as the photosphere (network) and has a filling factor 70-100% in the lower atmosphere. By analyzing the filling factors determined from intensity and circular polarization spectra, we have determined the complexity of stellar magnetic fields as the di ff erence between these two filling factors. It was found to increase toward lowermass stars.

In addition to 'surface' inhomogeneities, we have also identified di ff erences with the height in the atmosphere depending on the analyzed atomic or molecular species. In M dwarfs, FeH forms deeper and TiO higher in the atmosphere, while atomic lines form at intermediate heights. This has allowed us to establish that cool magnetic starspots in the lower atmosphere are smaller and less complex in earlier M dwarfs than in later M dwarfs, while in the upper atmosphere they are large and highly entangled in all M dwarfs. The network magnetic fields are more complex and entangled than starspots and fill practically the entire atmosphere of all M dwarfs.

More specifically, our findings are as follows:

The unspotted photosphere temperature ranges from 3800K for an M1 dwarf to 3100K for an M7 dwarf. These are somewhat higher than average e ff ective temperatures (photosphere + cool spots) of M dwarfs determined from observed colors or spectral fitting without accounting for cool spots.

Two kinds of magnetic regions with clearly di ff erent temperatures exist in M dwarf atmospheres: cold starspots and warm network (similar to the Sun). Their temperatures vary with the spectral class of stars. The network temperature is the same (within uncertainties) as the unspotted photosphere.

Fig. 6. A sketch of di ff erent magnetic regions in M dwarf atmospheres. Starspots expand toward higher layers and have a large filling factor in the upper atmosphere. Small-scale magnetic fields (network) dominate the lower atmospheres, also with a high filling factor. The geometrical scale given in axis on the left (km) belongs to the field lines. Atomic and molecular lines are formed at di ff erent optical depths (emphasized by di ff erent colors; the optical depth scale is given in the axis on the right) which allows us to probe changing magnetic strengths and structure size at di ff erent heights. The atmosphere parameters used for the sketch: e ff ective temperature of 3500K and magnetic field strength of 5 kG. Inside the magnetic spots, the temperature is about 2200K, which is not shown here.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00076v1-figure-0ffd2f13-eeb0-4e0f-a0e3-3b5714e28623.jpg)

The magnetic field strength is on average 5 kG in both starspot and network features, and this does not vary across the M1-M7 spectral class range. We conclude that this is an implication that magnetic structures become locally smaller in scale toward lower mass M dwarfs. Also, the plasma β parameter is smaller than 1 throughout the entire atmosphere of all Mdwarfs (above the optical depth 1 at 1.2 µ m), indicating that the plasma in M dwarf atmospheres is structured by the magnetic field, like in the solar chromosphere, corona and sunspot umbra.

This work is a major step toward a 3D reconstruction of magnetic fields in M dwarfs and a comparison with realistic 3D MHD models of sunspots and starspots. Analysing all Stokes IQUV parameters of both atomic and molecular lines in the next step of such a study is necessary for fully understanding magnetic structures on cool stars.

Network fields on M dwarfs are of mixed polarity and highly entangled (possibly small-scale, as on the Sun). They fill almost the entire lower and intermediate height M dwarf atmosphere. Polarimetric measurements alone may underestimate their magnetic flux (the complexity is about 70-99%).

Cool starspots on M dwarfs are also highly intermittent in polarity and appear as intermediate-scale features, similar to pores on the Sun. They are larger in a higher M dwarf atmosphere, where they can fill up to 80% surface area. Because of the mixed polarity, more than half of these features may remain undetected when only polarimetric data are analyzed (the complexity is about 10-70%).

The magnetic field complexity on M dwarfs (cf., entanglement) varies with the height in their atmospheres and depends on the spectral class. There is an indication that the magnetic flux emergence is larger in later M dwarfs (lower-mass), where highly complex, strong magnetic fields fill the atmosphere in all dimensions.

Acknowledgements. This work was supported by the European Research Council (ERC) Advanced Grant HotMol (ERC-2011-AdG 291659). SB thanks the Kavli Institute for Theoretical Physics (KITP), Santa Barbara, CA, for support during the program 'Better Stars, Better Planets' as well as program's organizers and participants for fruitful discussions about stellar magnetism. We thank the anonymous referee for the detailed review which helped to improve our manuscript.

## References

Afram, N. &amp; Berdyugina, S. V. 2015, A&amp;A, 576, 34 Afram, N., Berdyugina, S. V., Fluri, D. M., et al. 2007, A&amp;A,

- 473, L1

Afram, N., Berdyugina, S. V., Fluri, D. M., Solanki, S. K., &amp; Lagg, A. 2008, A&amp;A, 482, 387

Afram, N., Berdyugina, S. V., &amp; Reiners, A. 2009, in 'ASP Conf. Ser.', Vol. 405, Solar Polarization 5: In Honor of Jan Stenflo, ed. S. V . Berdyugina, K. N. Nagendra, &amp; R. Ramelli, 527

Allard, F., Hauschildt, P. H., Alexander, D. R., Tamanai, A., &amp; Schweitzer, A. 2001, ApJ, 556, 357

Anglada-Escud´ e, G., Amado, P. J., Barnes, J., et al. 2016, Nature, 536, 437 Barbuy, B., Schiavon, R. P., Gregorio-Hetem, J., Singh, P. D., &amp;Batalha, C. 1993, A&amp;A Suppl. Ser., 101, 409 Basri, G., Marcy, G. W., &amp; Valenti, J. A. 1992, ApJ, 390, 622 Bell, R. A., Edvardsson, B., &amp; B., G. 1985, Mon. Not. R. Astron. Soc., 212, 497 Berdyugina, S. V. 2005, Living Rev. Solar Phys., 2, no. 8 Berdyugina, S. V. 2009, in IAU Symposium, Vol. 259, Cosmic Magnetic Fields: From Planets, to Stars and Galaxies, ed. K. G. Strassmeier, A. G. Kosovichev, &amp; J. E. Beckman, 323332 Berdyugina, S. V. 2011, in Astronomical Society of the Pacific Conference Series, Vol. 437, Solar Polarization 6, ed. J. R. Kuhn, D. M. Harrington, H. Lin, S. V. Berdyugina, J. Trujillo-Bueno, S. L. Keil, &amp; T. Rimmele, 219 Berdyugina, S. V., Braun, P. A., Fluri, D. M., &amp; Solanki, S. K. 2005, A&amp;A, 444, 947 Berdyugina, S. V., Fluri, D. M., Ramelli, R., et al. 2006a, ApJ, 649, L49 Berdyugina, S. V., Frutiger, C., Solanki, S. K., &amp; Livingston, W. 2000, A&amp;A, 364, L101 Berdyugina, S. V., Harrington, D. M., Kuzmychov, O., et al. 2017, ApJ, 847, 61 Berdyugina, S. V., Petit, P., Fluri, D. M., Afram, N., &amp; Arnaud. 2006b, in ASP Conf. Ser., Vol. 358, Solar Polarization 4, ed. R. Casini &amp; B. W. Lites, 381 Berdyugina, S. V. &amp; Savanov, I. S. 1992, Astron. Zh., 69, 829 Berdyugina, S. V. &amp; Solanki, S. K. 2002, A&amp;A, 385, 701 Berdyugina, S. V., Solanki, S. K., &amp; Frutiger, C. 2001, in 'ASP Conf. Ser.', Vol. 248, 'Magnetic Fields Across the Hertzsprung-Russell Diagram', ed. G. Mathys, S. K. Solanki, &amp; D. T. Wickramasinghe, 9 Berdyugina, S. V., Solanki, S. K., &amp; Frutiger, C. 2003, A&amp;A, 412, 513 Burrows, A., Hubbard, W. B., Lunine, J. I., &amp; Liebert, J. 2001, Reviews of Modern Physics, 73, 719 Chabrier, G. &amp; Bara ff e, I. 1997, A&amp;A, 327, 1039 Chen, J., Bakker, J. M., Peters, A., et al. 2007, Phys. Chem. Chem. Phys., 9, 949 Crawford, F. H. 1934, Rev. Mod. Phys., 6, 90 Crozet, P., Ross, A., Alleq, N., et al. 2014, in Magnetic Fields throughout Stellar Evolution, Proceedings of the International Astronomical Union, ed. P. Petit, M. Jardine, &amp;H. C. Spruit, Vol. 302, 164 Crozet, P., Tourasse, G., Ross, A., Paletou, F., &amp; L´pez Ariste, o A. 2012, in ECLA - European Conference on Laboratory Astrophysics, EAS Publications Series, ed. C. Stehl´ e, C. Joblin, &amp; L. d'Hendecourt, Vol. 58, 63 Donati, J.-F. 2003a, in Astronomical Society of the Pacific Conference Series, Vol. 307, Solar Polarization, ed. J. Trujillo-Bueno &amp; J. Sanchez Almeida, 41 Donati, J.-F. 2003b, in EAS Publications Series, Vol. 9, EAS Publications Series, ed. J. Arnaud &amp; N. Meunier, 169 Donati, J.-F. &amp; Landstreet, J. 2009, Ann. Rev. A&amp;A, 47, 333 Donati, J.-F., Semel, M., Carter, B. D., Rees, D. E., &amp; Collier Cameron, A. 1997, MNRAS, 291, 658 Dulick, M., Bauschlicher Jr., C. W., Burrows, A., et al. 2003,

ApJ, 594, 651 Fouqu´ e, P., Moutou, C., Malo, L., et al. 2018, MNRAS, 475, 1960 Fowler, A. 1904, Proc. R. Soc. London, Ser. A, 73, 219 Frutiger, C., Solanki, S. K., Fligge, M., &amp; Bruls, J. H. M. J. 2000, A&amp;A, 358, 1109 Gillon, M., Triaud, A. H. M. J., Demory, B.-O., et al. 2017, Nature, 542, 456 Gizis, J. E., Reid, I. N., &amp; Hawley, S. L. 2002, AJ, 123, 3356 Herzberg, G. 1950, 'Molecular Spectra and Molecular Structure, I. Spectra of Diatomic Molecules', 2nd edn. (New York: Van Nostrand Company) Hill, E. L. 1929, Phys. Rep., 34, 1507 Johns-Krull, C. M. &amp; 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, &amp; S. Sciortino, 371 Kirkpatrick, J. D. 2005, ARA&amp;A, 43, 195 Kochukhov, O., Petit, P., Strassmeier, K. G., et al. 2017, Astronomische Nachrichten, 338, 428 Kronig, R. d. L. 1928, Phys. Rep., 31, 195 Kuzmychov, O., Berdyugina, S. V., &amp; Harrington, D. 2017, ApJ, 847, 60 Laborde, G. 1961, Ann. Astrophys., 24, 89 Mondrik, N., Newton, E., Charbonneau, D., &amp; Irwin, J. 2019, ApJ, 870, 10 Morin, J., Donati, J.-F., Petit, P., et al. 2010, Monthly Notices of the Royal Astronomical Society, 407, 2269 Mould, J. R. 1976, ApJ, 207, 535 Mould, J. R. &amp; Wallis, R. E. 1977, Mon. Not. R. Astron. Soc., 181, 625 Moutou, C., H´ ebrard, E. M., Morin, J., et al. 2017, MNRAS, 472, 4563 Ne ff , J. E., O'Neal, D., &amp; Saar, S. H. 1995, ApJ, 452, 879 Oehman, Y. 1934, ApJ, 80, 171 O'Neal, D., Ne ff , J. E., Saar, S. H., &amp; Cuntz, M. 2004, Astron. J., 128, 1802 O'Neal, D., Saar, S. H., &amp; Ne ff , J. E. 1996, ApJ, 463, 766 Rajpurohit, A. S., Reyl´ e, C., Allard, F., et al. 2013, A&amp;A, 556, A15 Ramsey, L. W. &amp; Nations, H. L. 1980, ApJ, 239, L121 Reiners, A. 2012, Living Rev. Solar Phys., 8, no. 1 Reiners, A. &amp; Basri, G. 2006, ApJ, 644, 497 Reiners, A. &amp; Basri, G. 2007, ApJ, 656, 1121 Rempel, M., Sch¨ ussler, M., &amp; Kn¨lker, M. 2009, ApJ, 691, 640 o Robinson, R. D. 1980, ApJ, 239, 961 Schmitt, J. H. M. M. &amp; Liefke, C. 2004, A&amp;A, 417, 651 Semel, M. 1989, A&amp;A, 225, 456 Semel, M., Donati, J.-F., &amp; Rees, D. E. 1993, A&amp;A, 278, 231 Sennhauser, C., Berdyugina, S. V., &amp; Fluri, D. M. 2009, A&amp;A, 507, 1711 Shulyak, D., Reiners, A., Seemann, U., &amp; Kochukhov, O. Piskunov, N. 2014, A&amp;A, 563, 35 Shulyak, D., Reiners, A., Wende, S., Kochukhov, O. Piskunov, N., &amp; Seifahrt, A. 2010, A&amp;A, 523, 37 Shulyak, D., Reiners, A., Wende, S., Kochukhov, O. Piskunov,

N., &amp; Seifahrt, A. 2011a, in 'ASP Conf. Ser.', Vol. 448, '16th Cambridge Workshop on Cool Stars, Stellar Systems,

and the Sun', ed. C. M. Johns-Krull, M. K. Browning, &amp;

A. A. West, 1263

- Shulyak, D., Seifahrt, A., Reiners, A., &amp; Kochukhov, O. Piskunov, N. 2011b, Mon. Not. R. Astron. Soc., 418, 2548

Silvestri, N. M., Hawley, S. L., &amp; Oswalt, T. D. 2005, AJ, 129, 2428

Solanki, S. K. 1987, PhD thesis, ETH, Zurich, Switzerland Solanki, S. K., R¨edi, I., &amp; Livingston, W. 1992, A&amp;A, 263, u 312

Strassmeier, K. G. 2009, A&amp;A, 17, 251

Strassmeier, K. G. 2011, Proc. of the IAU, 273, 174

Valenti, J. A., Johns-Krull, C. M., &amp; Piskunov, N. E. 2001, in ASP Conf. Ser., Vol. 223, 11th Cool Stars, Stellar Systems and the Sun, ed. R. J. Garc´ ıa L´pez, R. Rebolo, &amp; M. R. o Zapatero Osorio, 1579

- Valenti, J. A., Piskunov, N. E., &amp; Johns-Krull, C. M. 1998, ApJ, 498, 851
- Virgo, W. L., Steimle, T. C., &amp; Brown, J. M. 2005, ApJ, 628, 567

Vogt, S. S. 1979, Publ. Astron. Soc. Pac., 91, 616

Wallace, L., Livingston, W. C., Bernath, P. F., &amp; Ram, R. S. 1998, An Atlas of the Sunspot Umbral Spectrum in the Red and Infrared from 8900 to 15050 cm -1 (6642 to 11230 Å) (NOAO)

Wende, S., Reiners, A., Seifahrt, A., Shulyak, D., &amp; Kochukhov, O. 2011, in 'ASP Conf. Ser.', Vol. 448, '16th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun', ed. C. M. Johns-Krull, M. K. Browning, &amp; A. A. West, 1309

West, A. A., Hawley, S. L., Bochanski, J. J., et al. 2008, Astron.

J., 135, 785

Yang, H., Liu, J., Gao, Q., et al. 2017, ApJ, 849, 36 Zhang, R. &amp; Steimle, T. C. 2014, ApJ, 781, 51
N. Afram S. V. Berdyugina
Complexity of magnetic fields on red dwarfs
Magnetic fields in cool stars can be investigated by measuring Zeeman line
broadening and polarization in atomic and molecular lines. Similar to the Sun,
these fields are complex and height-dependent. Many molecular lines dominating
M-dwarf spectra (e.g., FeH, CaH, MgH, and TiO) are temperature -- and Zeeman --
sensitive and form at different atmospheric heights, which makes them excellent
probes of magnetic fields on M dwarfs. Our goal is to analyze the complexity of
magnetic fields in M dwarfs. We investigate how magnetic fields vary with the
stellar temperature and how "surface" inhomogeneities are distributed in height
-- the dimension that is usually neglected in stellar magnetic studies. We have
determined effective temperatures of the photosphere and of magnetic features,
magnetic field strengths and filling factors for nine M dwarfs (M1-M7). Our
chi^2 analysis is based on a comparison of observed and synthetic intensity and
circular polarization profiles. Stokes profiles were calculated by solving
polarized radiative transfer equations. Properties of magnetic structures
depend on the analyzed atomic or molecular species and their formation heights.
Two types of magnetic features similar to those on the Sun have been found: a
cooler (starspots) and a hotter (network) one. The magnetic field strength in
both starspots and network is within 3 kG to 6 kG, on average it is 5 kG. These
fields occupy a large fraction of M dwarf atmospheres at all heights, up to
100%. The plasma beta is less than one, implying highly magnetized stars. A
combination of molecular and atomic species and a simultaneous analysis of
intensity and circular polarization spectra have allowed us to better decipher
the complexity of magnetic fields on M dwarfs, including their dependence on
the atmospheric height. This work provides an opportunity to investigate a
larger sample of M dwarfs and L-type brown dwarfs.