Issue 
A&A
Volume 650, June 2021



Article Number  A53  
Number of page(s)  28  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/202039159  
Published online  07 June 2021 
Magnetic signatures on mixedmode frequencies
I. An axisymmetric fossil field inside the core of red giants
^{1}
AIM, CEA, CNRS, Université ParisSaclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 GifsurYvette, France
^{2}
Flatiron Institute, Simons Foundation, 162 Fifth Ave, New York, NY 10010, USA
email: lbugnet@flatironinstitute.org
^{3}
Instituto de Astrofísica de Canarias, 38200 La Laguna, Tenerife, Spain
^{4}
Universidad de La Laguna, Dpto. de Astrofísica, 38205 La Laguna, Tenerife, Spain
^{5}
University of Exeter, Department of Physics and Astronomy, Stoker Road, Devon, Exeter EX4 4GL, UK
^{6}
LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université de Paris, 5 Place Jules Janssen, 92195 Meudon, France
Received:
11
August
2020
Accepted:
2
March
2021
Context. The discovery of moderate differential rotation between the core and the envelope of evolved solarlike stars could be the signature of a strong magnetic field trapped inside the radiative interior. The population of intermediatemass red giants presenting surprisingly lowamplitude mixed modes (i.e. oscillation modes that behave as acoustic modes in their external envelope and as gravity modes in their core) could also arise from the effect of an internal magnetic field. Indeed, stars more massive than about 1.1 solar masses are known to develop a convective core during their main sequence. The field generated by the dynamo triggered by this convection could be the progenitor of a strong fossil magnetic field trapped inside the core of the star for the remainder of its evolution.
Aims. Observations of mixed modes can constitute an excellent probe of the deepest layers of evolved solarlike stars, and magnetic fields in those regions can impact their propagation. The magnetic perturbation on mixed modes may therefore be visible in asteroseismic data. To unravel which constraints can be obtained from observations, we theoretically investigate the effects of a plausible mixed axisymmetric magnetic field with various amplitudes on the mixedmode frequencies of evolved solarlike stars.
Methods. Firstorder frequency perturbations due to an axisymmetric magnetic field were computed for dipolar and quadrupolar mixed modes. These computations were carried out for a range of stellar ages, masses, and metallicities.
Conclusions. We show that typical fossilfield strengths of 0.1 − 1 MG, consistent with the presence of a dynamo in the convective core during the main sequence, provoke significant asymmetries on mixedmode frequency multiplets during the red giant branch. We provide constraints and methods for the detectability of such magnetic signatures. We show that these signatures may be detectable in asteroseismic data for field amplitudes small enough for the amplitude of the modes not to be affected by the conversion of gravity into Alfvén waves inside the magnetised interior. Finally, we infer an upper limit for the strength of the field and the associated lower limit for the timescale of its action in order to redistribute angular momentum in stellar interiors.
Key words: stars: oscillations / stars: magnetic field / stars: interiors / stars: evolution / stars: rotation
© L. Bugnet et al. 2021
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://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
Oscillations on the surface of solartype stars have been observed and studied across the HertzsprungRussel diagram thanks in large part to the data provided by CoRoT, Kepler, K2, and the TESS missions (e.g. Michel et al. 2008; Chaplin et al. 2010, 2020; Lund et al. 2017; Huber et al. 2019; García & Ballot 2019). While this list is hardly exhaustive, some key results are relevant to this paper and to the importance of the transport of angular momentum in the radiative regions of main sequence (MS) stars. Indeed, the radiative interior of the Sun seems to rotate as a solid body slightly slower than the equatorial surface rotation rate (García et al. 2007) until 0.25 solar radii (R_{⊙}; e.g. Thompson et al. 2003; Couvidat et al. 2003), and the nuclear core may rotate even faster (García et al. 2007). The study of some solarlike stars also shows that they present a nearly solidbody rotation (Benomar et al. 2015). Likewise, in subgiant (SG) stars and redgiant (RG) stars, some relevant works are the discoveries of the unexpectedly slow rotation rate of their cores (e.g. Deheuvels et al. 2012, 2014, 2016; Mosser et al. 2012, 2017; Gehan et al. 2018), and the surprisingly low amplitude of dipolar mixed modes in some red giants (RGs; García et al. 2014a; Mosser et al. 2012, 2017; Stello et al. 2016a). As of yet, there is no clear evolutionary model that yields internal rotation profiles akin to those observed (e.g. Eggenberger et al. 2012, 2017, 2019; Ceillier et al. 2013; Marques et al. 2013, for the loss of angular momentum on the giant branch), nor is there a robust explanation for the observed dipole mode amplitude suppression (Fuller et al. 2015; Cantiello et al. 2016; Lecoanet et al. 2017; Mosser et al. 2017; Loi & Papaloizou 2018).
Some magnetohydrodynamic fluid behaviours can strongly impact the rotation profile of the star both on secular timescales (Eggenberger et al. 2005; Cantiello et al. 2014; Fuller et al. 2014, 2019) and even on their dynamical timescales (Brun et al. 2005; Featherstone et al. 2009; Augustson et al. 2016). Moreover, the transport of chemical species, energy, and angular momentum by internal waves can also play a role (e.g. Belkacem et al. 2015; Pinçon et al. 2017).
However, none of these cited solutions fully explains the angular momentum transport from the subgiant stage until the end of the red giant branch (RGB). This raises the issue of the incompleteness or inaccuracy of the included set of physical processes driving the internal dynamics in stellar models at each step of the evolution. Magnetic fields are not considered for instance inside the radiative interior of solartype stars along their evolution in models (apart from a few studies that do include the effect of the TaylerSpruit dynamo, e.g. Cantiello et al. 2014), while we do know that at least weak fields must be present, resulting from the relaxation of past dynamo events (e.g. Braithwaite & Spruit 2004; Braithwaite 2008; Duez & Mathis 2010). Such internal magnetism may prevent differential rotation inside the radiative interior, a configuration that has been observed inside solarlike stars and the Sun (García et al. 2007; Benomar et al. 2015; Fossat et al. 2017).
The internal structure of SGs and RGs allows acoustic and gravity modes to couple to form mixed modes. As they probe the deepest layers of the star, they are of great interest for the understanding of physical processes taking place from the deepest layers of the radiative interior towards the surface of the star. For instance, they are known to provide estimates of the mean rotation rate of the core of SGs and RGs (e.g. Deheuvels et al. 2012, 2014, 2016; Mosser et al. 2012, 2017; Gehan et al. 2018). Fuller et al. (2015), Lecoanet et al. (2017), and Loi & Papaloizou (2017) suggested via different magnetic conversion mechanisms that the presence of a strong magnetic field inside the core of a RG can convert magnetogravity waves into Alfvén waves. This results in a loss of energy of the observed modes, which are no longer mixed modes and only present an acoustic signature. Their amplitude is therefore diminished in the power spectrum density (PSD). Even if these theories are controversial (e.g. Mosser et al. 2017), they point out the potentially large effect of magnetic fields on mixedmode amplitudes. The impact of buried pure toroidal or poloidal magnetic fields on acoustic mode frequencies was studied in the context of the SOlar and Heliospheric Observatory (SOHO, Domingo et al. 1995) mission (Solar Oscillation Imager/Michelson Doppler Imager; SOI/MDI, Scherrer et al. 1995) and the Global Oscillations at Low Frequencies (GOLF, Gabriel et al. 1995 instruments) for the observation of the Sun (e.g. Gough 1984; Gough & Thompson 1990; Dziembowski & Goode 1985, 1989; Takata & Shibahashi 1994; Kiefer & Roth 2018) but no trace of such fields has been found in the solar acoustic data. Rashba et al. (2007) developed the equivalent theory for the Sun’s g modes, whose detection is still controversial (García et al. 2007; Appourchaux et al. 2010, 2018; Fossat et al. 2017; Fossat & Schmider 2018; Schunker et al. 2018; Scherrer & Gough 2019). More recent studies (e.g. Hasan et al. 2005; Prat et al. 2019; van Beeck et al. 2020) focused on the impact of stable fields (with a magnetic configuration given by Duez & Mathis 2010 that consists in mixed poloidal and toroidal fields) on slowly pulsating B and γDoradus stars showing pure gravity modes. The effect of moderateamplitude magnetic fields on mixedmode frequencies of evolved stars was theoretically studied by Loi (2020) in the case of a nonrotating star. That study provides analytical (nonperturbative) expressions of the impact of magnetic fields on mixed mode frequencies that are complementary to the perturbative analysis presented in this paper.
Led by all these previous theoretical studies, we focus our efforts on understanding the core dynamics of SGs and RGs. We investigate the impact of a realistic axisymmetric fossil magnetic field buried inside the core of evolved lowmass (M_{⋆} ≲ 1.3 solar masses, M_{⊙}) and intermediatemass (1.3 M_{⊙} ≲ M_{⋆} ≲ 7 M_{⊙}) solarlike pulsators (i.e. SGSs and RGSs) on their observable mixedmode frequencies. In our study, as opposed to the work by Loi (2020), we consider the magnetic field amplitude to be small enough for its effects on the mixedmode frequencies to play as a firstorder perturbation, along with the firstorder effects resulting from the slow differential rotation of the star (e.g. Deheuvels et al. 2014; Gehan et al. 2018). This approximation is motivated by the nondiscovery to date of any magnetic signatures on the mixedmode frequency pattern of observed RGs.
For the magnetic field, we use a realistic mixed poloidal and toroidal configuration (Braithwaite & Spruit 2004; Braithwaite & Nordlund 2006; Duez et al. 2010a,b), and evaluate its impact on typical mixed modes computed using the stellar evolution code MESA (Paxton et al. 2011) and stellar oscillation code GYRE (Townsend & Teitler 2013). After describing the magneticfield configuration and its potential origin and evolution in Sect. 2, we provide in Sect. 3 the firstorder perturbative analysis leading to the magnetic and rotational shift of mixedmode frequencies. Section 4 investigates the shifting of the ℓ = 1 and ℓ = 2 mixed mode frequencies in the case of a M_{⋆} = 1.5 M_{⊙}, Z = 0.02 star along its evolution on the RGB. This star is massive enough to develop a convective core during the MS, which can lead to the intense production of magnetic energy. We provide values of the critical field strength associated to this axisymmetric topology above which the effect on mode frequencies should be visible in real data, and conclude on the validity of the perturbative approach depending on the magneticfield strength and the evolutionary stage of the star. In Sect. 5, we then follow the same approach as we do for the RG branch for the SG stage. During the SG branch (SGB), the nature of most mixed modes is transitioning from acoustic toward gravitydominated modes. Section 6 discusses the consequences of the presence of fossil magnetic fields inside evolved solarlike stars on angular momentum transport. In Sect. 7 we investigate the stellar mass and metallicity dependence of the magnetic splitting of mixed modes. Section 8 focuses on comparing the magnetic field amplitude needed for gmode conversion into Alfvén modes with those needed for magnetic splitting to be detectable in the data. Finally, we conclude on the large potential of this approach for the future detection of magnetic fields from the inversion of magneticfrequency splitting from real data.
2. Internal magnetic fields along the evolution
Magnetism in the depths of stars is very difficult to probe. Spectropolarimetry, which provides most measurements of stellar magnetic fields through the Zeeman effect, only provides information on magnetism emerging from the external layers of the star (e.g. Donati et al. 1997). Probing internal magnetism in evolved solarlike stars could be game changing, because magnetism is known to enhance angular momentum transport (Mestel & Weiss 1987; Charbonneau & MacGregor 1993; Gough & McIntyre 1998; Spruit 1999, 2002; Mathis & Zahn 2005; Fuller 2014; Fuller et al. 2019; Eggenberger et al. 2020). Thus, strong magnetic fields within subgiants and red giants could potentially explain their low coretoenvelope rotationrate ratio.
One powerful mechanism strengthening and sustaining magnetic fields is a convective dynamo, which can create strong magnetism from a weak initial field (e.g. Dikpati & Gilman 1999, for the Sun’s dynamo). Such a mechanism is likely to be the origin of the solar surface magnetic cycle (Brun et al. 2004; Brun & Browning 2017), and is probably active in any convective rotating region inside stars (e.g. Brun et al. 2005; Browning 2008; Brown et al. 2010, 2011; Augustson et al. 2012, 2015, 2016). However, internal layers of lowmass evolved stars are radiative during the SG and RG stages: no convective dynamo can take place in their depths. For a magnetic field to be present inside the radiative interior of evolved solarlike stars, either a dynamooriginated stochastic field has been preserved inside the radiative cavity following the last dynamo episode (fossil field scenario), or the field is transient and currently generated by magnetohydrodynamics instabilities (Spruit 2002; Zahn et al. 2007; Fuller et al. 2019). Spruit (2002) and Fuller et al. (2019) propose a solution to generate magnetic energy from the combination of the Tayler instability (Tayler 1973) and radial differential rotation. This rotational instability may result in largescale internal magnetic fields, but is out of the scope of this paper. We focus here on the socalled fossil field scenario, resulting from the relaxation of magnetic fields originating from a past convective episode, as described in Braithwaite (2008), Duez & Mathis (2010), and Mathis et al. (2010). In the case of giant stars, internal fields may form from two previous convective dynamo episodes, depending on the mass of the star. During the premain sequence (PMS), the star is fully convective: a similar dynamo process to the solar dynamo could take place inside the whole star. When the convective region reduces to a thin surface layer on the MS, the dynamooriginated magnetism may start to relax inside the radiative interior (Arlt et al. 2013; EmeriauViard & Brun 2017; Villebrun et al. 2019). As the star evolves on the MS, its radiative interior thus possesses a relaxed fossil field resulting from the convective dynamo taking place during the PMS. On the MS, either the star is of low mass (M_{⋆} ≤ 1.1 M_{⊙}) and the convective region simply reduces to a surface layer, or for M_{⋆} ≥ 1.1 M_{⊙} the convective region reduces to a surface layer while a convective core is formed due to the change in the hydrogen fusion mechanisms. Such a core likely also convectively generates a magnetic field via a dynamo (see Fig. 1). This field could therefore enhance the potentially already present relaxed magnetic field inside the radiative interior of the star. During the SG phase, the interior is completely radiative, and a fossil field is formed, which is also present during the RGB. The field is located inside the radiative region of the star, the extent of which varies with evolution.
Fig. 1. Magnetic configuration schemes (not to scale) following the evolution of low and intermediatemass stars from the main sequence. The represented configuration corresponds to the axisymmetric field topology considered in our study. It presents both a poloidal (blue lines) and a toroidal (red lines) component, with the axis of symmetry aligned with the rotation axis of the star. 
2.1. The fossil field scenario
In the fossil field scenario, the relaxing magnetic field may be preserved inside the radiative interior when the convection ends. If no disruptive processes occur inside the radiative interior from the MS towards the RGB, such as strong differential rotation (Aurière et al. 2015; Gaurat et al. 2015), the stochastic field will eventually stabilise into a fossil equilibrium configuration, reached when the Lorentz force balances hydrodynamics forces (Chandrasekhar 1958). In that case, and by assuming a constant magnetic diffusivity, the induction equation is written as:
where η is the magnetic diffusivity. If there is no turbulence in the region, the associated Ohmic relaxation time due to atomic processes is τ = R^{2}/η. For a ‘fossil field’, this is comparable to the lifetime of the star itself, being 10^{10} years in the case of a field inside the radiative interior of the Sun (e.g. Cowling 1945). Thus, if ever generated and never disrupted by hydrodynamical or other MHD processes, the dynamooriginated field should remain trapped for the remaining evolution of the star as a fossil field, and should still be present in its radiative interior during the RGB (Stello et al. 2016b). The intense magnetism discovered in some white dwarfs (Angel et al. 1981; Putney 1999; Schmidt et al. 2001) could then result from such a PMSMS dynamo field, surviving as a fossil field during the succeeding evolutionary stages. The primary challenges regarding these equilibrium magnetic fields are the determination of their 3D configuration and the estimation of their amplitude.
2.2. Fossil field topology
Purely toroidal and purely poloidal magnetic configurations are known to be unstable (e.g. Tayler 1973; Markey & Tayler 1973; Braithwaite 2006, 2007). The stability of mixed configurations with both poloidal and toroidal components is demonstrated in Tayler (1980), and the relative energy contained in the toroidal and poloidal components for the field to be stable has been evaluated numerically (Braithwaite 2008) and semianalytically (Akgün et al. 2013). In a similar vein, Braithwaite & Spruit (2004) simulated the relaxation of stochastic fields in a nonrotating radiative medium. Those simulations show that a stochastic field representing a dynamogenerated field inside a stably stratified region relaxes into a stable, largerscale, mixed poloidal and toroidal magnetic field. Duez & Mathis (2010) give the closest semianalytic description of such stable mixed toroidal and poloidal fossil fields. As opposed to previous studies (e.g. Gough & Thompson 1990) that considered purely toroidal or purely poloidal field topologies, we use this stable mixed formalism to investigate the effect on mixedmode frequencies of a hypothetical fossil field trapped inside the radiative interior of size R_{rad} of evolved stars (see Fig. 1). We do not consider in our study nonperturbative boundary conditions associated to more realistic magnetic field configurations (e.g. Roberts & Soward 1983; Campbell & Papaloizou 1986; Dziembowski & Goode 1996; Bigot et al. 2000). The mixed toroidal and poloidal expression from Duez & Mathis (2010) that minimises the energy of the system is dipolar. Such magnetism confined inside the radiative interior of the star is written as:
where ψ is the stream function satisfying
with μ_{0} the vacuum magnetic permeability expressed as 4π in cgs units, α a normalisation constant fixed by the chosen magneticfield amplitude, λ the eigenvalue of the problem to be determined, R_{rad} the radius of the radiative interior, and
with j_{1} (y_{1}) the firstorder spherical Bessel function of the first (second) kind (Abramowitz & Stegun 1972) and ρ = ρ(r) the density of the star.
In our study, we consider the field to be aligned with the rotation axis of the star; hence, nonaligned and nonfossil field configurations are out of the scope of this paper. Finally, the axisymmetric poloidal and toroidal magnetic field used in our analytical calculations can be expressed as
with the following correspondences to match the Duez & Mathis (2010) formalism:
with the prime symbol denoting the radial derivative.
In order to confine the field inside the radiative area we set λ as the smallest positive constant for B to vanish at the radiative–convective boundary (see Fig. 2). As opposed to massive stars, red giants present a strong internal density gradient. This property prevents from going to zero for any value of λ. As a result, we could not simultaneously cancel the b_{r} and b_{θ} components of the field at the radiative–convective boundary. To ensure that the field is confined to the radiative interior, we choose to eliminate only b_{r}. This leads to an azimuthal current sheet that might create instabilities (see Duez & Mathis 2010, for a detailed description). This is accomplished by searching for zeros of y_{1}(λ) instead (a more detailed description is given in Appendix B). We find that the corresponding smallest value for λ is about 2.80. This eigenvalue corresponds to the first zero of the function y_{1}, providing the most stable parametrisation that eliminates b_{r} at R_{rad}. This rather small value of λ leads to poloidally dominated fields as represented in Fig. 2.
Fig. 2. Mixed poloidal (black lines) and toroidal (colour scale) magnetic field modelled using the formalism of Duez & Mathis (2010); see Eq. (2). Values are normalised by the maximum radial field amplitude. The field is confined inside the radiative region of the star, the extent of which is given by the corresponding evolution model MESA and indicated by the black circle. Left to right: cases of a 1.5 M_{⊙}, Z = 0.02 star from the SG phase to RGB, with ages of 2.65, 2.75, and 2.80 Gyr, and effective temperatures of 5967, 5075, 4900 K, respectively. 
2.3. Estimating the evolving magnetic field amplitude
By considering an ideal fossilfield scenario, one can estimate the amplitude of the field along the evolution of the star by considering magneticflux conservation from the end of the more recent dynamo episode without considering any Ohmic loss of energy due to the reconnection of field lines during successive relaxations. Therefore, here we adopt the approach of giving an upperlimit of the field amplitude by first computing the rate of magnetic energy produced by stellar dynamos and transmitting its amplitude to the fossil field, a simple approach that is currently used in studies in asteroseismology (e.g. Fuller et al. 2015; Cantiello et al. 2016). However, we know that during the relaxation of fossil fields, loss of magnetic energy occurs at small scales. This has been observed in numerical simulations (e.g. Braithwaite 2008; Duez 2011; EmeriauViard & Brun 2017). How to properly quantify this Ohmic loss of energy and the corresponding turbulent scales all along the evolution of stars is an open and difficult question because a coherent theoretical modeling should be provided (e.g. Moffatt 2015; Hotta 2017) while our current MHD simulations assume magnetic diffusivities (and diffusivities in general) that are in general higher than in stellar interiors. This question will be examined in detail in a forthcoming article.
We consider the extent of the fossil field configuration to match the size of the radiative interior of the star during the SG and RG phases (see Fig. 1). This implies that no reconnections with the envelope field are considered in our simplified configuration. In our fossil field scenario, the last recorded internal dynamo is either the convective core dynamo for stars more massive than 1.1 M_{⊙}, as presented in Fig. 1, or the global convective dynamo during the PMS for lowmass stars.
As in Fuller et al. (2015) and Cantiello et al. (2016), we estimate the magnetic amplitude resulting from magnetic flux conservation from the end of the last convective dynamo episode towards the current giant stage, through:
with R_{init} the radius of the convective cavity of the last internal dynamo phase, and R_{rad, ⋆} the size of the current radiative sphere inside which the field is supposed to be confined. This flux conservation equation ensures that the entirety of the magnetic energy is conserved inside the radiative interior, with no exchanges with the surrounding convective zone, and with no disruption of the magnetic field by an external mechanism such as differential rotation. R_{init} and R_{rad, ⋆} are evaluated along the evolution using the MESA stellar evolution code (Paxton et al. 2011). The B_{0, init} field amplitude is estimated during the recession phase of the last internal dynamo by the use of rotating stellar models computed with STAREVOL (STAREVOL is used instead of MESA to take into account rotation during the early stages of the evolution of the star needed in order to estimate the magnetic field amplitude resulting from magnetostrophy and buoyancy regimes; see Amard et al. 2019; Astoul et al. 2019, for details). The field is evaluated following Augustson et al. (2019) and Astoul et al. (2019) by considering the magnetostrophic, equipartition, and buoyancydriven regimes (see Appendix A for details about how the different regimes are constructed). The turbulent equipartition regime corresponds to the state where the convective kinetic energy density of the fluid is fully converted into a magnetic energy density (superequipartition regime therefore refers to the fluid’s magnetic energy being greater than its kinetic energy). An equipartition or subequipartition regime is invoked for the Sun’s dynamo inside the convective envelope (Brun et al. 2017). The magnetostrophic regime is reached when Coriolis acceleration balances Lorentz force in the momentum equation. This corresponds to a maximumamplitude estimate of the magnetic field, and can be a superequipartition regime. The magnetostrophic regime is suspected by Augustson et al. (2019) to be applicable for stellar dynamos. The buoyancy dynamo regime corresponds to the case where Coriolis acceleration, buoyancy, and Lorentz forces all have the same order of magnitude. This regime is usually considered for rapidly rotating lowmass stars and planets with strong density gradient (Christensen et al. 2009).
2.3.1. Case of M_{⋆} ≲ 1.1 M_{⊙}
Lowmass stars do not develop a convective core during the MS, apart from a small convective interior that arises from the energy released during the CNO cycle before the temperature drops (down to ∼5600 K) leading to the pp fusion mechanism. However, such stars possess a convective envelope, schematically represented in Fig. 1. The top panel of Fig. 3 represents the amplitude of the magnetic field generated at the bottom of the convective envelope while the star evolves. When the radiative interior replaces the convection at the beginning of the MS (∼2 × 10^{7} years), the amplitude of the expected magnetic field that relaxes into the radiative interior at the boundary with the convective envelope is of the order of magnitude of respectively 10^{6}, 10^{5}, and 10^{4} G for the magnetostrophic, buoyancy, and equipartition regimes. Astoul et al. (2019) concluded that in order to reproduce the magnetic amplitude observed at the surface of rotating PMS and MS stars, the magnetostrophic regime should be considered. For the remainder of this paper we therefore use the amplitude of the dynamo field estimated by considering the magnetostrophic regime as an upper estimate of the magnetic amplitude inside the radiative interior. The equipartition approximation that is used by Cantiello et al. (2016) and Fuller et al. (2015) provides much smaller estimates of the field amplitude, as seen in Fig. 3: we consider these estimates as the lowest possible value of the fossil field in the remainder of the article.
Fig. 3. Top: magnetic field amplitude at the base of the convective envelope along the evolution of a M_{⋆} = 1 M_{⊙} star, with solar metallicity (Z = 0.0134), uniformly rotating with an initial rotation rate during the disklocking of the PMS of Ω = 2.5 μHz. The blue, red, and orange lines respectively represent the field estimated by considering the magnetostrophic, buoyancy, and equipartition regimes. The black dashed line represents the location of the limit radius between the internal radiative layer and the convective envelope. Bottom: same as top panel at the top of the convective interior along the evolution of the star. The black dashed line represents the location of the limit radius between the convective core (when existing) and the radiative layer. 
We shall also discuss the quick appearance of a convective core at the beginning of the MS. The characteristic timescale for a dynamo to emerge from convection and differential rotation is of the order of a year (by considering a typical solar internal rotation rate). Therefore, despite its comparatively short lifespan, this early convective core should still produce a dynamo field. Its amplitude would be of the order of magnitude represented in the bottom panel of Fig. 3. This newly generated field has an amplitude equivalent to that of the original field resulting from the dynamo in the convective envelope during the PMS: the presence of this convective core right after the end of the PMS should not modify the amplitude of the already relaxing field. In this way we ignore the interaction between the dynamo field and the surrounding fossil field as shown by Featherstone et al. (2009) in the case of Atype stars. Such interaction may enhance the dynamo process, and may lead to stronger core field amplitudes. For lowmass stars with no persistent convective core during the MS, we thus consider the relaxation and stabilisation of the field to occur at the end of the PMS inside the newly formed radiative interior, ignoring the early convective core, where B_{init} = B_{PMS} with values ranging from 10^{4} to 10^{6} G, depending on the chosen dynamo force balance regime.
2.3.2. Case of M_{⋆} ≳ 1.1 M_{⊙}
If the star is massive enough for a sustainable convective core to form during the MS, we then consider this convective core to be the location of the last internal dynamo episode from which the fossil field relaxed. In that case, R_{init} is equal to the radius corresponding to the maximum extent of the convective core during the MS (R_{MS}), which is estimated with evolution models from MESA (see Appendix E). Figure 4 represents the same diagrams as Fig. 3 for a M_{⋆} = 1.4 M_{⊙} star. The top panel indicates the relaxation of the PMS magnetic field with the same order of magnitude as the B_{PMS} field estimated for lowmass stars represented in Fig. 3. In the bottom panel of Fig. 4 we observe the formation of the convective core around 10 Myr, which lasts until the end of the MS. The corresponding field amplitudes are of the same order of magnitude as the different fields estimated from the relaxation of the magnetic field resulting from the PMS dynamo (see top panel of Fig. 4). In that case, we consider the relaxation to a fossil configuration of the field to occur inside the radiative interior at the end of the MS: B_{init} = B_{MS}, with values also ranging from 10^{4} to 10^{6} G. This new relaxing field eventually couples with the surrounding relaxed magnetic field from the PMS, which has an amplitude of the same order of magnitude.
As shown in Figs. 3 and 4, the expected magneticfield amplitude in the convective star during the PMS is very similar to that inside the convective core during the MS. Indeed, B_{init} is in both cases evaluated to be of the same order of magnitude (between 10^{4} and 10^{6} G). As a result, it appears from our study that there is no correlation between the development of a convective core during the MS and the fossil magnetic field amplitude on the RGB. This would imply that the discrepancy between low and intermediatemass star dipolar mode amplitudes on the RGB (Stello et al. 2016a) might not be provoked by the fossil magnetic field resulting from the core dynamo action on the MS. The main difference in field amplitude during the RGB in our scenario is therefore due to the contraction of the radiative interior, as the interiors of massive stars contract more efficiently than those of lowmass stars after the MS, resulting in larger magnetic field amplitudes during the RGB from Eq. (10).
We consider the dissipation of the fossil field to take place during the Ohmic timescale computed with its relaxed equilibrium global scale, and therefore we neglect the losses of energy induced by its relaxation from the previous dynamo field during the conversion of convective layers into stably stratified radiative ones throughout the evolution of stars. This is a significant approximation, as already pointed out before, as the reconnection of smallscale magnetic structures into a larger global structure would decrease the fossilfield initial amplitude (e.g. Braithwaite 2008; Duez 2011; Moffatt 2015; EmeriauViard & Brun 2017). This reconnection process could result into amplitude discrepancy between the fossil fields of low and intermediatemass stars. In addition, our estimation of magnetic field amplitudes throughout evolution is based on scaling laws and flux conservation only; we do not take into account any interaction between the fossil field in the radiative shell and the dynamo action inside the convective core of intermediatemass stars on the MS. Such interaction could enhance the dynamo action inside the convective core of intermediatemass stars (e.g. Featherstone et al. 2009) and thus result in stronger fossilfield amplitudes inside intermediatemass stars than inside lowmass stars on the RGB. Such interactions between stable and dynamo fields during the MS could therefore support the magnetic scenario as an explanation for the amplitude of dipolar modes discrepancy as presented by Stello et al. (2016a).
3. Firstorder frequency perturbation of mixed modes
In this section we develop the firstorder perturbation of the eigenfrequency of the oscillation modes due to the stellar rotation and internal magnetism.
3.1. Oscillations in nonrotating, nonmagnetised stars
The general linearised equation of motion of the fluid inside a nonrotating, nonmagnetised star can be written in the inertial frame as:
where ξ denotes eigenfunctions of the modes propagating inside the star. In the case of evolved solarlike oscillators, oscillation modes correspond to the superposition of mixed acoustic and gravity waves (Beck et al. 2011; Bedding et al. 2011). F represents all the applied forces in the inertial frame and is composed of the effect of pressure and density gradient for a nonrotating, nonmagnetised star (where the effects of the centrifugal acceleration and of the Lorentz force are neglected when computing the structure of the star).
We suppose a periodic Lagrangian displacement of the fluid inside the star due to oscillation modes , with ω the eigenfrequencies of the oscillations. The nonperturbative equilibrium state refers to the nonrotating, nonmagnetic oscillating star, for which the equation of motion can be written:
with ω_{0} the unperturbed eigenfrequencies, ξ_{0} the unperturbed eigenmodes, and F_{0} reduced to pressure and buoyancy forces. Unperturbed eigenmodes for slowly or nonrotating stars can be written as a function of the azimuthal order m and the degree ℓ of the mode (Unno et al. 1989):
with (r, θ, φ) the usual spherical coordinates, and the spherical harmonics of degree ℓ and azimuthal order m.
3.2. Firstorder perturbative forces
3.2.1. Rotational inertial forces
In order to connect our theory to observations, we approximate the star as a slowly rotating object with two layers, where each layer is rotating as a solid body. Let Ω_{core} be the angular rotation rate in the radiative interior and Ω_{env} be the angular rotation rate in the convective envelope. This approximation is justified because asteroseismology generally provides only one measurement inside the deep radiative interior of the star and one at the surface (because of the limitation of the number of measured modes and the spatial structure of the kernels or eigenfunctions). These rotation rates are considered small enough so that the resulting terms in Ω(r) only come into play as a firstorder perturbation (e.g. Gough & Thompson 1990). Observational studies of rotational splitting of mixed modes provide average values of surface and internal rotation for SGs and RGs. Specifically, we fix Ω_{core} ≃ 0.5 μHz and Ω_{env} ≃ Ω_{core}/10 (Deheuvels et al. 2014; Gehan et al. 2018). Finally, with R_{rad} being the size of the radiative interior,
The system is thus affected in the inertial frame at first order in Ω by the Coriolis acceleration (F_{c}, Eq. (16)) and by the advection linked to the rotating motion of the star relative to an observer in an inertial frame (F_{f}, Eq. (17)). In the following, for rotational effects, we make use of the decompositions ω = ω_{0} + ϵ_{Ω}ω_{1}, ξ = ξ_{0} + ϵ_{Ω}ξ_{1}, and F = F_{0} + ϵ_{Ω}F_{1}, with ϵ_{Ω} ≪ 1, which is the dimensionless parameter ϵ_{Ω} = Ω/Ω_{C}, with Ω_{C} = (GM/R^{3})^{−1/2}, which is the Keplerian critical angular velocity (see Gough & Thompson 1990). The Coriolis operator stands for:
where Ω = Ω(r)(cos θe_{r} − sin θe_{θ}) is the rotation vector (whose amplitude is small enough in the core and in the envelope for the effect of the Coriolis acceleration to be considered as a firstorder perturbation). Under the same approximation, the advection operator associated with the change of frame is written at first order as
3.2.2. The Lorentz force associated to magnetism
For RG with a moderate internal magnetic field, we consider the magnetic field energy to be weak enough for the effects of the unperturbed Lorentz force on the system to be negligible compared to the gravitational force (e.g. Augustson & Mathis 2018, and in prep.).
The system is thus also affected in the inertial frame at first order by the perturbed Lorentz force (δF_{L}, Eq. (18)). We again make use of decompositions associated with the perturbed Lorentz force: ω = ω_{0} + ϵ_{B}ω_{1}, ξ = ξ_{0} + ϵ_{B}ξ_{1} and F = F_{0} + ϵ_{B}F_{1}, with ϵ_{B} ≪ 1, the dimensionless parameter ϵ_{B} = B_{0}(μ_{0}GM^{2}/R^{4})^{−1/2} which compares the magnetic and gravitational forces (Gough & Thompson 1990). Finally, the perturbed linearised magnetic operator is given by
where δF_{L, j + t} represents the sum of the current and tension terms:
where the fluctuation of the magnetic field that comes from the linearised induction equation^{1} is
and δF_{L, c} is the compression term associated with the compressibility of the mode (Gough & Thompson 1990):
3.3. Firstorder perturbation equations
The firstorder momentum equation is written in the inertial frame as
with ξ_{1} being the perturbed eigenfunctions of the modes. We then apply the scalar product to Eq. (22). As F_{0} is Hermitian, Eq. (22) simplifies and leads to the expression of the frequency shift at first order in the inertial frame:
3.4. Analytical development of frequency shifts
We directly estimate the inertia of the modes, and the contribution of the Coriolis acceleration and of the advection term associated to the change of frame to the splitting in the inertial frame. Equation (24) is the inertia of the modes, independent of the azimuthal order m:
The effect of the Coriolis acceleration and of the change of frame are usually written together as the global rotational perturbation (Aerts et al. 2010). For clarity we chose to detail the two terms in Eqs. (25) and (26):
We verify that in the case of solid rotation Ω_{s}, the change of frame can as usually be written as
The calculation of ⟨ξ_{0}, δF_{L}(ξ_{0})/ρ⟩ involves 50 terms for which the angular integral does not cancel by geometrical combination of the spherical harmonics. These are reported in Appendix C. Depending on the nature of the mixed modes, dominant terms can be extracted from this massive expression (see e.g. works by Hasan et al. 2005; Rashba et al. 2007, for highradialorder pure gravity modes). In the case of subgiant stars, mixedmode patterns are strongly dominated by peaks located at the acoustic mode eigenfrequencies, and are denoted as p − m modes in the following study. On the contrary, modes located at internal gravitymode frequencies have very small acoustic signatures. They dominate the spectrum of red giants, and are denoted g − m modes. The dominant nature of the mode strongly influences the value of each of the 50 terms; we refer to Sects. 4.3 and 5.2 for the extraction of dominant terms and to Mathis et al. (2021) for asymptotic formulations in the case of gdominated modes and pdominated modes.
3.5. Validity of the perturbative analysis
In our scenario the magnetic field perturbation on mixedmode frequencies plays a role at first order. For this approximation to be valid, the magnetic field amplitude should be small when compared to the zerothorder processes. Thus, the frequency of the g − m modes should be larger than the characteristic Alfvén frequency ω_{A}, usually written as
with k the wave vector. For highorder radial modes, B.k ≃ B_{0}b_{r}k_{r}, with b_{r} close to unity, k_{r} the vertical wave number scaling for g − m modes as , the critical field amplitude value for the perturbative study to be valid can be expressed as:
with ρ_{c} the central density, and N_{max} the maximum value of the BruntVäisälä angular frequency profile defined as:
where g is the gravity of the hydrostatic background. In order to ensure that the frequency of the mode is much larger than the Alfvén frequency, we set the limit of validity of the perturbative analysis for g − m modes to B_{c, g}/1000.
In the case of highradialorder p − m modes, the Alfvén speed
should be compared to the sound speed of the gas (c_{s}) through S_{l} the Lamb frequency, defined as:
For highradialorder p − m modes, the minimum field above which the perturbative analysis may not be valid, considering the minimum value of the Lamb angular frequency S_{l, min} inside the radiative core, is written as
In the remainder of the present study, we set the limit of validity for the perturbative analysis for p − m modes to B_{c, p}/1000.
4. Investigation of the magnetic impact on mixedmode frequencies of red giants
To evaluate the different contributions from Eq. (23) to the frequency perturbation of mixed modes, we model stellar structures with the MESA evolutionary code (Paxton et al. 2011), and we compute the frequency of the associated oscillation modes using the GYRE pulsation code (Townsend & Teitler 2013). In this section, we estimate frequency perturbations from rotation and magnetism on mixed dipolar (ℓ = 1) and quadrupolar (ℓ = 2) modes for the reference star with M_{⋆} = 1.5 M_{⊙} and Z = 0.02 along the RGB. This star is characteristic of intermediatemass stars observed on the RGB (e.g. Yu et al. 2016). In Sect. 5 the effect of internal magnetism on mixed modes during the SG stage is also investigated.
4.1. Expected fossil magnetic field signature on mixedmode frequencies during the RGB
From the conservation of the magnetic flux from the end of the MS, we evaluate the magnetic amplitude range expected for a typical red giant with M_{⋆} = 1.5 M_{⊙} and Z = 0.02. The star has oscillation modes centred around the frequency of maximum mode power ν_{max} of about 172 μHz. For such a star, the initial radiative radius R_{init} is about 92% of the star’s total radius at the beginning of the subgiant phase (R_{⋆} = 1.8 × 10^{11} cm, R_{init} = 1.6 × 10^{11} cm). This reduces to about 13% of the total radius for the considered evolutionary stage star at the middle of the RGB (R_{⋆} = 3.7 × 10^{11} cm, R_{rad} = 5.0 × 10^{10} cm). The equipartition regime, which is used to estimate field amplitudes in the study of Cantiello et al. (2016) leads to B_{0} ∼ 0.7 MG, while the magnetostrophic regime sets the upper boundary of the field amplitude to B_{0} ∼ 7 MG (these boundaries will be extended along the evolutionary path of the reference star for the analysis of the impact of magnetism along the evolution of the star in Sect. 4.5). In the following paragraph, we investigate the effect of the change of internal structure of the star along its evolutionary path on the RGB on the effect of magnetism at mixedmode frequencies. We choose to use a magneticfield amplitude value of 1 MG, consistent with the presence of a convective core dynamo during the MS.
As a first step, we verify that our perturbative approach can be used for this typical M_{⋆} = 1.5 M_{⊙}, Z = 0.02, ν_{max} ≃ 172 μHz red giant. By fixing the validity limit of the perturbative analysis to a reasonable value of B_{0} = B_{c, g}/1000, the critical value for the magneticfield amplitude at ν_{max} is B_{0} ≃ 3 MG from Eq. (29). Regarding the relatively large value of this limit, we conclude that the typical magnetic frequencyshift amplitudes of g − m modes around ν_{max} represented in Fig. 5 belong to the valid frequency range of the perturbative analysis. However, a similar calculation for a frequency of 100 μHz leads to a critical field value (divided by 1000) of B_{0} ∼ 1 MG, equal to the value used to build Fig. 5. The MG order of magnitude is too high to ensure that the perturbative analysis is valid for evolved red giants with ν_{max} lower than ∼100 μHz. This upper boundary will further be discussed along the evolution of the reference star in Sect. 4.5. The reasonable limiting value of B_{0} = B_{c, p}/1000 for the pertubative analysis to be valid for p − m modes leads to field amplitudes of about B_{0} ≃ 10 MG for typical S_{l, min} ≃ 200 μHz for the considered red giant by using Eq. (33). The following study with a field amplitude of 1 MG is therefore consistent with a perturbative regime for p − m modes.
Fig. 5. Frequency splittings (δν) calculated for mixed modes frequencies of a M_{⋆} = 1.5 M_{⊙}, Z = 0.02, ν_{max} ≃ 172 μHz red giant, corresponding to (a) magnetic splitting, (b) the Coriolis effect, (c) the change from corotating to inertial frame, and (d) total frequency shifts. Red lines and points represent splittings for the m = 0 component. The blue and green lines represent the magnetic splitting for the m = 1 and m = −1 components respectively, that overlap on panel a. The grey dashed line represents the ζ function as defined by Goupil et al. (2013), the beige vertical line indicates the position of ν_{max}, and the grey vertical line the frequency of the ℓ = 1 gdominated mode closest to ν_{max}, ν_{1, g, max}, with squares indicating the corresponding splitting values associated with the different orders m. 
4.2. Results for a typical M_{⋆} = 1.5 M_{⊙}, Z = 0.02, ν_{max} ∼ 172 μHz star on the RGB
The different terms constituting Eq. (23) are evaluated for mixed modes ℓ = 1 and ℓ = 2 frequencies within the [ν_{max} − 7Δν : ν_{max} + 7Δν] interval, which ensures that detectable solarlikeoscillation modes are inside this region for most stars. The unperturbed frequencies of mixed modes computed with GYRE (Townsend & Teitler 2013) are degenerated because of the choice of a nonrotating and nonmagnetic equilibrium. Each unperturbed mode is thus expected to be split into 2ℓ+1 components because of both rotational and magnetic effects when applying Eq. (23). Both magnetic firstorder terms (reported in Appendix C) and the mode inertia (see Eq. (23)) depend on m^{2} or are independent of m; thus, we expect the magneticfrequency splittings of the ±m modes to be equal.
Panel a of Fig. 5 represents magnetic shifts on ℓ = 1 mixedmode frequencies (δν_{mag, core, m}=−⟨ξ_{0}, F_{L}(ξ_{0})⟩/(2ω_{0}⟨ξ_{0},ξ_{0}⟩)), resulting from a field of amplitude B_{0} = 1 MG affecting the reference star on the RGB. The magnetic field components b_{r}, b_{θ}, and b_{φ} expressed by Eqs. (7)–(9) describe the stable axisymmetric fossil field from the configuration given by Eq. (2). The magnetic splittings of the m = −1 and m = 1 modes are equal as expected (C_{1, [1, −1]}), and are larger than the splitting for m = 0 modes (C_{1, 0}) by a factor C_{1, [1, −1]}/C_{1, 0} = 2, as calculated by Hasan et al. (2005). A characteristic pattern in the magnetic splittings can be seen: the global trend follows a power law in ν, with deviations toward lower splitting values equivalently spaced in frequency. Such deviations are associated with the ζ function, a measure of the g nature of the mode through the ratio of the mode inertia in the g cavity over the total mode inertia (Deheuvels et al. 2012; Goupil et al. 2013):
This characteristic pattern is very similar to the pattern of mixed modes split by rotation as studied by Goupil et al. (2013) and is represented by the grey dashed line in Fig. 5. Panel b of Fig. 5 represents the frequency splittings due to the Coriolis acceleration only on ℓ = 1 mixedmode frequencies, whereas panel c shows the effect of the change of frame on mixed mode frequencies due to the rotation of the star. These two rotational components are evaluated by considering typical redgiant rotation rates from Gehan et al. (2018) of Ω_{core}/(2π) = 0.5 μHz, and Ω_{env}/(2π) = 0.05 μHz. Combined together, these terms are well known from previous studies (e.g. Aerts et al. 2010; Mosser et al. 2015; Vrard et al. 2015), and can also be expressed as a function of the ζ function.
Dips in the ζ function (grey dashed line) indicate the location of p − m modes whereas the value of ζ tends towards 1 for g − m modes. As the rotation rate of evolved solarlike stars is much higher deeper below the surface than at the surface, the gmode cavity is more affected by rotation. This leads to a global rotation splitting (including both b and c components) that is larger for g − m than for p − m modes, which also probe the slowly rotating envelope. As a consequence of the confinement of the magnetic field inside the gmode cavity, gdominated mixed modes are also more affected by the magnetic field than pdominated mixed modes. Indeed, the minimum values of ζ correlate with the minimum in the magnetic splitting in panel a.
Finally, the bottom right panel provides the global splitting due to rotation and magnetism under the previously detailed rotational and magnetic configurations. Unperturbed frequencies are generally shifted towards higher frequencies for modes around ν_{max}, with a greater shift for prograde (m = 1) and retrograde (m = −1) components than for the zonal modes (m = 0), which are only affected by magnetism. As a consequence of this global frequency change, magnetism leads to a break of symmetry of the usually studied rotational gdominated triplet, as seen by the nonregular vertical spacing between two consecutive grey squares. Rotation plays against magnetism for the shifting of m = −1, p − m modes (at frequencies corresponding to dips of the ζ function), resulting in upward peaks in the bottom right diagram at high frequencies.
In Fig. 6 the same approach is taken for quadrupolar ℓ = 2 mixed modes. In panel a, magnetic splits are represented as a function of the frequency of the modes, depending on the value of m. As opposed to the ℓ = 1 case, zonal modes are more affected by magnetism than the m = 1, −1 and even more so than the m = 2, −2 mixed modes, by factors of C_{2, 0}/C_{2, [1, −1]} = 1.125 and C_{2, 0}/C_{2, [2, −2]} = 1.8 respectively, as in Hasan et al. (2005). Panels c and d represent the contribution from the Coriolis acceleration and from the change of frame. Panel d shows the global frequency shifts for m ∈ { − 2, −1, 0, 1, 2} in the inertial frame due to magnetism and rotation. The ζ function and the location of ν_{1, g, max} are reported in each panel for reference. As for ℓ = 1 mixed modes, p − m modes are less affected by both rotation and magnetism than g − m modes. Each perturbed quintuplet corresponding to one given unperturbed mixed mode is expected to have nonconstant spacings between its components, as shown by the nonconstant vertical spacing between the quintuplet components in panel d. The asymmetry of ℓ = 1 and ℓ = 2 mixedmode multiplets are schematically represented in Fig. 7, in which we can observe the simplified behaviour of dipolar and quadrupolarmode frequency patterns due to internal magnetism.
Fig. 6. As in Fig. 5 but for ℓ = 2 mixed modes. The red line represents the m = 0 mode, the green line with blue points the m = −1, 1 components, and the green line with purple points the m = −2, 2 components. 
Fig. 7. Sketch representing the effect of rotation and magnetism on a single unperturbed mixed mode frequency of degree ℓ = 1 and ℓ = 2. 
4.3. Asymptotic expressions of the magnetic perturbations on g–m mode frequencies
In order to interpret the previously presented asymptotic patterns of magnetic splittings, we investigate the frequency dependency of the magnetic splitting at low (ω_{0} ≪ N) and high (ω_{0} ≫ S_{l}) frequency, with N the angular BruntVäisälä frequency (Eq. (30)) and S_{l} the Lamb angular frequency (Eq. (32)) in rad s^{−1}.
In both cases, the oscillation vertical wavelength is much smaller than the characteristic distance of change of the equilibrium state describing the surrounding fluid. In the Cowling approximation, and for highorder modes, the equation of nonradial oscillations can be written (e.g. Aerts et al. 2010):
In the approximation of highorder mixed modes, the local radial displacement is dominated by the horizontal one (ξ_{r} ≪ ξ_{h}) in the radiative region. The highradialorder modes have a small vertical wavelength, allowing us to perform an asymptotic JeffreysWentzelKramersBrillouin (JWKB) analysis (Hasan et al. 2005; Prat et al. 2019). Terms with highorder ξ_{h} derivatives dominate because , and k_{r} ≫ 1.
These approximations allow us to estimate the dominant terms composing Eq. (23) in the case of lowfrequency g − m modes:
From Eq. (23), the magnetic contribution to the frequency perturbation is
with, in the case of highorder gravity modes,
where we used Eqs. (36) and (39). In the lowfrequency regime, where ω_{0} ≪ N ≪ S_{l}, Eq. (35) becomes
In the Cowling approximation, we obtain a set of radial equations of momentum (see for instance Alvan et al. 2013), with p_{1} the perturbation associated with the gas pressure:
and
Finally, we have that
We write N = N_{max} × f(r), with N_{max} the maximum of the BruntVäisälä frequency inside the radiative interior and f(r) containing all the radial dependence of the N profile. As we look for the frequency dependence of ℐ_{g}, and by the use of the JWKB solution (Froman & Froman 2005), we obtain
with F(r) containing all factors that are independent of ω. Considering highorder g − m modes, we evaluate
and
This leads to
The magnetic splitting of g − m modes is thus proportional to to leading order. In Fig. 8, we represent the magnetic splittings δν_{mag} = δν_{mag, core, m} normalised to one of m = −1, 0, 1 mixed modes for the reference red giant from panel a of Fig. 5. The black line indicates the normalised frequency power law from Eq. (49). As a result, g − m modes perfectly follow the frequency power law, as inferred from the JWKB analysis. However, p − m modes that are less affected by magnetism are not reproduced by this analysis. For a complete theoretical description of asymptotic behaviours, including the origin of the p − m dips through the ζ function and the theoretical estimation of angular integrals setting the amplitude of δν_{mag}, we refer to our paper (Mathis et al. 2021).
Fig. 8. Normalised magnetic splittings versus mixedmode frequencies computed for a simulated M_{⋆} = 1.5 M_{⊙}, Z = 0.02, ν_{max} ≃ 172 μHz red giant. The normalised power law following 1/ν^{3} is superimposed in black. 
4.4. Characterisation of the multiplet asymmetry
Regardless of the degree ℓ of the mixed modes, the g − m multiplet formed by the lifting of degeneracy by firstorder perturbations linked to rotation and magnetism is asymmetric. This asymmetry can be quantified for each m through the calculation of an asymmetry measure δ_{ℓ, m} (similar formalism to that in Deheuvels et al. 2017, adapted for the consideration of magnetism in addition to rotation):
This measure is equal to zero when the m and −m components of the multiplet are perfectly symmetric around the m = 0 component, and reaches 1 or −1 when one of the m components overlaps with the m = 0 mode. For the reference red giant studied in this section, with the considered Ω_{core} = 10 × Ω_{env} = 0.5 μHz, B_{0} = 1 MG, total frequency perturbations are estimated at ν_{1, g, max}: δν_{1, 1} ≃ 0.41 μHz, δν_{1, −1} ≃ −0.076 μHz, δν_{1, 0} ≃ −0.11 μHz, δν_{2, 2} ≃ 0.93 μHz, δν_{2, −2} ≃ −0.74 μHz, δν_{2, 1} ≃ 0.56 μHz, δν_{2, −1} ≃ −0.27 μHz, and δν_{2, 0} ≃ 0.17 μHz from d panels of Figs. 5 and 6. These values lead to asymmetry measures of δ_{1, 1} ≃ 0.28, δ_{2, 2} ≃ −0.09, and δ_{2, 1} ≃ −0.06. We notice the change of sign of the measure from ℓ = 1 to ℓ = 2 modes, where the m = 0 mode is closer to the m = −1 mode for ℓ = 1, and closer to the m = 1 mode for ℓ = 2. The absolute value of the asymmetry measure is higher for ℓ = 1 modes than for ℓ = 2, meaning that the impact of magnetism compared to rotation is stronger for ℓ = 1 modes. This is not surprising considering that the acoustic dipolar modes couple with internal gravity modes, which are sensitive to the magnetic field, deeper inside the radiative interior than ℓ = 2 modes for a given frequency.
Figure 9 depicts the changes in the magnetic splittings along the evolutionary track of the reference star from the base of the RGB, during which mixed modes can be approximated as g − m modes. Squares indicate for each evolutionary stage the ν_{1, g, max} frequency defined in the previous section, for each of the colourcoded evolutionary stages reported on the HertzsprungRussell diagram in the righthand panel of Fig. 9. The middle panel is composed of the superposition of the magnetic splittings (equivalent to panel a in Fig. 5) corresponding to a field of amplitude B_{0} = 1 MG for the m = 1 mixed modes at the given evolutionary stages, as a function of their unperturbed frequencies. We notice that for stars ascending the RGB, the value of ν_{1, g, max} is globally monotonously decreasing. As a consequence, the magnetic splitting value of g − m modes increases as the star evolves on the RGB (the base of the RGB corresponds to ν_{1, g, max} ≲ 550 μHz for the considered reference star). The p − m modes, which are less affected by magnetism generally, are also more split as they evolve along the RGB as seen by the minima on the curve of the middle panel. The increase of δν_{mag, core, m} with the evolution of the star on the RGB at a given magnetic field amplitude should make the detection of the magnetic effect on mixedmode (g − m and p − m) frequencies easier for evolved red giants than younger ones on the RGB.
Fig. 9. Left: mixedmode frequencies at different evolutionary stages along the RGB indicated in the right panel. Purple dots represent the individual frequencies of mixed modes, with the shade of the dots allowing to visualise the evolution of the mixedmode order n_{pg} (the order is higher at high frequencies). Squares represent for each evolutionary stage the frequency of the ℓ = 1 gdominated mode closest to ν_{max} (see Fig. 10 for the determination of this frequency). Middle: m = 1 magnetic splitting of ℓ = 1 mixed modes located inside [ν_{max} − 7Δν : ν_{max} + 7Δν] for each of the evolutionary stages denoted in the right panel. Squares represent for each evolutionary stage the frequency of the ℓ = 1 gdominated mode closest to ν_{max} (ν_{1, g, max}, see Fig. 5) Right: HertzsprungRussell diagram of the M_{⋆} = 1.5 M_{⊙}, Z = 0.02 star ascending the RGB. Coloured stars indicate the different evolutionary stages represented in the other panels. 
Similar results are obtained for different magnetic field amplitudes B_{0}, ranging from 0.1 to 10 MG, which may arise from past dynamo events (see Sect. 2.3). As a summary, the colour map in the left panel of Fig. 10 represents values of m = 1, g − mmodefrequency splittings due to magnetism only, as a function of the frequency ν_{1, g, max}, which is a proxy for the evolutionary stage along the RGB, and of the applied magnetic field amplitude. The critical value of the magnetic field B_{c}/1000 estimated by Eq. (29) for g − m modes at all points along a star’s evolution delimits the top of the colour map. Specifically, results above this line are beyond the domain of applicability of the perturbative analysis and have been removed. The magnetic signature increases as the star evolves, as deduced from Fig. 9. The magnetic signature also increases with the magnetic field amplitude, as expected from the dependency in the magnetic splitting expression, and as already deduced from Fig. 9.
Fig. 10. Left: m = 1 mixedmode mean frequency perturbation due to magnetism only (δν_{mag} in μHz) as a function of mixedmode frequencies (ν_{1, g, max} in μHz) and the magnetic field amplitude in Gauss for a M_{⋆} = 1.5 M_{⊙}, Z = 0.02. Typical limiting frequencies are indicated: black dashed lines represent the frequency resolution in the PSD corresponding to 4 years of Kepler continuous observations (∼8 nHz), 3 years of PLATO two main fields (∼11 nHz), 1 year of TESS continuous viewing zone data (∼30 nHz), a typical line width of radial modes, which are an upper limit for the line width of dipolar mixed modes (∼0.12 μHz), 2 months of PLATO stareandstep observations (0.19 μHz), and 27 days of TESS individual sectors (0.37 μHz). The coloured area is limited by the critical magneticfield amplitude range for the validity of the perturbative analysis. We choose to limit the results at B_{c, g}/1000, as defined in Sect. 3.5. Red lines represent the expected magneticfield amplitude by considering magneticflux conservation from the end of the MS, with original fossil field amplitudes of 0.1 and 1 MG, as detailed in Sect. 2.3. Right: same as left panel, but simplified with the colour map replaced by typical limiting frequency values corresponding to Kepler data. The yellow area corresponds to the combination of ν_{1, g, max} and magnetic amplitude for which magnetic splittings cannot be visible in Kepler observations. The blue area corresponds to small magnetic shifts that should be visible in Kepler observations. The purple area corresponds to very large magnetic shifts, larger than the typical line width of radial modes. 
4.5. Detectability of the magnetic signature
For comparison, we represent the position in the diagram of different typical observational frequencies: the frequency resolutions of the Kepler, TESS, and PLATO data, and the typical value of ℓ = 0 mode line widths as estimated by Vrard et al. (2017) and Mosser et al. (2018). The line width of ℓ = 1 mixed modes of normal amplitude was shown by Benomar et al. (2014) to be linked to that of ℓ = 0 modes by
Here, this radial mode line width is considered to be an upper limit for mixedmode line width values. All these typical frequencies provide constraints on the minimum field amplitude needed at each evolutionary stage in order for magnetic effects to be visible in asteroseismic observations. When considering a star from the continuous viewing zone (CVZ) of the TESS satellite, meaning that the star has been observed for about 1 year, with ν_{max} ≃ 172 μHz, the lower bound value for the magnetic field amplitude to have detectable signatures is about 0.4 MG. In the righthand panel, the same diagram is shown but simplified, emphasising the (B_{0}, ν_{1, g, max}) combinations for which the magnetic signature should be either easily detectable (purple area), detectable (blue area), or undetectable (yellow area) based on 4 years of observations, such as those by Kepler.
For simplification, we consider the magnetic signature to be observable when δν_{mag} is larger than the data frequency resolution (δf). However, this threshold should be discussed and considered carefully as described by the following criteria:
1. It usually takes a few δf for a signal to be detectable within observational data because of the spreading of the signal over several frequency bins. However, typical line widths of g − m mixed modes are very small, as measured from Eq. (51) where ζ → 1. For the typical g − m mode located at ν_{1, g, max}, we estimate ζ ≃ 0.98 from Fig. 5, which leads to an estimation of g − m mode line widths of about 2.4 nHz by considering the typical line width of radial modes of 120 nHz (e.g. Mosser et al. 2018). The minimal frequency resolution of the Kepler satellite being greater than ∼7.9 nHz, g − m mixed modes are not resolved, and a dilution factor must be considered (Dupret et al. 2009; Mosser et al. 2018). This results in a lower limit for g − m mode line widths of 2δf/π, evaluated at 5 nHz for the Kepler fouryear observations (see also Mosser et al. 2018, for the estimation of ℓ = 1mode line widths for a typical red giant). We conclude that the line width of g − m modes is of the order of magnitude of the resolution of the data. Therefore, the detection limit of the magnetic signature at δν_{mag, core, m} ≃ δf is consistent. In contrast, even though p − m modes have a larger line width of about 70 nHz, the magnetic effect has a much smaller amplitude, and these are therefore much more difficult to detect. We keep the detection limit at δν_{mag, core, m} ≳ δf which is pertinent for g − m modes, but one should use this lower boundary with care especially when looking at pdominated mixed modes.
2. In addition to the magnetic splitting, one should also add the shift due to the rotation of the star from Eqs. (16) and (17), whose measurements come with their own uncertainties. Thereby, the magnetic effect should be large enough for its signature to be easily distinguishable from rotational pattern adjustment errors. These errors have been estimated by Mosser et al. (2018) at about Δν/200, corresponding to about 1.3 times the Kepler typical frequency resolution. Once again, the lower limit of detectability evaluated at δf may be too small, especially in the case of noisy data.
3. When adding rotational perturbations, the magnetic signature can only be measured through the use of the asymmetry measure δ_{ℓ, m} defined in Eq. (50), as magnetic and rotational effects add up. For the effect of magnetism to be detectable, the criterion is changed from δν_{mag, core, m} ≳ δf to
in the presence of rotation. From a perfectly symmetric rotational ℓ = 1 triplet, at given frequencies [ − δν_{rot, 1}, 0, δν_{rot, 1}], perturbations by magnetism produce the shifts [δν_{mag, 1}, δν_{mag, 1}/2, δν_{mag, 1}] (see Sect. 4.1). The criterion from Eq. (52) may then be rewritten as
with ν_{0, ℓ, 0} the unperturbed frequency of the m = 0 component. Simplifying Eq. (53) leads back to the equation δν_{mag, core, m} ≳ δf. The chosen minimum limit of detection is thus conserved when the star rotates.
4. For a few red giants, a departure of the symmetric rotational triplet due to buoyancy glitches can be observed, which is caused by strong chemical gradients generated by the first dredgeup and left behind by the retreating envelope (Cunha et al. 2015, 2019; Jiang et al. 2020). Mosser et al. (2018) shows that KIC 3216736 is the only red giant among the 200 studied that exhibits buoyancy glitches, with only its m = 0 component visible. As deduced analytically by Cunha et al. (2015, 2019), buoyancy glitches are very rare on the RGB, and glitchinduced oscillation variation occur only at the luminosity bump. As a consequence, glitches will be neglected in the rest of our study.
5. If the star rotates fast enough, secondorder and higherorder asymmetric rotational perturbations of the centrifugal and Coriolis accelerations can affect the symmetric rotational frequency pattern (Dziembowski & Goode 1992; Suárez et al. 2006). Such secondorder perturbations should affect both ℓ = 1 and ℓ = 2 mixed modes. However, Deheuvels et al. (2017) emphasise that a measured core rotation rate of Ω_{c} ≃ 710 nHz along with a fivetimessmaller envelope rotation rate are much too low for secondorder rotational effects to significantly contribute to rotational splittings. Therefore, secondorder rotational effects should not produce significant asymmetric perturbations in the spectrum of red giants or subgiants. Moreover, we consider only rotation as a perturbation in the case of red giants as they are considered as slow rotators. Indeed, Ouazzani et al. (2013) show that the effect of rotation can be expressed by perturbative calculations when Ω/(2π)/(ΔP/P^{2})≲2, which is the case for red giants according to core rotation measurements by Gehan et al. (2018). For the study of rapid rotators with nonperturbative developments, we refer to Prat et al. (2019) and van Beeck et al. (2020) (for gravity modes), and to Reese et al. (2006) (for acoustic modes).
6. Latitudinal differential rotation may also induce asymmetries, but as evaluated through nonperturbative calculations with the Adiabatic Code of Oscillation including Rotation (ACOR) by Deheuvels et al. (2017), typical latitudinal differential rotation profiles lead to a very small asymmetry measure of ∼10^{−3}. This corresponds to δν_{mag, core, m} ≃ 2 × 10^{−3}δν_{rot, core, m}, evaluated at δν_{mag, core, m} ≃ 0.5 nHz for the rotational splittings associated with core rotation of 0.5 μHz. Given the frequency resolution of asteroseismic data ranging from ∼7.9 (Kepler 4 years data) to ∼380 nHz (TESS 1 month data), and considering typical latitudinal differential rotation inside solarlike stars, the effect of latitudinal differential rotation on the symmetry of the mixedmode pattern is therefore negligible.
7. Neardegeneracy effects occur by the combination of rotation and mode mixing: when two mixed modes with the same (ℓ, m) combination have frequencies that are too close to each other (i.e. the frequency spacing between the two mixed modes is smaller than the rotation rate), their frequencies are perturbed (Dziembowski & Goode 1992; Suárez et al. 2006). We refer to the complete study of Deheuvels et al. (2017) for the theoretical development of neardegeneracy effects on the asymmetry of rotational multiplets. We emphasise the fact that neardegeneracy effects produce increasing asymmetry measure δ_{ℓ, m} when ℓ increases (δν_{degeneracy, ℓ = 1} ≪ δν_{degeneracy, ℓ = 2}). This can be interpreted as being due to the fact that the frequency separation between two consecutive ℓ = 1 mixed modes is much larger than the separation between two consecutive ℓ = 2 mixed modes. In the case of KIC 7341231 studied by Deheuvels et al. (2017), no asymmetries were found in the ℓ = 1 triplet, whereas ℓ = 2 multiplet asymmetries are δ_{2, 2}≳0.14. In order to disentangle neardegeneracy from magnetic effects, measures of ℓ = 1 and ℓ = 2 mixed mode asymmetries are effective. More precisely, an asymmetry measure such that δν_{ℓ = 1} ≳ δν_{ℓ = 2} is a clear indicator that the magnetic effects are larger than those of neardegeneracy.
In conclusion, (B_{0}, ν_{1, g, max}) areas delimited in the right panel of Fig. 10 should be used with care. In addition to the characteristic frequency positions, expected magneticfield amplitudes at given evolutionary stages from Eq. (29) are represented by red lines in each panel, considering amplitudes at the end of the MS of B_{ms} = 0.1 and 1 MG. The maximum expected magnetic amplitude is represented by the B_{ms} = 1 MG upper red line on the diagram (magnetic field amplitude along the evolution; we recall that this results from the conservation of the magnetic field flux from the end of the MS as explained in Sect. 2.3). We note that the value of B_{ms} = 0.01 MG corresponding to the equipartition regime leads to magnetic field amplitude values that are too low on the RGB to be represented on the diagram. From Fig. 10 we conclude that the more evolved the star the easier it is to detect its internal magnetic signature. More specifically, magnetic fields of B_{MS} ≃ 1 MG at the end of the MS produce frequency perturbations that are too small to be detected inside young RG with ν_{max} ≳ 475 μHz and are too large for the perturbative analysis to be valid for older RG with ν_{max} ≲ 190 μHz. Globally, the expected range for the magneticfield amplitude to reach sufficiently large values for the magnetic signature to be observable depends upon the duration of the observing campaign; for example, the red lines crossing the different Kepler boundaries are shown in the right panel of Fig. 10.
4.6. Stretched spectrum
As a way of summarising previous results concerning the asymmetry of the perturbed rotational triplet and the amplitude of this perturbation by magnetism, we illustrate the stretched spectrum corresponding to perturbed frequencies of mixed modes in Fig. 11. This visualisation technique was initially developed by Vrard et al. (2015) to estimate the observational period spacing (ΔΠ_{1}) of mixed modes, and it is quite often employed to estimate the internal rotation rate of red giants (e.g. Gehan et al. 2018). In these diagrams (see Fig. 11), ℓ = 1 mixedmode frequencies are represented by coloured dots (green: m = −1, red: m = 0, blue: m = 1) as a function of the corrected period τ modulo ΔΠ_{1}, defined via the differential equation:
Fig. 11. Synthetic stretched period echelle diagram built from Eq. (54), showing modelled mixedmode frequencies affected by rotation with Ω_{core} = 0.5 μHz and magneticfield strengths increasing left to right with B_{0} = 0.1 MG, B_{0} = 0.2 MG, B_{0} = 0.5 MG, and B_{0} = 1 MG in the radiative interior. Purely rotationally split components are given for reference by background crosses on the first two panels. Coloured dots indicate the positions of the ℓ = 1 mixed modes of azimuthal order m ∈ {1, 0, −1}, respectively blue, red, and green dots. The horizontal black line delimits domains for which the rotational splitting is small compared to (see Gehan et al. 2018, for more details). Synthetic patterns above this line are typically representative of patterns of mixed modes for subgiants and early giants, and those below are associated with evolved red giants. The magnetic effect is small compared to rotational effects in white areas, is of the same order of magnitude as in orange areas, and is larger than δν_{rot, core, m} in red areas. 
with ζ defined by Eq. (34), and ν the observational frequency of the mode. For a rotating star without magnetism, the period separation between two mixedmodes of the same azimuthal order m is given by integrating Eq. (54):
with δν_{rot, core, m} the rotational perturbation of the mode due to the core rotation. We show in Appendix D that this stretched period spacing can be rewritten in the presence of magnetism as
with δν_{mag, core, m} the magnetic perturbation in the core of the star. The method described by Gehan et al. (2018) to estimate the rotation period of the star can therefore still be applied to magnetised stars due to the similarities between Eqs. (55) and (56). However, the resulting value is no longer an estimate of δν_{rot, core} but rather of (δν_{rot, core}+δν_{mag, core}). A second analysis step is necessary in order to separate the rotational and magnetism signatures. This is accomplished through the measurement of the asymmetry of the multiplet described in Sect. 4.4.
In Fig. 11 we represent four stretched spectra. From left to right, we increase the core magnetic field amplitude (B_{0} ∈ {0.1, 0.2, 0.5, 1} MG). In the first two panels, the positions of nonmagnetised modes are represented by faded crosses for comparison. They are constructed following the method described in Gehan et al. (2018). As in the case of nonmagnetised rotating stars, we obtain three ridges, corresponding to the m ∈ { − 1, 0, 1} components of the mixed modes. We can see in the left panel that for a weak magnetic field the rotational ridges are nearly identical to the magnetically influenced ones, with the triplet being slightly shifted towards higher frequencies as expected given our discussion in Sect. 4.4. The horizontal black line separates two rotational regimes. Below this line the effect of rotation is moderate, where the resulting rotational splitting of the modes is of the order of . This regime is associated with crossing of the three ridges in the stretched spectrum, wherein the rotational splitting leads to the crossing of mixed modes with different n_{g}. In contrast, above the black line, the core rotation rate is small enough so that the individual multiplets do not overlap with each other.
We now construct criteria similar to those in Sect. 4.5 characterising the magnetic effect on the rotational triplet. By studying the m = −1 component of the triplet, one can compare the effect of rotation and magnetism on the frequency of the modes. Indeed, magnetism and rotation have opposite effects on the frequency of the m = −1 component. We conclude that when the stretched period of the m = −1 mode increases with the frequency of the mode, magnetic effects dominate rotation effects. We identify three different regimes depending on the curvature of the m = −1 ridge of the stretched spectrum:
1. When the green ridge representing m = −1 mode components is concave, the effect of magnetism is negligible relative to rotational effects. This scenario is indicated by the white areas in Fig. 11.
2. When δν_{mag} approaches δν_{rot}, the curvature of the m = −1 ridge reverses and becomes convex. At that point, magnetism and rotation have comparable effects on mixedmode frequencies. This scenario corresponds to the orange areas in Fig. 11.
3. When the magnetic field is even larger, its effects can dominate the rotational effects. Such instances are shown as red areas in Fig. 11. One must take caution in this regime, as the magnetic field amplitudes may be large enough to violate the firstorder approximation considered in this paper.
We choose not to indicate the frequency values on the yaxes because the position of the ridges of the stretched spectrum depends on the choice of the integration constant, and thus on the minimum frequency we consider, where we set all τ values to zero.
5. Magnetic effect on mixed mode frequencies of subgiants
Studying magnetic effects on mixed modes inside subgiants is more complicated than studying those occurring on the RGB because of the transition from p − m to g − m modes that occurs during this evolutionary stage. As a consequence, the driving terms listed in Appendix C are no longer fully represented by Eq. (39), and the mode inertia is no longer simply Eq. (36) during the subgiant stage (Hekker & ChristensenDalsgaard 2017), leading to variations in the magneticsplitting patterns.
5.1. Evolution of the magnetic splitting during the subgiant stage
The subgiant phase is much shorter than the RGB, especially for intermediatemass stars; it lasts ∼0.1 Gyr for a M_{⋆} = 1.5 M_{⊙}, Z = 0.02 star, during which the frequency of maximum power varies between approximately 700 and 500 μHz. Figure 12 represents the same three panels as Fig. 9 at this earlier stage, which is indicated by the position of the considered stars in the right panel. In the left panel, we see that the value of ν_{1, g, max} is no longer monotonously decreasing as it was during the RGB. In the middle panel, the amplitude of the magnetic signature for modes contained in the [ν_{max} − 7Δν : ν_{max} + 7Δν] range are represented for all the considered evolutionary states by the blackdotted lines. Green squares indicate the position of ν_{1, g, max}. We observe that for evolved subgiants (and early red giants, e.g. stars older than ∼2.675 Gyr), the magnetic signature is consistent with the asymptotic theoretical pattern detailed in Sect. 4.3. However, younger red giants have larger magnetic frequency splittings at ν_{1, g, max}, and even more substantial signatures for p − m modes, as seen in the middle panel of Fig. 12. We thus investigate the asymptotic regime where highradialorder p − m modes dominate the frequency spectrum.
Fig. 12. Left: a range of mixed mode frequencies is shown for different evolutionary stages along the subgiant phase (and early RGB). The individual stages are indicated on the evolutionary track plotted in the right panel. Purple dots represent the individual frequencies of mixed modes, with the shade of the dots allowing to visualise the evolution of the mixedmode order n_{pg}. Squares represent the frequency of the ℓ = 1 gdominated mode closest to ν_{max} for each evolutionary stage. We refer to Fig. 5 for the determination of this frequency. Middle: magnetic splitting of ℓ = 1, m = 1 mixed modes located inside the range [ν_{max} − 7Δν : ν_{max} + 7Δν] for each of the evolutionary stages denoted in the right panel. Squares represent the frequency of the ℓ = 1 gdominated mode closest to ν_{max} (see Fig. 5 for further reference). Right: HertzsprungRussell diagram of a M_{⋆} = 1.5 M_{⊙}, Z = 0.02 star evolving on the subgiant phase. Coloured stars indicate the different evolutionary stages represented in the other panels. 
5.2. Analytic expression of the magnetic perturbations for highradialorder p–m modes
In the approximation of acoustic modes, the local radial displacement is much larger than the horizontal one (ξ_{h} ≪ ξ_{r}). The highradialorder modes have a small wavelength, allowing us to perform an asymptotic JWKB analysis, where terms with highorder ξ_{r} derivatives dominate. This approximation is valid for highfrequency, highorder mixed modes, for which the splitting components can be expressed as:
In the case of highorder acoustic modes, the dominant magnetic ratio becomes:
At high frequency, where ω_{0} ≫ S_{l} ≫ N, Eq. (35) may be written as
Taking S_{l} = S_{l, min} × g(r), with S_{l, min} the minimum value of the Lamb frequency inside the acoustic cavity and g(r) containing all radial dependencies of S_{l}, the solution can be expressed as
We thus arrive at the following expression for the amplitude scaling of the radial displacement:
Therefore, for a rapidly oscillating radial eigenfunction (), and i ∈ {θ, φ},
By employing the scaling for ξ_{r} in Eq. (61), we can see that
and
For highfrequency p − m modes, the magnetic splitting is therefore proportional to the unperturbed frequency of the mode at first order. This asymptotic behaviour explains the average rise in the measured δν_{mag, core, m} with ν_{1, g, max} at high frequencies (early subgiants; see Fig. 12). As for red giants, we check this power law for a typical M_{⋆} = 1.5 M_{⊙}, Z = 0.02, ν_{max} = 750 μHz. In Fig. 13 the normalised magnetic splitting of the ℓ = 1, m ∈ {1, 0, −1} modes within the [ν_{max} − 7Δν : ν_{max} + 7Δν] range is represented by the green dotted line. The subgiant mixedmode frequency space is less populated, as p − m modes completely dominate the spectrum. We observe a g − m mode among p − m modes at ∼565 μHz. The black line corresponds to the normalised power law describing the p − mmode magnetic signature from Eq. (67), proportional to the mixedmode unperturbed frequencies ν_{0} = ω_{0}/(2π). We confirm that the modelled p − m modes follow this ν_{0} asymptotic trend well at high frequency. We refer to Mathis et al. (2021) for a complete study of the asymptotic behaviour of p − mmode frequencies in the presence of magnetism.
Fig. 13. Normalised magnetic splittings of the p − m modes versus mixedmode frequencies computed for a simulated M_{⋆} = 1.5 M_{⊙}, Z = 0.02, ν_{max} ≃ 750 μHz subgiant. The normalised power law following ν is superimposed in black. 
5.3. Detectability of magnetic signature along the subgiant stage
As done for stars on the RGB (see Sect. 4.5), we now investigate the detectability of such magnetic signatures within the power spectrum density of subgiants. As this evolutionary stage is a transition period between p − m and g − m modes, we choose to represent the detectability of the two asymptotic regimes p − m and g − m in Fig. 14. In contrast with Fig. 10, in which we follow the ascension of the RGB using the ν_{1, g, max} proxy, we choose to study the star during the SG stage by following its age, as ν_{1, g, max} does not evolve monotonously from consecutive subgiant stages. The top panels of Fig. 14 represent the amplitude of the magnetic splitting (not including rotational effects) affecting ℓ = 1, m = 1 g − m modes along the evolution of the star during the SGB. The abscissa are reversed in order to retain the same orientation as used in Fig. 10. From right to left in the diagrams, we observe a decrease in the magnetic effect on g − m modes as the star evolves along on the SGB, followed by an increase of the magnetic signature as the star reaches the RGB. The topright panel delimits regions of detectability in the (Age, B_{0}) space. For a similar field amplitude, the magnetic signatures of the g − m modes are harder to detect in the power density spectra of subgiants than in those of a red giant. We observe the same change of slope in the colour map for p − m modes (see the bottom panels of Fig. 14), which are the majority among the mixed modes in subgiants. This change is due to the switch of asymptotic regime. In all panels, the red lines delineate the expected magnetic field amplitude evaluated by conservation of the magnetic flux inside the radiative interior from the end of the MS. By comparing the relative position of these red lines with the frequencies characterising asteroseismic data, if such magnetic field amplitudes (B ≲ 1 MG) are indeed present in the subgiants’ radiative interiors, their effect on mixedmode frequencies is not detectable with data from any current satellite. In order to reach the frequency resolution needed to be able to identify magnetic effects on mixed mode frequencies during the SGB (δf ≲ 10^{−3} μHz), the observation duration must be longer than about 30 years. Therefore, we do not present further analysis of the magnetic signature on mixedmode frequencies for subgiants.
Fig. 14. Top panels: same as Fig. 5 but for younger stars (ν_{max} ≳ 350 μHz) with the abscissa changed to the age of the star (reversed to be consistent with the convention used in Fig. 5). Bottom panels: same as top panels but considering p − m modes instead of g − m modes. 
6. Angular momentum transport by fossil magnetic fields in evolved solarlike stars
The transport of angular momentum inside stars is a consequence of internal dynamical mechanisms. Understanding the transport of angular momentum inside stars therefore allows us to understand their global dynamical evolution. Current observational constraints on the transport inside stars begin with an estimation of the profile of their internal rotation rate (e.g. Beck et al. 2012). Constraining the core rotation rate also helps to estimate the surface rotation rate (Gallet & Bouvier 2013; Spada & Lanzafame 2020), and the measurement of the surface rotation rate provides an estimate of the age of the star during the MS through gyrochronology (e.g. Barnes 2003, 2010; Angus et al. 2015). Stars with magnetic activity may possess signatures of their surface rotation rate in their PSD through the periodic reduction of their brightness by a few percent due to magnetised dark spots at the surface (e.g. McQuillan et al. 2013; van Saders & Pinsonneault 2013; Mathur et al. 2014; García et al. 2014b; Ceillier et al. 2017; Santos et al. 2019). However, internal rotation rate measurements are more difficult to obtain in general; for instance even the rotation rate of the core of the Sun has yet to be measured, as it requires the presence of g modes to efficiently probe the deepest layers of the radiative interior (e.g. García et al. 2008; Mathur et al. 2008). In the case of subgiants and red giants that possess mixed modes, the internal rotation rate of a restricted sample of evolved stars has recently been measured (e.g. Deheuvels et al. 2012, 2020; Mosser et al. 2012; Gehan et al. 2018). As detailed in Sect. 2.3.2, the observed rotationrate ratio (Ω_{core}/Ω_{env} ≃ 5 − 10) between the core and the envelope of evolved solarlike stars (Deheuvels et al. 2012, 2014; Mosser et al. 2015; Vrard et al. 2015; Gehan et al. 2018) is not consistent with the strong contraction of the core after the MS that should lead to much higher ratio values (e.g. Eggenberger et al. 2012, 2019; Ceillier et al. 2013; Cantiello et al. 2014). A dynamical process must be identified to transport angular momentum from the contracting core towards the envelope in order for this observation to be understood. In the recent literature that tackles this important question, most of the attention has been given to the potential effects of Maxwell stresses triggered by unstable fields in stably stratified radiative regions (e.g. Fuller et al. 2019; Eggenberger et al. 2020; Den Hartogh et al. 2020; Jouve et al. 2020). Here we recall the potentially significant efficiency with which a stable axisymmetric field redistributes angular momentum along poloidal field lines, a wellknown result since Ferraro (1937), Mestel & Weiss (1987). We aim to provide a quantitative estimate of the characteristic timescale needed for such a field to redistribute angular momentum for field amplitude that could be detected in seismic data. As such, we are closing the loop: we are looking for seismic signatures of a potential axisymmetric fossil field as a potential candidate to explain the strong angular momentum transport revealed by the observed weak surfacetocore rotation contrast. Given the amplitude of the field that we predict to be able to detect, we compute the characteristic timescale on which it redistributes angular momentum. If it is very short compared to evolution timescales, as expected, this confirms that it is an excellent candidate for the observed strong extraction of angular momentum.
By considering a largescale axisymmetric fossil magnetic field trapped inside the radiative interior of evolved solarlike pulsators, we investigate the impact of such magnetism on the internal rotation profile of the star (see e.g. Mestel & Weiss 1987; Charbonneau & MacGregor 1993, for previous studies on the Sun). We consider the field lines to be closed, without any reconnection with the dynamogenerated field in its convection zone. Such configurations are depicted in Fig. 2, and following Mestel & Weiss (1987), consider the following poloidal current due to the toroidal component of the axisymmetric magnetic field:
The toroidal component of the equation of induction can be written
where B_{p} is the poloidal magnetic field vector. The current has a component perpendicular to B_{p}, and so it exerts a torque that changes the rotation profile according to
The rotation profile is therefore affected in time according to the toroidal component of the momentum equation
Assuming no turbulence, and given that the Ohmic diffusion timescale due to atomic processes is very large, changes to the poloidal component of the magnetic field can be neglected at the leading order. This leads to the partial differential equation that describes the redistribution of angular momentum along the poloidal field lines by Alfvén waves:
If the variation of Ω along the magnetic field lines of B_{p}, which define the coordinate s, is small with r and with the scale of variation of B_{p}, then
From this equation, Alfvén waves transport angular momentum leading to the Ferraro isorotation law (B_{p}⋅∇Ω) = 0 where rotation becomes constant along the poloidal field lines because they are considered to be fixed in time. Given the defintion of our axisymmetric magnetic field, B_{p} scales with B_{0}, and as the magnetic field extent covers all the radiative interior delimited by R_{rad}, the characteristic time for angular momentum transport derived from Eq. (73) is equal to the Alfvén time:
with v_{A} the Alfvén speed defined by Eq. (31). As shown in Ferraro (1937), Mestel (1953), and Mestel & Weiss (1987), this characteristic time to flatten the rotational profile of the radiative interior is very short.
Figure 15 evaluates the characteristic time for magnetism to flatten the rotational profile of the radiative interior of the reference star along the RGB (left panel) and along the SGB (right panel) as a function of ν_{1, g, max} and age, respectively, and as a function of the magnetic field amplitude ranging from 0.1 to 10 MG. As expected, a greater magnetic amplitude implies a shorter characteristic transport time. We do not observe any significant impact of the evolutionary stage of the star on the transport characteristic time during the SGB and RGB. In any case, the characteristic time to flatten the radiative interior of the star is of the order of 1 year and thus such considered magneticfield amplitudes (as represented by red lines in Fig. 15) are very efficient at transporting angular momentum inside the radiative region. With this orderofmagnitude analysis we reconfirm that magnetism is a very good candidate for rapidly transporting angular momentum after the MS. It may actually be too powerful to maintain the slight differential rotation rate observed inside evolved solarlike stars (Eggenberger et al. 2012, 2017; Cantiello et al. 2014).
Fig. 15. Characteristic time of angular momentum transport τ in years as expressed by Eq. (74), depending on the given evolutionary stage and the magnetic field amplitude. The magnetic field amplitude expected at each stage from the conservation of flux from the end of the MS is indicated by red lines for B_{MS} = 0.1 and 1 MG. Left: transport along the RGB with ν_{1, g, max} as the abscissa. Right: transport along the subgiant branch with age as the abscissa. 
The isolatedfield scenario used for these calculations is however a strong assumption; this should be discussed. It has been proposed that the primordial magnetic field buried inside the radiative zone of the Sun inhibits the spread of the tachocline (e.g. Rüdiger & Kitchatinov 1997; Gough & McIntyre 1998; Barnes et al. 1999). Such confined magnetism can explain the quasiuniform rotation rate of the radiative interior of the Sun. However, Brun & Zahn (2006) and Strugarek et al. (2011) showed that we may expect the confined magnetic field to spread by Ohmic diffusion towards the envelope, and to eventually reconnect with the convection zone dynamo field. Due to Ferraro’s law of isorotation, the latitudinal differential rotation of the convection zone would then imprint on the radiation zone. Such a phenomenon is however not observed in the radiative interior of the Sun. On the other hand, from a more theoretical standpoint, it has been argued that with a proper ordering of timescales of the dynamical processes, this spread can be prevented (e.g. Wood et al. 2011; AcevedoArreguin et al. 2013). Indeed, this question of tachocline confinement and its role in global dynamics is still an open question; it is nevertheless still relevant for the radiative interior of red giants. If the fossil field reconnects with the surface dynamo field, we expect the radiative interior to undergo a small amount of differential rotation despite the strong magnetic field.
7. Dependency of the magnetic splitting on stellar parameters
In the above sections, we derive our study from mixed modes of a typical M_{⋆} = 1.5 M_{⊙}, Z = 0.02 star during its evolution from the base of the subgiant stage towards the top of the RGB. In the following section, we explore the impact of a change of mass and metallicity on the previous results. From theoretical studies, stellarevolution models, and observations, we know that massive stars evolve on shorter timescales than less massive ones. It is also well known that metalrich stars live longer. We follow the evolution of stars with different initial masses (M_{⋆} ∈ {1, 1.5, 2, 3} M_{⊙}) and different metallicities (Z ∈ {0.0002, 0.002, 0.02, 0.04}) in the HertzsprungRussell diagram in right panel of Figs. 16 and 17. In the left panels, for each mass we observe the evolution of the magnetic splitting δν_{mag, core, m} of the central g − m mode as the star evolves.
Fig. 16. Left: magnetic frequency splitting calculated at the central mode frequency ν_{1, g, max} versus the age of the different evolutionary stages reported in the right panel for Z = 0.02, M_{⋆} ∈ {1, 1.5, 2, 3} M_{⊙}. Middle: magnetic frequency splitting calculated at the central mode frequency ν_{1, g, max} for the different evolutionary stages and masses reported in the right panel. Right: HertzsprungRussell diagram corresponding to the colourcoded Z = 0.02, M_{⋆} ∈ {1, 1.5, 2, 3} M_{⊙} stars. The star symbols indicate the position of the star at which measurements are reported in the left and middle panels. 
Fig. 17. Same as Fig. 16 but with varying metallicity; we represent M_{⋆} = 1.5 M_{⊙}, Z ∈ {0.04, 0.02, 0.002, 0.0002} stars. 
7.1. Mass dependency
To investigate the effect of stellar mass on magnetic splitting, we consider low and intermediatemass solartype oscillators from the subgiant towards the red giant phases, for which the convective envelope is thick enough for acoustic modes to be excited at the surface. From the left panel of Fig. 16, it takes around 12.5 Gyr for the M_{⋆} = 1 M_{⊙} star to ascend the RGB, whereas the M_{⋆} = 3 M_{⊙} reaches this stage in less than 1 Gyr. As a first result, significant magnetic signatures arise on frequency spectra at an earlier time for more intermediatemass stars than for lighter ones. This first effect is due to the fact that intermediatemass stars take less time to evolve from the SG to the RG stage than lowmass stars.
We also represent the value of the magnetic splittings as a function of the dominant g − m mode frequencies ν_{1, g, max} in the middle panel of Fig. 16. Intermediatemass stars show the transition between the p − m and g − mdominated mixed modes (see Fig. 12) at lower frequencies than lowmass stars. We interpret this mass dependency as follows: the frequency range inside which we can detect mixedmode form is set by the frequency range of acoustic modes that propagate all over the star, strongly correlated with the size of the star. As a result, acoustic modes propagating inside intermediatemass stars have lower eigenfrequencies than those inside lowmass stars for a given evolutionary stage. In other words, the frequency associated with the base of the RGB where g − m modes dominate is lower for intermediatemass stars than for lowmass stars. However, the mass difference does not significantly impact the magnetic frequencysplitting values for a given g − m mode frequency. The small mass dependency seen at low frequency in the middle panel comes from the mixedmode nature of the mode, which is more dominated by its acoustic nature for massive stars than for lowmass stars at a given frequency on the RGB. The consequence is a slightly smaller effect of magnetism on mixedmode frequencies for more massive stars during the RGB.
7.2. Metallicity dependency
In Fig. 17 the same three panels as in Fig. 16 are represented, with the change in mass replaced by a change in metallicity from Z = 0.0002 to Z = 0.04 for a fixed M_{⋆} = 1.5 M_{⊙}. As expected, the more metallic the star, the slower it evolves as seen in the left panel. The effect of metallicity on the magnetic signature at a given mixedmode frequency of the stars is negligible during the RGB. However, the transition from SG to RG arises at lower frequencies for lowmetallicity stars. This leads to a large spread in magneticsignature values at a given frequency (see the middle panel of Fig. 17), similar to what is observed when varying the mass of the star shown in the middle panel of Fig. 16.
8. Comparison with lowamplitude dipolar mixedmode critical field
The magnetic greenhouse effect, as proposed by Fuller et al. (2015) and supported by the study of Lecoanet et al. (2017), consists in a complete transfer of energy of magnetised gravity waves towards Alfvén waves, geometrically trapped inside the core of the star. The resulting oscillation mode is then purely acoustic, with a great loss of power compared to the corresponding mixed mode due to the loss of the gravity wave energy trapped inside the core. Such a mechanism implies the total disappearance of g components in the ℓ = {1, 2} mode regions in the PSD. For this reason, Mosser et al. (2017) led an observational study looking for signatures of mixedmode residuals inside lowamplitude ℓ = 1 regions. The complete rotational mixedmode frequency pattern is adjusted: it reconstructs the observed ℓ = 1 pattern well. These latter authors conclude that even for stars showing lowamplitude inside the ℓ = 1 region, the remaining oscillation modes have a mixed p and g nature. However, Loi & Papaloizou (2020) set a warning about the physical conditions under which the suppression of mixedmode amplitude may arise along with the disappearance of the g components: the authors support that interactions between gravity modes and magnetic fields may result in various behaviours depending on the configuration and strength of the magnetic field, and of the star’s structure and stratification. This latter study brings perspectives concerning the validity of the theory proposed by Fuller et al. (2015), which may be adapted for mixedmode amplitude suppression by magnetism without loosing all the mixed nature of the modes.
This topic being very controversial, we do not intend to form any conclusions as to whether or not the greenhouse effect is the key to mode suppression. However, we provide the comparison of the critical magneticfield amplitude needed for mode suppression –as fixed by Fuller et al. (2015)– with the minimum field amplitudes needed for the magnetic signature of mixedmode frequencies to be detectable in observational data. Table 1 contains the approximated values of the magnetic field corresponding to the minimum detection threshold inside data from Sect. 4.5 (B_{min}) and of the critical field needed for mode suppression to occur according to the theory by Fuller et al. (2015). Results depend on the duration of observation, and on the considered frequency of the modes. In any case, we observe that as the star evolves (corresponding to a decrease of ν_{max}) the two critical magneticfield amplitudes are of the same order of magnitude, whereas the critical field for mode suppression is much larger for lessevolved red giants. This means that if magnetic suppression by the greenhouse mechanism is at work inside RGs, one should also observe magnetic signatures on mixedmode frequencies, especially for young RGs (see also the work of Rieutord 2017).
Comparison of the critical field B_{c, low − amplitude} needed for mode suppression by the greenhouse effect (Fuller et al. 2015) and the minimum magnetic field amplitude leading to frequency shifts of about the frequency resolution inside data.
9. Discussion and perspectives
In this study, we investigate the effects of an axisymmetric mixed poloidal and toroidal fossil field aligned with the rotation axis of the star on mixedmode frequencies during the SG and RG evolutionary stages. We concentrate on the fossil field scenario, in which a buried magnetism is present inside the radiative interior of the stars. This field results from the relaxation of a magnetic field generated from previously active dynamos. By using different force and energy balance regimes, we show that the expected fossil magneticfield amplitude stabilised at the end of the MS should be in the range of [0.1 − 1] MG. This range is independent of the mass of the star, and of the presence or not of a sustainable convective core during the MS. Such field amplitudes are sufficiently large to be pertinent for our study of the effects of magnetism on mixedmode frequencies. Indeed, our perturbative analysis of mixedmode frequencies of evolved solarlike stars provides constraints for the expected field amplitudes inside evolved stars. Such field amplitudes during the RGB may permit the detection of magnetic signatures on g − m mixedmode frequencies. Specifically, the amplitude of magnetically induced frequency shifts is of about ∼0.1 μHz during the RGB when considering field amplitudes between 0.1 and 1 MG, and this shift is larger than the frequency resolution of longterm observations by the Kepler, TESS, and PLATO missions. Depending on the nature (p − m or g − m) of the mixed mode and of its frequency, the magnetically induced frequency shift may also exceed the typical line width of the mode. Along the RGB, the effect of magnetism on g − m modes, which propagate inside the magnetised radiative interior, is larger than on p − m modes, which also propagate inside the convective envelope. This difference between g − m and p − m modes is also visible when evaluating the effect of the rotation of the star, where p − m modes probing the more slowly rotating envelope are less affected by rotation than g − m modes that primarily provide information about the more rapidly rotating core. Our results therefore mostly concern the effect of magnetism on g − m modes, despite the fact that they are more difficult to observe than p − m modes because of their relatively low amplitude in current data sets.
However, during the SG stage, p − m modes dominate the observed spectra and are more affected by magnetism than during the RGB. This is reflected in the large extent of the radiative interior: R_{rad} → R_{⋆} during the SG phase while R_{rad} → 0 during the RGB. As a result, p − m modes probe more highly magnetised plasma when the star is young. Fossil magnetism resulting from active dynamos before or during the MS is however not strong enough throughout the SGB to produce detectable changes either on their p − m or g − m mixedmode frequencies. We do not exclude the possibility of having stronger fields than that obtained by the balanced regimes. For example, Fuller et al. (2015) proposed a mechanism that uses the frequency of magnetogravity waves to estimate field amplitudes as large as 10 MG in KIC 8561221 during the SG phase.
The two g and p dominated mode regimes correspond to two different asymptotic cases: ω ≪ N ≪ S_{l} for g − m modes, and N ≪ S_{l} ≪ ω for p − m modes. Asymptotic power laws corresponding to the magnetic effect on p − m and g − m mode frequencies are provided, where lowfrequency g − mmode magnetic splittings behave as δω_{mag, g} ∼ 1/ω^{−3} whereas p − mmode magnetic splittings behave as δω_{mag, p} ∼ ω. These asymptotic power laws possess amplitude scaling factors, and correlate with the ζ modecoupling function. Complete expressions for these are derived in Mathis et al. (2021). Considering the typical balanced field amplitudes, we argue that, as for rotational perturbations, these firstorder expressions are sufficient to ascertain the effect of magnetism on mixedmode frequencies. As the star evolves on the RGB, the effect of magnetism on mixedmode frequencies gets larger, where very evolved red giants or AGB stars with ν_{max} ≲ 100 μHz may be nonperturbatively impacted by magnetism. For such evolved stars, we refer to the study of Loi (2020).
We show that the considered axisymmetric magnetism, aligned with the rotation axis of the star, acts as a new perturbation of the already present rotational multiplet. The magnetised mixed multiplet with order n_{pg} is made asymmetric by the presence of a magnetic field: all ℓ = 1 and ℓ = 2 multiplet components are shifted towards higher frequencies, with the amplitude of the shifts depending on the azimuthal order m of the mode. We argue that this asymmetry can be distinguished from other asymmetry sources such as nondegenerate effects. Adjustment methods such as described in Vrard et al. (2015), Mosser et al. (2015), and Gehan et al. (2018) may be adapted to investigate magnetic asymmetries. If the magnetic effect is small compared to the rotational impact on mixedmode frequencies, the multiplet is simply shifted towards higher frequencies. However, we may observe crossings of the components between successive n_{pg} mixedmode multiplets if the magnetic field amplitude is sufficiently large, yet a strong core rotation also yields similar crossings.
If such isolated fossil fields exist inside the radiative interior of evolved stars, they would also affect the rotational profile of the radiative zone: we show that within the assumptions of Mestel & Weiss (1987), the rotation is frozen to poloidal field lines within a few years, eventually leading to an almost flat rotational profile inside the radiative interior. A reconnection of the fossil and dynamogenerated magnetic fields in the convection zone at the tachocline may however permit a small amount of differential rotation to persist inside the radiative interior.
Under the hypothesis that g modes may be trapped inside the core of the star as described by Fuller et al. (2015), we also show that magnetic signatures may be detectable in the frequency pattern of a star before the complete suppression of the mode amplitudes, especially in the case of young red giants. This result is of great interest because the magnetic signature within mixedmode frequencies may appear in the PSD with ℓ = 1 modes of normal amplitude.
By combining the Maxwell Faraday and the Maxwell Ampére equations with the Ohm’s law and by considering an ideal plasma with infinite electric conductivity (low Ohmic diffusivity η, see Eq. (1), usually verified inside stars), the induction equation reduces to ∂B/∂t = ∇ ∧ (u∧B), with u the velocity field.
Acknowledgments
We thank the referee for very useful and detailed comments that allow to improve the quality of our study and the article. L. Bugnet, V. Prat, S. Mathis, A. A. Stoul, and K. Augustson acknowledge support from the European Research Council through ERC grant SPIRE 647383. All CEA members acknowledge support from GOLF and PLATO CNES grants of the Astrophysics Division at CEA. S. Mathur acknowledges support by the Ramon y Cajal fellowship number RYC201517697. L. Amard acknowledges funding from the European Research Council (grant agreement No. 682393 AWESoMeStars). We made great use of the megyr python package for interfacing MESA and GYRE codes. In honour of our dear friend and colleague Michael J. Thompson.
References
 Abramowitz, M., & Stegun, I. 1972, Handbook of Mathematical Functions (New York: Dover) [Google Scholar]
 AcevedoArreguin, L. A., Garaud, P., & Wood, T. S. 2013, MNRAS, 434, 720 [Google Scholar]
 Aerts, C., ChristensenDalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology, Astronomy and Astrophysics Library (Dordrecht: Springer Netherlands) [Google Scholar]
 Akgün, T., Reisenegger, A., Mastrano, A., & Marchant, P. 2013, MNRAS, 433, 2445 [Google Scholar]
 Alvan, L., Mathis, S., & Decressin, T. 2013, A&A, 553, A86 [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]
 Angel, J. R. P., Borra, E. F., & Landstreet, J. D. 1981, ApJ, 45, 457 [Google Scholar]
 Angus, R., Aigrain, S., ForemanMackey, D., & McQuillan, A. 2015, MNRAS, 450, 1787 [Google Scholar]
 Appourchaux, T., Belkacem, K., Broomhall, A. M., et al. 2010, A&ARv, 18, 197 [Google Scholar]
 Appourchaux, T., Boumier, P., Leibacher, J. W., & Corbard, T. 2018, A&A, 607, A108 [Google Scholar]
 Arlt, R., Leussu, R., Giese, N., Mursula, K., & Usoskin, I. G. 2013, MNRAS, 433, 3165 [Google Scholar]
 Astoul, A., Mathis, S., Baruteau, C., et al. 2019, A&A, 631, A111 [CrossRef] [EDP Sciences] [Google Scholar]
 Augustson, K., & Mathis, S. 2018, SF2A2018: Proceedings of the Annual Meeting of the French Society of Astronomy and Astrophysics, 113 [Google Scholar]
 Augustson, K. C., Brown, B. P., Brun, A. S., Miesch, M. S., & Toomre, J. 2012, ApJ, 756, 169 [NASA ADS] [CrossRef] [Google Scholar]
 Augustson, K., Brun, A. S., Miesch, M., & Toomre, J. 2015, ApJ, 809, 149 [Google Scholar]
 Augustson, K. C., Brun, A. S., & Toomre, J. 2016, ApJ, 829, 92 [Google Scholar]
 Augustson, K. C., Brun, A. S., & Toomre, J. 2019, ApJ, 876, 83 [Google Scholar]
 Aurière, M., KonstantinovaAntova, R., Charbonnel, C., et al. 2015, A&A, 574, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barnes, S. A. 2003, ApJ, 586, 464 [NASA ADS] [CrossRef] [Google Scholar]
 Barnes, S. A. 2010, ApJ, 722, 222 [NASA ADS] [CrossRef] [Google Scholar]
 Barnes, G., Charbonneau, P., & MacGregor, K. B. 1999, ApJ, 511, 466 [Google Scholar]
 Beck, P. G., Bedding, T. R., Mosser, B., et al. 2011, Science, 332, 205 [Google Scholar]
 Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Bedding, T. R., Mosser, B., Huber, D., et al. 2011, Nature, 471, 608 [Google Scholar]
 Belkacem, K., Marques, J. P., Goupil, M. J., et al. 2015, A&A, 579, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Benomar, O., Belkacem, K., Bedding, T. R., et al. 2014, ApJ, 781, L29 [Google Scholar]
 Benomar, O., Takata, M., Shibahashi, H., Ceillier, T., & García, R. A. 2015, MNRAS, 452, 2654 [NASA ADS] [CrossRef] [Google Scholar]
 Bigot, L., Provost, J., Berthomieu, G., Dziembowski, W. A., & Goode, P. R. 2000, A&A, 356, 218 [NASA ADS] [Google Scholar]
 Braithwaite, J. 2006, A&A, 453, 687 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Braithwaite, J. 2007, A&A, 469, 275 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Braithwaite, J. 2008, MNRAS, 386, 1947 [Google Scholar]
 Braithwaite, J., & Nordlund, A. 2006, A&A, 450, 1077 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Braithwaite, J., & Spruit, H. C. 2004, Nature, 431, 819 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Brown, B. P., Browning, M. K., Brun, A. S., Miesch, M. S., & Toomre, J. 2010, ApJ, 711, 424 [Google Scholar]
 Brown, B. P., Miesch, M. S., Browning, M. K., Brun, A. S., & Toomre, J. 2011, ApJ, 731, 69 [Google Scholar]
 Browning, M. K. 2008, ApJ, 676, 1262 [Google Scholar]
 Brun, A. S., & Browning, M. K. 2017, Liv. Rev. Sol. Phys., 14, A4 [Google Scholar]
 Brun, A. S., & Zahn, J. P. 2006, A&A, 457, 665 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Brun, A. S., Miesch, M. S., & Toomre, J. 2004, ApJ, 614, 1073 [Google Scholar]
 Brun, A. S., Browning, M. K., & Toomre, J. 2005, ApJ, 629, 461 [Google Scholar]
 Brun, A. S., Strugarek, A., Varela, J., et al. 2017, ApJ, 836, 192 [NASA ADS] [CrossRef] [Google Scholar]
 Campbell, C. G., & Papaloizou, J. C. B. 1986, MNRAS, 220, 577 [Google Scholar]
 Cantiello, M., Mankovich, C., Bildsten, L., ChristensenDalsgaard, J., & Paxton, B. 2014, ApJ, 788, 93 [Google Scholar]
 Cantiello, M., Fuller, J., & Bildsten, L. 2016, ApJ, 824, 14 [CrossRef] [Google Scholar]
 Ceillier, T., Eggenberger, P., García, R. A., & Mathis, S. 2013, A&A, 555, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ceillier, T., Tayar, J., Mathur, S., et al. 2017, A&A, 605, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chandrasekhar, S. 1958, Proc. Natl. Acad. Sci., 44, 842 [Google Scholar]
 Chaplin, W. J., Appourchaux, T., Elsworth, Y., et al. 2010, ApJ, 713, 169 [Google Scholar]
 Chaplin, W. J., Serenelli, A. M., Miglio, A., et al. 2020, Nat. Astron., 4, 383 [Google Scholar]
 Charbonneau, P., & MacGregor, K. B. 1993, ASP Conf. Ser., 40, 464 [Google Scholar]
 Christensen, U. R., Holzwarth, V., & Reiners, A. 2009, Nature, 457, 167 [Google Scholar]
 Couvidat, S., Garca, R. A., TurckChize, S., et al. 2003, ApJ, 597, L77 [Google Scholar]
 Cowling, T. G. 1945, MNRAS, 38, 105 [Google Scholar]
 Cunha, M. S., Stello, D., Avelino, P. P., ChristensenDalsgaard, J., & Townsend, R. H. 2015, ApJ, 805, 127 [Google Scholar]
 Cunha, M. S., Avelino, P. P., ChristensenDalsgaard, J., et al. 2019, MNRAS, 490, 909 [Google Scholar]
 Deheuvels, S., García, R. A., Chaplin, W. J., et al. 2012, ApJ, 756, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Deheuvels, S., Doǧan, G., Goupil, M. J., et al. 2014, A&A, 564, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Brandão, I., Silva Aguirre, V., et al. 2016, A&A, 589, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Ouazzani, R. M., & Basu, S. 2017, A&A, 605, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Ballot, J., Eggenberger, P., et al. 2020, A&A, 641, A117 [EDP Sciences] [Google Scholar]
 Den Hartogh, J. W., Eggenberger, P., & Deheuvels, S. 2020, A&A, 634, L16 [CrossRef] [EDP Sciences] [Google Scholar]
 Dikpati, M., & Gilman, P. A. 1999, ApJ, 512, 417 [Google Scholar]
 Domingo, V., Fleck, B., & Poland, A. I. 1995, Sol. Phys., 162, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Donati, J.F., Semel, M., Carter, B. D., Rees, D. E., & Cameron, A. C. 1997, MNRAS, 291, 658 [Google Scholar]
 Duez, V. 2011, Astron. Nachr., 332, 983 [Google Scholar]
 Duez, V., & Mathis, S. 2010, A&A, 517, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Duez, V., Braithwaite, J., & Mathis, S. 2010a, ApJ, 724, L34 [Google Scholar]
 Duez, V., Mathis, S., & TurckChièze, S. 2010b, MNRAS, 402, 271 [Google Scholar]
 Dupret, M. A., Belkacem, K., Samadi, R., et al. 2009, A&A, 506, 57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dziembowski, W. A., & Goode, P. R. 1985, ApJ, 296, 27 [Google Scholar]
 Dziembowski, W. A., & Goode, P. R. 1989, ApJ, 347, 540 [Google Scholar]
 Dziembowski, W. A., & Goode, P. R. 1992, ApJ, 394, 670 [Google Scholar]
 Dziembowski, W. A., & Goode, P. R. 1996, ApJ, 458, 338 [Google Scholar]
 Eggenberger, P., Maeder, A., & Meynet, G. 2005, A&A, 440, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eggenberger, P., Montalbán, J., & Miglio, A. 2012, A&A, 544, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eggenberger, P., Lagarde, N., Miglio, A., et al. 2017, A&A, 599, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eggenberger, P., Deheuvels, S., Miglio, A., et al. 2019, A&A, 621, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eggenberger, P., Den Hartogh, J. W., Buldgen, G., et al. 2020, A&A, 631, L6 [Google Scholar]
 EmeriauViard, C., & Brun, A. S. 2017, ApJ, 846, 8 [Google Scholar]
 Featherstone, N. A., Browning, M. K., Brun, A. S., & Toomre, J. 2009, ApJ, 705, 1000 [Google Scholar]
 Ferraro, V. C. A. 1937, MNRAS, 97, 458 [Google Scholar]
 Fossat, E., & Schmider, F. X. 2018, A&A, 612, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fossat, E., Boumier, P., Corbard, T., et al. 2017, A&A, 604, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Froman, N., & Froman, P. O. 2005, Physical Problems Solved by the Phaseintegral Method (Cambridge, UK: Cambridge University Press), 228 [Google Scholar]
 Fuller, J. 2014, Icarus, 242, 283 [Google Scholar]
 Fuller, J., Lecoanet, D., Cantiello, M., & Brown, B. 2014, ApJ, 796, 17 [Google Scholar]
 Fuller, J., Cantiello, M., Stello, D., Garcia, R. A., & Bildsten, L. 2015, Science, 350, 423 [Google Scholar]
 Fuller, J., Piro, A. L., & Jermyn, A. S. 2019, MNRAS, 485, 3661 [NASA ADS] [Google Scholar]
 Gabriel, A. H., Grec, G., Charra, J., et al. 1995, Sol. Phys., 162, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Gallet, F., & Bouvier, J. 2013, A&A, 556, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 García, R. A., & Ballot, J. 2019, Liv. Rev. Sol. Phys., 16, A4 [Google Scholar]
 García, R. A., TurckChièze, S., JiménezReyes, S., et al. 2007, Science, 316, 1591 [Google Scholar]
 García, R. A., Mathur, S., Ballot, J., et al. 2008, Sol. Phys., 251, 119 [Google Scholar]
 García, R. A., Mathur, S., Pires, S., et al. 2014a, A&A, 568, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 García, R. A., Ceillier, T., Salabert, D., et al. 2014b, A&A, 572, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gaurat, M., Jouve, L., Lignières, F., & Gastine, T. 2015, A&A, 580, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gehan, C., Mosser, B., Michel, E., Samadi, R., & Kallinger, T. 2018, A&A, 616, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gough, D. 1984, Interdiscip. Sci. Rev., 9, 240 [Google Scholar]
 Gough, D. O., & McIntyre, M. E. 1998, Nature, 394, 755 [Google Scholar]
 Gough, D. O., & Thompson, M. J. 1990, MNRAS, 242, 25 [Google Scholar]
 Goupil, M. J., Mosser, B., Marques, J. P., et al. 2013, A&A, 549, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hasan, S. S., Zahn, J. P., & ChristensenDalsgaard, J. 2005, A&A, 444, 29 [Google Scholar]
 Hekker, S., & ChristensenDalsgaard, J. 2017, A&ARv, 25, A1 [Google Scholar]
 Hotta, H. 2017, ApJ, 843, 52 [Google Scholar]
 Huber, D., Chaplin, W. J., Chontos, A., et al. 2019, AJ, 157, 245 [Google Scholar]
 Jiang, C., Cunha, M., ChristensenDalsgaard, J., & Zhang, Q. 2020, MNRAS, 495, 621 [Google Scholar]
 Jouve, L., Lignières, F., & Gaurat, M. 2020, A&A, 641, A13 [EDP Sciences] [Google Scholar]
 Kiefer, R., & Roth, M. 2018, ApJ, 854, 74 [CrossRef] [Google Scholar]
 Lecoanet, D., Vasil, G. M., Fuller, J., Cantiello, M., & Burns, K. J. 2017, MNRAS, 466, 2181 [Google Scholar]
 Loi, S. T. 2020, MNRAS, 496, 3829 [Google Scholar]
 Loi, S. T., & Papaloizou, J. C. B. 2017, MNRAS, 3225, 3212 [Google Scholar]
 Loi, S. T., & Papaloizou, J. C. B. 2018, MNRAS, 477, 5338 [Google Scholar]
 Loi, S. T., & Papaloizou, J. C. B. 2020, MNRAS, 491, 708 [Google Scholar]
 Lund, M. N., Aguirre, V. S., Davies, G. R., et al. 2017, ApJ, 835, 172 [Google Scholar]
 Markey, P., & Tayler, R. 1973, MNRAS, 163, 77 [Google Scholar]
 Marques, J. P., Goupil, M. J., Lebreton, Y., et al. 2013, A&A, 549, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., & Zahn, J. P. 2005, A&A, 440, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., Duez, V., & Braithwaite, J. 2010, Proc. Int. Astron. Union, 6, 270 [Google Scholar]
 Mathis, S., Bugnet, L., Prat, V., et al. 2021, A&A, 647, A122 [EDP Sciences] [Google Scholar]
 Mathur, S., EffDarwich, A., García, R. A., & TurckChièze, S. 2008, A&A, 484, 517 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathur, S., García, R. A., Ballot, J., et al. 2014, A&A, 562, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 McQuillan, A., Aigrain, S., & Mazeh, T. 2013, MNRAS, 432, 1203 [Google Scholar]
 Mestel, L. 1953, MNRAS, 113, 716 [Google Scholar]
 Mestel, L., & Weiss, N. O. 1987, MNRAS, 226, 123 [NASA ADS] [CrossRef] [Google Scholar]
 Michel, E., Baglin, A., Auvergne, M., et al. 2008, Science, 322, 558 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Moffatt, H. K. 2015, J. Plasma Phys., 81, 905810608 [Google Scholar]
 Mosser, B. 2015, EAS Publ. Ser., 73, 3 [Google Scholar]
 Mosser, B., Elsworth, Y., Hekker, S., et al. 2012, A&A, 537, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Vrard, M., Belkacem, K., Deheuvels, S., & Goupil, M. J. 2015, A&A, 584, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Belkacem, K., Pinçon, C., et al. 2017, A&A, 598, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Gehan, C., Belkacem, K., et al. 2018, A&A, 618, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ouazzani, R.M., Goupil, M. J., Dupret, M.A., & Marques, J. P. 2013, A&A, 554, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [Google Scholar]
 Pinçon, C., Belkacem, K., Goupil, M. J., & Marques, J. P. 2017, A&A, 605, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Prat, V., Mathis, S., Buysschaert, B., et al. 2019, A&A, 627, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Putney, A. 1999, 11th European Workshop on White Dwarfs (San Francisco: Astronomical Society of the Pacific) [Google Scholar]
 Rashba, T. I., Semikoz, V. B., TurckChièze, S., & Valle, J. W. 2007, MNRAS, 377, 453 [Google Scholar]
 Reese, D., Lignières, F., & Rieutord, M. 2006, A&A, 455, 621 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rieutord, M. 2017, EPJ Web Conf., 160, 02011 [Google Scholar]
 Roberts, P. H., & Soward, A. M. 1983, MNRAS, 205, 1171 [NASA ADS] [Google Scholar]
 Rüdiger, G., & Kitchatinov, L. L. 1997, Astron. Nachr., 318, 273 [Google Scholar]
 Santos, A. R. G., García, R. A., Mathur, S., et al. 2019, ApJS, 244, 21 [Google Scholar]
 Scherrer, P. H., & Gough, D. O. 2019, ApJ, 877, 42 [Google Scholar]
 Scherrer, P. H., Bogart, R. S., Bush, R. I., et al. 1995, Sol. Phys., 162, 129 [Google Scholar]
 Schmidt, G. D., Vennes, S., Wickramasinghe, D. T., & Ferrario, L. 2001, MNRAS, 328, 203 [Google Scholar]
 Schunker, H., Schou, J., Gaulme, P., & Gizon, L. 2018, Sol. Phys., 293, 95 [Google Scholar]
 Spada, F., & Lanzafame, A. C. 2020, A&A, 636, A76 [CrossRef] [EDP Sciences] [Google Scholar]
 Spruit, H. C. 1999, A&A, 349, 189 [NASA ADS] [Google Scholar]
 Spruit, H. C. 2002, A&A, 381, 923 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Stello, D., Cantiello, M., Fuller, J., Garcia, R. A., & Huber, D. 2016a, PASA, 33, e011 [Google Scholar]
 Stello, D., Cantiello, M., Fuller, J., et al. 2016b, Nature, 529, 364 [Google Scholar]
 Strugarek, A., Brun, A. S., & Zahn, J. P. 2011, A&A, 532, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Suárez, J. C., Goupil, M. J., & Morel, P. 2006, A&A, 449, 673 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Takata, M., & Shibahashi, H. 1994, PASJ, 46, 301 [Google Scholar]
 Tayler, R. 1973, MNRAS, 161, 365 [Google Scholar]
 Tayler, R. 1980, MNRAS, 191, 151 [Google Scholar]
 Thompson, M., ChristensenDalsgaard, J., Miesch, M. S., & Toomre, J. 2003, ARA&A, 41, 81 [Google Scholar]
 Townsend, R. H., & Teitler, S. A. 2013, MNRAS, 435, 3406 [Google Scholar]
 Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial Oscillations of Stars (Tokyo: University of Tokyo Press) [Google Scholar]
 van Beeck, J., Prat, V., Van Reeth, T., et al. 2020, A&A, 638, A149 [EDP Sciences] [Google Scholar]
 van Saders, J. L., & Pinsonneault, M. H. 2013, ApJ, 776, 67 [Google Scholar]
 Villebrun, F., Alecian, E., Hussain, G., et al. 2019, A&A, 622, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vrard, M., Mosser, B., & Barban, C. 2015, EPJ Web Conf., 101, 06067 [Google Scholar]
 Vrard, M., Mosser, B., & Barban, C. 2017, EPJ Web Conf., 160, 04012 [Google Scholar]
 Woltjer, L. 1959, ApJ, 130, 405 [Google Scholar]
 Wood, T. S., McCaslin, J. O., & Garaud, P. 2011, ApJ, 738, 47 [Google Scholar]
 Yu, J., Huber, D., Bedding, T. R., et al. 2016, MNRAS, 463, 1297 [Google Scholar]
 Zahn, J. P., Brun, A. S., & Mathis, S. 2007, A&A, 474, 145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Magnetic scalinglaw regimes
