From stellar coronae to gyrochronology: A theoretical and observational exploration

Context. Stellar spin down is the result of a complex process involving rotation, dynamo, wind, and magnetism. Multiwavelength surveys of solar-like stars have revealed the likely existence of relationships between their rotation, X-ray luminosity, mass losses, and magnetism. They impose strong constraints on the corona and wind of cool stars. Aims. We aim to provide power-law prescriptions of the mass loss of stars, of their magnetic ﬁeld, and of their base coronal density and temperature that are compatible with their observationally-constrained spin down. Methods. We link the magnetic ﬁeld and the mass-loss rate from a wind torque formulation, which is in agreement with the distribution of stellar rotation periods in open clusters and the Skumanich law. Given a wind model and an expression of the X-ray luminosity from radiative losses, we constrained the coronal properties by assuming di ﬀ erent physical scenarios linking closed loops to coronal holes. Results. We ﬁnd that the magnetic ﬁeld and the mass loss are involved in a one-to-one correspondence that is constrained from spin down considerations. We show that a magnetic ﬁeld, depending on both the Rossby number and the stellar mass, is required to keep a consistent spin down model. The estimates of the magnetic ﬁeld and the mass-loss rate obtained from our formalism are consistent with statistical studies as well as individual observations and they give new leads to constrain the magnetic ﬁeld-rotation relation. The set of scaling-laws we derived can be broadly applied to cool stars from the pre-main sequence to the end of the main sequence (MS), and they allow for stellar wind modeling that is consistent with all of the observational constraints available to date.


