Issue 
A&A
Volume 647, March 2021



Article Number  A122  
Number of page(s)  11  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/202039180  
Published online  18 March 2021 
Probing the internal magnetism of stars using asymptotic magnetoasteroseismology
^{1}
AIM, CEA, CNRS, Université ParisSaclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 GifsurYvette Cedex, France
email: stephane.mathis@cea.fr, lisa.bugnet@cea.fr
^{2}
Instituto de Astrofísica de Canarias, 38200 La Laguna, Tenerife, Spain
^{3}
Universidad de La Laguna, Dpto. de Astrofísica, 38205 La Laguna, Tenerife, Spain
Received:
13
August
2020
Accepted:
23
November
2020
Context. Our knowledge of the dynamics of stars has undergone a revolution through the simultaneous large amount of highquality photometric observations collected by spacebased asteroseismology and groundbased highprecision spectropolarimetry. They allowed us to probe the internal rotation of stars and their surface magnetism in the whole HertzsprungRussell diagram. However, new methods should still be developed to probe the deep magnetic fields in these stars.
Aims. Our goal is to provide seismic diagnoses that allow us to probe the internal magnetism of stars.
Methods. We focused on asymptotic lowfrequency gravity modes and highfrequency acoustic modes. Using a firstorder perturbative theory, we derived magnetic splittings of their frequencies as explicit functions of stellar parameters.
Results. As in the case of rotation, we show that asymptotic gravity and acoustic modes can allow us to probe the different components of the magnetic field in the cavities in which they propagate. This again demonstrates the high potential of using mixedmodes when this is possible.
Key words: asteroseismology / waves / stars: magnetic field / stars: oscillations / methods: analytical
© S. Mathis 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
Along with rotation, magnetism is one of the two corner stones of stellar dynamics. Magnetism deeply modifies the interactions of every type of stars with their surrounding environment and their evolution, in particular through the transport of angular momentum and chemicals that they potentially trigger in their interiors (e.g. Brun & Browning 2017; Braithwaite & Spruit 2017). In this framework, two revolutions have occurred in stellar physics in the past two decades.
The first revolution is spacebased helio and asteroseismology (e.g. ChristensenDalsgaard 2015; García & Ballot 2019; Aerts 2021). These techniques have allowed us to probe the structure, the internal rotation, and the magnetic activity of the Sun and of stars with high precision. Helioseismology analyses with spacedbased instruments (i.e. the Michelson Doppler Imager (MDI; Scherrer et al. 1995) and Global Oscillations at Low Frequencies (GOLF; Gabriel et al. 1995) instruments onboard the SOlar and Heliospheric Observatory (SOHO; Domingo et al. 1995) spacecraft) revealed the solar rotation profile down to 0.25 R_{⊙} (R_{⊙} is the radius of the Sun) with an outer radiative core that rotates as a solid body. Strong efforts are still ongoing to constrain the rotation of the nuclear core (e.g. García et al. 2007; Mathur et al. 2008). Spacebased asteroseismology with the Kepler space mission then allowed us to probe the rotation of stars from their surface to their core in the case of mainsequence and evolved low and intermediatemass stars (Aerts et al. 2019, and references therein). First, acoustic modes have been used to study the internal rotation of mainsequence lowmass stars (e.g. Benomar et al. 2015). Next, gravitoinertial modes allowed us to determine the rotation near the interface of the radiative envelope with the convective core and in the convective core of fastrotating intermediatemass stars (e.g. Van Reeth et al. 2016; Ouazzani et al. 2020, respectively). Finally, mixed gravitoacoustic modes in evolved low and intermediatemass stars provide precise constraints on the rotation of the internal radiative zone down to the core and the averaged rotation of the envelope (e.g. Mosser et al. 2012; Gehan et al. 2018; Deheuvels et al. 2014, 2020). These observations have revealed coretosurface rotation ratios that are weaker by several orders of magnitude than those predicted by the most modern rotating stellar models for the whole HertzsprungRussel diagram (e.g. Eggenberger et al. 2012; Ceillier et al. 2013; Marques et al. 2013; Cantiello et al. 2014; Ouazzani et al. 2019). This demonstrates that a powerful mechanism acts throughout the evolution of stars, including our Sun, and extracts angular momentum from the deep interior. The transport induced by stable (e.g. Mestel & Weiss 1987; Mathis & Zahn 2005) or unstable magnetic fields in stably stratified stellar radiation zones is one of the potential mechanisms (e.g. Spruit 2002; Fuller et al. 2019). These discoveries have triggered the development of a strong theoretical basis from which we interpret seismic observations. For oscillation modes, with frequencies far higher than the rotation frequency, the perturbative theory allows us to compute frequency splittings due to the Coriolis and centrifugal accelerations (e.g. Aerts et al. 2010; Goupil et al. 2013, and references therein). It also allows us to compute splittings due to the change of reference frame from the stellar corotating one to the observer frame. These splittings are then used to constrain the internal stellar rotation. This method has been intensively used to probe the solar rotation profile (e.g. Thompson et al. 1996; Couvidat et al. 2003) and the coretosurface rotation ratio in evolved stars (e.g. Beck et al. 2012; Deheuvels et al. 2012, 2014, 2015, 2020; Spada et al. 2016). For lowfrequency gravitoinertial modes with frequencies of the same order of magnitude as the rotation frequency, the traditional approximation of rotation (TAR), which is not perturbative, can be used (e.g. Lee & Saio 1997; Townsend 2003; Van Reeth et al. 2018; Mathis & Prat 2019). The traditional approximation of rotation can be used when the buoyancyrestoring force dominates the Coriolis force in the direction of the stable entropy and chemical stratification (Mathis 2009). When the variation in the period spacing (i.e. the difference between the periods of two consecutive modes) is considered as a function of the mode period, the (differential) rotation of the oscillations propagation cavity can be constrained (Bouabid et al. 2013; Van Reeth et al. 2018). This method is currently used to provide key information on the rotation rate of the radiative layers surrounding the convective core in fastrotating intermediatemass stars (e.g. Van Reeth et al. 2016; Christophe et al. 2018; Aerts et al. 2019; Li et al. 2019, 2020). For each case, the study of the effect of rotation on asymptotic gravity (and gravitoinertial), acoustic, and mixed gravitoacoustic modes, which rapidly oscillate in the radial direction, has provided a deep understanding of their modification by the Coriolis (and centrifugal) acceleration and key tools for seismic modelling and inversions.
The second revolution is groundbased highprecision spectropolarimetry that probes the strength and geometric configuration of magnetic fields at the surface of stars in the whole HertzsprungRussell diagram from the premain sequence to the late stages of their evolution (e.g. Donati & Landstreet 2009; Aurière et al. 2015; Wade et al. 2016). Large spectropolarimetric surveys thus led our knowledge of dynamogenerated and fossil magnetic fields to a new level of understanding, in particular through the simultaneous development of theoretical models and threedimensional (3D), global, nonlinear magnetohydrodynamical (MHD) simulations (e.g. Brun & Browning 2017; Braithwaite & Spruit 2017). However, spectropolarimetry does not give access to the distribution of the magnetic field inside stars. Once again, asteroseismology is the answer. The method is the same as for rotation, that is, specific signatures of the effect of magnetism on the excitation, the propagation, and the damping of stellar oscillations are looked for. First studies have considered the effects of simple magnetic configurations such as purely dipolar poloidal fields (aligned or inclined to the rotation axis) or purely toroidal fields (e.g. Goode & Thompson 1992; Takata & Shibahashi 1994). However, our knowledge of magnetic topologies and their stability has strongly progressed in the past two decades through the simultaneous development of stellar spectropolarimetry and of 3D MHD simulations. For instance, we now have a better understanding of the formation and the topologies of fossil fields in stellar radiation zones (Braithwaite & Spruit 2004). They should be a combination of dipolar poloidal and toroidal fields to be in their lowest energy stable equilibrium state (Duez & Mathis 2010). Therefore it becomes mandatory to study the seismic signatures of such complex topologies. This has recently been undertaken for the fossil fields in the radiative envelope of upper mainsequence stars (Prat et al. 2019) and the radiative core of lowmass and intermediatemass evolved stars (Bugnet et al. 2021) as well as for the dynamo fields generated in the convective envelope of solartype stars (Kiefer & Roth 2018).
In this framework, asymptotic analyses have been very useful to provide powerful seismic diagnoses for probing the rotation of stellar interiors (Goupil et al. 2013). In this article, we show how a method like this, which has been developed for rotation, can be applied to magnetism. In Sect. 2 we first recall the asymptotic theory for acoustic, gravity, and mixed gravitoacoustic modes. In Sect. 3 we identify the dominant terms in the magnetic splittings for each class of modes, and we derive asymptotic expressions that are the counterpart of those known for stellar rotation. In Sect. 4 we present a proofofconcept application to a typical intermediatemass red giant star. Finally, in Sect. 5, we make a synthesis of the results and discuss how asymptotic predictions can allow us to probe the distribution of the magnetic field at different depths in stars, and the applications of these predictions to different stellar types that can be foreseen.
2. Asymptotic theory for (mixed) stellar oscillation modes
We first recall the basics of the theory of the oscillations of nonrotating and nonmagnetic stars, which is necessary to compute firstorder magnetic splittings in the asymptotic limits of lowfrequency gravity (g) modes and highfrequency pressure (p) modes. The Lagrangian displacement of an oscillation eigenmode is expanded as
where ξ_{r} and ξ_{h} are the radial functions of its vertical and horizontal components, ω is the angular frequency, are the spherical harmonics, l is its degree, and m the azimuthal order, ∇_{h} = ∂_{θ}(⋅)e_{θ} + 1/sin θ ∂_{φ}(⋅)e_{φ} is the horizontal gradient, {e_{r},e_{θ},e_{φ}} are the spherical unit vectors, and Re is the real part of a complex number. Following Gough (1993) and Hekker & ChristensenDalsgaard (2017), who used the adiabatic linearised hydrodynamical equations assuming the Cowling approximation, we derive the Schrödingerlike wave equation
where X = c^{2}ρ^{1/2}div ξ. The sound speed c is defined as , where ρ and P are the density and the pressure of the hydrostatic background, respectively, and Γ_{1} = (∂lnP/∂lnρ)_{S} is the first adiabatic exponent, with S being the macroscopic entropy. The vertical wave number k_{r} is given by
We have identified the Lamb frequency S_{l} defined by
is the horizontal wave number, and the BruntVäisälä frequency N defined by
where g is the gravity of the hydrostatic background. Finally, ω_{c} is the acoustic cutoff frequency, defined by
where H = −(dlnρ/dr)^{−1} is the density scale height.
2.1. Asymptotic g modes
In the lowfrequency regime for which ω < < {S_{l},N}, the Lagrangian displacement becomes mostly horizontal (i.e. ξ_{r}< < ξ_{h}) and the vertical wave number reduces to (e.g. Hekker & ChristensenDalsgaard 2017)
Using the reduced canonical variable
which has been introduced by Shibahashi (1979), and applying the JWKB (for Jeffreys, Wentzel, Kramers, and Brillouin) method (we refer to Appendix C and to Erdelyi 1956; Fröman & Fröman 2005), we derive the asymptotic expression for ξ_{h}(r)
where A_{W} is the amplitude fixed by the excitation and damping mechanisms and the inertia of the mode (Samadi et al. 2015), is a phase, and r_{t; i} (r_{t; e}) is the internal (external) turning point for which k_{r} vanishes.
2.2. Asymptotic p modes
In the highfrequency regime for which ω > > S_{l}, the Lagrangian displacement becomes mostly vertical (i.e. ξ_{h} < < ξ_{r}) and the vertical wave number reduces to
Using the reduced canonical variable
which has been introduced by Shibahashi (1979), and applying again the JWKB method (Appendix C), we derive the asymptotic expression for ξ_{r}(r)
where A_{V} is the mode amplitude, is a phase, and R is the radius of the star.
3. Asymptotic magnetic splittings
3.1. Studied magnetic configuration
For a sufficiently moderate magnetic field, we can perturb the wave equation to the firstorder by the linearised Lorentz force
where B is the stellar magnetic field, δB = ∇ ∧ (ξ∧B) is its fluctuation, and μ_{0} = 4π (in c.g.s.) is the magnetic permeability of vacuum. We obtain the general expression for the magnetic frequency splitting
where (where the asterisk is the complex conjugate) and ω_{0} is the angular frequency of the unperturbed mode; we refer to Prat et al. (2019) for details on the perturbation formalism. A supplementary term, ∇ ⋅ (ρξ)(∇∧B) ∧ B/(μ_{0}ρ), should be taken into account in the general case (Goode & Thompson 1992). However, we can demonstrate again using the JWKB approximation that it is not the dominant term for the asymptotic modes. It scales as k_{r}, while the dominant terms scales as while k_{r} H > > 1 in the JWKB limit. Therefore, we assume that the characteristic length scale of the wave is shorter than the length scales of the variation in the hydrostatic background (H) and in the magnetic field configuration.
The mode inertia is expressed as
In this work, we chose to focus on a dipolar mixed (i.e. with both a poloidal and a toroidal components) axisymmetric configuration
where B_{0} is the field amplitude and the functions b_{i}(r) (with i = {r,θ,φ}) describe the radial dependence of each component. On the one hand, it can model a stable fossil field that has been formed in a stably stratified radiation zone (e.g. Duez & Mathis 2010). On the other hand, it can model the first largescale component of a dynamogenerated magnetic field (e.g. Brun & Browning 2017).
3.2. Asymptotic g modes
We first focus on lowfrequency asymptotic g modes, which probe radiative zones, that is, cores of lowmass stars. For these modes, we have ω_{0} < < {S_{l},N} and thus k_{r} ≈ (N/ω_{0})k_{h} > > k_{h} with . Because of their weak compressibility, k ⋅ ξ ≈ 0, and they are thus quasihorizontal with ξ_{v} < < ξ_{h}. Therefore we selected the terms that involve the products to isolate the dominant terms in the general expression of the magnetic splitting provided by Eqs. (B.1) and (B.2). Then, using the results provided in Eq. (9), we have ξ_{h} ∝ (k_{r})^{−1/2}exp[i∫k_{r}dr′] if the studied lowfrequency mode varies rapidly in space in comparison to the characteristic length scale of the variation in the stellar hydrostatic structure and in the magnetic configuration. We thus have and . Because k_{r}H > > 1 (because N > > ω, and where H is the characteristic length of the variation in the hydrostatic background) in the asymptotic regime, the dominant terms are thus those , that is, . As in Hasan et al. (2005) and Rashba et al. (2007), we can therefore identify the dominant terms in the complete expression of the magnetic splitting given in Eq. (B.1). After a last integration by parts, we obtain
where
Because the angular coupling coefficient C_{l, m} is independent of the sign of m in the case of the considered dipolar magnetic configuration, we can already predict that the magnetic frequency splittings for pro and retrograde modes (m > 0 and m < 0, respectively) will be the same. This will allow us to potentially distinguish them from the rotationinduced frequency splittings which depend on the sign of m with opposite values when computed to the firstorder in Ω (Ω is the rotation of the star); this is discussed in details in Sect. 4. In the case of mainsequence lowmass stars and of evolved stars, we can assume that r_{t; i} ≈ 0 and r_{t; e} ≈ R_{core} for lowfrequency g modes, where R_{core} is the radius of the radiative core. In the case of mainsequence intermediatemass stars, we have r_{t; i} ≈ R_{CC} and r_{t; e} ≈ R, where R_{CC} is the radius of the convective core.
In addition, the mostly horizontal lowfrequency g modes allows us to probe the orthogonal radial component of the field. This can be understood from the form of the linearised Lorentz force, which implies terms of the form ξ ∧ B (cf. Eq. (13)). These terms thus couple the Lagrangian displacement of the mostly horizontal mode with the orthogonal radial component of the field.
Using Eq. (9), this can be written as
where and .
Using the properties of rapidly oscillating integrals (see Appendix C), this reduces to
where we have introduced the dimensionless radius x = r/R (with x_{t; i} = r_{t; i}/R and x_{t; e} = r_{t; e}/R), the maximum of the BruntVäisälä frequency N_{max} and its dimensionless radial profile defined such that . The form of this integral is really interesting because we recover the one known for asymptotic rotational splittings, but for the radial component of the magnetic field. We can then define a kernel as Goupil et al. (2013) (see the end of Sect. 3.4.), which opens the path for potential inversions when magnetic seismic signatures can be detected (we refer to Bugnet et al. 2021 for the corresponding values of the critical amplitudes).
This result can be generalised for g modes with frequencies of the same order of magnitude as the inertial frequency 2Ω, where Ω is the angular velocity of the star. This is the case for instance in fastrotating stars such as γ Doradus stars (e.g. Van Beeck et al. 2020) and SPB stars (e.g. Prat et al. 2019). When the buoyancy force exceeds the Coriolis force in the direction of the stable entropy or chemical stratification, we can assume the traditional approximation of rotation (TAR; e.g. Eckart 1960; Lee & Saio 1997; Townsend 2003; Bouabid et al. 2013; Mathis 2009). In this approximation, the horizontal projection of the rotation vector is neglected that allows us to separate variables when the wave propagation equation is solved as in the nonrotating case. Then, we obtain for g modes modified by rotation, in other words, for gravitoinertial modes (hereafter gi modes),
with
where ν = 2Ω/ω_{0} is the spin parameter and Λ_{k, m}(ν) is the horizontal eigenvalues of the Hough functions (e.g. Hough 1898; Lee & Saio 1997; Townsend 2003) {H_{r}(cos θ),H_{θ}(cos θ),H_{φ}(cos θ)} that generalise the spherical harmonics when the rotation within the TAR is taken into account (we refer to Appendix A for the details of their definition).
3.3. Asymptotic p modes
We can apply the same method to the case of highfrequency acoustic modes, which probe the external layers of stars. In the same way, we isolate using Eq. (12) the dominant terms in Eq. (B.2), which are
where is computed in Appendix D.3. As for the angular coupling coefficients C_{l, m} computed for asymptotic g modes, we also identify here that D_{l, −m} = D_{l, m}. We have introduced the internal turning point for which S_{l}(r_{t}) = ω_{0}. In the case of lowdegree highfrequency acoustic modes, we have r_{t} → 0. Using the asymptotic JWKB solution given in Eq. (12) and the theory of rapidly oscillating integrals (see Appendix C), we finally obtain
where . We again recover exactly the same form as for the asymptotic rotational splittings, but for the horizontal components of the magnetic field. As in the case of lowfrequency g modes, highfrequency acoustic modes, which have mostly vertical Lagrangian displacements, allow us to probe the orthogonal horizontal (latitudinal and azimuthal) components of the field.
3.4. Mixed modes
As in Goupil et al. (2013), we can express the magnetic splittings for mixed (gravitoacoustic) modes using the expressions of the splittings of asymptotic gravity and acoustic modes:
where we have introduced the
function, with I_{g} and I being the inertia of their gdominated component (the socalled g − m modes) and their total inertia, respectively. It quantifies their behaviour as a gdominated mode when ζ = 1 or as a pdominated mode (the socalled p − m modes) when ζ → 0.
It is interesting here to make the analogy with the case of rotation studied in Goupil et al. (2013). They showed that for the dipolar modes (l = 1)
where
with Ω_{c} = (GM/R^{3})^{1/2} the critical Keplerian velocity, and
Using Eqs. (20) and (24), we identify the analogy between Ω in the core and in the case of gdominated modes and between Ω in the envelope and for pdominated modes. This shows the potential power of mixed modes to probe the magnetism of stars where they propagate from their surfaces to their cores.
4. Key application: evolved lowmass and intermediatemass stars
Because of all the key results they have provided within the theory of stellar evolution and the study of the internal dynamics of stars based on spacebased asteroseismology (e.g. Hekker & ChristensenDalsgaard 2017; Aerts et al. 2019), we chose here to apply our theoretical results to a typical 1.5 M_{⊙} 2.8 Gyr old intermediatemass red giant star (with a solar metallicity Z = 0.02). The profiles of its BruntVaïsälä and Lamb frequencies for dipolar and quadrupolar modes (N, S_{l = 1}, S_{l = 2}, respectively) are represented in Fig. 1 (left panel) where the position of ν_{max} = 172.2 μHz, the frequency of maximum oscillation power, is reported by the thick horizontal orange line.
Fig. 1.
Left: profiles of the Lamb frequencies (S_{l}; red lines) for l = {1,2} and of the BruntVäisälä frequency (N; blue line) for the studied 1.5M_{⊙} 2.8 Gyr old intermediatemass red giant star (with a solar metallicity Z = 0.02); the position of ν_{max} is also reported (orange line). Right: topology of the fossil magnetic field in the stably stratified radiative core of the studied red giant (its limit is given by the thick black line). The thin black lines are the poloidal field lines, while the intensity of the toroidal field (normalised by the maximum values of the radial field) is given in colours. 
First, we computed the topology of a mixed (with both poloidal and toroidal components) relaxed fossil magnetic field in its radiative core following the method introduced in Prat et al. (2019) and the theoretical modelling for this field by Duez & Mathis (2010). As in Bugnet et al. (2021), the chosen boundary conditions are b_{r} = 0 at the boundary of the radiative core and convective envelope. Its amplitude was fixed to B_{0} = 10^{6} G and its geometry is reported in Fig. 1 (right panel). For this possible field amplitude (Cantiello et al. 2016), magnetic frequency splittings should be detectable within Kepler data (Bugnet et al. 2021). The thin black lines represent the poloidal field lines, and the colours give the amplitude of the toroidal component of the field (the boundary of the radiative core is given by the thick black line).
Then, we computed the magnetic splittings with their complete expression (δω_{mag}/ω_{0})_{num} given in Appendix B using the combination of the GYRE and MESA codes (Townsend & Teitler 2013; Paxton et al. 2011) for the dipolar (l = 1) and quadrupolar (l = 2) mixed modes (represented by the dots in the left panel (right panel) of Fig. 2). We computed the corresponding asymptotic prediction (δω_{mag}/ω_{0})_{asympt} using Eq. (20) (solid lines in Fig. 2). The complete and asymptotic expressions were compared based on the relative error defined as [(δω_{mag}/ω_{0})_{asympt}−(δω_{mag}/ω_{0})_{num}]/(δω_{mag}/ω_{0})_{num} (Fig. 2, bottom panels). In the lowfrequency regime, the agreement is excellent for gdominated mixed modes with a relative error lower than 4% for dipolar modes and 3% for quadrupolar modes outside the dips of the ζ function. These results are coherent with the high radial orders of the obtained mixed modes, that is, n ∈ [−138,−24] for the dipolar modes and n ∈ [−244,−56] for the quadrupolar modes, for which the JWKB approximation works well. This demonstrates the high potential of using the asymptotic magnetic splittings for intensive seismic modelling where they can be very useful to explore a broad space of stellar parameters (see e.g. the work by Van Beeck et al. 2020, for mainsequence intermediatemass stars).
Fig. 2.
Firstorder magnetic splitting for dipolar (l = 1, left top panel) and quadrupolar (l = 2, right top panel) mixed modes in the studied 1.5 M_{⊙} red giant star. The dots correspond to the full expression for the splitting (Eqs. (B.1) and (B.2)) computed numerically using the GYRE and MESA codes, and the solid lines are the asymptotic solutions given in Eq. (20). The ζ function related to the mixed character of modes propagating in evolved lowmass stars is plotted as a dotted grey line. Finally, the vertical thick orange line indicates the position of ω_{max} = 2πν_{max}. Bottom panels: relative errors between the asymptotic magnetic splittings as derived in Eq. (20) and their full expression given in Appendix B for dipolar (left panel) and quadrupolar (right panel) modes, respectively. 
As this has already been identified in Eq. (18), the effects of the studied dipolar mixed (poloidal + toroidal) fossil magnetic field on the frequencies of gdominated mixed modes are the same for pro and retrograde modes. This strongly differs from the case of the frequency splittings induced by the stellar rotation, which are due to the combination of the Coriolis acceleration and of the change in reference frame (e.g. Aerts et al. 2010),
where F_{c}(ξ_{0}) = 2iω_{0}Ω ∧ ξ_{0} is the Coriolis acceleration operator and F_{f}(ξ_{0}) = − 2mω_{0}Ωξ_{0} is the advection operator due to the change of reference frame. In the case of a radial differential rotation, we have
and
Therefore, δω_{rot} is proportional to m, leading to opposite rotational frequency splittings for pro and retrograde modes. This is illustrated in Fig. 3, where we represent in the middle column this rotational splitting when a typical twozone differential rotation with a solidbody rotation Ω_{core}/(2π) = 0.5 μHz in the radiative core and Ω_{envelope}/(2π) = (Ω_{core}/10)/2π = 0.05 μHz in the convective envelope is assumed (Gehan et al. 2018). When added to the magnetic frequency splitting, which is the same for pro and retrograde modes, this leads to an asymmetry of the total frequency splitting that scales with the squared amplitude of the field (cf. Eqs. (17) and (20)). This asymmetry will allow asteroseismologists to probe the potential presence of an (axisymmetric) stable fossil magnetic field in the core of evolved lowmass and intermediatemass stars. We refer to Bugnet et al. (2021) for a detailed characterisation of this asymmetry in seismic data. We did not consider the case of pdominated modes during the subgiant phase because Bugnet et al. (2021) demonstrated that magnetic signatures would not be detectable for them.
Fig. 3.
Leadingorder perturbations on the frequencies of dipolar (l = 1, top panels) and quadripolar (l = 2, bottom panels) mixed modes in the studied 1.5 M_{⊙} red giant star. Left panels: magnetic angular frequency splitting δω_{mag} (for B_{0} = 1MG) of the m ∈ [ − l, l] modes. Middle panels: angular frequency splitting δω_{rot} for a plausible twolayer differential rotation (with Ω_{core} = 0.5 μHz and Ω_{env} = Ω_{core}/10) of the same mixed modes is represented. Right panels: combination of magnetic and rotational splittings of mixedmode angular frequencies δω_{tot} = δω_{mag} + δω_{rot}. The vertical orange line indicates in each panel the expected frequency of maximum power of the modes ν_{max}. The ζ function is represented in each case by the grey dashed line. 
5. Conclusion and perspectives
We have established the asymptotic values of the magnetic splittings of lowfrequency gravity and gravitoinertial (gi) modes and highfrequency acoustic modes. These results were applied to the case of mixed gravitoacoustic modes such as those that propagate in evolved low and intermediatemass red giant stars. They were derived in the case of a largescale axisymmetric mixed (poloidal + toroidal) dipolar field that can be representative of a stable fossil field or of the largescale axisymmetric dipolar component of a dynamogenerated field. The main results that we obtained are listed below.
The expressions obtained for the magnetic splittings in the asymptotic limits of lowfrequency g (gi) modes and highfrequency p modes are very similar in their form to those of rotational splittings in the same limits. This opens the path to potential inversions of the internal distribution of magnetic fields in stellar interiors when these signatures can be detected in stellar oscillation frequency spectra.
For asymptotic g (gi) modes, the splittings scale as the normalised squared ratio of the local Alfvén frequency to the mode frequency integrated along the buoyancy radius. For asymptotic p modes, they scale as the normalised squared ratio of the Alfvén velocity to the sound speed integrated along the acoustic path.
In the case of the dipolar configuration we studied, magnetic frequency splittings are the same for pro and retrograde modes (m > 0 and m < 0, respectively). This allows us to distinguish them from rotational frequency splittings that have opposite values for pro and retrograde modes. This leads to an asymmetry of the total frequency splittings when the two effects are added. It could be detectable within highprecision asteroseismic data for sufficient field amplitude. We refer to Bugnet et al. (2021) for a complete seismic characterisation of this asymmetry and for the study of this required amplitude in the evolution of lowmass and intermediatemass stars.
In the case of the largescale mixed dipolar configuration, which is treated here, lowfrequency g (gi) modes and highfrequency p modes explore different components of the magnetic field because of the form of the linearised Lorentz force, which involves ξ ∧ B terms. On the one hand, mostly horizontal lowfrequency g (gi) modes, which propagate in the central regions of stars, probe the orthogonal radial component of the field. On the other hand, mostly vertical highfrequency p modes, which propagate in the external regions of stars, probe the orthogonal horizontal (latitudinal and azimuthal) components of the field. As in the case of the rotation, where mixed modes have allowed us to probe the rotation of the core and of the external envelope with their gdominated (g − m) and pdominated (p − m) components, respectively, this opens the path to probe the stellar magnetic field from the surface to the core of stars. Because high frequencies of asymptotic p modes will be more distant than low frequencies of asymptotic g modes from the Alfvén frequency, the seismic signatures of the external field can be more difficult to detect in the case of a lowamplitude field. This can be fixed in the case of bright stars such as those observed by the Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015), and in the future by the PLAnetary Transits and Oscillations of stars (PLATO) mission (Rauer et al. 2014), by combining spectropolarimetic observations to determine the configuration and the strength of the field at the surface and seismic signatures to probe the core magnetism. This strategy is summarised in Fig. 4.
Fig. 4.
Global strategy for probing the internal magnetism of stars using asteroseismology (in potential synergy with spectropolarimetry). Here, the scheme is designed for evolved low and intermediatemass stars. 
In a near future, the derived asymptotic relations should be tested and systematically used for all classes of pulsators. In addition, more complex (nondipolar and nonaxisymmetric) magnetic topologies must be studied (e.g. Prat et al. 2020). This should allow us to build a complete knowledge of the magnetism of stars from their surfaces to their cores, step by step.
Acknowledgments
The authors thank the referee and Pr. J. ChristensenDalsgaard for their very constructive comments and remarks that allowed us to improve the article. St. M., L. B., V. P., and K. A. acknowledge support from the European Research Council through ERC grant SPIRE 647383. All the members from CEA 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. We made great use of the megyr python package for interfacing MESA and GYRE codes.
References
 Abramowitz, M., & Stegun, I. A. 1972, Handbook of Mathematical Functions (New York: Dover Publications, INC.) [Google Scholar]
 Aerts, C. 2021, Rev. Mod. Phys., 93, 015001 [Google Scholar]
 Aerts, C., ChristensenDalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology (Springer) [Google Scholar]
 Aerts, C., Mathis, S., & Rogers, T. M. 2019, ARA&A, 57, 35 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Aurière, M., KonstantinovaAntova, R., Charbonnel, C., et al. 2015, A&A, 574, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Benomar, O., Takata, M., Shibahashi, H., Ceillier, T., & García, R. A. 2015, MNRAS, 452, 2654 [NASA ADS] [CrossRef] [Google Scholar]
 Bouabid, M. P., Dupret, M. A., Salmon, S., et al. 2013, MNRAS, 429, 2500 [NASA ADS] [CrossRef] [Google Scholar]
 Bugnet, L., Prat, V., Mathis, S., et al. 2021, A&A, submitted [arXiv:2102.01216] [Google Scholar]
 Braithwaite, J., & Spruit, H. C. 2004, Nature, 431, 819 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Braithwaite, J., & Spruit, H. C. 2017, R. Soc. Open Sci., 4, 160271 [Google Scholar]
 Brun, A. S., & Browning, M. K. 2017, Liv. Rev. Sol. Phys., 14, 4 [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]
 ChristensenDalsgaard, J. 2015, From the Sun to the distant stars, and back Future research directions (Cambridge University Press), 276 [Google Scholar]
 Christophe, S., Ballot, J., Ouazzani, R. M., Antoci, V., & Salmon, S. J. A. J. 2018, A&A, 618, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Couvidat, S., García, R. A., TurckChièze, S., et al. 2003, ApJ, 597, L77 [NASA ADS] [CrossRef] [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 [Google Scholar]
 Deheuvels, S., Ballot, J., Beck, P. G., et al. 2015, A&A, 580, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Ballot, J., Eggenberger, P., et al. 2020, A&A, 641, A117 [EDP Sciences] [Google Scholar]
 Domingo, V., Fleck, B., & Poland, A. I. 1995, Sol. Phys., 162, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Donati, J. F., & Landstreet, J. D. 2009, ARA&A, 47, 333 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Duez, V., & Mathis, S. 2010, A&A, 517, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eckart, C. 1960, Hydrodynamics of Oceans and Atmospheres (Pergamon Press) [Google Scholar]
 Eggenberger, P., Montalbán, J., & Miglio, A. 2012, A&A, 544, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Erdelyi, A. 1956, Asymptotic Expansions (Dover Publications, Inc.) [Google Scholar]
 Fröman, N., & Fröman, P. O. 2005, Physical Problems Solved by the PhaseIntegral Method (NorthHolland Publishing Company) [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]
 García, R. A., & Ballot, J. 2019, Liv. Rev. Sol. Phys., 16, 4 [Google Scholar]
 García, R. A., TurckChièze, S., JiménezReyes, S. J., et al. 2007, Science, 316, 1591 [Google Scholar]
 Gehan, C., Mosser, B., Michel, E., Samadi, R., & Kallinger, T. 2018, A&A, 616, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Goode, P. R., & Thompson, M. J. 1992, ApJ, 395, 307 [NASA ADS] [CrossRef] [Google Scholar]
 Gough, D. O. 1993, Astrophysical Fluid Dynamics  Les Houches 1987, 399 [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, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hekker, S., & ChristensenDalsgaard, J. 2017, A&ARv, 25, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hough, S. S. 1898, Philos. Trans. R. Soc. London Ser. A, 191, 139 [Google Scholar]
 Kiefer, R., & Roth, M. 2018, ApJ, 854, 74 [CrossRef] [Google Scholar]
 Lee, U., & Saio, H. 1997, ApJ, 491, 839 [NASA ADS] [CrossRef] [Google Scholar]
 Li, G., Bedding, T. R., Murphy, S. J., et al. 2019, MNRAS, 482, 1757 [NASA ADS] [CrossRef] [Google Scholar]
 Li, G., Van Reeth, T., Bedding, T. R., et al. 2020, MNRAS, 491, 3586 [NASA ADS] [CrossRef] [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. 2009, A&A, 506, 811 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., & Prat, V. 2019, A&A, 631, A26 [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., & Zahn, J. P. 2005, A&A, 440, 653 [NASA ADS] [CrossRef] [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]
 Mestel, L., & Weiss, N. O. 1987, MNRAS, 226, 123 [NASA ADS] [CrossRef] [Google Scholar]
 Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012, A&A, 548, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ouazzani, R. M., Marques, J. P., Goupil, M. J., et al. 2019, A&A, 626, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ouazzani, R. M., Lignières, F., Dupret, M. A., et al. 2020, A&A, 640, A49 [EDP Sciences] [Google Scholar]
 Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [Google Scholar]
 Prat, V., Mathis, S., Buysschaert, B., et al. 2019, A&A, 627, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Prat, V., Mathis, S., Neiner, C., et al. 2020, A&A, 636, A100 [CrossRef] [EDP Sciences] [Google Scholar]
 Rashba, T. I., Semikoz, V. B., TurckChièze, S., & Valle, J. W. F. 2007, MNRAS, 377, 453 [NASA ADS] [CrossRef] [Google Scholar]
 Rauer, H., Catala, C., Aerts, C., et al. 2014, Exp. Astron., 38, 249 [Google Scholar]
 Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telescopes Instrum. Syst., 1, 014003 [Google Scholar]
 Samadi, R., Belkacem, K., & Sonoi, T. 2015, EAS Publ. Ser., 73–74, 111 [CrossRef] [Google Scholar]
 Scherrer, P. H., Bogart, R. S., Bush, R. I., et al. 1995, Sol. Phys., 162, 129 [Google Scholar]
 Shibahashi, H. 1979, PASJ, 31, 87 [NASA ADS] [Google Scholar]
 Spada, F., Gellert, M., Arlt, R., & Deheuvels, S. 2016, A&A, 589, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Spruit, H. C. 2002, A&A, 381, 923 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Takata, M., & Shibahashi, H. 1994, PASJ, 46, 301 [Google Scholar]
 Thompson, M. J., Toomre, J., Anderson, E. R., et al. 1996, Science, 272, 1300 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Townsend, R. H. D. 2003, MNRAS, 340, 1020 [Google Scholar]
 Townsend, R. H. D., & Teitler, S. A. 2013, MNRAS, 435, 3406 [NASA ADS] [CrossRef] [Google Scholar]
 Van Beeck, J., Prat, V., Van Reeth, T., et al. 2020, A&A, 638, A149 [EDP Sciences] [Google Scholar]
 Van Reeth, T., Tkachenko, A., & Aerts, C. 2016, A&A, 593, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Van Reeth, T., Mombarg, J. S. G., Mathis, S., et al. 2018, A&A, 618, A24 [NASA ADS] [EDP Sciences] [Google Scholar]
 Varshalovich, D. A., Moskalev, A. N., & Khersonskii, V. 1988, Quantum Theory of Angular Momentum (World Scientific) [CrossRef] [Google Scholar]
 Wade, G. A., Neiner, C., Alecian, E., et al. 2016, MNRAS, 456, 2 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Hough functions
The radial Hough functions (Hough 1898) are defined by H_{r}(cos θ) = f(μ = cos θ), where f is the solution of the socalled tidal Laplace equation,
where ν = 2Ω/ω is the spin factor. In the nonrotating case (ν = 0), the eigenvalue Λ (they are generally indexed by an integer k for a given m; see e.g. Lee & Saio 1997) reduces to l(l + 1), where l is the angular degree of the mode, and H_{r} simplifies to the classical associated Legendre polynomial .
The latitudinal and azimuthal Hough functions are respectively derived,
and
where ′ is the total latitudinal derivative. H_{φ} has the same parity as H_{r} with respect to θ = π/2, while H_{θ} has the opposite parity. In the nonrotating case (ν = 0), H_{θ} and H_{φ} reduce to and , respectively.
Appendix B: General expression of magnetic splittings
We recall all the nonzeroaverage terms of the work of the Lorentz force μ_{0}δF_{L} ⋅ ξ^{*}. All these terms involve either only poloidal components of the magnetic field or only the toroidal component, and they are grouped accordingly as in Prat et al. (2019). All terms are composed of a purely radial part multiplied by a purely latitudinal part, and the prime symbol (′) is a total radial or latitudinal derivative, depending on the considered part.
B.1. Poloidal terms
We define A = [(rb_{θ})′ + b_{r}]. The terms that involve poloidal components are obtained,
B.2. Toroidal terms
The terms involving the toroidal component are derived,
Appendix C: JWKB method for ordinary differential equations and rapidly oscillating integrals
On the one hand, studying the dynamics of asymptotic lowfrequency g and highfrequency p modes implies solving Schrödingerlike equations of the form
where λ is a large parameter and is a normalised vertical wave number. Applying the JWKB (for JeffreysWentzelKrammerBrillouin) method leads for λ > > 1 to the following solution (e.g. Erdelyi 1956; Fröman & Fröman 2005):
On the other hand, we have to compute integrals of the form
to calculate their magnetic splittings. In the cases where λ > > 1, they can be approximated by
where the x_{i} are the points such that f′(x_{i}) = 0 (Erdelyi 1956). If , the x_{i} are the turning points such that .
Appendix D: Algebra for computing horizontal integrals
D.1. Spherical harmonics
D.1.1. Definitions
The spherical harmonics are defined by (e.g. Varshalovich et al. 1988)
where is the associated Legendre function, and the normalisation coefficient
They obey the orthogonality relation,
where dΩ = sin θ dθ dφ and where the complex conjugate spherical harmonic is given by
D.1.2. Linear recurrence relations
We introduce useful linear recurrence relations that should be used to analytically compute the C_{l, m} and D_{l, m} coefficients. Using again Varshalovich et al. (1988), we have on the one hand
where
and
On the other hand, we obtain
where
and
D.1.3. Products of spherical functions and specific integrals
Integrals involving products of spherical harmonics should also be computed to calculate C_{l, m} and D_{l, m}. Using the normalisation and the orthogonality of spherical harmonics (Eqs. (D.2)–(D.3)) and their complex conjugate (Eq. (D.4)), we can write
where we define the integral like in Varshalovich et al. (1988),
with the 3jWigner coefficients (⋅⋅⋅) that are related to the classical ClebschGordan coefficients by
Then, using the initial definition of spherical harmonics (cf. Eqs. (D.1)–(D.2)), we deduce the expansion for the product of two associated Legendre functions,
where
The last relevant integral is
where Γ are the usual gamma functions (Abramowitz & Stegun 1972).
D.2. Application to C_{l, m} coefficients
It is possible to provide a fully analytical expression of C_{l, m} by using Eqs. (D.8), (D.9), (D.16), (D.17). However, it is very cumbersome, and the only interesting tractable expression is obtained in the case when l = m, where
The values of C_{l, m} for the dipolar (l = 1) and the quadrupolar (l = 2) modes are
D.3. Application to D_{l, m} coefficients
We have to compute
Using Eq. (D.4) and the identity
where is obtained with Eq. (D.2), we obtain
where ℐ has been defined in Eq. (D.13).
The values of D_{l, m} for the dipolar (l = 1) and the quadrupolar (l = 2) modes are
All Figures
Fig. 1.
Left: profiles of the Lamb frequencies (S_{l}; red lines) for l = {1,2} and of the BruntVäisälä frequency (N; blue line) for the studied 1.5M_{⊙} 2.8 Gyr old intermediatemass red giant star (with a solar metallicity Z = 0.02); the position of ν_{max} is also reported (orange line). Right: topology of the fossil magnetic field in the stably stratified radiative core of the studied red giant (its limit is given by the thick black line). The thin black lines are the poloidal field lines, while the intensity of the toroidal field (normalised by the maximum values of the radial field) is given in colours. 

In the text 
Fig. 2.
Firstorder magnetic splitting for dipolar (l = 1, left top panel) and quadrupolar (l = 2, right top panel) mixed modes in the studied 1.5 M_{⊙} red giant star. The dots correspond to the full expression for the splitting (Eqs. (B.1) and (B.2)) computed numerically using the GYRE and MESA codes, and the solid lines are the asymptotic solutions given in Eq. (20). The ζ function related to the mixed character of modes propagating in evolved lowmass stars is plotted as a dotted grey line. Finally, the vertical thick orange line indicates the position of ω_{max} = 2πν_{max}. Bottom panels: relative errors between the asymptotic magnetic splittings as derived in Eq. (20) and their full expression given in Appendix B for dipolar (left panel) and quadrupolar (right panel) modes, respectively. 

In the text 
Fig. 3.
Leadingorder perturbations on the frequencies of dipolar (l = 1, top panels) and quadripolar (l = 2, bottom panels) mixed modes in the studied 1.5 M_{⊙} red giant star. Left panels: magnetic angular frequency splitting δω_{mag} (for B_{0} = 1MG) of the m ∈ [ − l, l] modes. Middle panels: angular frequency splitting δω_{rot} for a plausible twolayer differential rotation (with Ω_{core} = 0.5 μHz and Ω_{env} = Ω_{core}/10) of the same mixed modes is represented. Right panels: combination of magnetic and rotational splittings of mixedmode angular frequencies δω_{tot} = δω_{mag} + δω_{rot}. The vertical orange line indicates in each panel the expected frequency of maximum power of the modes ν_{max}. The ζ function is represented in each case by the grey dashed line. 

In the text 
Fig. 4.
Global strategy for probing the internal magnetism of stars using asteroseismology (in potential synergy with spectropolarimetry). Here, the scheme is designed for evolved low and intermediatemass stars. 

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.