By considering the fluid in the rotating frame under the effect of magnetism (via the Lorentz force) and rotation (via the Coriolis acceleration), assuming stationarity and neglecting viscosity and the centrifugal acceleration, the NavierStokes equation in a convection zone may be written as:
with g the gravitational acceleration.
For this equation, several force balances allow us to evaluate the amplitude of the magnetic field in different regimes. The balance between forces is a partitioning that depends upon the convective Rossby number. It is more subtle than a typical asymptotic analysis because it is meant to encompass three convective Rossby number regimes smoothly:

At high Rossby number (Ro ≫ 1) flows are weakly rotationally constrained but possess a strong dynamo;

for Ro close to one, flows are modestly rotationally constrained, again with a strong dynamo;

for low Ro (Ro ≪ 1), the magnetostrophic regime is reached, corresponding to rotationally constrained and magnetised flows.
It is found in Augustson et al. (2019) that the partitioning of forces is dominated by combinations of the inertial terms (Reynolds stresses), the Coriolis acceleration, and the Lorentz force (or Maxwell stresses) for all Rossby numbers. However, depending upon Rossby number regime, the relative influence of each of these three strongest forces changes.
The force balance described in Augustson et al. (2019) fits the available simulation data well when the force balance is described as
where I is the magnitude of the inertial forces, C that of the Coriolis force, and L that of the Lorentz force. Using a scaling argument about the dominant length scales, velocity, and magnetic field amplitudes at that scale, one can show that in this case
where Ro ≈ v_{l}/(2Ω_{0}l) and l is that dominant length scale, a describes the efficiency of a smallscale local dynamo that depends only on the helical turbulence generated by convection, and b refers to how close to perfect magnetostrophy the system is. In this framework, three specific regimes can be reached:

The magnetostrophic regime (low Rossby number) is reached when the Lorentz force balances the Coriolis acceleration and other forces are much smaller in magnitude. This regime can be achieved for a sufficiently high rotation rate. From Eq. (A.1), we deduce that
with R_{conv} the size of the convective zone, and v_{conv} the convective velocity. This can be rewritten in terms of convective Rossby number Ro = v_{conv}/(2ΩR_{conv}) and kinetic energy density associated with convective motions as:

Dynamo action converts a fraction of the kinetic energy of the convective motions into magnetic energy. In the equipartition regime (high Rossby number regime where the Coriolis acceleration can be neglected), the assumption is thus that the dynamo is efficient enough such that the magnetic energy density ME = B^{2}/(8π) balances the convective kinetic energy density KE of the fluid:

In the modest Rossby number regime, all three forces have roughly the same magnitude. So, ultimately, in the Rossby number regime that is close to unity, all three forces play nearly equal roles. By neglecting inertial forces, considering constant density, and with ∇.v = 0, the curl of Eq. (A.1) can be written
Comparing the three terms composing Eq. (A.7), the balance is reached for
We refer to the Augustson et al. (2019) and Astoul et al. (2019) studies for additional details on the estimation of the buoyancy dynamo field strength. This field has an intermediate value between the equipartition and magnetostrophy regimes, which set the upper and lower boundaries of the expected magnetic amplitude during episodes of internal convective dynamo action.
Appendix B: Magnetic field stable topology
We recall the stable configuration theoretical expression (Duez & Mathis 2010) that is used in our study to represent fossil fields aligned with the rotation axis of the star:
where ψ is the stream function:
with μ_{0} the vacuum magnetic permeability, α a normalisation constant, λ the eigenvalue of the problem that fixes the shape of the magnetic configuration, R_{rad} the radius of the radiative cavity, and
with j_{1} (y_{1}) the firstorder spherical Bessel function of the first (second) kind (Abramowitz & Stegun 1972).
In order for the field to be confined inside the radiative interior of evolved solarlike stars, ψ(r, θ) (and thus A(r)) should go to zero at the radiative–convective boundary located by R_{rad}. There are two options to cancel ψ(R_{rad}, θ) for any value of θ:

Cancelling allows us to confine b_{r} and b_{θ} inside the radiative interior and to keep the field in its more stable configuration (Woltjer 1959; Duez et al. 2010a).

Cancelling y_{1}(λ) allows only b_{r} to go to zero at the radiative boundary. Nonzero B_{θ} generates an azimuthal current sheet that potentially creates instabilities (Duez & Mathis 2010).
We observe in panel a of Fig. B.1 that the Ψ function fails to reach zero at the edge of the convective interior for any value of λ_{1} when trying to suppress . The first minimum of the function ψ(R_{rad}, θ) is found at λ_{1} ≃ 7.62, and leads to a ψ_{λ1}(r, θ) function proportional to the blue curve represented in panel b and to the magnetic field topology represented in panels c and d of Fig. B.1, for which the field is not trapped inside the radiative interior. Such magneticfield configuration is not the most stable one according to the studies of Braithwaite (2008) and Duez & Mathis (2010).
Fig. B.1. Attempt to find the first zero of at the radiative–convective boundary. We consider the typical red giant with M = 1.5 M_{⊙}, Z = 0.02 on the RGB. Panel a: trend of the evolution of the ψ function at r = R_{rad} with the eigenvalue λ. The value λ_{1} closest to the eigenvalue is indicated by the red vertical line. Panel b: shape of the functions , and ψ_{λ1}(r, θ) with λ_{1} ≃ 7.62, the vertical dashed line indicates the surface. Panel c: magnetic field components with λ_{1} ≃ 7.62. Panel d: resulting magnetic field topology with λ_{1} ≃ 7.62. 
By performing instead the analysis with a constant density profile (i.e. searching for zeros of the function ) allows the search for λ_{1} to converge. We therefore conclude that the steep density profile inside red giants prevents the integral from converging to zero easily.
The only remaining option to find the eigenvalue λ_{1} that cancels ψ at the radiative boundary for all θ is therefore to search for the zeros of y_{1}(λ) instead of those of the integral in Eq. (4). Results are represented in Fig. B.2, with λ_{1} ≃ 2.80 the eigenvalue used in this study. This method allows b_{r} to cancel at the radiative boundary, but does not constrain the horizontal component of the field. As shown by panel c of Fig. B.2, angular components are still very close to zero at the radiative–convective boundary. As a result, this magnetic field topology is stable (Duez & Mathis 2010) and can be considered as trapped inside the radiative interior along the evolution.
Fig. B.2. Same legend as Fig. B.1 for the first zero of at the radiative–convective boundary instead. In this case we find the first eigenvalue λ_{1} ≃ 2.80. 
Appendix C: Nonzero average terms of the Lorentz work
In this section, we present all the nonzero terms composing the work of . They either involve only poloidal components of the magnetic field, or toroidal components only. The prime symbol (′) indicates a total derivative, either radial or latitudinal depending on the considered variable.
C.1. Poloidal terms
By using A = [(rb_{θ})′ + b_{r}], the terms of the work of involving poloidal components write:
and the work of is composed of:
C.2. Toroidal terms
The terms that involve the toroidal component of the field are for the work of :
and:
for the work of .
Appendix D: Stretched spectrum spacing Δτ_{m} in the presence of magnetism
In order to investigate the internal rotation rate of evolved solarlike stars possessing mixed modes, Mosser (2015) concentrates on the rotational splitting on g − m modes. We assume that solarlike stars are slow rotators, and therefore neglect the centrifugal acceleration and other secondorder rotational effects which scale as Ω^{2}. From the study of Goupil et al. (2013) the global rotational splitting of a mixed mode is written as
For g − m modes the ζ function is very close to 1, so δν_{rot, m} ≃ δν_{rot, g}ζ ≃ δν_{rot, core}ζ, providing a measure of the core rotation rate of the star. Therefore, the unperturbed frequency of g − m modes can be written as ν_{g} = ν − δν_{rot, g}ζ with ν the measured frequency in the presence of rotation.
For a rotating and magnetised star, we rewrite the unperturbed frequency as a function of the observed frequency through:
with δν_{mag, g} the frequency perturbation due to magnetism only. The validity of this expression comes from the asymptotic study contained in the paper (Mathis et al. 2021). As a consequence, the period of mixed modes of azimuthal component m varies as
leading to
With the use of a limiting case where (δν+δν_{mag, core, m})/ν → 0, one obtains:
The period spacing between two consecutive mixed modes is therefore written as:
with ΔP the period spacing between two consecutive m = 0 axisymmetric mixed modes. Noting that Δτ_{m} = ΔP_{m}/ζ and by using the fact ΔP can be written as ΔP = ζΔΠ_{1}, Eq. (D.6) becomes:
Appendix E: MESA inlist
In this appendix we report the MESA inlist used to calculate the stellar evolution models of the 1.5 M_{⊙}, Z = 0.02 star:
&star_job ! begin with a premain sequence model create_pre_main_sequence_model = .true. / !end of star_job namelist
&controls ! starting specifications initial_mass = 1.5 ! in Msun units initial_z = 0.02 use_Type2_opacities = .true. Zbase = 0.02 max_model_number = 900 !____________ MISC profile_interval = 10 history_interval = 10 max_num_profile_models = 1000 calculate_Brunt_N2 = .true. star_history_name = 'history.data' profile_data_prefix = 'profile' profiles_index_name = 'profiles.index' set_min_D_mix = .true. min_D_mix = 1d1 !____________ Output pulse files for GYRE pulse_data_format = 'GYRE' write_pulse_data_with_profile = .true. add_center_point_to_pulse_data = . true. add_double_points_to_pulse_data = . true. !____________ WIND cool_wind_RGB_scheme = 'Reimers' cool_wind_AGB_scheme = 'Blocker' RGB_to_AGB_wind_switch = 1d4 Reimers_scaling_factor = 0.2 Blocker_scaling_factor = 0.5 use_accreted_material_j = .true. accreted_material_j = 0 !_____________ OVERSHOOTING overshoot_scheme(1) = 'exponential' overshoot_zone_type(1) = 'any' overshoot_zone_loc(1) = 'any' overshoot_bdy_loc(1) = 'any' overshoot_f(1) = 0.015 overshoot_f0(1) = 0.004 !_____________ MESH mesh_delta_coeff = 0.7 varcontrol_target = 0.7d3 predictive_mix(1) = .true. predictive_superad_thresh(1) = 0.005 predictive_avoid_reversal(1) = 'he4' predictive_zone_type(1) = 'any' predictive_zone_loc(1) = 'core' predictive_bdy_loc(1) = 'top' dX_div_X_limit_min_X = 1d4 dX_div_X_limit = 5d1 dX_nuc_drop_min_X_limit = 1d4 dX_nuc_drop_limit = 1d2 / ! end of controls namelist
All Tables
Comparison of the critical field B_{c, low − amplitude} needed for mode suppression by the greenhouse effect (Fuller et al. 2015) and the minimum magnetic field amplitude leading to frequency shifts of about the frequency resolution inside data.
All Figures
Fig. 1. Magnetic configuration schemes (not to scale) following the evolution of low and intermediatemass stars from the main sequence. The represented configuration corresponds to the axisymmetric field topology considered in our study. It presents both a poloidal (blue lines) and a toroidal (red lines) component, with the axis of symmetry aligned with the rotation axis of the star. 