Introduction
The rotation of stars is subject to complex evolution during their whole life.During the early stage of their lifetime, solar-type stars spin up as they contract during the pre-main sequence.Once the zero-age main sequence (ZAMS) is reached, their moment of inertia remains relatively constant (Armitage & Clarke 1996) while they lose mass and angular momentum through the flow of a magnetized stellar wind (Schatzman 1962;Weber & Davis 1967;Mestel 1968).It results in a slow-down of their rotation as they age, which approximately follows the empirical Skumanich's law: Ω ∝ t −0.5 (Skumanich 1972).This spin down also depends on the stellar mass (Weber & Davis 1967;Matt et al. 2015).Indeed, during most of the pre-main sequence, lower mass stars tend to remain fast rotators for a longer time than the higher mass stars.Then, after hundreds of millions of years, the slowest rotators converge toward a sequence in which the rotation rate increases with mass.These phenomena make gyrochronology possible (Barnes 2003), thereby allowing for the estimation of stellar ages through measurements of rotation periods and masses.
Understanding the feedback loop between rotation, dynamo action, magnetism, and wind is a key issue in predicting the behavior of the solar-like stars as they evolve.It is also important to understand the evolution of star-planet systems and to follow potential planetary migrations (Zhang & Penev 2014;Benbakoura et al. 2019).Our ability to track the stellar rotation evolution strongly relies on the wind braking modeling.
Most angular momentum evolution models rely on the Kawaler (1988) prescription when assessing the wind torque, which is expressed in this model as a power-law depending on the magnetic field, the mass-loss rate, the mass, and the radius of the star.Several modifications have since been added to this formulation (Krishnamurthi et al. 1997;Bouvier et al. 1997;Reiners & Mohanty 2012).For instance, to account for fast rotators on the ZAMS, a saturation of the braking torque is required (Barnes & Sofia 1996).Magnetohydrodynamical (MHD) simulations can also be used to assess of the angular momentum loss due to the wind.For instance, Matt et al. (2012), following Matt & Pudritz (2008), simulated the flow of a stellar wind along the opened field lines of a dipolar configuration to estimate the torque by taking the influence of stellar rotation on the wind acceleration into account (Sakurai 1985).More recently, modified versions of this formulation have been presented to take the influence of the magnetic topology into account.Réville et al. (2015a) considered the magnetic flux through the open field lines to build a topology-independant wind torque, while Finley & Matt (2017, 2018) relied on a broken powerlaw behavior to deal with combined geometries.Garraffo et al. (2016Garraffo et al. ( , 2018) ) accounted for magnetic topology by means of a modulating factor to the angular momentum loss estimated for a dipolar configuration.In general, in most of the prescriptions of spin-evolution torques, it is assumed that the wind carries away angular momentum at a rate proportional to Ω 3 during the main sequence in order to follow the Skumanich law.It is important to note that the Skumanich law is currently questioned for evolved stars.Some recent studies have shown a substantial decrease of the wind braking efficiency for evolved stars around the solar value of the Rossby number (van Saders et al. 2016) even if this alternative scenario seems to be in disagreement with solar twin studies (Lorenzo-Oliveira et al. 2019).
To estimate the angular momentum loss, we need to know the precise stellar magnetic field and the accurate properties of its wind (such as the mass loss induced).The latter is bound to physical parameters in the corona, such as the plasma temperature and density.Several observational trends constrain and correlate those different quantities: the mass loss (Wood et al. 2002(Wood et al. , 2005;;Jardine & Collier 2019), the X-ray activity (Pizzolato et al. 2003;Wright et al. 2011;Reiners et al. 2014), the magnetic field (Vidotto et al. 2014;See et al. 2017), and the rotation rate of the star for different ages (Agüeros et al. 2011;McQuillan et al. 2014;Gallet & Bouvier 2015).Coupling all those quantities is therefore necessary when designing a consistent model of stellar spin down.
A preliminary exploration was carried out by Blackman & Owen (2016), who present a simplified model for the coupled time evolution of the relevant quantities on the basis of a pressure-driven isothermal wind (Parker 1958).In this framework, a dynamo-induced magnetic field dictates the behavior of the stellar wind and the X-ray luminosity through a coronal equilibrium.More recently, Skumanich (2019) has focused on the connections between the physical parameters regulating the stellar spin down by considering their rotational evolution.By assuming a direct correlation between the mass loss and the magnetic field, he studied the influence of the Skumanich law on all of the relevant quantities by means on a variety of observational trends.Such a study suggests that the magnetic field and the rotation of a cool star should scale linearly, while the mass loss should scale quadratically with stellar rotation.
More generally, stellar wind models require knowledge on the coronal temperature and density.To be consistent with the aforementioned correlations, one needs to constrain those two quantities from tracers of the coronal activity, such as the soft X-ray emission of the star.In this spirit, by relying on a 1D polytropic and magnetocentrifugal wind, Holzwarth & Jardine (2007) provided scaling laws that are in accordance with the rotational evolution of the X-ray luminosity (Ivanova & Taam 2003) and the empirical mass loss-X-ray flux correlation from Wood et al. (2005).
In building on those previous studies, the main goal of this paper is to infer from stellar spin down considerations some prescriptions of the magnetic field, the mass loss, the coronal temperature, and the coronal density as a function of fundamental stellar parameters (such as mass, radius, and rotation rate) in order to be consistent with all of the observational trends.Each of those constraints are successively introduced to eliminate the largest number of free parameters involved in the spin down process.In our attempt to extract the most important interdependencies between the various physical mechanisms involved in the stellar spin down theory, we had to make some simplifying assumptions; however, we were careful to retain all of the key mechanisms.Our scaling laws provide a novel and systematic way to connect all of these mechanisms together and they should be seen as a first approach to systematically characterize these complex relationships between the stellar parameters.In Sect. 2 we introduce the theoretical framework we used for the torque parametrization.In Sect.3, observational constraints are leveraged to unveil the inter-dependency of the magnetic field and the mass-loss rate of cool stars.We further derived the associated prescriptions for the coronal temperature and the coronal density for a given wind model.In Sect.4, we summarize our prescriptions and give a practical application in the case of Eridani.All of those results are then summarized, discussed, and put in perspective in Sect. 5.

Stellar wind torques of solar-type stars:
Theoretical approach 2.1.Fundamental stellar parameters and architecture of the model First of all, designing a consistent model for the stellar spin down requires taking an inventory of the coupling between the various physical mechanisms, and their associated control parameters, involved in the process.Stellar rotation and magnetism strengthen the supersonic flow of a stellar wind (Weber & Davis 1967), which itself entails a mass loss of Ṁ.The wind then carries away angular momentum, leading to a braking torque Γ wind that spins down the star.The logical sequence from fundamental stellar parameters to the wind braking torque and the role of the X-ray luminosity inside the architecture of the model is summarized in Fig. 1, where we show all of the interdependencies between the various physical mechanisms and their control parameters.Theoretical assumptions are represented in blue and are essentially related to the wind torque parametrization, dynamo scalings, as well as the choice of a wind model.Observational trends, which are in red, are considered to be constraints for the physical models (see Appendix F for an extensive view of those ingredients and their caveats).Solar-type stars generate a magnetic field through a dynamo in their envelope (Brun et al. 2004(Brun et al. , 2015) ) that is triggered by turbulent convective movements, which can be influenced by stellar rotation (Durney & Latour 1978).This effect can be quantified by the Rossby number Ro, which has been shown to be a relevant quantity in characterizing the magnetic activity of the stars (Noyes et al. 1984).Direct spectropolarimetric studies also exhibit scaling laws linking the stellar magnetic field and this dimensionless number (Vidotto et al. 2014;See et al. 2017).However, multiple definitions of the Rossby number are available in the literature, which has forced the community to clarify this aspect (Landin et al. 2010;Brun et al. 2017;Amard et al. 2019).In this work, the Rossby number is normalized to the solar value.This way, different prescriptions of the Rossby number can be equally used as long as they only differ by a proportionality factor.This is generally the case of main-sequence stars, which we focus on here.However, to deal with individual stars, we used the stellar Rossby number for simplicity purposes, which is expressed as where P rot is the stellar rotation period and τ c is the convective turnover time.The latter itself can be assessed at different locations or computed in a variety of ways, which are often linked to a given stellar evolution model (Landin et al. 2010;Cranmer & Saar 2011;Sadeghi Ardestani et al. 2017) or a set of observations (Wright et al. 2011).
Coronal activity, which is measured with the X-ray emission of the star, has been correlated with the stellar magnetic field (Pevtsov et al. 2003;Vidotto et al. 2014) and the mass loss (Wood et al. 2002(Wood et al. , 2005)).Observations have also exhibited a relationship between the stellar rotation and the  1.Sketch of the multiple couplings between the inputs of the model (stellar mass, radius, and rotation rate), the magnetic field, the X-ray luminosity, the coronal temperature and density, the mass loss, and the wind braking torque.In blue: model assumptions.In red: observational constraints.
X-ray luminosity for main-sequence stars (Noyes et al. 1984;Pizzolato et al. 2003;Wright et al. 2011;Vidotto et al. 2014;Reiners et al. 2014), allowing us to link the coronal activity to fundamental stellar parameters.This gives us the opportunity to closely link those physical quantities together through stellar activity considerations.
Furthermore, the corona of the stars that we consider is thought to be heated up through magnetic processes (for a review, see Mathioudakis et al. 2013;Cranmer et al. 2017) to reach a typical coronal temperature T c of around 10 6 K (Suzuki & Inutsuka 2006).Such a hot and diluted medium (with a density n c ) expands and generates a transsonic wind (Parker 1958).Observationally, scaling laws have been established between the coronal temperature and the X-ray emission of the star (Preibisch 1997;Johnstone & Güdel 2015;Wood et al. 2018).From a theoretical point of view, the latter can be deduced from coronal properties by estimating the radiative losses in the X-ray emitting area of the corona (Ivanova & Taam 2003;Aschwanden 2004;Blackman & Owen 2016).Therefore, it could be possible to provide a prescription for the coronal temperature and density in terms of the stellar mass, radius, and rotation rate.

Torque parametrization
Solar-like stars spin down due to the angular momentum extraction by the stellar wind (Schatzman 1962).Weber & Davis (1967) showed with a 1D model in the stellar equatorial plane that the angular momentum loss Γ wind can be estimated at the Alfvén radius r A as in which Ω is the stellar rotation rate and Ṁ is the mass loss.
The Alfvén radius can be expressed as a function of the stellar magnetic field, the mass loss, and other fundamental stellar parameters, such as the stellar mass, the stellar radius, and the stellar rotation rate (Kawaler 1988).For instance, Matt et al. (2012) present the following expression for this characteristic distance by assuming a dipolar magnetic field: where M is the stellar mass, R is the stellar radius, and K 1 , K 2 , and m are constants, which were set using 2D MHD simulations.More precisely, K 1 is used to calibrate the solar wind torque, K 2 is the efficiency of the magnetocentrifugal acceleration (Sakurai 1985), and m stands for a magnetic topology parameter (Réville et al. 2015a).For a dipolar field, we use the value proposed in Matt et al. (2012), that is, m = 0.2177.We note that f = Ω / GM /R 3 is the break-up ratio, which is obtained by dividing the stellar rotation rate at the equator of the star by the Keplerian angular velocity.The magnetization parameter Υ (Matt & Pudritz 2008) is defined as where B is the magnetic field strength at the stellar equator and v esc = √ 2GM /R is the escape velocity.With this formulation, the wind braking torque becomes A170, page 3 of 18 A&A 635, A170 (2020) We rewrite this power-law expansion of the torque to adopt the following generic expression where K = √ 2K 2 is a constant.

Stellar magnetic field prescription
Computing the wind braking torque requires estimating the stellar magnetic field and the mass-loss rate.Here, we focus on the surface magnetic field strength at the stellar equator B .For the sake of simplicity, we assume a power-law expression for the magnetic field as a function of the stellar mass, the stellar radius, and the Rossby number as where p B , m B , and r B are unspecified exponents for the time being.
Since we present a significant number of power-law prescriptions in the rest of this work, we have to fix a generic notation for the different exponents.More precisely, for each estimate of a physical quantity, we write a given exponent with a lowercase letter indicating the variable of the power law (p for the rotation period, r for the stellar radius, and m for the stellar mass).Each letter has a subscript in capital letters representing the quantity for which we give a prescription.For example, to express a quantity A as a function of the stellar mass M , we write (8)

Mass loss prescription
Because of the wind torque parametrization and the stellar magnetic field prescription, we consider a power-law expression of the stellar mass loss for the sake of consistency, as follows: where p Ṁ , m Ṁ , and r Ṁ are constrained later.Those exponents follow the nomenclature presented in Sect.2.3.However, the mass loss can be obtained from coronal quantities through a wind model, for example, by assuming a radial polytropic pressure-driven outflow, with an index γ (cf.Appendix A): where n , T are respectively the solar values of the density and temperature at the base of the corona.We note that T min, = (1 − 1/γ) Gm p M /2k B R ≈ 11(1 − 1/γ) MK is the minimal temperature needed at the base of the corona for the Sun to obtain a transsonic wind.For instance, for γ = 1.05, which is the value commonly used in the literature (Matt et al. 2012;Réville et al. 2015a;Finley & Matt 2017), we have T min, ≈ 0.52 MK.Such an expression is used in Sect.3.5.

General formulation of the torque
The parametrization of the torque with Eq. ( 6), together with the power-law expressions of the stellar magnetic field and the mass loss (Eqs.( 7) and ( 9)) lead to the following formulation of the torque as a function of fundamental stellar parameters: From this generic formulation, it is now possible to constrain the different exponents by several observational trends into account.

Observational constraints
3.1.Relationships between stellar parameters Some observational trends are based on a set of main-sequence stars in the unsaturated rotation regime (for which the Rossby number is greater than a certain threshold) and, therefore, they take scaling laws between stellar parameters into account, depending on a specific stellar evolution model.In particular, we need to consider mass-radius and mass-luminosity relationships during the main sequence.Here, we assume the general correlations where η and ξ are constants that depend on which stellar model is being considered.Therefore, the upcoming formulations, which are only valid during the main sequence, can accommodate any evolutionary model.As a default, we use the typical values η = 4 and ξ = 0.9 (Kippenhahn & Weigert 1994).

Constraints form stellar rotational evolution
As already introduced, stellar spin down studies in open clusters, which were first performed on the Pleiades, Ursa Major, and the Hyades, show that the rotation rate of evolved main-sequence stars tends to converge at the solar rate on a sequence where Skumanich 1972).By assuming that the moment of inertia of the star is constant, the wind braking torque is constrained in the unsaturated regime to scale as the cube of the stellar rotation rate In the following sections, we assume that gyrochonology is valid thanks to the Skumanich law.This way, if a decrease of the wind braking efficiency is genuinely occurring for evolved stars (van Saders et al. 2016), then we only consider solar-type stars that are younger than the Sun to ensure the validity of Eq. ( 14).A stalling of the magnetic braking could be modeled in our formalism with a "resaturation" regime where the Rossby number is greater than a certain threshold Ro break .Matt et al. (2015) studied the different dependencies of the wind braking torque in more detail to explain some characteristic features of the distribution of stellar rotation periods in open clusters and Kepler stars as a function of their mass.To this end, they focus on two kinds of stellar populations: the slow rotators in an unsaturated regime and the fast saturated rotators.The saturation threshold is given by a value of the Rossby number Ro sat , which is assumed to be independent of any stellar parameters, at least at zeroth-order.To explain the mass dependency of the stellar spin down, they took the Skumanich law in the unsaturated regime and a linear saturation for the wind braking torque into account, leading to the following prescription 15) with Γ = 6.3 × 10 30 erg, a = 3.1, and b = 0.5.This prescription can be compared with our formulation, in Eq. ( 11), by assuming f 1 and by neglecting secular changes of the stellar parameters.It is thus possible to link the magnetic field of the star and the mass loss through the following conditions Equation ( 17) is similar to the condition obtained by Skumanich (2019) if we use p B = 1/β and p Ṁ = α/β, according to his notation.However, no correlation between B and Ṁ is assumed here.The wind torque parametrization therefore gives us the opportunity to infer the mass loss prescription from the magnetic field prescription and reciprocally.
In Fig. 2 we illustrate the interdependencies in this first set of exponents.More precisely, the p B exponent is expressed as a function of p Ṁ from Eq. ( 17) for a dipolar field (solid black line), a quadrupolar field (dashed black line) and an octupolar field (dotted black line) in the case of an unsaturated rotation regime.The dashed blue lines correspond to magnetic field prescriptions from See et al. (2017;S17) and Vidotto et al. (2014;V14).The blue crosses represent the magnetic field and mass loss prescriptions from Johnstone et al. (2015a), Tu et al. (2015) and Sadeghi Ardestani et al. (2017), resp.J15, T15, SA17.Those scaling laws were derived to reproduce rotation rates of open clusters and are in good agreement with a square root spin down law.As already pointed out by Skumanich (2019), a wide range of exponents is admissible from the different prescriptions considered.
The magnetic topology also has a significant influence on the Ṁ − B prescriptions.Indeed, the value of the m exponent decreases with an increasing complexity of the topology, corresponding to higher-order multipoles (Réville et al. 2015a).Here, we use m = 0.15 as a quadrupolar field (see the dashed black line in Fig. 2) and m = 0.11 as an octupolar field (see the dotted black line in Fig. 2).More complex magnetic fields lead to steeper slopes in Fig. 2 and, therefore, to a less constrained magnetic field.In regards to the mass loss, it tends to be proportional to Ro −2 .Only one prescription of B and Ṁ is compatible with the Skumanich law for any magnetic topology (see the dark red dashed lines in Fig. 2) and corresponds to the one highlighted by Skumanich (2019).This configuration leads to a linear magnetic field-rotation relation and a quadratic mass loss-rotation relation.This way, such a prescription may be used as a first estimate of B and Ṁ based on spin down considerations without having to assume a particular magnetic topology.However, in the following sections, we use unspecified prescriptions for the magnetic field and the mass-loss rate in order to study the influence of additional observational constraints.It is important to note that complex fields can significantly modify the wind braking torque itself (Réville et al. 2015a;Garraffo et al. 2016).Furthermore, from those considerations, Garraffo et al. (2018) were able to reproduce the bimodal distribution of slow and fast rotators.
Since the dipole component tends to dominate the wind braking torque for mixed geometries (Finley & Matt 2018), such a topology is assumed in the following sections by default.

