Issue 
A&A
Volume 631, November 2019



Article Number  A111  
Number of page(s)  13  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201936477  
Published online  05 November 2019 
Does magnetic field impact tidal dynamics inside the convective zone of lowmass stars along their evolution?
^{1}
AIM, CEA, CNRS, Université ParisSaclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 GifsurYvette, France
email: aurelie.astoul@cea.fr
^{2}
IRAP, Observatoire MidiPyrénées, Université de Toulouse, 14 avenue Edouard Belin, 31400 Toulouse, France
^{3}
Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France
^{4}
Observatoire de Genève, Université de Genève, 51 chemin des Maillettes, 1290 Sauverny, Switzerland
Received:
7
August
2019
Accepted:
21
September
2019
Context. The dissipation of the kinetic energy of wavelike tidal flows within the convective envelope of lowmass stars is one of the key physical mechanisms that shapes the orbital and rotational dynamics of shortperiod exoplanetary systems. Although lowmass stars are magnetically active objects, the question of how the star’s magnetic field impacts largescale tidal flows and the excitation, propagation and dissipation of tidal waves still remains open.
Aims. Our goal is to investigate the impact of stellar magnetism on the forcing of tidal waves, and their propagation and dissipation in the convective envelope of lowmass stars as they evolve.
Methods. We have estimated the amplitude of the magnetic contribution to the forcing and dissipation of tidally induced magnetoinertial waves throughout the structural and rotational evolution of lowmass stars (from M to Ftype). For this purpose, we have used detailed grids of rotating stellar models computed with the stellar evolution code STAREVOL. The amplitude of dynamogenerated magnetic fields is estimated via physical scaling laws at the base and the top of the convective envelope.
Results. We find that the largescale magnetic field of the star has little influence on the excitation of tidal waves in the case of nearlycircular orbits and coplanar hotJupiter planetary systems, but that it has a major impact on the way waves are dissipated. Our results therefore indicate that a full magnetohydrodynamical treatment of the propagation and dissipation of tidal waves is needed to properly assess the impact of starplanet tidal interactions throughout the evolutionary history of lowmass stars hosting shortperiod massive planets.
Key words: magnetohydrodynamics (MHD) / waves / planetstar interactions / stars: evolution / stars: magnetic field / stars: rotation
© A. Astoul et al. 2019
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Over the last two decades, a large variety of exoplanetary systems has been discovered, primarily through photometric transit and radial velocity observations (e.g. Perryman 2018). Among these systems, several populations of exoplanets orbit very close to their host stars (with orbital periods of a few days), such as hotJupiters, superEarths, and miniNeptunes (Mayor et al. 1997; Schlichting 2014). In these short period systems as well as in tight binary stars, tides induced by each other’s body drive the rotational and orbital evolutions of the system through dissipation mechanisms (the socalled tidal dissipation, see e.g. Zahn 1977; Hut 1980; Zahn & Bouchet 1989; Goldreich & Nicholson 1989; Witte & Savonije 2002).
Usually, tides are split into two components following the work of Zahn (1966a,b,c, 1975, 1977). First, the equilibrium tide is the quasihydrostatic response of the main body to tidal perturbations induced by the companion. It is materialised by a largescale deformation in the perturbed body, that is displayed as a nearequatorial tidal bulge in the direction of the companion. In addition, waves are excited in its interior as the equilibrium tide is not an exact solution to the equations of motion, providing an additional driving force (Ogilvie 2014). The tidallyforced waves correspond to the dynamical tides. The convective zone (CZ) of a perturbed rotating body is the seat of inertial waves that get dissipated by turbulent friction (e.g. Ogilvie & Lin 2004, 2007), while the radiative zone supports gravitoinertial waves which are dissipated by thermal damping and turbulent friction (e.g. Zahn 1975; Terquem et al. 1998; Goodman & Dickson 1998; Barker & Ogilvie 2010).
The dissipation of the dynamical tides in the CZ is very efficient in young stars (premain sequence and early main sequence) whereas for older stars the dissipation of the equilibrium tide dominates mainly because of the slower stellar rotation, as demonstrated for instance by Bolmont & Mathis (2016) (see also Strugarek et al. 2017; Gallet et al. 2017; Beck et al. 2018). The efficiency of the tidal excitation and viscous dissipation of inertial and gravitoinertial waves can be measured with the tidal quality factor Q. This quantity reflects the fact that the perturbed body undergoes a forced oscillation and dissipates a fraction of the associated energy during each oscillation period. It has been evaluated by Ogilvie & Lin (2004, 2007), with and without differential rotation (see also Guenel et al. 2016a,b, for this particular topic), in the case of giant planets and solartype stars. Recently, Mathis (2015) and Gallet et al. (2017) explored the influence of mass, age and rotation based on the frequencyaveraged dissipation estimates of Ogilvie (2013) to understand the behaviour of tidal dissipation along the evolution of stars. They emphasised that the variation of these parameters could drastically modify the strength of tidal dissipation with a higher frequencyaveraged tidal dissipation for lowmass stars. Similarly, Bolmont et al. (2017) highlighted the importance of the stellar metallicity. Within this context, the variations of tidal dissipation along stellar evolution have a strong impact on the orbital architecture of compact planetary systems and the planet survival (Bolmont & Mathis 2016; Benbakoura et al. 2019).
In the aforementioned studies of starplanet tidal interactions, an important ingredient is missing though: stellar magnetism. In the Sun and solarlike stars, magnetism is revealed by external magnetic features such as sunspots, prominences or flares (Donati & Landstreet 2009). The magnetic fields of solarlike stars originates from a powerful dynamo mechanism, sustained by turbulent convection and differential rotation in the convective envelope of the star (Brun & Browning 2017). Recent endeavours have been carried out to assess the effects of magnetism on tidallyexcited inertial waves in stars (Wei 2016, 2018; Lin & Ogilvie 2018). In the presence of a magnetic field, tidal waves excited in the CZ become magnetoinertial waves. Moreover, they feel the magnetic tension of the largescale magnetic field which affects their propagation and dissipation (Finlay 2008). These magnetically modified inertial waves have a broader range of propagation frequencies, compared to the hydrodynamical case, and can be dissipated through both viscous and Ohmic processes. Specifically, the transition between hydrodynamical and magnetohydrodynamical (MHD) regimes have been explored in a shearing box model with a uniform magnetic field (Wei 2016) and in spherical geometry with both a uniform field directed along the zaxis and a dipolar magnetic field (Lin & Ogilvie 2018). The authors of these studies both stressed that the Lehnert number Le (Lehnert 1954) determines how important magnetism is to tidal dissipation. This dimensionless number compares the Alfvén velocity to the rotation speed of the body. Additionally, how magnetism influences the effective tidal force to excite magnetoinertial waves remains to be addressed, as mentioned in Lin & Ogilvie (2018, Appendix B). In practice, given the equation of motion tidal waves are excited by an effective body force driven by the Coriolis acceleration of the equilibrium tidal flow in the nonmagnetised case (Ogilvie 2005). In the presence of a magnetic field, the Lorentz force acting on the equilibrium tide is likely to play a role in the excitation of tidal waves. Thus, this motivates the study of the impact of stellar magnetism on both dissipation and excitation of dynamical tides inside the CZ and along the evolution of lowmass stars.
The paper is organised as follows. In Sect. 2.1 we work out the contribution of magnetism to tidal forcing and derive a criterion to assess its importance relative to nonmagnetised forces. This criterion features the Lehnert number and thus the magnetic field and the rotation speed of the body, as well as the tidal forcing frequency. Simple scaling laws are applied in Sect. 2.2 to estimate the dynamodriven magnetic field in the convective envelope of lowmass stars. Thanks to the stellar evolution code STAREVOL (Sect. 2.3, see also Amard et al. 2019, and references therein), the strength of a mean magnetic field is given, in Sect. 2.4, at the base and top of the CZ for various lowmass stars. Thus, we evaluate in Sect. 2.5 the Lehnert number as a function of age, mass, initial rotation, and radius in the CZ of these stars. We then estimate in Sect. 2.6 the Lehnert number along with the rotation and tidal frequencies in several observed shortperiod exoplanetary systems to assess the importance of the star’s magnetic field on the tidal forcing. Based on our estimates of the Lehnert number, we compare in Sect. 3 the relative importance of Ohmic over viscous dissipations of tidally induced magnetoinertial waves throughout the evolution of lowmass stars. In Sect. 4, we examine how smallscale magnetic fields impact tidal forcing. Finally, we present in Sect. 5 the conclusions and perspectives of this work.
2. The Lorentz force influence on tidal forcing
The purpose of this section is to quantify the contribution of the stellar magnetic field to the tidal excitation of magnetoinertial waves in the CZ of lowmass stars.
2.1. A criterion to settle the importance of magnetism
We aim for an approach as general as possible, yet we restrict our model adopting a solidbody rotation with an angular frequency Ω. In particular, we do not assume a specific geometry for the magnetic field. We linearise the momentum and induction equations to derive the magnitude of the magnetic tidal forcing (the effective tidal force arising from the Lorentz force) and compare it with the classical hydrodynamical tidal forcing. We introduce the selfgravitational potential Φ_{0}, as well as the gravitational potential perturbation Φ, and the external tidal potential Ψ (see e.g. Zahn 1966a; Ogilvie 2013). The continuity and entropy equations are as given in Zahn’s paper, with the addition of Ohmic heating to the entropy equation. The momentum equation for tidal perturbations in the corotating frame can be written as:
where ρ_{0} is the mean density and u, p, and ρ are the perturbed velocity, pressure, and density, respectively. We include the volumetric viscous force F_{ν} = ρ_{0}ν∇^{2}u, which represents the effective action of turbulent convection on tidal flows with ν the effective socalled eddyviscosity (e.g. Zahn 1966b, 1989; Ogilvie & Lesur 2012; Mathis et al. 2016) which we assume to be constant in the CZ. Moreover,
is the linearised Lorentz force, where B_{0} and b are the largescale and perturbed magnetic fields, respectively. Note also that we have no background flow (u_{0} = 0) as we work in the rotating frame, the action of the convective flows is parametrised as a diffusion, and any differential rotation and associated meridional flows are neglected. Furthermore, the linearised induction equation is:
with η the magnetic turbulent diffusivity, related to the eddyviscosity by the relationship η = ν/Pm, where Pm is the turbulent magnetic Prandtl number often chosen close to unity (e.g. Camargo & Tasso 1992; Jurčišinová et al. 2013; Käpylä et al. 2019).
Following Ogilvie (2005, 2013), we decompose all physical perturbed quantities X into a nonwave like part associated with the equilibrium tide denoted as X_{e}, and a wavelike part X_{d} related to the dynamical tides. The equilibrium tidal flow u_{e} is defined as the velocity resulting from the hydrostatic adjustment of the primary due to the perturbation induced by the companion in the rotating frame of the tidal bulge. This frame rotates at the corresponding tidal frequency nΩ_{o}/2 (see Remus et al. 2012), where n labels the temporal harmonic of the orbital motion of the perturber (when projecting the tidal potential on the spherical harmonics basis) and Ω_{o} is the associated orbital frequency. In the adiabatic case, this hydrostatic equilibrium leads to (Zahn 1966a):
This equation comes from Eq. (1) without the left hand side and dissipative terms. As a first step, we also neglect the deformation of the stellar structure induced by rotation and magnetic field. In the case of magnetic fields this is a reasonable assumption except in the low density region near the stellar surface (Duez et al. 2010). For the centrifugal acceleration, this is a fair hypothesis for slow and median rotators while potentially strong deformation should be taken into account for young rapid rotators (see e.g. Gallet & Bouvier 2013, Fig. 7). We split the equation of induction in the corotating frame, accounting for the equilibrium and dynamical tides decomposition:
where ξ_{e} is the equilibrium tide displacement, defined by u_{e} = ∂_{t}ξ_{e}, given a mean static magnetic field B_{0}. We also introduce u_{d}, the perturbed flow of the dynamical tide. We assume here that B_{0} does not vary over the tidal timescale which is a few days. This is corroborated by the fact that the largescale magnetic field varies very little (far below one order of magnitude) within several years for the majority of observed stars (Vidotto et al. 2014). Moreover, we choose to define the magnetic field associated with the equilibrium tide as the field it induces by the advection of B_{0}. We neglect the Ohmic diffusion acting on the equilibrium tide because its time scale R^{2}/η, where R is the radius of the star which is also the length of variation of this flow, is much larger than its typical time of variation (in the range of a few days for Hot Jupiter), even when considering a turbulent magnetic diffusivity. One can not do the same assumption for dynamical tides since they involve potentially smaller length scales, for example along waves’ attractors. The equation of induction for dynamical tides (Eq. (4b)) follows from this definition when writing the equation of induction for the sum of the equilibrium and dynamical tide perturbations (Eq. (2)) since it is a linear equation.
In the momentum equation, we use the Cowling approximation (Cowling 1941) for the dynamical tides i.e. we neglect their perturbed gravitational potential Φ_{d}. In addition, as our model applies to a convective region (i.e. adiabatically stratified), the term ρ_{d}∇Φ_{0} is neglected because it is related to the buoyancy force associated with gravity waves. Using Eq. (3), the momentum equation for tidallyforced magnetoinertial waves becomes:
where the wavelike part encompassing the propagation of tidal waves (on the left hand side of the equation), is excited by an effective force driven by the equilibrium tidal flow (on the right hand side). This force can be written as:
It is worth noting that the action of turbulent friction on the hydrostatic flow has been neglected for the same reasons we ignored the ohmic diffusion in Eq. (4a). The term f_{hydro} comprises the driving inertial force and the Coriolis acceleration (see the Appendix B in Ogilvie 2005) while f_{mag} embodies the action of the Lorentz force on the hydrostatic displacement and has been derived by Lin & Ogilvie (2018, see Appendix B). These authors studied the propagation and dissipation of magnetoinertial waves excited by the effective forcing induced solely by the Coriolis acceleration of the equilibrium tide (in short f(u_{e}) = f_{hydro}). However, they also suggest that a largescale magnetic field can potentially interact with the equilibrium tide for sufficiently large Lehnert numbers (typically Le > 0.1). For this reason, we propose to examine the relative importance of both forcings f_{mag} and f_{hydro}. We use R as the typical length scale of the largescale magnetic field and of the equilibrium tide, which involves largescale flows. Henceforth, we can give the order of magnitude of the different forcings:
and their ratio:
where we define the Lehnert number as . We also introduce the Dopplershifted tidal Rossby number Ro_{t} = σ_{t}/(2Ω) with σ_{t} = u_{e}/ξ_{e} the related tidal frequency, and a dimensionless frequency ratio . According to Eq. (6), magnetism needs to be taken into account in the tidal forcing whenever .
2.2. Scaling laws to estimate stellar magnetic fields
The determination of the Lehnert number inside the convective envelope of lowmass stars requires knowledge of the internal magnetic field of these stars. However, we are currently only able to constrain the magnetic field of stars at their surface. Indeed, thanks to ZeemanDoppler Imaging, one can reconstruct the topology and strength of largescale, stellar magnetic fields (Donati et al. 2006, 2007).
Regarding the Sun, internal magnetic fields can be assessed indirectly from surface tracers like sunspots, as a manifestation of flux ropes emerging through the surface (Charbonneau 2013). This approach is based on the socalled interface dynamo theory where the magnetic field is generated by a convective dynamo and pumped into the tachocline (the interface between the radiative and convective zones), where it eventually becomes strong enough to be buoyantly unstable, rise through the convection zone, and emerge at the surface as sunspots (Spiegel & Zahn 1992; Charbonneau 2014; Brun & Browning 2017). This picture has been recently questioned, especially by Wright & Drake (2016) by studying fully convective stars. It is indeed possible that the spotforming magnetic fields can be generated in the bulk (Strugarek et al. 2017) or in the shallow layers (Basu 2016) of the CZ instead of overshoot layers beneath the coreenvelope interface. Through helioseismology, Gough & Thompson (1990) and Antia et al. (2000) have placed an upper bound of 30 T on the toroidal magnetic field strength at the base of the CZ. On the contrary, the mean magnetic field at the Sun’s surface is significantly smaller, about a few Gauss (∼10^{−4} T), even though sunspots are the seats of local intense magnetic fields (several tenths of a Tesla).
As far as younger or smaller stars than the Sun are concerned, Vidotto et al. (2014) gives an overview of known largescale surface magnetic fields, that appear to vary from the globally weak fields of the order of the Gauss for solarlike stars to the strong Teslastrength fields of M dwarf and T Tauri stars. They establish a relationship between surface magnetic field and the convective Rossby number Ro = u_{c}/(2Ωl_{c}), where u_{c} and l_{c} are the convective velocity and length, respectively. This dimensionless quantity is generally calculated at half the mixinglength αH_{p}/2 (Landin et al. 2010; Gilliland 1986) above the base of the CZ, with H_{p} the pressure scale height and α the mixinglength theory coefficient.
Nonetheless, as indicated above, it is difficult to estimate the magnetic field inside stars. In this context, 3D global nonlinear simulations (Strugarek et al. 2017; EmeriauViard & Brun 2017) and scaling relationships for stellar dynamos (Augustson et al. 2019, and references therein) can help us to estimate the internal magnetic field strength in the convective envelope of lowmass stars. Three scaling laws are described hereafter and the derivation of the related dynamoinduced magnetic field is detailed in Appendix A. We have:
– the (turbulent) equipartition, that is often used to give an averaged, rough estimate of the magnetic field’s amplitude in the bulk of the CZ (Brun & Browning 2017). It assumes that the dynamo is efficient i.e. that the system is equally good at generating magnetic field as it is at generating flows. This balance is also used in moderately active plages at the solar surface, while in the active sunspots superequipartition (i.e. magnetic energy is greater than kinetic energy) can be fairly common (Donati & Landstreet 2009).
– the buoyancy dynamo regime, in which the Coriolis, buoyancy and Lorentz forces are taken to have the same order of magnitude assuming a low atomic magnetic Prandtl number, that is the ratio of atomic viscosity to magnetic diffusivity (Davidson 2013; Augustson et al. 2019). This assumption is well verified for fast rotating giant planets, young contracting stars like TTauri and rapidlyrotating lowmass stars (Christensen et al. 2009).
– the magnetostrophic regime, for which the force balance is realised between the Coriolis and Lorentz forces. This balance, also called magnetostrophy, gives an upper estimate of the magnetic field. The magnetostrophic regime generally gives fields in superequipartition (Brun et al. 2015; Augustson et al. 2019).
We summarise in Table 1 the order of magnitude estimation of the magnetic field (named hereafter B_{dyn}) within the convective envelope of a lowmass star as obtained by the three aforementioned scaling laws. These relationships involve the kinetic and magnetic energy densities of the convective flow and respectively, along with the convective Rossby number Ro = u_{c}/(2Ωl_{c}) at the base of the CZ.
Magnetic fields derived from simple balances (forces or energies) written in the second column.
The main objective of this paper is to determine and quantify the impact of a largescale magnetic field on the excitation and dissipation of the dynamical tides. We still discuss in Sect. 4 to what extent the smallscale component of the stellar magnetic field can influence the tidal flows. To allow comparison with the amplitude of observed surface magnetic fields of lowmass stars, we extrapolate a largescale, dipolar, surface magnetic field from the scaling laws at the base of the CZ (Table 1). First, we suppose that the dipolar component of the star’s magnetic field at the base of the CZ is a fraction γ of the dynamoinduced magnetic field at that location:
The factor γ encapsulates both the ratio of largescale to smallscale magnetic fields (similarly to the filling factor in Reiners 2012; See et al. 2019) and the part of the total energy that is available in the dipolar component of the magnetic field.
Then, we infer the dipolar component of the surface magnetic field simply as:
where R is the radius of the star, the subscripts “base” and “top” are the position inside the CZ, and “dip” and “dyn” refer to the dipolar and dynamoinduced magnetic fields, respectively. As the dipolar magnetic field at the surface of the Sun is well known, typically 4 G (DeRosa et al. 2012), γ can be estimated for the Sun by using Eqs. (7) and (8):
with R_{⊙}, B_{dyn, ⊙}(r_{base}) and r_{base, ⊙} obtained from the grid models of the STAREVOL evolution code (see the Sect. 2.3 below) for a 1 M_{⊙} star of the age of the Sun. In the following we assume that γ in Eq. (7) is independent of the mass and the age of the star, and that it takes the Sun’s current value as in Eq. (9). This is a strong assumption since γ is close to a filling factor that depends on the Rossby number and therefore on the angular frequency of the star (See et al. 2019). Nevertheless, refining the expression of this factor would not change the final conclusions of this paper which are robust to several orders of magnitude, as we see later. The factor γ however depends on the scaling law used to estimate B_{dyn, ⊙}.
Finally, the dipolar component of a star’s surface magnetic field will be estimated as
2.3. The stellar evolution code STAREVOL
To estimate the magnetic field and then the Lehnert number in the convective envelope of lowmass stars of different ages via the scaling laws described in Sect. 2.2, we use the 1D stellar evolution code STAREVOL (Amard et al. 2019). Initial masses of the stars range from 0.4 to 1.4 M_{⊙}, given a solar metallicity Z = 0.0134 (Asplund et al. 2009), and a mixing length parameter α = 1.9730. This latter is defined by the calibration of the standard solar model and used to model convective regions according to the mixing length theory. Basic input microphysics like the equation of state, nuclear reactions or opacities, are described in Amard et al. (2016) and Lagarde et al. (2012). The initial rotation periods are fixed using the calibration for fast (1.6 days), median (4.5 days) and slow (9 days) rotators from Amard et al. (2019). The rotation is assumed to be uniform inside the CZ but varies dramatically with stellar ages (Gallet & Bouvier 2013). As a result, the evolution of the surface angular velocity dictates that of the Lehnert number. During the first few Myr of the premain sequence (PMS), the surface angular velocity of the stars remains stationary as the result of stardisk magnetic interactions (Zanni & Ferreira 2013; Gallet & Bouvier 2015; Amard et al. 2016). This holds over the disk lifetime, typically a few Myr (Rebull et al. 2004; Bell et al. 2013; Gallet & Bouvier 2015). After the dissipation of the disk, the gravitational contraction of the star leads to an increase in the angular velocity, in order to conserve angular momentum, until the star begins hydrogen fusion at the zeroage main sequence (ZAMS). From this stage onward, magnetised stellar winds apply a torque on the star, spinning it down throughout its main sequence (MS) lifetime. In the STAREVOL code, the effects of a stellar wind acting from the early PMS to the tip of the MS are implemented using the prescription given by Matt et al. (2015). In the following, the base of the CZ refers to a height that is located 0.002R above the bottom of the CZ as computed in the STAREVOL model. This convention avoids erratic numerical behaviour of the mixing length convective velocity at the interface between the radiative and convective zones. Furthermore, the top of the CZ refers to the radius where the convective velocity vanishes in the STAREVOL models.
2.4. Estimation of the dipolar magnetic field at the base and the top of a CZ
Figure 1 shows the evolution of the magnetic field of a 0.9 M_{⊙} star along its lifetime for the three different initial rotation rates stated in Sect. 2.3. It should be specified that we chose a 0.9 M_{⊙} star rather than 1 M_{⊙} to add measurements of the mean dipolar magnetic fields at the surface of stars in the early MS (see next paragraph). At the base of the CZ (left panel), the magnetic field is calculated with the scaling laws listed in Table 1. At the top of the CZ (right panel), we use these dynamoinduced magnetic fields to extrapolate dipolar magnetic fields near the surface by using Eq. (10). The results obtained with fast (1.6 days), median (4.5 days), and slow (9 days) initial rotations are plotted with solid, dashdotted, and dashed curves, respectively. We note that the magnetic field decreases in time after about 50 Myr in the two panels for both the magnetostrophic and buoyancy dynamo regimes. It is due to the fact that these regimes depend on a positive power of the angular velocity (in the denominator of the Rossby number, see Table 1) that decreases after the ZAMS as a result of the stellar wind action on the star’s surface. After 1 Gyr, this decline is well described by the empirical Skumanich relationship (Weber & Davis 1967; Skumanich 1972). At the age of the present Sun (∼4.6 Gyr), the buoyancy dynamo and magnetostrophic regime at the base of the CZ give the order of magnitude of the toroidal magnetic field strength expected in the Sun at the tachocline, typically a few to a few tens of Tesla (Charbonneau 2013).
Fig. 1. Magnetic field versus the age of a 0.9 M_{⊙} star for three different initial rotation periods. At the base of the CZ (left panel), the magnetic field is estimated from the scaling laws listed in Table 1. At the top of the CZ (right panel), the same scaling laws are used to extrapolate a dipolar magnetic field near the surface (see Eq. (10)). The “stars” symbols () that represent mean dipolar magnetic fields at the surface of 0.9 M_{⊙} stars are taken from See et al. (2017), and the age of the stars are from Ge et al. (2006), Janson et al. (2008), Line et al. (2015), and Folsom et al. (2016). Multiple observations of the dipolar magnetic field of a star are joined by a line. 
We have added on the plot for the top of the CZ (righthand panel) the average unsigned measured dipolar field strength of 0.9 M_{⊙} stars (See et al. 2017). We have adopted a conservative error estimate of 0.434 dex in log B_{dip} for all stars. Note that HD 22049 displays several values of the amplitude of the dipolar magnetic field measured at different times. The ages of the stars and their errors are taken from Folsom et al. (2016) for the early MS stars and Ge et al. (2006), Janson et al. (2008), Line et al. (2015) for the three oldest stars. It is interesting to note that the strength of the observed dipolar magnetic fields seems to be steady for end PMS and early MS stars. This finding is consistent with the saturation levels observed by Vidotto et al. (2014), See et al. (2017) for the largescale and dipolar surface magnetic fields at low Rossby numbers. As a matter of fact, low Rossby numbers means high angular frequencies and therefore stars in the vicinity of the ZAMS. However, whether or not this observed saturation level is an atmospheric effect or a dynamo related phenomenon has yet to be determined (Vilhu 1984; Jardine & Unruh 1999).
Dipolar fields of young stars seem to fit quite well with the magnetostrophic regime along with fast initial rotation. We must nevertheless specify that we realised the same plot with more observed dipolar magnetic fields (See et al., in prep., priv. comm.) featuring weaker magnetic fields for PMS and early MS stars, and so consistent with median and slow initial rotation. That being said, we have to bear in mind that several strong assumptions have been made to derive the dipolar magnetic field near the surface (Eq. (10), see also Appendix B). These prescriptions are nonetheless sufficient for this work given the robustness of our results as we see in Sect. 2.6.
We have plotted in Fig. 2 the ratio of the observed and estimated dipolar magnetic fields versus the mass of various lowmass stars distributed from PMS to MS stages. The magnetic field B_{sim} is calculated with the magnetostrophic regime and a median initial rotation. The magnetic field B_{obs} is again taken from See et al. (2017) and ages are from Vidotto et al. (2014) based on different methods (see the last quoted paper for more details). We note that the surface dipolar magnetic field of stars is well reproduced by B_{sim} within an order of magnitude. The estimate tends to deviate from for relatively massive (M_{⋆} ≳ 1.2 M_{⊙}) or very lowmass (M_{⋆} ≲ 0.6 M_{⊙}) stars, which is not surprising considering that the ratio γ of dipolar/non dipolar magnetic field is determined from solar parameters (see Eq. (9)).
Fig. 2. Ratio of the mean unsigned observed dipolar magnetic field (See et al. 2017; Vidotto et al. 2014) over estimated dipolar magnetic field as a function of the mass of the star. The age of the stars is indicated in colour and ranges from 21 Myr to 17.2 Gyr. Magnetostrophic regime and median initial rotation have been used to plot the ratio . Multiple observations of the dipolar magnetic field of a star are again joined by a line. A typical error bar for M_{⋆} and is indicated in grey. We have adopted a conservative error estimate of 0.434 dex in log B_{dip} and 0.1 M_{⊙} in M_{⋆}. 
2.5. Lehnert number for a lowmass star along its evolution
Using the estimates of the dynamoinduced magnetic field in Table 1, we can express the related Lehnert number at the base of the CZ (see Table 2). Within each regime, Le depends on the convective Rossby number, and the ratio between the convective length scale and the stellar radius l_{c}/R. This ratio results from the different length scales used in the magnetic scaling laws and in the definition of the Lehnert number (the length scales l_{c} and R respectively). Specifically, we choose a definition of the Lehnert number that is consistent with the previous works of Lin & Ogilvie (2018) and Wei (2018).
At the top of the CZ, we use the dipolar magnetic field derived in Eq. (10) to estimate the Lehnert number near the surface:
where we recall that “top” and “base” refer to the position in the CZ and that Le_{base} is taken from Table 2.
In Fig. 3, the Lehnert number squared of a 1 M_{⊙} star is shown as a function of age, at the base (left panel, see Table 2) and top (right panel, Eq. (11)) of the CZ. The different scaling laws, with fast and slow initial rotations, are shown with the same layout as Fig. 1. First, we note that Le^{2} remains always smaller than unity, consistent with the previous works of Lin & Ogilvie (2018) and Wei (2018). Whether we look at the base or the top, Le^{2} increases with decreasing initial rotation speed from the PMS until about 1 Gyr. This is the opposite behaviour to the magnetic field strength, since the Lehnert number decreases with the angular velocity (see Table 2 and Eq. (11)). Finally, Le^{2} is greater at the base than at the top of the CZ as it scales with B^{2}.
Fig. 3. Lehnert number squared at the base (left panel) and the top (right panel) of the CZ as a function of age for a 1 M_{⊙} star, and for different magnetic scaling laws and initial rotations (see the legends in the left panel). 
The mass dependence of the Lehnert number is illustrated in Fig. 4, which displays Le^{2} against age for stars ranging from 0.7 to 1.4 M_{⊙}. We have used slow initial rotation and the magnetostrophic regime to plot Le^{2} at both the base and the top of the CZ. At the base of the CZ, we obtain the same features as in Fig. 3 i.e. Le^{2} reaches a minimum at the ZAMS and increases during the MS. We highlight a quite small dispersion of the curves, though Le^{2} still grows slightly with mass for M_{⋆} ≲ 1.2 M_{⊙} stars. This trend is driven mostly by the decay of the convective turnover time τ_{c} = l_{c}/u_{c} when the mass grows (note that Ro = (τ_{c}2Ω)^{−1} in Table 2). When M_{⋆} ≳ 1.2 M_{⊙}, the decrease in l_{c}/R at the base of the CZ for growing masses helps to reverse this behaviour. The drop in Le^{2} with mass is even more pronounced at the top of the CZ for M_{⋆} ≳ 1 M_{⊙}. These changes simply reflects the fact that the convective envelope shrinks with mass, which naturally leads to an increase in convective length and a decrease in convective turnover time and density at the base of the CZ (see Table 2 and Eq. (11)).
Fig. 4. Evolution of the Lehnert number squared over time, at the base (solid lines) and the top (dotted dashed lines) of the CZ, for various lowmass stars sorted from 0.7 to 1.4 M_{⊙}. Slow initial rotation has been used here and a magnetostrophic balance has been assumed. 
2.6. The influence of magnetism on tidal forcing for observed starplanet systems
The ratio of the magnitude of the Lorentz forcing to the hydrodynamical forcing (Eq. (6)) does not only depend on the Lehnert number, but also on the Dopplershifted Rossby number Ro_{t} and the ratio of frequencies which need to be estimated from a twobody system. Ogilvie (2014) gives the expression of the tidal frequency in the fluid frame, σ_{t} in our notations, using the integers l, m, n coming from the spherical harmonics functions on which the gravitational potential is projected (we refer the reader to Sect. 2.1 in Ogilvie 2014). The tidal frequency is σ_{t} = nΩ_{o} − mΩ, where m is the azimuthal order of the spherical harmonics, and n and Ω_{o} have been introduced in Sect. 2.1. In the previous section, it has been shown that Le^{2} does not go much beyond 10^{−3}, regardless of the age and the mass of the star at the base and top of the CZ (Figs. 3 and 4). This means that the ratio f_{mag}/f_{hydro} (Eq. (6)) is likely to be small compared to unity, unless the rotation frequency is much greater (by at least a factor of a thousand) than the tidal frequency. Given the definition of the tidal frequency introduced above, the closer we get to a resonance between orbital and rotation frequencies, the more important the ratio f_{mag}/f_{hydro} will be.
For the sake of simplicity we consider here systems with an almost circular and coplanar orbit. This allows us to reduce the number of pairs (m, n) because the tidal potential components are labelled by these integers and depend on eccentricity and stellar obliquity. Moreover, the dominant term in the tidal potential is the quadrupolar component as long as the planet and its host star are well separated, namely l = 2 with l the degree of the spherical harmonics (Mathis & Le PoncinLafitte 2009; Ogilvie 2014). Within this assumption and the limits of low eccentricity and obliquity, (m, n)∈{(2, 2),(0, 1),(2, 1),(2, 3),(1, 0),(1, 2)} (see Ogilvie 2014, for more mathematical details). When the orbit is strictly circularised and coplanar, the asynchronous tide acts alone and the only matching pair of integers is (2, 2). For the other pairs, the eccentricity or obliquity tides can be dominant.
In Table B.1, we have thus listed known starplanet systems satisfying the following conditions:
– a nearcircular orbit: we choose the eccentricity such as e < 0.1.
– a low skyprojected obliquity: λ < 30°. λ is the skyprojected angle between the stellar spin axis and the axis perpendicular to the planet orbit. Ideally, we should use the true obliquity ψ but this quantity is more difficult to determine than λ and too few measurements exist. However, when both values of ψ and λ exist for the selected systems, they are quite similar and far from the threshold of 30°.
– a planet orbiting close to its host star: P_{o} < 10 days, where P_{o} is the orbital period of the planet, so that stellar tidal effects are important.
Under these conditions the chosen systems are mostly hotJupiterlike systems, with host stars ranging from 0.7 to 1.4 M_{⊙} in order to have a similar structure to that of the Sun (namely with a convective envelope and a radiative region below it during the MS). They have been picked out using the Extrasolar Planets Encyclopaedia^{1} (e.g. Schneider et al. 2011) which presents the orbital period and eccentricity of the planet, along with the TEPcat^{2} (Southworth 2011) database to find the skyprojected obliquity. Then, we removed the systems for which the age of the star was not known. For the remaining systems, stellar period has been found in the literature: the related references are reported in the last column of Table B.1.
In Fig. 5, we present the ratio of the magnetic to hydrodynamic forcings as a function of the mass of the host star. The quantity has been calculated using orbital and stellar rotation periods (see Table B.1) for the pair (m, n) that minimises this quantity and thus maximises f_{mag}/f_{hydro}. To calculate the Lehnert number squared, we have used the star’s rotation period and radius displayed in Table B.1, coupled with density, convective length and velocity given by our grid of STAREVOL models at the closest age and mass of the host star. For each system, Le^{2} is evaluated at the base and the top of the CZ. Moreover, the magnetostrophic regime has been selected because it best reproduces the observed surface magnetic fields (see Figs. 1 and 2). As expected, the ratio f_{mag}/f_{hydro} at the base is greater than at the top of the CZ, consistent with the relative magnitude of the largescale magnetic field inside the convective envelope. Moreover, at the base of the CZ the higher the mass, the greater this ratio. On the contrary, we observe a drop in f_{mag}/f_{hydro} at the top of the CZ for stars more massive than 1.2 M_{⊙}, similar to what we notice in Fig. 4. More importantly, we point out that f_{mag}/f_{hydro} is always smaller than unity, regardless of the mass of star and the position inside the stellar envelope. Only HATP24 (b) may feature the ratio of forcings around unity within the error bar for (m, n) = (2, 1) at the base of the CZ. Indeed, the rotation period of the star HATP24 is nearly twice the orbital period of the planet HATP24 b, which implies that Ro_{t} is very close to zero, thereby making f_{mag}/f_{hydro} close to unity, despite the small value of Le^{2}. It is worth noticing that Ogilvie (2009) and Lin & Ogilvie (2018) mentioned a significant impact of the magnetic field on nonwave like motions as soon as Le > 0.1. In our study, some starplanet systems (like HATP13 and HD 149026) feature a Lehnert number greater than 10^{−1} and yet a ratio f_{mag}/f_{hydro} far below unity. This stresses the fact that the Lehnert number is not the only quantity to come into play and that the relative amplitude of the tidal and rotation frequencies has to be considered to conclude on the impact of magnetism on the tidal forcing.
Fig. 5. Dependence of the ratio of the Lorentz forcing relative to the hydrodynamic forcing with respect to the mass of the host star for the selected starplanet systems (see Table B.1). The magnetic forcing term f_{mag} is estimated with the magnetostrophic balance at the base and the top of the CZ and a stellar model evolved with a slow initial rotation. The tidal frequency has been calculated with the pair (m, n) that minimises and thus maximises f_{mag}/f_{hydro}. 
From the results of this section, we conclude that the tidal forcing arising from the Lorentz force remains small in comparison to a pure hydrodynamical forcing. This conclusion is important, as it stresses that adopting a Coriolisdriven tidal forcing is justified to study the propagation and dissipation of tidallyforced magnetoinertial waves in the convective envelope of lowmass stars, despite the presence of a largescale, dynamo generated magnetic field, as was done in Wei (2016, 2018) and Lin & Ogilvie (2018).
3. The relative importance of viscous over Ohmic dissipation for magnetoinertial waves
In the previous section, we have demonstrated that the Lorentz force has a weak contribution to the tidal forcing of (magneto) inertial waves. It can, however, affect their propagation (hence the name magnetoinertial waves) and dissipation as the Lorentz force acts on the wavelike part of the equation of motion and we have introduced the Ohmic diffusion in the induction equation (see Eqs. (5) and (4b), respectively). Wei (2016) studied the dissipation of these waves by turbulent friction and magnetic diffusion processes using a local Cartesian model of an isentropic convective region. He compared the importance of Ohmic versus viscous dissipations, especially at resonances. Varying the Lehnert number, Wei (2016) found that the transition from a viscousdominated regime to a regime dominated by Ohmic dissipation occurs when the Lehnert number is greater than ∼10^{−4}–10^{−3}, for average atomic Ekman and Prandtl numbers close to those expected in the Sun or in Jupiter. This led him to the conclusion that when Le is larger than 10^{−3}, magnetic effects on tidal dissipation should be taken into account.
This work has been taken up by Lin & Ogilvie (2018) in which they studied the propagation and the kinetic and magnetic energy dissipations of tidally forced magnetoinertial waves in a spherical shell. They showed that, at high Lehnert numbers, dissipation is no longer focused along the shear layers that are shaped by rotation and viscosity following attractors of characteristics as in the pure hydrodynamical case. Once again, they identified a critical Lehnert number that separates the regimes dominated by viscosity or by Ohmic dissipation. They expressed this critical Lehnert number with the help of characteristic timescales of ray tracing. Specifically, the width of an inertial wave beam ℓ is deduced by equating the magnetic diffusion timescale τ_{η} = ℓ^{2}/η and the inertial wave propagation time τ_{i} = R/V_{g}, where V_{g}∼ℓΩ is the group velocity (see Fig. 6). Indeed, the higher the Ohmic diffusion (or the viscosity), the larger the spread of the inertial wave beam. Furthermore, hydrodynamical terms prevail over MHD ones when Alfvén waves (produced by the deformation of the magnetic field by the inertial flow) do not have time to distort the wave beam. In other words the hydrodynamical terms dominate when τ_{a} > τ_{i} where
is the typical time for an Alfvén wave to transversely cross the wave beam (see Fig. 6). Using these heuristic considerations, Lin & Ogilvie (2018) have shown that the propagation of inertial waves is little influenced by a magnetic field as long as
with Em = η/(2ΩR^{2}) the magnetic Ekman number. This prediction has been inferred in the context of a low Pm. This condition is generally satisfied in solarlike stars. In the Sun, the atomic magnetic Prandtl number varies from 10^{−6} at the surface to 10^{−1} at the base of the CZ (Zeldovich 1983). When viscosity dominates, for instance in the core of massive stars (see Fig. 2 in Augustson et al. 2019), the same relationship holds with the viscous Ekman number Ek instead of Em.
Fig. 6. Sketch of an inertial wave beam with a magnetic field (red arrow), an Alfvén wave (dashed red line) and an inertial wave propagating from the left to the right (blue arrows). The time for an Alfvén wave to transversely cross the wave beam of length ℓ(τ_{i} = τ_{η}), is greater than the time for an inertial wave to go through the wave beam of length L ∼ R. 
In Fig. 7, the Lehnert number is illustrated as a function of stellar age, where the threshold defined in Lin & Ogilvie (2018) is included with atomic and turbulent magnetic diffusivities: and being the atomic and turbulent magnetic Ekman numbers, respectively. The parameter Em_{ato} has been computed thanks to the Braginskii prescription for plasma diffusivities (Appendix B, Braginskii 1965; Augustson et al. 2019) using the grid of STAREVOL models. The turbulent magnetic Ekman number (Em_{tur}) is derived assuming that the eddymagnetic diffusivity takes the simple form: η_{tur} = u_{c}l_{c}/3. In this approach, η_{tur} is equivalent to the eddyviscosity that means a magnetic turbulent Prandtl number of the order of unity (see e.g. Camargo & Tasso 1992). We highlight that in both panels of Fig. 7, the threshold calculated with an atomic magnetic diffusivity is much lower than Le derived in the various regimes, by at least three order of magnitudes near the surface and six order of magnitude at the base of the CZ. Consequently, the Ohmic dissipation largely outbalances the viscous dissipation along the lifetime of a 1 M_{⊙} when taking Lin and Ogilvie’s threshold assessed with a magnetic atomic diffusivity η_{ato}. In contrast, the limit estimated with the turbulent magnetic Ekman number is of the same order of magnitude as the Lehnert number in the equipartition regime at the base of the CZ, whereas the threshold is slightly smaller than the Lehnert number curves at the top. Note that the turbulent magnetic Ekman number is close to the value of the magnetic Ekman number used in Lin & Ogilvie’s (previously quoted) paper for their simulations. This explains why the transition from an hydrodynamical to a fully magnetic regime is carried out at similar Lehnert number in their case and in ours when using Em_{tur} here (see Figs. 2–4 in Lin & Ogilvie 2018). By choosing an eddymagnetic diffusivity, both Ohmic and viscous dissipations have to be taken into account in the dissipation calculation.
Fig. 7. Lehnert number at the base (left panel) and top (right panel) of the CZ versus the age of a 1 M_{⊙} star. Slow initial rotation has been used in both panels. The thresholds for which magnetic field impacts the propagation of inertial waves is estimated with different Ekman numbers and drawn in orange (see the legend). 
4. The impact of a smallerscale magnetic field
The stellar dynamo is a multiscale process. Indeed, largescale and smallscale dynamos coexist inside a star’s CZ and produce magnetic fields at different length scales (Brandenburg & Subramanian 2005; Strugarek et al. 2013, and references therein). At the base of the CZ, the convective length scale computed by STAREVOL is about one tenth of the radius of the star and decreases drastically towards the top of the CZ (∼10^{−4} R in our model). In Sect. 2.2, we have made the assumption that a dynamolike magnetic field is the result of turbulent convective motions featured by the convective velocity at the base of the CZ, associated to a relatively largescale convective length. Then, we have used this dynamoinduced magnetic field to extrapolate a dipolar magnetic field near the surface of the star (see Eq. (10)). However, one can question the role of the smallscale dynamo fields on tidal excitation and dissipation throughout the convective envelope.
The dynamical tides will possibly interact with the smallerscale magnetic field. Indeed, the scale of variation of dynamical tides along the inertial wave beam is of order R (Sect. 3), but the transverse length scale l of this beam is smaller. The balance between magnetic diffusion and inertial timescales (see also Sect. 3) leads to i.e. the transverse length scale of the beam varies from one hundredth to one tenth of the stellar radius when using typical values of the turbulent Ekman number (see Fig. 7). This orthogonal length scale is thus comparable to the length scale of the convection and of the corresponding magnetic field.
The scale of variation of the equilibrium tide is also of order R (see e.g. Remus et al. 2012). The smallscale components of the star’s magnetic field can collectively affect the largescale flow of the equilibrium tide provided their correlations sustain a largescale effective Lorentz force. We plot Le^{2} associated with a smallscale magnetic field in Fig. 8 as a function of the normalised radius in the whole CZ. The typical length to be used in the expression of f_{mag} (Eq. (6)) is no longer R but l_{c} which represents the characteristic length of the fluctuating component of the dynamogenerated magnetic field. We redefine the Lehnert number here as , where B_{dyn} (Table 1) depends upon l_{c}, which in turn varies inside the CZ. In the three regimes listed in Table 2, Le now depends only on the Rossby number (i.e. each scaling has to be multiplied by R/l_{c}). The profile of Le^{2} is shown for magnetostrophic and equipartition regimes at different evolutionary stages. In both regimes, the Lehnert number follows the Rossby number tendency obtained in Mathis et al. (2016). Indeed, the authors pointed out that Ro always increases with radius, regardless of the changes in mass and stellar evolutionary phases. In the early stages of the 1 M_{⊙} star’s evolution (PMS and ZAMS), the probability of having a significant magnetic forcing e.g. f_{mag}/f_{hydro} > 1 is weak at the base of the CZ (Le^{2} < 10^{−3}), but becomes high near the surface (Le^{2} > 10) based on the analyses of Sects. 2.5 and 2.6. The chances are even higher when the star is older (from the midpoint to the end of the MS) since Le^{2} is already greater than 10^{−3} at the base of the CZ in both regimes.
Fig. 8. Lehnert number squared against the normalised radius in the magnetostrophic (dark red) and equipartition (green) regimes, for a 1 M_{⊙} star. The reader is referred to the body of Sect. 4 for the new definition of Le that involves a smallscale magnetic field. The Lehnert number squared is plotted at different epochs: amid the premain sequence (midPMS; ∼30 Myr), at zeroage main sequence (ZAMS; ∼60 Myr), amid the main sequence (midMS; ∼5.3 Gyr), and toward the end of the main sequence (TAMS; ∼10.5 Gyr). 
From this analysis, we should keep in mind that smallscale magnetic fields can affect the equilibrium and dynamical tides. Indeed, Le^{2} can be significantly enhanced near the surface of the star, a fortiori for all three regimes described in Table 1. This may change the results of Sect. 2.6, namely boost the effect of magnetism on tidal forcing so that f_{mag} is no longer negligible in front of f_{hydro}.
5. Conclusions and discussions
The influence of magnetism on tidal interaction along the evolution of lowmass stars has been investigated through its impact on tidal excitation and dissipation. We have first derived an analytic criterion to quantify the Lorentz contribution to the tidal forcing. The relative importance of Coriolislike versus Lorentzlike forcings relies on the Lehnert number and characteristic frequencies of starplanet systems, i.e. the tidal frequency and the rotation frequency of the main body. With the help of simple scaling laws, we have estimated the magnitude of a dynamogenerated magnetic field near the tachocline. Then, a largescale dipolar magnetic field has been inferred near the surface of the star from these scaling relationships. For this purpose, we have used grids of lowmass stars computed by the stellar evolution code STAREVOL (Amard et al. 2019). This allowed us to compare our magnetic field estimates at the surface with the observations of the surface dipolar magnetic field. We find that the magnetostrophic regime (resulting from the balance between Coriolis and Lorentz forces) gives a good estimate, within an order of magnitude, of the magnetic fields of PMS and MS stars when assuming a median or fast initial rotation. Nevertheless, improvements still need to be made in terms of magnetic scaling laws in order to better reproduce the observed magnetic fields of Vidotto et al. (2014) and See et al. (2017) (see Appendix B).
Subsequently, we have estimated the Lehnert number with these prescriptions for the stellar dynamo at the base and the top of the CZ. We have explored how it varies for different stellar masses and initial rotation speeds as a function of time. In all cases, the Lehnert number is small compared to unity. This means that magnetism can play a role in tidal forcing only for small tidal to rotational frequency ratios. In practice, this condition is verified neither at the base nor at the top of the CZ when using the largescale magnetic fields in the primary component of quasicircular and coplanar starplanet systems. Although the ratio f_{mag}/f_{hydro} increases with the mass of the host star until 1.2 M_{⊙} (see Fig. 5) and is greater at the base of the CZ, it remains below 1 for all the studied systems. Indeed, only specific and sharp enough resonances between the tidal and the rotational frequency can increase the f_{mag}/f_{hydro} ratio. This statement is also valid for noncircular and noncoplanar systems, the tidal frequency being just a different combination of orbital and rotation frequencies. Eventually, one should note that all the selected systems (Table B.1) have a host star on the mid and late MS. But it should not change the fact that largescale magnetic fields have little influence on the tidal excitation because Le^{2} is even weaker for younger stars (as we have seen in Fig. 4).
In contrast, the dissipation of the dynamical tide is strongly affected by stellar magnetism. Indeed, we have shown that viscous dissipation is no longer the main process of energy dissipation. Ohmic dissipation is at least as important as viscous dissipation (see Fig. 7) in the whole CZ for a 1 M_{⊙} star as well as other lowmass stars. This means that a full MHD treatment is needed to analyse the propagation and dissipation of tidallyforced (magneto) inertial waves inside the CZ of a lowmass star.
It can be added that whatever the energy mix distributed in the toroidal and poloidal components of the largescale magnetic field at the surface (which is also variable over time Kochukhov et al. 2017), it does not change the main conclusions of our paper because these energies are comparable within an order of magnitude (See et al. 2015).
In Sect. 4, we have addressed the question of the effect of smallscale components of the dynamoinduced magnetic field on tidal interactions. We have pointed out that they are likely to play a role on both dynamical and equilibrium tides and consequently on the tidal forcing and dissipation. We first demonstrated that the transverse component of the typical length scale of the dynamical tide is commensurate with the convective mixing length in the bulk of the CZ. In this respect, it would also mean that the modelling of the friction applied by turbulent convection on the dynamical tides should go beyond an eddy viscosity model (Ogilvie & Lesur 2012) to incorporate magnetic effects. Then, we showed that the Lehnert number squared estimated with a smallscale magnetic field can be much greater than Le^{2} calculated with a largescale field. As a result, the impact of a smallscale dynamo induced magnetic field on tidal forcing could be not as negligible as the effect of a largerscale magnetic field. Quantifying this impact now requires abinitio modelling of tidal forcing in turbulent convective envelopes sustaining a smallscale dynamo.
Acknowledgments
We would like to thank the anonymous referee for the helpful comments and suggestions regarding our work. A. Astoul, K. Augustson, E. Bolmont, and S. Mathis acknowledge funding by the European Research Council through the ERC grant SPIRE 647383. The authors acknowledge the PLATO CNES funding at CEA/IRFU/DAp and IRAP. The authors further thank V. See for fruitful discussions and the use of his data. F. Gallet acknowledges financial support from a CNES fellowship. A. S. Brun acknowledges funding by ERC WHOLESUN 810218 grant, INSU/PNST, and CNES Solar Orbiter. This work has been carried out within the framework of the NCCR PlanetS supported by the Swiss National Science Foundation. This research has made use of NASA’s Astrophysics Data System.
References
 Amard, L., Palacios, A., Charbonnel, C., Gallet, F., & Bouvier, J. 2016, A&A, 587, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Amard, L., Palacios, A., Charbonnel, C., et al. 2019, A&A, 631, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Anderson, D. R., Barros, S. C. C., Boisse, I., et al. 2011, PASP, 123, 555 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, D. R., Brown, D. J. A., Collier Cameron, A., et al. 2014, MNRAS, submitted [arXiv:1410.3449] [Google Scholar]
 Antia, H. M., Chitre, S. M., & Thompson, M. J. 2000, A&A, 360, 335 [NASA ADS] [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Augustson, K. C., Brun, A. S., & Toomre, J. 2019, ApJ, 876, 83 [NASA ADS] [CrossRef] [Google Scholar]
 Barker, A. J., & Ogilvie, G. I. 2010, MNRAS, 404, 1849 [NASA ADS] [Google Scholar]
 Basu, S. 2016, Liv. Rev. Sol. Phys., 13, 2 [Google Scholar]
 Beck, P. G., Mathis, S., Gallet, F., et al. 2018, MNRAS, 479, L123 [NASA ADS] [CrossRef] [Google Scholar]
 Bell, C. P. M., Naylor, T., Mayne, N. J., Jeffries, R. D., & Littlefair, S. P. 2013, MNRAS, 434, 806 [NASA ADS] [CrossRef] [Google Scholar]
 Benbakoura, M., Réville, V., Brun, A. S., Le PoncinLafitte, C., & Mathis, S. 2019, A&A, 621, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bolmont, E., Gallet, F., Mathis, S., et al. 2017, A&A, 604, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bolmont, E., & Mathis, S. 2016, Celest. Mech. Dyn. Astron., 126, 275 [Google Scholar]
 Braginskii, S. I. 1965, Rev. Plasma Phys., 1, 205 [Google Scholar]
 Brandenburg, A., & Subramanian, K. 2005, Phys. Rep., 417, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Brown, D. J. A. 2014, MNRAS, 442, 1844 [NASA ADS] [CrossRef] [Google Scholar]
 Brun, A. S., & Browning, M. K. 2017, Liv. Rev. Sol. Phys., 14, 4 [Google Scholar]
 Brun, A. S., García, R. A., Houdek, G., Nandy, D., & Pinsonneault, M. 2015, Space Sci. Rev., 196, 303 [NASA ADS] [CrossRef] [Google Scholar]
 Camargo, S. J., & Tasso, H. 1992, Phys. Fluids B, 4, 1199 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Charbonneau, P. 2013, Solar and Stellar Dynamos: SaasFee Advanced Course 39 Swiss Society for Astrophysics and Astronomy (SpringerVerlag: Berlin, Heidelberg), 39 [Google Scholar]
 Charbonneau, P. 2014, ARA&A, 52, 251 [Google Scholar]
 Christensen, U. R., Holzwarth, V., & Reiners, A. 2009, Nature, 457, 167 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Cowling, T. G. 1941, MNRAS, 101, 367 [NASA ADS] [Google Scholar]
 Czesla, S., Salz, M., Schneider, P. C., Mittag, M., & Schmitt, J. H. M. M. 2017, A&A, 607, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Davidson, P. A. 2013, Geophys. J. Int., 195, 67 [NASA ADS] [CrossRef] [Google Scholar]
 DeRosa, M. L., Brun, A. S., & Hoeksema, J. T. 2012, ApJ, 757, 96 [Google Scholar]
 Donati, J. F., Catala, C., Landstreet, J. D., & Petit, P. 2006, in ASP Conf. Ser., ed. R. Casini, 358, 362 [Google Scholar]
 Donati, J.F., & Landstreet, J. D. 2009, ARA&A, 47, 333 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Donati, J.F., Jardine, M. M., Gregory, S. G., et al. 2007, MNRAS, 380, 1297 [NASA ADS] [CrossRef] [Google Scholar]
 Duez, V., Mathis, S., & TurckChièze, S. 2010, MNRAS, 402, 271 [NASA ADS] [CrossRef] [Google Scholar]
 EmeriauViard, C., & Brun, A. S. 2017, ApJ, 846, 8 [NASA ADS] [CrossRef] [Google Scholar]
 Esposito, M., Covino, E., Desidera, S., et al. 2017, A&A, 601, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 FerrazMello, S., Tadeu dos Santos, M., Folonier, H., et al. 2015, ApJ, 807, 78 [NASA ADS] [CrossRef] [Google Scholar]
 Finlay, C. C. 2008, in Dynamos, eds. P. Cardin, & L. Cugliandolo (Amsterdam: Elsevier), Les Houches, 88, 403 [CrossRef] [Google Scholar]
 Folsom, C. P., Petit, P., Bouvier, J., et al. 2016, MNRAS, 457, 580 [NASA ADS] [CrossRef] [Google Scholar]
 Gallet, F., & Bouvier, J. 2013, A&A, 556, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gallet, F., & Bouvier, J. 2015, A&A, 577, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gallet, F., Bolmont, E., Mathis, S., Charbonnel, C., & Amard, L. 2017, A&A, 604, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gandolfi, D., Hébrard, G., Alonso, R., et al. 2010, A&A, 524, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ge, J., van Eyken, J., Mahadevan, S., et al. 2006, ApJ, 648, 683 [NASA ADS] [CrossRef] [Google Scholar]
 Gilliland, R. L. 1986, ApJ, 300, 339 [NASA ADS] [CrossRef] [Google Scholar]
 Goldreich, P., & Nicholson, P. D. 1989, ApJ, 342, 1079 [NASA ADS] [CrossRef] [Google Scholar]
 Goodman, J., & Dickson, E. S. 1998, ApJ, 507, 938 [NASA ADS] [CrossRef] [Google Scholar]
 Gough, D. O., & Thompson, M. J. 1990, MNRAS, 242, 25 [Google Scholar]
 Granata, V., Nascimbeni, V., Piotto, G., et al. 2014, Astron. Nachr., 335, 797 [NASA ADS] [CrossRef] [Google Scholar]
 Guenel, M., Baruteau, C., Mathis, S., & Rieutord, M. 2016a, A&A, 589, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Guenel, M., Mathis, S., Baruteau, C., & Rieutord, M. 2016b, ArXiv eprints [arXiv:1612.05071] [Google Scholar]
 Hébrard, G., Collier Cameron, A., Brown, D. J. A., et al. 2013, A&A, 549, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hirano, T., SanchisOjeda, R., Takeda, Y., et al. 2012, ApJ, 756, 66 [NASA ADS] [CrossRef] [Google Scholar]
 Hut, P. 1980, A&A, 92, 167 [NASA ADS] [Google Scholar]
 Janson, M., Reffert, S., Brandner, W., et al. 2008, A&A, 488, 771 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jardine, M., & Unruh, Y. C. 1999, A&A, 346, 883 [NASA ADS] [Google Scholar]
 Jurčišinová, E., Jurčişin, M., Remecký, R., & Zalom, P. 2013, Phys. Rev. E, 87, 043010 [NASA ADS] [CrossRef] [Google Scholar]
 Käpylä, P. J., Rheinhardt, M., Brandenburg, A., & Käpylä, M. J. 2019, A&A, submitted [arXiv:1901.00787] [Google Scholar]
 Kochukhov, O., Petit, P., Strassmeier, K. G., et al. 2017, Astron. Nachr., 338, 428 [NASA ADS] [CrossRef] [Google Scholar]
 Lagarde, N., Decressin, T., Charbonnel, C., et al. 2012, A&A, 543, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Landin, N. R., Mendes, L. T. S., & Vaz, L. P. R. 2010, A&A, 510, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lanza, A. F., Damiani, C., & Gandolfi, D. 2011, A&A, 529, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lehnert, B. 1954, ApJ, 119, 647 [NASA ADS] [CrossRef] [Google Scholar]
 Lin, Y., & Ogilvie, G. I. 2018, MNRAS, 474, 1644 [NASA ADS] [CrossRef] [Google Scholar]
 Line, M. R., Teske, J., Burningham, B., Fortney, J. J., & Marley, M. S. 2015, ApJ, 807, 183 [NASA ADS] [CrossRef] [Google Scholar]
 Mancini, L., Esposito, M., Covino, E., et al. 2015, A&A, 579, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mancini, L., Esposito, M., Covino, E., et al. 2018, A&A, 613, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S. 2015, A&A, 580, L3 [Google Scholar]
 Mathis, S., & Le PoncinLafitte, C. 2009, A&A, 497, 889 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., AuclairDesrotour, P., Guenel, M., Gallet, F., & Le PoncinLafitte, C. 2016, A&A, 592, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Matt, S. P., Brun, A. S., Baraffe, I., Bouvier, J., & Chabrier, G. 2015, ApJ, 799, L23 [NASA ADS] [CrossRef] [Google Scholar]
 Maxted, P. F. L., Serenelli, A. M., & Southworth, J. 2015, A&A, 577, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mayor, M., Queloz, D., Udry, S., & Halbwachs, J. L. 1997, in IAU Colloq. 161: Astronomical and Biochemical Origins and the Search for Life in the Universe, eds. C. Batalli Cosmovici, S. Bowyer, & D. Werthimer, 313 [Google Scholar]
 MohlerFischer, M., Mancini, L., Hartman, J. D., et al. 2013, A&A, 558, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ogilvie, G. I. 2005, J. Fluid Mech., 543, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Ogilvie, G. I. 2009, MNRAS, 396, 794 [NASA ADS] [CrossRef] [Google Scholar]
 Ogilvie, G. I. 2013, MNRAS, 429, 613 [NASA ADS] [CrossRef] [Google Scholar]
 Ogilvie, G. I. 2014, ARA&A, 52, 171 [Google Scholar]
 Ogilvie, G. I., & Lesur, G. 2012, MNRAS, 422, 1975 [NASA ADS] [CrossRef] [Google Scholar]
 Ogilvie, G. I., & Lin, D. N. C. 2004, ApJ, 610, 477 [NASA ADS] [CrossRef] [Google Scholar]
 Ogilvie, G. I., & Lin, D. N. C. 2007, ApJ, 661, 1180 [NASA ADS] [CrossRef] [Google Scholar]
 Perryman, M. 2018, The Exoplanet Handbook (Cambridge: Cambridge University Press) [CrossRef] [Google Scholar]
 Rebull, L. M., Wolff, S. C., & Strom, S. E. 2004, AJ, 127, 1029 [NASA ADS] [CrossRef] [Google Scholar]
 Reiners, A. 2012, Liv. Rev. Sol. Phys., 9, 1 [Google Scholar]
 Remus, F., Mathis, S., & Zahn, J. P. 2012, A&A, 544, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Schlichting, H. E. 2014, ApJ, 795, L15 [NASA ADS] [CrossRef] [Google Scholar]
 Schneider, J., Dedieu, C., Le Sidaner, P., Savalle, R., & Zolotukhin, I. 2011, A&A, 532, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 See, V., Jardine, M., Vidotto, A. A., et al. 2015, MNRAS, 453, 4301 [NASA ADS] [CrossRef] [Google Scholar]
 See, V., Jardine, M., Vidotto, A. A., et al. 2017, MNRAS, 466, 1542 [NASA ADS] [CrossRef] [Google Scholar]
 See, V., Matt, S. P., Folsom, C. P., et al. 2019, ApJ, 876, 118 [NASA ADS] [CrossRef] [Google Scholar]
 Skumanich, A. 1972, ApJ, 171, 565 [NASA ADS] [CrossRef] [Google Scholar]
 Southworth, J. 2011, MNRAS, 417, 2166 [NASA ADS] [CrossRef] [Google Scholar]
 Spiegel, E. A., & Zahn, J. P. 1992, A&A, 265, 106 [NASA ADS] [Google Scholar]
 Strugarek, A., Brun, A. S., Mathis, S., & Sarazin, Y. 2013, ApJ, 764, 189 [NASA ADS] [CrossRef] [Google Scholar]
 Strugarek, A., Beaudoin, P., Charbonneau, P., Brun, A. S., & Nascimento, J. D. 2017, Science, 357, 185 [NASA ADS] [CrossRef] [Google Scholar]
 Terquem, C., Papaloizou, J. C. B., Nelson, R. P., & Lin, D. N. C. 1998, ApJ, 502, 788 [NASA ADS] [CrossRef] [Google Scholar]
 TregloanReed, J., Southworth, J., & Tappert, C. 2013, MNRAS, 428, 3671 [NASA ADS] [CrossRef] [Google Scholar]
 TregloanReed, J., Southworth, J., Burgdorf, M., et al. 2015, MNRAS, 450, 1760 [NASA ADS] [CrossRef] [Google Scholar]
 Vidotto, A. A., Gregory, S. G., Jardine, M., et al. 2014, MNRAS, 441, 2361 [NASA ADS] [CrossRef] [Google Scholar]
 Vilhu, O. 1984, A&A, 133, 117 [NASA ADS] [Google Scholar]
 Weber, E. J., & Davis, Jr., L. 1967, ApJ, 148, 217 [NASA ADS] [CrossRef] [Google Scholar]
 Wei, X. 2016, ApJ, 828, 30 [Google Scholar]
 Wei, X. 2018, ApJ, 854, 34 [Google Scholar]
 Winn, J. N., Noyes, R. W., Holman, M. J., et al. 2005, ApJ, 631, 1215 [NASA ADS] [CrossRef] [Google Scholar]
 Witte, M. G., & Savonije, G. J. 2002, A&A, 386, 222 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wright, N. J., & Drake, J. J. 2016, Nature, 535, 526 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Zahn, J. P. 1966a, Ann. Astrophys., 29, 313 [NASA ADS] [Google Scholar]
 Zahn, J. P. 1966b, Ann. Astrophys., 29, 489 [NASA ADS] [Google Scholar]
 Zahn, J. P. 1966c, Ann. Astrophys., 29, 565 [NASA ADS] [Google Scholar]
 Zahn, J. P. 1975, A&A, 41, 329 [NASA ADS] [Google Scholar]
 Zahn, J. P. 1977, A&A, 500, 121 [NASA ADS] [Google Scholar]
 Zahn, J. P. 1989, A&A, 220, 112 [Google Scholar]
 Zahn, J. P., & Bouchet, L. 1989, A&A, 223, 112 [NASA ADS] [Google Scholar]
 Zanni, C., & Ferreira, J. 2013, A&A, 550, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zeldovich, Y. B. 1983, Magnetic Fields in Astrophysics (New York: Gordon and Breach) [Google Scholar]
Appendix A: Magnetic field scaling laws
Several estimates of the magnetic field of stars can be derived thanks to energetic or force balances in the equations for heat transport (items ii and iii below) or for momentum (items i, and also ii and iii). Based on the MHD equations, and assuming stationarity, the NavierStokes equation of convective motions is
where we have introduced g the gravitational acceleration of the star. At the base of the CZ, we assume that the gradient and the velocity scale like and u ≈ u_{c}.

i
Magnetostrophy
In this regime, the Lorentz force balances the Coriolis force i.e. (3)≈(1) in Eq. (A.1). Thus,
where is the kinetic energy density associated with convective motions, and Ro = u_{c}/(2Ωl_{c}) the convective Rossby number.

ii
Equipartition regime
In this regime we assume that the dynamo is efficient enough such that the convective kinetic energy density KE of the fluid is equivalent to the magnetic energy density ME = B^{2}/(2μ_{0}) i.e. (3)≈(2) in Eq. (A.1):

iii
Buoyancy dynamo
Taking the curl of Eq. (A.1) and neglecting viscous and inertial forces we obtain
assuming that the velocity is divergencefree. The buoyancy dynamo lies on the fact that each term in Eq. (A.4) has the same order of magnitude (Davidson 2013; Augustson et al. 2019). If the Rossby number is small enough, namely the rotation frequency of the body is larger than the convective frequency (the inverse of the convective turnover time), the convective flows become almost columnar, and aligned along the rotation axis. In the stellar evolution grids we use, this assumption is at least verified in the first half of the CZ from the base. As a consequence, two characteristic length scales appear, l_{∥} and l_{⊥}: these lengths are parallel and perpendicular to the rotation axis, respectively. A direct result is that l_{⊥} < l_{∥}, and we assume that l_{c} ≃ l_{∥} so that the convective Rossby number is now Ro = u_{c}/(2Ωl_{∥}). It seems consistent with the high values of the convective length scale we find at the base of the envelope for different stars (several tenths as the CZ thickness). Thus, taking the characteristic scales of each quantities in Eq. (A.4):
where by definition of l_{∥}, otherwise because l_{⊥} < l_{∥}. Then, we set up p_{g} = gu_{c} the buoyancy power (or rate of buoyancy work) per unit of mass. We can also express the magnetic energy in term of buoyancy power because of the equivalence of forces in Eq. (A.4). Therefore, in order to respect dimensional homogeneity, the magnetic energy per unit of mass satisfies the following relationship:
where we are only using l_{∥}, again because l_{∥} > l_{⊥} since we are looking for the largescale dynamo magnetic field. From the equivalences in Eq. (A.5), we can express g as a function of Ω, u_{c}, l_{∥} and the energy by unit of mass B^{2}/(ρμ_{0}):
Injecting Eq. (A.7) in , and all of this in Eq. (A.6), one can find:
Considering that Ro = u_{c}/(2Ωl_{∥}), Eq. (A.8) can be written finally in the form:
Appendix B: Limitations of the best scaling for
In Figs. 1 and 2, we observe that B_{dip} calculated with the magnetostrophic regime and a fast initial rotation at the top of the CZ is in adequacy with the observed dipolar magnetic field at the surface. We present here some limitations of this comparison. In Fig. B.1, we display the angular velocity (and the period in the vertical right axis) of the 0.9 M_{⊙} stars featured in Fig. 1 against their age, compared to the different rotational evolutions given by STAREVOL. The angular velocities are found in See et al. (2017) and references therein. They are obtained mainly through photometry while the age of the stars are mostly determined by the comparison between isochrones and colourmagnitude diagrams. Note that the error bars of the angular velocities are too small to be visible on the graph. We point out that the rotation periods of the stars seem to be in agreement with the slow or median rotation profiles provided by STAREVOL, unlike what we observed for the surface dipolar magnetic fields in Fig. 1. To some extent, the same assessment is made when displaying rotation versus time diagram for other lowmass stars (not presented here). There are several possible explanations for this result. The discrepancy could be first attributed to the simplicity of our dynamoinduced magnetic fields models. Then, it is possible that the metallicity affects the rotation evolution models as it has been demonstrated by Amard et al. (in prep.) (we remind the reader that the metallicity is fixed in Sect. 2.3 to the solar value). Furthermore, as we said in Sect. 2.4, this observation tends to be more nuanced with increasing data. As a result, we have to cautiously interpret the fact that B_{dip} estimated with the magnetostrophic regime and a fast initial rotation match the observed values of the dipolar magnetic field at the surface.
Fig. B.1. Rotation frequency on the vertical left axis (and period on the vertical right axis) versus the age of 0.9 M_{⊙} stars (), given fast, median, and slow initial rotations. Ages, stellar periods and their associated error bars can be found in the papers quoted in Fig. 1. 
Starplanet systems with a quasicircular and quasicoplanar orbit, e < 0.1 and λ < 30° with e and λ the eccentricity and the skyprojected spinorbit misalignment angle, respectively.
All Tables
Magnetic fields derived from simple balances (forces or energies) written in the second column.
Starplanet systems with a quasicircular and quasicoplanar orbit, e < 0.1 and λ < 30° with e and λ the eccentricity and the skyprojected spinorbit misalignment angle, respectively.
All Figures
Fig. 1. Magnetic field versus the age of a 0.9 M_{⊙} star for three different initial rotation periods. At the base of the CZ (left panel), the magnetic field is estimated from the scaling laws listed in Table 1. At the top of the CZ (right panel), the same scaling laws are used to extrapolate a dipolar magnetic field near the surface (see Eq. (10)). The “stars” symbols () that represent mean dipolar magnetic fields at the surface of 0.9 M_{⊙} stars are taken from See et al. (2017), and the age of the stars are from Ge et al. (2006), Janson et al. (2008), Line et al. (2015), and Folsom et al. (2016). Multiple observations of the dipolar magnetic field of a star are joined by a line. 

In the text 
Fig. 2. Ratio of the mean unsigned observed dipolar magnetic field (See et al. 2017; Vidotto et al. 2014) over estimated dipolar magnetic field as a function of the mass of the star. The age of the stars is indicated in colour and ranges from 21 Myr to 17.2 Gyr. Magnetostrophic regime and median initial rotation have been used to plot the ratio . Multiple observations of the dipolar magnetic field of a star are again joined by a line. A typical error bar for M_{⋆} and is indicated in grey. We have adopted a conservative error estimate of 0.434 dex in log B_{dip} and 0.1 M_{⊙} in M_{⋆}. 

In the text 
Fig. 3. Lehnert number squared at the base (left panel) and the top (right panel) of the CZ as a function of age for a 1 M_{⊙} star, and for different magnetic scaling laws and initial rotations (see the legends in the left panel). 

In the text 
Fig. 4. Evolution of the Lehnert number squared over time, at the base (solid lines) and the top (dotted dashed lines) of the CZ, for various lowmass stars sorted from 0.7 to 1.4 M_{⊙}. Slow initial rotation has been used here and a magnetostrophic balance has been assumed. 

In the text 
Fig. 5. Dependence of the ratio of the Lorentz forcing relative to the hydrodynamic forcing with respect to the mass of the host star for the selected starplanet systems (see Table B.1). The magnetic forcing term f_{mag} is estimated with the magnetostrophic balance at the base and the top of the CZ and a stellar model evolved with a slow initial rotation. The tidal frequency has been calculated with the pair (m, n) that minimises and thus maximises f_{mag}/f_{hydro}. 

In the text 
Fig. 6. Sketch of an inertial wave beam with a magnetic field (red arrow), an Alfvén wave (dashed red line) and an inertial wave propagating from the left to the right (blue arrows). The time for an Alfvén wave to transversely cross the wave beam of length ℓ(τ_{i} = τ_{η}), is greater than the time for an inertial wave to go through the wave beam of length L ∼ R. 

In the text 
Fig. 7. Lehnert number at the base (left panel) and top (right panel) of the CZ versus the age of a 1 M_{⊙} star. Slow initial rotation has been used in both panels. The thresholds for which magnetic field impacts the propagation of inertial waves is estimated with different Ekman numbers and drawn in orange (see the legend). 

In the text 
Fig. 8. Lehnert number squared against the normalised radius in the magnetostrophic (dark red) and equipartition (green) regimes, for a 1 M_{⊙} star. The reader is referred to the body of Sect. 4 for the new definition of Le that involves a smallscale magnetic field. The Lehnert number squared is plotted at different epochs: amid the premain sequence (midPMS; ∼30 Myr), at zeroage main sequence (ZAMS; ∼60 Myr), amid the main sequence (midMS; ∼5.3 Gyr), and toward the end of the main sequence (TAMS; ∼10.5 Gyr). 

In the text 
Fig. B.1. Rotation frequency on the vertical left axis (and period on the vertical right axis) versus the age of 0.9 M_{⊙} stars (), given fast, median, and slow initial rotations. Ages, stellar periods and their associated error bars can be found in the papers quoted in Fig. 1. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.