In the text 
Fig. 2. Mixed poloidal (black lines) and toroidal (colour scale) magnetic field modelled using the formalism of Duez & Mathis (2010); see Eq. (2). Values are normalised by the maximum radial field amplitude. The field is confined inside the radiative region of the star, the extent of which is given by the corresponding evolution model MESA and indicated by the black circle. Left to right: cases of a 1.5 M_{⊙}, Z = 0.02 star from the SG phase to RGB, with ages of 2.65, 2.75, and 2.80 Gyr, and effective temperatures of 5967, 5075, 4900 K, respectively. 

In the text 
Fig. 3. Top: magnetic field amplitude at the base of the convective envelope along the evolution of a M_{⋆} = 1 M_{⊙} star, with solar metallicity (Z = 0.0134), uniformly rotating with an initial rotation rate during the disklocking of the PMS of Ω = 2.5 μHz. The blue, red, and orange lines respectively represent the field estimated by considering the magnetostrophic, buoyancy, and equipartition regimes. The black dashed line represents the location of the limit radius between the internal radiative layer and the convective envelope. Bottom: same as top panel at the top of the convective interior along the evolution of the star. The black dashed line represents the location of the limit radius between the convective core (when existing) and the radiative layer. 

In the text 
Fig. 4. Same as Fig. 3 but for a M_{⋆} = 1.4 M_{⊙} star. 