Constraints on the mass-loss rate
Given the one-to-one correspondence between the B and Ṁ exponents, constraints on the mass loss affect the magnetic field and reciprocally.First, the mass loss is enhanced by stellar rotation (Wood et al. 2005;Suzuki et al. 2013;Holzwarth & Jardine 2007), leading to p Ṁ ≥ 0. Therefore, we obtain Furthermore, Wood et al. (2002Wood et al. ( , 2005) ) showed a correlation between the mass loss and the X-ray stellar flux F X for unsaturated main-sequence stars, which is expressed as where w is a constant between w min ≈ 0.3 and w max ≈ 1.9, which is consistent with the obervations from Wood et al. (2005).To convert this correlation into a Ṁ(Ro, R , M ) prescription, we consider a relationship between the coronal activity and the stellar rotation for unsaturated main-sequence stars: where L X is the X-ray stellar luminosity, L is the luminosity of the star, and p L is an exponent between p L,min = 2 and p L,max = 3 (Pizzolato et al. 2003;Wright et al. 2011;Reiners et al. 2014).If we take the scaling laws from Sect.3.1 into account, it is possible to estimate the mass loss as a function of the Rossby number and the stellar mass from stellar activity considerations.Indeed, the X-ray flux becomes which leads to the following mass loss: For identification purposes, the exponents of Ṁ can be inferred from the set {η, ξ, p L , w}: To be consistent with the data from Wood et al. (2005), the most flexible constraint on the p Ṁ exponent can be obtained from Eq. ( 27) as According to Eq. ( 29), the values of p Ṁ , leading to (F X , Ṁ) outside the envelope of the data from Wood et al. (2005), define the exclusion red regions in Figs. 2 and 3. Equations ( 27) and ( 28) then give an additional constraint on the different exponents of the mass loss prescription This alone dictates the mass and radius dependency of the stellar magnetic field, which is easier to observe, through Eqs. ( 17)-( 20), such that 31) Equation ( 31) shows that the mass-radius dependency of the magnetic field can be inferred from its Rossby number dependency, as shown in Fig. 3, where the green area represents the exponents of the magnetic field prescription, which is compatible with the different constraints we considered.The negative values of ξr B + m B show that the magnetic field should decrease with the stellar mass, which seems to be in agreement with dynamo models and observations (Johns-Krull & Valenti 2000; Brun et al. 2015).
One can also notice in Fig. 3 that a stronger Rossby number dependency, for high values of p B , leads to a weaker but non negligible dependency on stellar mass through the Ṁ−F X correlation.Futhermore, by considering the different exclusion regions (see the hatched areas in Fig. 3), we can infer that a magnetic field, which is solely a function of the Rossby number (i.e., r B = m B = 0), cannot be consistent with the observed rotational dependency of L X and Ṁ as well as with the Ṁ−F X correlation.To fulfill this set of conditions, we need to add an explicit mass dependency to the stellar magnetic field.
Only a linear combination of r B and m B is constrained here because of the mass-radius relationship.However, it is possible to discriminate the mass and radius dependencies by considering the magnetic topology.