In the text 
Fig. 5. Frequency splittings (δν) calculated for mixed modes frequencies of a M_{⋆} = 1.5 M_{⊙}, Z = 0.02, ν_{max} ≃ 172 μHz red giant, corresponding to (a) magnetic splitting, (b) the Coriolis effect, (c) the change from corotating to inertial frame, and (d) total frequency shifts. Red lines and points represent splittings for the m = 0 component. The blue and green lines represent the magnetic splitting for the m = 1 and m = −1 components respectively, that overlap on panel a. The grey dashed line represents the ζ function as defined by Goupil et al. (2013), the beige vertical line indicates the position of ν_{max}, and the grey vertical line the frequency of the ℓ = 1 gdominated mode closest to ν_{max}, ν_{1, g, max}, with squares indicating the corresponding splitting values associated with the different orders m. 

In the text 
Fig. 6. As in Fig. 5 but for ℓ = 2 mixed modes. The red line represents the m = 0 mode, the green line with blue points the m = −1, 1 components, and the green line with purple points the m = −2, 2 components. 

In the text 
Fig. 7. Sketch representing the effect of rotation and magnetism on a single unperturbed mixed mode frequency of degree ℓ = 1 and ℓ = 2. 

In the text 
Fig. 8. Normalised magnetic splittings versus mixedmode frequencies computed for a simulated M_{⋆} = 1.5 M_{⊙}, Z = 0.02, ν_{max} ≃ 172 μHz red giant. The normalised power law following 1/ν^{3} is superimposed in black. 