Constraints on the magnetic field
The stellar magnetic field can be constrained in a variety of ways.Spectropolarimetric studies have exhibited correlations between the large scale magnetic field and other stellar parameters.For example, in Vidotto et al. (2014), B is shown to scale as Ro −1.38±0.14 , which leads to a mass loss of p Ṁ = −1.41± 0.22.The X-ray luminosity of the star can also be considered to infer the magnetic field.Here, we assume a correlation of L X −B such that (Pevtsov et al. 2003;Vidotto et al. 2014) which leads to p B = p L /v.For instance, in assuming that p B = −1.38 ± 0.14 and p L = 2.5 ± 0.5, this leads to v = 1.81 ± 0.6, which can be compared with the value v = 1.61 ± 0.15, which was obtained in Vidotto et al. (2014) through a power-law fit.
Since the magnetic field is enhanced by the stellar rotation (Noyes et al. 1984;Brandenburg & Saar 2000;Petit et al. 2008), the exponent p B has to be nonnegative, leading to the following constraint on the mass-loss rate in the unsaturated regime At saturation, with Eq. ( 22), such a condition gives p B = 0 and p Ṁ = 0, which means a mass loss and a magnetic field are independent of the stellar rotation rate.
In the following sections, we consider lower and upper bounds for the values of p B and p Ṁ in order to be consistent with observational trends.The measured stellar magnetic field from Zeeman broadening and Zeeman-Doppler Imaging (ZDI; see Montesinos & Jordan 1993;Vidotto et al. 2014;See et al. 2017) have only exhibited linear or super-linear dependencies between the large-scale magnetic field and the Rossby number.In order to take this into account, we thus adopt p B ≥ 1.The See et al. (2017) prescription, that is, p B = 1.65, is used as an upper bound.The values of all of the exponents associated with these two scenarios are given in Table 1.It is a common knowledge that ZDI maps, representing the large-scale unsigned magnetic flux, do not give any information about the small-scale magnetic field, which might be dominant for young fast rotating stars.However, a correlation between measurements from Zeeman-Doppler imaging and Zeeman broadening (See et al. 2019) gives some confidence in the general trends found in the literature.
The bounds we assume above are greater than the exponents predicted by some scaling laws for stellar dynamos (p B = 0 for the equipartition, p B = 1/4 for the buoyancy regime, and p B = 1/2 for the magnetostrophy regime; for a review, see Augustson et al. 2017).This discrepancy likely comes from the fact that dynamo scaling laws and ZDI observations do not relate to the same magnetic field.The measured magnetic field from ZDI studies corresponds to the average unsigned photospheric flux B V , which is an estimate of the large scale magnetic field at the stellar surface.Dynamo-based scaling laws aim to estimate the stellar magnetic field over a wide spectral range.The latter can be linked to the average unsigned surface field strength, B I , obtained from Zeeman broadening, by means of a filling factor f that represents the fraction of the stellar surface which is magnetized (see Reiners 2012 for more details).By exhibiting a correlation between B V and B I , See et al. ( 2019) estimated a filling factor from a dynamo-produced magnetic field in the equipartition regime and show a strong Rossby number dependency in their f estimate, which could be a means of explaining this difference.
In a nutshell, it is possible to rely on a wind torque parametrization to provide B and Ṁ estimates that are consistent with spin down and X-ray emission constraints (cf.Eqs. ( 17)-( 20) and ( 31)), given for example the B −Ro relation (i.e., the p B exponent).In a quest to better understand the link between stellar and wind properties, we expand on our set of prescriptions to probe the coronal properties of the star by relying on a wind model.

Constraints on the coronal properties
3.5.1.Probing the coronal properties: The role of magnetic topology As seen in the previous section, the mass loss behavior can be constrained by the stellar magnetic field and the wind braking torque.This way, the knowledge of the Ṁ expression allows us to go further by inferring prescriptions for the coronal temperature and density by means of the X-ray stellar emission.Indeed, their rotational dependency has often been constrained in the literature through the high-energy activity of the star (Mestel & Spruit 1987;Ivanova & Taam 2003;Holzwarth & Jardine 2007).
However, X-ray emission and wind acceleration, at least for the fast wind, are believed to arise from different regions in the corona, that is, from dead zones and coronal holes, respectively.Therefore, here, we assume that the coronal temperature T c and the base density n c in the open-field regions are the quantities ruling the mass loss for a given wind model (cf.Eq. ( 10) for example), whereas the X-ray luminosity can be inferred from radiative losses by knowing the temperature T l and the density n l in closed-field regions (cf.Appendix B).In other words, the Ṁ prescription from the previous sections allows us to link T c and n c , while the Rossby dependency of L X , which has to be consistent with Eq. ( 24), correlates T l and n l .All of those connexions are summarized in Fig. 4. Since T l can be obtained from observations (Preibisch 1997;Johnstone & Güdel 2015;Wood et al. 2018), we need an additional constraint to relate (T c , n c ) to (T l , n l ).To this end, we consider three scenarios to couple the openfield and the closed-field regions.In scenario 1 (in red in Fig. 4), we consider a common scaling law for T c and T l (Johnstone et al. 2015b;Ó Fionnagain & Vidotto 2018), which leads to different scalings for the densities.In scenario 2 (in purple in Fig. 4), we assume an entropy equilibrium between the closed loops and the coronal holes.In scenario 3 (in blue in Fig. 4), we assume a common scaling for n c and n l (Ivanova & Taam 2003;Holzwarth & Jardine 2007;See et al. 2014), which leads to different scalings for the temperatures.
It is worth noticing that scenarios 1 and 3 lead to drastically different trends for the temperature T c .Therefore, a scenario where both the density and the temperature possess the same scaling in the open and closed field regions is unrealistic.
Such an approach allows us to deal with the local magnetic field distribution in an admittedly simplified way by distinguishing quiet open-field regions from closed loops associated with active regions.This distinction translates into different coronal temperature and density prescriptions, which are always chosen to be compatible with our derived wind massloss rate in the polytropic formalism.Of course, at the surface of the star, the plasma dynamics and heating mechanisms are much more involved than in our simplified approach (Wedemeyer-Böhm et al. 2009).Still this is a first step in characterizing the general properties of the coronae of cool stars, and we intend to consider more realistic modeling for the detailed coronal heating mechanism in the future.