In the text 
Fig. 9. Left: mixedmode frequencies at different evolutionary stages along the RGB indicated in the right panel. Purple dots represent the individual frequencies of mixed modes, with the shade of the dots allowing to visualise the evolution of the mixedmode order n_{pg} (the order is higher at high frequencies). Squares represent for each evolutionary stage the frequency of the ℓ = 1 gdominated mode closest to ν_{max} (see Fig. 10 for the determination of this frequency). Middle: m = 1 magnetic splitting of ℓ = 1 mixed modes located inside [ν_{max} − 7Δν : ν_{max} + 7Δν] for each of the evolutionary stages denoted in the right panel. Squares represent for each evolutionary stage the frequency of the ℓ = 1 gdominated mode closest to ν_{max} (ν_{1, g, max}, see Fig. 5) Right: HertzsprungRussell diagram of the M_{⋆} = 1.5 M_{⊙}, Z = 0.02 star ascending the RGB. Coloured stars indicate the different evolutionary stages represented in the other panels. 

In the text 
Fig. 10. Left: m = 1 mixedmode mean frequency perturbation due to magnetism only (δν_{mag} in μHz) as a function of mixedmode frequencies (ν_{1, g, max} in μHz) and the magnetic field amplitude in Gauss for a M_{⋆} = 1.5 M_{⊙}, Z = 0.02. Typical limiting frequencies are indicated: black dashed lines represent the frequency resolution in the PSD corresponding to 4 years of Kepler continuous observations (∼8 nHz), 3 years of PLATO two main fields (∼11 nHz), 1 year of TESS continuous viewing zone data (∼30 nHz), a typical line width of radial modes, which are an upper limit for the line width of dipolar mixed modes (∼0.12 μHz), 2 months of PLATO stareandstep observations (0.19 μHz), and 27 days of TESS individual sectors (0.37 μHz). The coloured area is limited by the critical magneticfield amplitude range for the validity of the perturbative analysis. We choose to limit the results at B_{c, g}/1000, as defined in Sect. 3.5. Red lines represent the expected magneticfield amplitude by considering magneticflux conservation from the end of the MS, with original fossil field amplitudes of 0.1 and 1 MG, as detailed in Sect. 2.3. Right: same as left panel, but simplified with the colour map replaced by typical limiting frequency values corresponding to Kepler data. The yellow area corresponds to the combination of ν_{1, g, max} and magnetic amplitude for which magnetic splittings cannot be visible in Kepler observations. The blue area corresponds to small magnetic shifts that should be visible in Kepler observations. The purple area corresponds to very large magnetic shifts, larger than the typical line width of radial modes. 