Coronal temperature and density: X-ray luminosity consistency
The temperature and the density in closed loops are connected to the X-ray luminosity, which must be consistent with the prescriptions we adopt in Eq. ( 24).From radiative losses considerations, the X-ray luminosity can be expressed as (cf.Appendix B) A170, page 7 of 18 In order to standardize the different prescriptions, we consider a power-law expression for the temperature and the density in the dead zones, that is, Furthermore, for main-sequence stars in the unsaturated regime, we assume a correlation between the X-ray flux and the coronal temperature in the closed loops, such that (Preibisch 1997;Johnstone & Güdel 2015;Wood et al. 2018) with an exponent j l = 0.26 (Johnstone & Güdel 2015).Equation (25), with relationships from Sect.3.1, then result in a T l (Ro, M ) formulation This leads by identification to The rotational dependency of the density can therefore be deduced from the magnetic field prescription for a given set {p L , j l } by ensuring the consistency with Eq. ( 24): The X-ray luminosity therefore gives us the possibility to constrain the Rossby dependency of the temperature and the density in closed loops.
We now need to study the connection between the X-ray luminosity and the temperature in coronal holes, which will allow us to infer the mass-radius dependency of T c from its rotational dependency, as we did for the mass-loss rate in Sect.3.3.X-ray emission and wind acceleration have similar sources that are closely linked to the heating of the corona, which is probably due to the transport of energy from the photosphere through weakly dissipative Alfvén waves.These phenomena involve steep density gradients (Heyvaerts & Priest 1983) and nonlinear interactions between inward and outward perturbations (Velli et al. 1989), among others (for a review, see Mathioudakis et al. 2013;Cranmer et al. 2017).While this process is efficient in closed loops thanks to the magnetic topology, it requires wave reflections, for example through the parametric instabilities in coronal holes (Réville et al. 2018), thus reducing the heating efficiency in those regions.To take this behavior into account, we assume a correlation between the X-ray flux and the coronal temperature T c , such as in the closed loops case, with an unspecified exponent that is different from j l because of the possible difference of heating efficiency between open-field and closed-field regions.

Coronal temperature and density: Polytropic wind model considerations
The temperature and the density in coronal holes need to be consistent with our mass loss prescription, which is constrained in Sects.3.1-3.4.We consider the following power-law expression for the coronal temperature: Such a prescription has to fulfill the condition that Ṁ > 0, which defines a maximal value for the p T exponent (see Appendix C for more details).
In knowing the value of T c for a given wind model, it is then possible to infer the coronal density to obtain a consistent massloss rate.As we see in Sect.2.4 with the Eq. ( 10), the mass loss can be expressed as a function of stellar parameters and coronal properties.To be consistent with the power-law prescription from Eq. ( 9), the coronal density n c has to be expressed as It is important to keep in mind that this equation is valid if c s /v esc 1, where c s is the speed of sound in the stellar corona and v esc is the escape velocity at the stellar surface (cf.Appendix A).In light of this condition, high coronal temperatures could invalidate this analytical expression of n c .
If we assume that the coronal temperature varies slightly for the stellar parameters we consider, we can approximate this expression to a power law on the main sequence, such that (see details in Appendix D) with: and Therefore, for a given mass-loss rate prescription, a oneto-one correspondence between {p T , r T , m T } and {p n , r n , m n } occurs by considering a pressure-driven polytropic wind.Furthermore, we are able to infer ξr T + m T from the p T exponent thanks to a T c −F X correlation (cf.Appendix C).Given that the Rossby dependency of T l and n l is already known through Eqs. ( 39) and ( 40), respectively, we only need one additional constraint to fully determine the expression of T c and n c .Here, we chose to connect the rotational dependency of the temperature and the density in open and closed regions by means of the three scenarios we present in Sect.3.5.1.Now, we detail the implications of these scenarios in Sects.3.5.4-3.5.6.

Scenario 1: Single temperature scaling
To connect closed loops to coronal holes, we can adopt a single scaling law for the temperature (i.e., T l ∝ T c ), leading to the same X-ray flux-temperature correlation in both regions.This way, the Johnstone & Güdel (2015) prescription, coupled with the L X −Ro relation, and the results from Sect.3.1 provide a complete expression for T c .One cannot assume weak variations of the coronal temperature, which means that a power-law expression for n c may be a loose approximation in this scenario (cf.Appendix D).However, the coronal density can be directly inferred from our wind model through Eq. ( 42) in order to keep a consistent mass loss.Those considerations then allow us to estimate T c and n c by means of the following prescriptions In a saturated rotation regime, it is impossible from Eq. ( 50) to keep a constant value for both T c and n c .For simplicity purposes, we assume a single temperature at saturation, the corresponding coronal density being inferred from the Ṁ prescription.
The robustness of the single temperature scaling hypothesis can, nevertheless, be questioned.In closed loops, available estimates of T l (Preibisch 1997;Johnstone & Güdel 2015;Wood et al. 2018) fall back on an emission measure weighted average coronal temperature, based on heavy ions emission (Güdel 2007).From the wind model point of view, which is based on the modeling of a perfectly ionized hydrogen gas, the electron temperature, which is similar to the proton temperature in this context, is required to compute the mass-loss rate.Thus, a single scaling law for both temperatures may appear as a loose assumption.Indeed, if different radial profiles of temperature are observed in the solar wind for these two populations (Cranmer et al. 2017), with a higher temperature for heavy ions, their dependency on stellar parameters is still unknown.
Furthermore, if we assume T c ∝ T l , slower rotators may have a colder corona, leading to extremely strong densities to keep a consistent mass loss.For example, for a solar twin with a rotation period of 56 days, the single temperature scaling scenario leads to n c ≈ 19.58 n , which has to be compared to n l ≈ 0.48 n , according to Eq. ( 40).Coronal holes of slow rotators would be abnormally dense, which may suggest that a single scaling law for the temperature might be inconsistent with the other hypothesis of our formalism, especially the choice of our wind model.Because of those points, we consider other scenarios in what follows.

Scenario 2: Entropy equilibrium
A second possibility is to assume entropy equilibrium between closed and open field regions, allowing both density and temperature to vary simultaneously.This allows us to derive the following relation: where γ ad = c p /c v = 5/3 is the standard adiabatic exponent for an ideal gas.If we assume a power-law expression for n c as in Eq. ( 43), this balance then gives the temperatures and densities for the Rossby dependency as follows: Along with Eq. ( 44), which dictates the rotational dependency of the mass-loss rate in our wind model, one can express the p T and p n exponents as follows: The T c −F X correlation (cf.Appendix C), along with the wind model through Eqs. ( 45) and ( 46), dictate the mass-radius dependency of T c and n c as (56) Equations ( 53)-( 56) then define the power-law expressions of T c and n c for this scenario.At high rotation rates, as in Sect.3.3, we see that the linear saturation of the wind braking torque leads to a magnetic field and a mass loss that is independent of the Rossby number, that is, p B = p Ṁ = 0. Similar behavior has been observed for the X-ray activity of stars (Pizzolato et al. 2003;Wright et al. 2011).Here, we assume for simplicity that p L = 0 in the saturated regime.By relying on calculations that are similar to those presented in Sects.3.5.2and 3.5.3,we can show that p n = p T = 0 is an acceptable choice if we assume a similar behavior in closed loops.

Scenario 3: Single density scaling
Let us finally consider that the coronal density has a similar behavior in the closed-field and open-field regions.This way, two different scaling laws arise for the temperatures in closedfield or open-field regions.In this case, the coronal density n c , which launches the wind in the open-field regions, is proportional to the density in closed loops n l .Therefore, it can be expressed as a power law with a Rossby exponent given by Eq. ( 40) to ensure a consistent X-ray luminosity.Assuming a power-law expression for n c is a reasonable hypothesis if the temperature in coronal holes only presents small variations with respect to stellar parameters (cf.Appendix D).Such a supposition seems to be consistent with Suzuki et al. (2013), who performed simulations of flux tubes that were heated by Alfvén wave dissipation in coronal holes.Indeed, they predict a weak dependency of the coronal temperature on the stellar magnetic field, which can be extrapolated to more fundamental stellar parameters through a dynamo relationship (see Eq. ( 7)).We can now assess the Rossby-dependency of the coronal temperature (i.e., the p T exponent) through Eq. ( 44).All of those aspects allow us to determine the rotational dependency of T c and n c thanks to the following prescriptions As for the entropy equilibrium scenario, the mass-radius dependency of T c and n c is determined through the wind model Equations ( 57)-( 60) therefore give us the possibility to estimate the coronal temperature and density in open regions.As in Sect.3.5.5, the assumption of a coronal temperature and density, independent of the Rossby number in the saturated rotation regime, is consistent with the different observational constraints we imposed.
Figure 5 shows the T c and n c estimates for a solar twin with a rotation period of 15 days (corresponding to Ro = 0.6 in the Sadeghi Ardestani et al. 2017 prescription) in the three physical scenarios we analyzed.We set η = 4, ξ = 0.9 (Kippenhahn & Weigert 1994), γ = 1.05, and T = 1.5 MK (Réville et al. 2016).One can notice that a stronger dependency of the magnetic field on the Rossby number, that is, higher values of p B , leads to weaker values of T c and n c .However, this is apart from the coronal temperature in Scenario 1, which remains constant (see the light red band in panel A of Fig. 5).Furthermore, coronal temperatures obtained in scenarios 2 and 3 (see the blue and purple bands in panel A of Fig. 5) only vary a little around the solar value, which is consistent with the assumptions we made to express n c as a power law.Since the star we consider here is a solar twin that rotates more rapidly (i.e., with Ro/Ro < 1), the positiveness criterion of the mass-loss rate (cf.Appendix C) sets a maximal coronal temperature (see the dark gray line in panel A of Fig. 5).From this condition, one can affirm that scenarios 2 and 3 are compatible with a transsonic wind.However, in scenario 1, the condition that Ṁ > 0 is only ensured for p L ≈ 2; the coronal temperature is too high otherwise given a reasonable choice of (γ, T ).
Considering our lower and upper bounds for p B as described in Sect.3.4, single density scaling and entropy equilibrium lead to coronal temperatures close to the Holzwarth & Jardine (2007) prescription (see the black horizontal line in Fig. 5).Scenario 3 leads to temperatures that are weaker than the predicted value of Holzwarth & Jardine (2007), resulting in slightly higher densities for a given mass-loss rate (see the blue band in panel B of Fig. 5).On the contrary, marginally higher temperatures and lower densities can be observed in scenario 2 (see the purple band in panel B of Fig. 5).Scenario 3, with a T c prescription that is independent of the magnetic field and the mass-loss scaling laws, results in the highest values of the coronal temperature and the lowest values of the coronal density shown in Fig. 5 (in light red in Fig. 5).In light of those different aspects, scenarios 2 and 3 seem more likely to account for a consistent stellar spin down with a polytropic pressure-driven wind.

Scaling laws and observations of individual systems
4.1.Using our scaling laws: A practical guide

Rossby number and rotation regime
To deal with individual systems, we need to know their Rossby number to rely on our prescriptions.For the sake of simplicity, we use the stellar Rossby number here (cf.Eq. ( 1)) and more precisely the Sadeghi Ardestani et al. ( 2017) prescription for the convective turnover time.Following their approach, this characteristic time is defined as the ratio between the pressure scale height and the convective velocity that is estimated with the mixing length theory.They computed the relevant quantities at half a pressure scale height over the base of the convective zone and followed their evolution by means of the CESAM stellar evolution code (Morel & Lebreton 2008).The formulation they obtained has the advantage of being valid during the premain sequence and the main sequence for metallicities ranging from [Fe/H] = −0.5 to 0.5 by falling back on the stellar convective mass as the control parameter.Given the age of the system we consider, such an expression can be simplified for mainsequence stars in order to obtain a power law that only depends on the stellar mass and the stellar radius, which leads to the following formulation for the convective turnover time: This prescription leads to a solar value of Ro = 1.113 and a saturation value of Ro sat = 0.09.

Scaling laws: Numerical values of the exponents
In the previous sections, we have been able to constrain the magnetic field, the mass loss, the coronal temperature, and the coronal density from wind braking considerations by assuming the  Statistical ZDI studies, through the bounds of the p B exponents and prescriptions from spin down considerations, allow us to predict the admissible trends for the Ṁ−F X correlation (cf.Eq. ( 23)).This is shown in Fig. 6 as a light gray area along with the observed Ṁ ∝ R 2 F w X relationship, which was initially published by Wood et al. (2005).One can notice that the predicted values of the w exponent are compatible with the observations from Wood et al. (2005) and Jardine & Collier (2019).Furthermore, they present a bias toward active stars with low mass-loss rates, leading to weaker slopes than the Wood et al. (2005) prescription of w = 1.34 (see the gray dotted line in Fig. 6).Nevertheless, this behavior is in agreement with the results from Suzuki et al. (2013), which predicts an exponent of w = 0.82 (see the red dashed line in Fig. 6).

Normalization factors
If the exponents are fixed thanks to the previous calculations, the normalization of the scaling laws is still an open issue.By default, in the above analytical development, all of the relevant quantities have been normalized to the solar values.However, we need to take additional constraints into account for each physical parameter in a solar configuration.
Regarding the magnetic field, ZDI studies show a significant scatter in the dataset is used to exhibit correlations between the large scale magnetic field and other stellar parameters.Such a dispersion is taken into account here in the magnetic field normalization by considering the See et al. (2017) dataset for which the average dipolar field strength at their value for the solar Rossby number lies between B = 0.6 and 4 G.
for m = 0.22, we have 7.9 × 10 −15 ≤ Ṁ [M yr −1 ] ≤ 1.47 × 10 −13 .The normalization of the coronal temperature and density is determined by performing 1D simulations of a pressure-driven polytropic wind with γ = 1.05, using the starAML routine (Réville et al. 2015b).The value of T is tuned to provide an average solar wind velocity at 1 AU of 444 km s −1 , leading to T = 1.5 MK (Réville et al. 2016).The density at the base of the solar corona is then computed to be consistent with T and Ṁ , which results in 2.49 × 10 7 ≤ n [cm −3 ] ≤ 4.63 × 10 8 .
Since the exponents of the scaling laws and their normalization factors are well-defined, we can now compare the different prescriptions to observations of individual systems.We focus on ZDI studies that constrain the stellar magnetic field and mass loss measurements from astrospheres' Lyα absorption.