In the text 
Fig. 11. Synthetic stretched period echelle diagram built from Eq. (54), showing modelled mixedmode frequencies affected by rotation with Ω_{core} = 0.5 μHz and magneticfield strengths increasing left to right with B_{0} = 0.1 MG, B_{0} = 0.2 MG, B_{0} = 0.5 MG, and B_{0} = 1 MG in the radiative interior. Purely rotationally split components are given for reference by background crosses on the first two panels. Coloured dots indicate the positions of the ℓ = 1 mixed modes of azimuthal order m ∈ {1, 0, −1}, respectively blue, red, and green dots. The horizontal black line delimits domains for which the rotational splitting is small compared to (see Gehan et al. 2018, for more details). Synthetic patterns above this line are typically representative of patterns of mixed modes for subgiants and early giants, and those below are associated with evolved red giants. The magnetic effect is small compared to rotational effects in white areas, is of the same order of magnitude as in orange areas, and is larger than δν_{rot, core, m} in red areas. 

In the text 
Fig. 12. Left: a range of mixed mode frequencies is shown for different evolutionary stages along the subgiant phase (and early RGB). The individual stages are indicated on the evolutionary track plotted in the right panel. Purple dots represent the individual frequencies of mixed modes, with the shade of the dots allowing to visualise the evolution of the mixedmode order n_{pg}. Squares represent the frequency of the ℓ = 1 gdominated mode closest to ν_{max} for each evolutionary stage. We refer to Fig. 5 for the determination of this frequency. Middle: magnetic splitting of ℓ = 1, m = 1 mixed modes located inside the range [ν_{max} − 7Δν : ν_{max} + 7Δν] for each of the evolutionary stages denoted in the right panel. Squares represent the frequency of the ℓ = 1 gdominated mode closest to ν_{max} (see Fig. 5 for further reference). Right: HertzsprungRussell diagram of a M_{⋆} = 1.5 M_{⊙}, Z = 0.02 star evolving on the subgiant phase. Coloured stars indicate the different evolutionary stages represented in the other panels. 

In the text 
Fig. 13. Normalised magnetic splittings of the p − m modes versus mixedmode frequencies computed for a simulated M_{⋆} = 1.5 M_{⊙}, Z = 0.02, ν_{max} ≃ 750 μHz subgiant. The normalised power law following ν is superimposed in black. 

In the text 
Fig. 14. Top panels: same as Fig. 5 but for younger stars (ν_{max} ≳ 350 μHz) with the abscissa changed to the age of the star (reversed to be consistent with the convention used in Fig. 5). Bottom panels: same as top panels but considering p − m modes instead of g − m modes. 

In the text 
Fig. 15. Characteristic time of angular momentum transport τ in years as expressed by Eq. (74), depending on the given evolutionary stage and the magnetic field amplitude. The magnetic field amplitude expected at each stage from the conservation of flux from the end of the MS is indicated by red lines for B_{MS} = 0.1 and 1 MG. Left: transport along the RGB with ν_{1, g, max} as the abscissa. Right: transport along the subgiant branch with age as the abscissa. 

In the text 
Fig. 16. Left: magnetic frequency splitting calculated at the central mode frequency ν_{1, g, max} versus the age of the different evolutionary stages reported in the right panel for Z = 0.02, M_{⋆} ∈ {1, 1.5, 2, 3} M_{⊙}. Middle: magnetic frequency splitting calculated at the central mode frequency ν_{1, g, max} for the different evolutionary stages and masses reported in the right panel. Right: HertzsprungRussell diagram corresponding to the colourcoded Z = 0.02, M_{⋆} ∈ {1, 1.5, 2, 3} M_{⊙} stars. The star symbols indicate the position of the star at which measurements are reported in the left and middle panels. 

In the text 
Fig. 17. Same as Fig. 16 but with varying metallicity; we represent M_{⋆} = 1.5 M_{⊙}, Z ∈ {0.04, 0.02, 0.002, 0.0002} stars. 

In the text 
Fig. B.1. Attempt to find the first zero of at the radiative–convective boundary. We consider the typical red giant with M = 1.5 M_{⊙}, Z = 0.02 on the RGB. Panel a: trend of the evolution of the ψ function at r = R_{rad} with the eigenvalue λ. The value λ_{1} closest to the eigenvalue is indicated by the red vertical line. Panel b: shape of the functions , and ψ_{λ1}(r, θ) with λ_{1} ≃ 7.62, the vertical dashed line indicates the surface. Panel c: magnetic field components with λ_{1} ≃ 7.62. Panel d: resulting magnetic field topology with λ_{1} ≃ 7.62. 

In the text 
Fig. B.2. Same legend as Fig. B.1 for the first zero of at the radiative–convective boundary instead. In this case we find the first eigenvalue λ_{1} ≃ 2.80. 

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.