A star studied through astrosphere's Lyα absorption and
Zeeman-Doppler Imaging: Eridani We applied our formalism to Eridani, a young active K2V dwarf, which hosts an exoplanet and a debris disk.We study this individual system by only using the minimal information required in our formalism in order to test the different scaling laws.In practical terms, we only rely on the stellar mass and the rotation period, which are essential when using the different prescriptions (see Table 2 for numerical values of those stellar parameters).An estimate of the age indicates that the star is in the main sequence.The stellar luminosity, the stellar radius, and the X-ray luminosity can be estimated through Eqs. ( 12), (13), and (24) (see Table 1 for the associated exponents).The values obtained from those correlations, which are shown in Table 2, are in good agreement with their observed analogs.In practice, those relationships are crucial to estimate B and Ṁ since inaccuracy in the determination of L X would lead to an erroneous value for the mass-loss rate through the Ṁ−F X correlation.
All of those considerations allow us to estimate the stellar magnetic field and the mass loss with our formalism.We illustrate those predictions by the blue bands in Fig. 7. Without further assumptions, our model predicts a broad range of possible B (panel A) and Ṁ (panel B) values.Since we are performing a systematic study of Eridani, we first assume that the star follows the observational trends inferred from ZDI studies.Only the exponents situated in the light gray area that are  (Jeffers et al. 2014;See et al. 2017) and the massloss rate (Wood et al. 2002).The light gray band corresponds to the upper and lower bounds of p B from a large sample of ZDI studies.The black contours delimit the region of B and Ṁ, which is consistent with statistical ZDI studies and individual measurements of the large-scale magnetic field.In beige: estimates consistent with all of the constraints considered.defined in Fig. 7 are considered (1 ≤ p B ≤ 1.65).A modification in the B −Ro relationship then corresponds to a horizontal shift in Fig. 7.A change of normalization between B = 0.6 G (dashed blue line in Fig. 7) and B = 4 G (dotted blue line in Fig. 7), which is linked to the scatter in the ZDI studies' dataset, results in a vertical shift in Fig. 7. Browsing this parameter space then gives 3 ≤ B [G] ≤ 47 and 1.9 The large-scale magnetic field of the star has been monitored for almost seven years, between January 2007 and October 2013, by Jeffers et al. (2014).We only consider the dipolar component of their ZDI maps, leading to field values between 6.15 and 19.8 G (See et al. 2017), thus constraining our estimates of the magnetic field (see the red band at the top of Fig. 7).By taking those observations into account, we can constrain the normalization factors and derive a narrower range for Ṁ: 7.59 × 10 −14 ≤ Ṁ [M yr −1 ] ≤ 4.59 × 10 −13 (see the black contours in panel B of Fig. 7).
Furthermore, Wood et al. (2002) measured the mass-loss rate of the star through its astrosphere's Lyα absorption, leading to ).Therefore, we can assume a mass loss measured through Eridani's astrosphere between 3 × 10 −13 M yr −1 and 1.2 × 10 −12 M yr −1 (see the red band in panel B of Fig. 7).This constraint further refines the acceptable parameters of our model, which leads to B between 6.15 and 8 G (see the beige areas in Fig. 7).It is worth noticing that our B and Ṁ prescriptions are consistent with both statistical approaches and individual measurements, which is evident in the small beige areas.Furthermore, the wind of Eridani has been modeled in 3D with an MHD model by Alvarado-Gómez et al. (2016).They find a mass-loss ranging from 2.77 × 10 −14 to 10.2 × 10 −14 M yr −1 for January 2010.Our approach is also compatible with these values (see panel B in Fig. 2).The range of Ṁ obtained (between 3 × 10 −13 and 4.59 × 10 −13 M yr −1 ) allowed us to estimate the coronal temperature T c and the coronal density n c for this particular star.Depending on the scenario considered (single density scaling, entropy equilibrium, or single temperature scaling), since the normalization factor T was fixed in Sect.4.1.3,the range of p B available provides directly upper and lower bounds of T c with our prescriptions (see Table 1 and Appendix E for numerical values).The range of density is then determined by connecting the lower and upper bounds of Ṁ and T c .The corresponding numerical values for the coronal properties are shown in Table 3.
One can see that the entropy equilibrium scenario yields slightly higher temperatures compared to the single density scaling hypothesis, therefore resulting in lower densities for a given mass-loss rate.In the single temperature scaling scenario, T c is independent of the B and Ṁ prescriptions because of the T l −F X correlation (cf.Eq. ( 37)).Moreover, here, the star is characterized by a low Rossby number in the unsaturated regime, leading to a higher coronal temperature than those derived from the two other scenarios.Hence, the coronal density, in the case of a single temperature scaling, reaches the lowest values in order to keep a consistent mass loss.
To sum up, we show with Eridani that it is possible to have an analytical prediction for the large-scale magnetic field and the mass-loss rate that is in agreement with all of the observational constraints available and, additionally, that it is feasible to infer from the values obtained a range for the coronal properties of the star, according to different scenarios.Furthermore, with a systematic approach, relying on scaling laws and statistical considerations gives a quite large range for B and Ṁ, compared to well-constrained quantities coming from individual studies.In the case of Eridani, such a guess can deviate from the measured value by a factor of 6 at most in B and about one order of magnitude in Ṁ.If for the study of an individual star, additional measurements are required to reduce the interval of confidence, we see how powerful our approach is to guess key trends along with stellar properties in an ensemble approach.Given our minimal set of hypothesis, our formalism provides a good estimate of the relevant quantities from general scaling laws and statistical trends, which is compatible with individual studies.those prescriptions are based on the magnetic energy content in the stellar interior over a wide spectral range, whereas in this work we consider the large scale magnetic field at the stellar surface.To link the two approaches, one can fall back on a filling factor, which may depend on stellar rotation (See et al. 2019).
Stellar metallicity was not directly taken into account in our study.However, it could significantly affect the coronal density and the stellar mass loss (Suzuki 2018), thus influencing the wind braking torque.In our work, the effect of metallicity on the stellar structure is included in the Rossby dependency.Therefore, studying the influence of metallicity on the rotational evolution of solar-type stars may be a promising avenue to test different torque prescriptions.
Furthermore, coronal temperature and density are directly linked to the choice of a wind model.In this work, to connect the mass-loss rate to the coronal properties, we relied on an expression for Ṁ valid in the case of a nonmagnetized outflow of a pressure-driven polytropic wind.Assuming a more realistic model, such as a magnetized wind, would introduce corrections in the mass loss expression due to the magnetocentrifugal effect (Preusse et al. 2005;Johnstone 2017) involving, for instance, the stellar rotation rate and the Alfvén radius.This may lead to implicit relations between the different prescriptions.Given that our formalism relies on the Matt et al. (2015) wind braking torque, which does not take such an effect into account, a pressure-driven hydrodynamic polytropic wind may be more suitable in this context to keep a consistent model.Modifications to the wind torque and the mass-loss rate may be required to deal with very fast rotators.More complex models could also be investigated, such as a polytropic gas with a spatially varying polytropic index (Johnstone et al. 2015b).
In this paper, we do not consider a slow and a fast wind, which would be a way of improving our model ever more.Indeed, for a pressure-driven wind, higher temperatures lead to faster winds.However, an anticorrelation is observed between the terminal speed of the two components of the solar wind and the coronal temperature of the source region, which may be due to a difference in the altitude of the heating region between the fast and the slow component (Geiss et al. 1995;Schwadron & McComas 2003).Therefore, more realistic wind acceleration processes, including the influence of coronal heating, have to be taken into account to deal with a fast and slow wind (Réville et al. 2020;Riley et al. 2019).
We show that the p L exponent (L X ∝ L Ro −p L ) is one of the most important parameters allowing us to constrain all of the different prescriptions efficiently, especially the B and Ṁ massradius dependency as well as the expression of the coronal properties.This way, uncertainty as to this exponent (Pizzolato et al. 2003;Wright et al. 2011) may lead to a significant scatter in our scaling laws.Therefore, our prescriptions could be significantly tightened if the interval of confidence of the p L exponent could be reduced.
In the case of evolved stars, the decrease of the wind braking efficiency (van Saders et al. 2016) was not studied in this work.However, the influence of this phenomenon on the different wind parameters and its eventual inconsistency with other observational constraints may be an application of our formalism.It could be possible to introduce a "resaturation" regime at high Rossby numbers (Sadeghi Ardestani et al. 2017) and to look for hints of breaking of gyrochonology in the physical quantities involved in stellar spin down.This prescription is assumed to be accurate for Tl ≈ 2−10 MK and acceptable as an average down to 0.4 MK.
To estimate the radius of the dead zone, we assume a dipolar magnetic field in the closed-field region, where the pressure gradient and the centrifugal force are not strong enough to distort the field lines.The closed loops then trap hot gas and prevent the flow of a stellar wind, allowing us to neglect the ram pressure of the gas in this area.Here, the edge of the dead zone is assimilated to a limit of confinement of the plasma and therefore can be estimated with the following pressure equilibrium: p th (r dz ) = p mag (r dz ), (B.3)where p th (r dz ), p th (r dz ) are the thermal and the magnetic pressures of the gas that are estimated at the edge of the dead zone, respectively.Since we consider a dipolar magnetic field, Eq. (B.3) becomes: If we consider that the dead zone is filled with an ideal gas of constant temperature Tl = T l and a density at the edge of the dead zone that evolves similarly to the density at the base of the corona n l , that is, n l (r dz ) = Kn l in which K is a constant independent of any stellar parameter, then p th (r dz ) = n l (r dz )k B T l = Kn l k B T l .The radius of the dead zone then becomes: where: If we assume that nl varies the same way as the density at the base of the corona, we can express the X-ray luminosity as a function of more explicit stellar parameters as follows: By assuming r 3 dz /R 3 1, the X-ray luminosity becomes: Fig.1.Sketch of the multiple couplings between the inputs of the model (stellar mass, radius, and rotation rate), the magnetic field, the X-ray luminosity, the coronal temperature and density, the mass loss, and the wind braking torque.In blue: model assumptions.In red: observational constraints.

Fig. 2 .
Fig. 2. Correlation p B − p Ṁ for a dipolar (solid black line), a quadrupolar (dashed black line), and an octupolar (dotted black line) topology in the unsaturated regime.Dashed blue lines: magnetic field prescriptions from See et al. (2017, S17) and Vidotto et al. (2014, V14).Blue crosses: magnetic field and mass loss prescriptions from Johnstone et al. (2015a, J15), Tu et al. (2015, T15), and Sadeghi Ardestani et al. (2017, SA17).In dark gray: excluded region in the dipolar case.Dashed dark red lines: prescriptions compatible with the Skumanich law for any magnetic topology.In light red: exponents corresponding to mass losses outside the envelope of Wood et al. (2005) data.The light gray band corresponds to the upper and lower bounds from Zeeman-Doppler imaging (ZDI) statistical studies (see Sect. 3.4 for more details).

Fig. 3 .
Fig. 3. Correlation between p B and ξr B + m B in the unsaturated regime for p L = 2 (black line) and p L = 3 (dotted).Here, η = 4 and ξ = 0.9.In dark gray: excluded region.In red: exponents corresponding to mass losses outside the envelope of Wood et al. (2005) data.In green: exponents consistent with the mass-loss constraints.The light gray band corresponds to the upper and lower bounds from ZDI statistical studies (see Sect. 3.4 for more details).

Fig. 4 .
Fig. 4. Sketch of the coupling between the coronal temperature and density in the open-field and the closed-field regions, the X-ray luminosity, and the mass-loss rate.Several scenarios are considered to connect the closed loops to the coronal holes.Scenario 1, in red: single temperature scaling.Scenario 2, in purple: entropy equilibrium.Scenario 3, in blue: single density scaling.

Fig. 5 .
Fig. 5. Evolution of T c (panel A) and n c (panel B) as a function of p B for a solar twin with a rotation period of 15 days.Solid line: p L = 2. Dashed lines: p L = 3.In black: Holzwarth & Jardine (2007) prescription.In light red: single temperature scaling scenario.In purple: entropy equilibrium scenario.In blue: single density scaling scenario.In hatched gray: values of coronal temperatures leading to Ṁ < 0. The light gray band corresponds to the upper and lower bounds of p B , which we consider to be consistent with ZDI statistical studies.

Fig. 6 .
Fig.6.Mass loss per unit area as a function of the X-ray flux.The red dots come fromWood et al. (2005; W05)  and the blue dots are fromJardine & Collier (2019; JC19).In green: measurements from the atmospheric evaporation of GJ 436b(Vidotto & Bourrier 2017, VB17).Gray dotted line:Wood et al. (2005) scaling law corresponding to w = 1.34.Red dashed line:Suzuki et al. (2013) scaling law corresponding to w = 0.82.The light gray region corresponds to the upper and lower bounds we considered to be consistent with ZDI statistical studies.

Fig. 7 .
Fig. 7. Panel A: estimate of B for Eridani as a function of p B .Panel B: estimate of Ṁ as a function of p B .In blue: predictions of our prescriptions.Dashed blue lines: normalizations corresponding to B = 0.6 G. Dotted blue lines: normalizations corresponding to B = 4 G. Red bands: observational constraints on the dipolar component of the magnetic field(Jeffers et al. 2014;See et al. 2017) and the massloss rate(Wood et al. 2002).The light gray band corresponds to the upper and lower bounds of p B from a large sample of ZDI studies.The black contours delimit the region of B and Ṁ, which is consistent with statistical ZDI studies and individual measurements of the large-scale magnetic field.In beige: estimates consistent with all of the constraints considered.

Table 1 .
Parameters defining the Γ wind , B , Ṁ, T c , and n c prescriptions for p B = 1 and p B = 1.65 in the single density scaling scenario.