Perturbative analysis of the effect of a magnetic field on gravito-inertial modes

,


Introduction
The transport of angular momentum and chemical elements in stellar radiative zones is a key issue in modern stellar evolution theory.Purely hydrodynamical models fail to account for all the observational constraints; in particular, they are not efficient enough to reproduce the core rotation rates of red giant or intermediate-mass stars measured by seismology (Marques et al. 2013;Ceillier et al. 2013;Ouazzani et al. 2019).Models involving magnetic fields are a promising way of solving this angular momentum transport problem (Gough & McIntyre 1998;Fuller et al. 2019;Gouhier et al. 2022).However, observational constraints on the strength and topology of radiative zone magnetic fields are rare, limiting our ability to test these models.
Spectropolarimetry enables magnetic field measurements at the surface of stars with a radiative envelope.It was found that 5-10% of intermediate-mass and massive main-sequence stars have detectable magnetic fields.The vast majority of them are strong (>100 Gauss) and stable in time, and have a simple large scale topology (Aurière et al. 2007;Donati & Landstreet 2009;Wade et al. 2016), whereas very weak (∼1 Gauss) fields have been detected on a few very bright intermediate-mass stars (Lignières et al. 2014;Blazère et al. 2016).Asteroseismology has recently provided constraints on internal magnetic fields.First, the absence of the dipolar mixed gravito-acoustic modes in the oscillation spectrum of a fraction (∼20%) of red giant stars has been interpreted as due to strong magnetic fields in their radiative core (Fuller et al. 2015;Stello et al. 2016).Accordingly, if magnetic fields exceed a certain limit, gravity waves transform into short-wavelength Alfvèn waves that soon get dis-sipated (Lecoanet et al. 2017;Rui & Fuller 2023).It has nevertheless been argued that this interpretation is incompatible with the seismic properties of stars with weakly depressed dipolar mixed modes (Mosser et al. 2017).Second, a clear signature of magnetic fields, the asymmetry of the = 1 triplet, has been discovered in the oscillation spectrum of 13 red giant stars.This data, together with a theoretical description of the effects of an arbitrary magnetic field on these oscillation modes, enabled the measurement of the strength of the radial component of the magnetic field in the vicinity of the H-burning shell of these stars, and constraints on their topology (Li et al. 2022(Li et al. , 2023)).A related seismic signature, the deviation from the uniform period spacing of gravity modes (hereafter g modes), led to the detection of core magnetic fields in 11 more red giant stars (Deheuvels et al. 2023).
Beyond red giants, seismic signatures of magnetic fields can of course be searched in other pulsating stars.γ Doradus (γ Dor) stars are good candidates because (i) they oscillate in the lowfrequency range, which is most affected by magnetic fields, (ii) their modes are gravity and Rossby modes, which are sensitive to the radiative layers just outside the convective core where a magnetic field is most probably generated by a convective dynamo, and (iii) these gravity and Rossby modes are identified in more than 600 γ Dor stars (Li et al. 2020) observed by Kepler (Borucki et al. 2010) and in about 60 γ Dor stars (Garcia et al. 2022) observed by TESS (Transiting Exoplanet Survey Satellite, Ricker et al. 2015).However, the two seismic diagnostics available to search for magnetic fields in slowly rotating red giants, namely the asymmetry of the = 1 triplet and the deviation from the g-mode uniform period spacing, are not relevant for rapidly rotating stars such as γ Dor stars.A theoretical analysis of the effects of an arbitrary magnetic field in a rapidly rotating star is therefore necessary to search for and interpret potential magnetic seismic signatures in γ Dors.
The effects of a magnetic field on oscillation modes were first considered for acoustic modes in Unno et al. (1989), Gough & Thompson (1990).Both the magnetic field and the rotation were assumed to be small enough to be treated perturbatively.This approach was later applied to high-order g modes, including mixed gravito-acoustic modes.In these studies, the poloidal magnetic field has been assumed to take the simple form of a dipolar field either aligned (Hasan et al. 2005;Gomes & Lopes 2020;Bugnet et al. 2021;Mathis et al. 2021) or inclined (Loi 2021) with respect to the rotation axis.In Li et al. (2022), we extended to an arbitrary field the perturbative analysis of magnetic and rotation effects on g modes.This extension is significant as we do not want seismic measurements of a magnetic field to depend on an a priori assumption about the magnetic field topology.
In γ Dor stars, the rotation cannot be treated with a perturbative method, as the Coriolis force significantly affects the g modes and induces new modes such as the Rossby modes (hereafter the r modes).A well-known approximation, the traditional approximation of rotation (TAR), allows Coriolis effects to be taken into account while keeping the eigenvalue problem separable in the spherical coordinates.Although the TAR has its limitations (Friedlander 1987;Gerkema & Shrira 2005;Ballot et al. 2012;Ouazzani et al. 2017Ouazzani et al. , 2020)), it is most often accurate and has been successful in identifying oscillation modes and measuring core rotation rates in γ Dor stars (Bedding et al. 2015;Van Reeth et al. 2015, 2016;Christophe et al. 2018;Li et al. 2019Li et al. , 2020)).The magnetic effects on TAR gravito-inertial modes have already been studied in the case of dipolar poloidal magnetic fields either aligned (Prat et al. 2019) or inclined (Prat et al. 2020) with respect to the rotation axis.In both studies, frequency shifts induced by the dipolar magnetic field have been computed for one particular stellar model.In the present paper, we extend these works (i) by considering an arbitrary magnetic field and (ii) by deriving approximate analytical formulae relating the magnetic frequency shifts to a weighted average of the internal magnetic field, the rotation rate, and stellar structure properties for the modes that are most frequently identified in γ Dor stars.We then propose simple methods to search for magnetic seismic signatures and potentially measure magnetic fields in the frequency spectra of γ Dor or SPB stars.
The paper is organised as follow: the formula of the frequency shift induced by an arbitrary magnetic field on high-order gravito-inertial modes is derived in Sect. 2. In Sect.3, properties of the magnetic frequency shift are studied for the modes observed in γ Dor: approximate analytical formulas are derived and tested (Sect.3.1), the magnetic shift of the different modes are compared to help detect magnetic signature in observed spectra (Sect.3.2), and methods to detect magnetic signatures are tested on synthetic spectra (Sect.3.3).In Sect.4, the results are summarised and discussed.

Derivation of the magnetic shift
Time-harmonic, ∝ exp(iωt), small perturbations of a magnetic and uniformly rotating star are solutions of the following eigenvalue problem: where ξ is the displacement vector, 2iωΩ∧ξ is the Coriolis force term, and L 0 is a linear operator, that has the same form as in a non-rotating, non-magnetic star, with ρ and ψ the equilibrium density and gravitational potential, and p , ρ , and ψ the perturbations to pressure, density, and gravitational potential.L L (ξ) is the linear operator associated with the Lorentz force, with B the star magnetic field, B = ∇ × (ξ × B) the Eulerian perturbation to B, and µ 0 the magnetic permeability (e.g.Li et al. 2022).
We performed a first-order perturbative analysis of the effects of a magnetic field, assuming L L (ξ) L 0 (ξ)+2iωΩ∧ξ .The oscillation frequency and the displacement vector were written as the sum of unperturbed and perturbed quantities, that is, ω = ω 0 + ω 1 , ξ = ξ 0 + ξ 1 , with ω 1 ω 0 and ξ 1 ξ 0 , and inserted into Eq.( 1).At zero order, the system reduces to the eigenvalue problem governing the oscillations of a rotating non-magnetic star: In the low-frequency domain, this problem can simplified using the TAR (Eckart 1960;Lee & Saio 1997;Townsend 2003).In this framework, the centrifugal force, the perturbations of the gravitational potential, and the latitudinal component of the rotation vector, Ω, are neglected.The eigenvalue problem then becomes separable in the spherical coordinates so that the eigenfunctions, ξ 0 , can be written as where the variations in latitude are given by the three functions, H r , H θ , H φ , which are known as Hough functions.The radial Hough function, H r , is the solution of Laplace's tidal equation, where µ = cos θ, θ is the colatitude, and s = 2Ω ω 0 is the spin parameter.For a non-rotating star, the eigenvalues, Λ, are equal to ( + 1) for all azimuthal numbers, − ≤ m ≤ , and the eigenfunctions are a linear combination of the 2 + 1 spherical harmonics Y m , denoting the degree of the spherical harmonic.The m-degeneracy is lifted in a rotating star.The eigenvalues, Λ ,m (s), and eigenfunctions, H ,m r (θ, s), of Laplace's tidal equation that correspond to g modes at a vanishing rotation are labelled with and m.For the r modes that only exist in rotating stars, a different label, k, is used instead of .In the following, the , m, s or k, m, s dependence will most often be implicit to A2, page 2 of 18 ease the reading.The horizontal Hough functions, H θ , H φ , can be derived from H r through Eqs.(B.1) and (B.2).
In the radial direction, the eigenvalue problem is identical to the radial eigenvalue problem of the non-rotating case except that ( + 1) is replaced by Λ.Since γ Dor stars oscillate in highorder modes, the Wentzel-Kramers-Brillouin (WKB) solution of the radial eigenvalue problem provides a good approximate solution as long as the mode radial wavelength is much shorter than the scales of variation in the star structure.It reads: with r i and r o are the inner and outer boundaries of the g-mode cavity, n is the mode radial order, and N(r) is the Brunt-Väisälä frequency.As Λ depends on s, Eq. ( 7) is an implicit equation for s.
For a given ( , m), it relates s with n.In the following, we shall also use the WKB solution for the horizontal component of the displacement vector, r 2 represents the square of the local radial wavelength.
At first order, the eigenvalue problem, Eq. ( 1), corresponds to the following eigenvalue problem for the perturbed frequency, ω 1 , and the perturbed eigenfunctions, ξ 1 : As discussed in Gough & Thompson (1990) and Mathis & Bugnet (2023), the modification of the hydrostatic equilibrium by the magnetic field is neglected as its effect on low-frequency g modes is negligible compared to the effect of the perturbed Lorentz force L L (ξ 0 ).Taking the scalar product of this equation with ξ 0 and using the self-adjoint character of the operator, L 0 + 2iωΩ∧, we obtain where the inner product is with V being the star volume.This general expression of the magnetic shift was already derived in Prat et al. (2019).
In the case of slowly rotating stars, Hasan et al. (2005) showed that the expression of the magnetic shift can be significantly simplified for low-frequency g modes provided the magnetic field varies over scales larger than the mode radial wavelength.As detailed in Appendix A, this simplification also holds for the low-frequency gravito-inertial modes as long as B φ /B r (N/ω 0 ) and B θ /B r (N/2Ω) 1/2 (N/ω 0 ) 1/2 .In this regime, the magnetic frequency shift only involves the horizontal displacements and the radial field component, where This formula generalises to an arbitrary magnetic field the expression obtained by Prat et al. (2019Prat et al. ( , 2020) ) for an oblique dipolar field.It is valid for low-frequency (or equivalently highorder) TAR gravito-inertial modes in a magnetic field that is not strongly dominated by its horizontal components.
Using the WKB solution for ξ h (r) and the stationary phase approximation, the integral involving ξ h (r) can be simplified as where f (r) is an arbitrary smooth function.The magnetic frequency shift then becomes where The term I depends on the stellar structure within the oscillating cavity, and B 2 r φ is the azimuthal average of B 2 r , while K r and K θ are, respectively, radial and latitudinal weight functions, as they verify that In a non-rotating star, this expression reduces to where Ŷm (θ) is the latitudinal part of the spherical harmonic of degree, , and the azimuthal number, m, defined as Y m (θ, φ) = Ŷm (θ)e imφ .As was expected, we recovered the diagonal terms of the magnetic matrix derived by Li et al. (2022).We note that the non-diagonal terms of the magnetic matrix are not relevant for the rapidly rotating stars considered in the present study.Indeed, A2, page 3 of 18 they arise in slowly rotating from the coupling of the nearly degenerate components of the multiplet.
To analyse the properties of the magnetic frequency shift in a rotating star, we found it convenient to rewrite Eq. ( 12) as where we introduced two dimensionless factors, and and two different averages of B 2 r , the average of B 2 r in the azimuthal and radial directions, the average in the radial direction being weighted by K r (r) and performed between the inner and outer radii of the mode cavity, and the average of B 2 r performed over the whole volume of the mode cavity, with a K r weight in the radial direction.
The factor, F, depends on the mode, that is on (or k), m, and s, but not on the magnetic field or the stellar structure.The factor, T B , is a weighted average of B 2 r scaled by B 2 r .It depends on the field topology but also on the mode through the latitudinal weight function, K θ (θ).The product, I B 2 r , depends on the stellar structure and on the mode considered through the radii of the mode cavity.

Analysis of the magnetic frequency shift for the modes identified in γ Dor stars
In the following, we focus our analysis on the four types of mode that are most frequently identified in γ Dor stars: the two sectoral prograde g modes (also named Kelvin modes), ( = 1, m = −1) and ( = 2, m = −2), the r mode labelled (k = −2, m = 1), and the zonal dipolar g mode ( = 1, m = 0).What is actually identified in these stars are frequency patterns formed by different radial orders, n, of these modes.Among the 611 γ Dors analysed by Li et al. (2020), 594 stars show a frequency pattern of ( , m) = (1, −1) modes, 172 a pattern of ( , m) = (2, −2) modes, and 110 a pattern of (k = −2, m = 1) modes, while 29 stars show a pattern of ( , m) = (1, 0) modes, although only 11 stars among them have a typical γ Dor fast rotation (P rot ∼ 1 day).Retrograde dipolar modes have also been observed but only in slowly rotating stars and will not be considered here.The radial orders are rather high and the patterns cover a large interval of radial orders.Indeed, for the two prograde sectoral modes, the radial order distribution has a median value of n ∼ 50 and a mean pattern length of 30 radial orders.For the (k = −2, m = 1) r modes, the median is lower, n ∼ 36, but the mean pattern length is also 30.
Similar radial orders translate into different spin parameter intervals for the different modes.Indeed, for the prograde dipolar mode ( = 1, m = −1), the spin parameter distribution peaks at s = 5 and goes up to ∼20, while it peaks at s = 2.5 and goes up to ∼10 for the quadrupolar sectoral modes ( = 2, m = −2).The interval of (k = −2, m = 1) r mode spin parameters range from 6 to ∼35, with a peak at s ∼ 9.And among the 11 fast rotating stars with zonal dipolar g modes ( = 1, m = 0), the spin parameter ranges from ∼1 to ∼2.5.
For a given frequency pattern, that is, for a given star and a fixed ( , m) or (k, m), the magnetic frequency shift varies with the spin parameter, s (or equivalently the radial order n).In Sect.3.1, we study these variations and derive and test approximate analytical formulas for the frequency shift of the selected modes.In Sect.3.2, we compare the magnetic shifts of these different modes and study a simple seismic diagnostic of the presence of magnetic fields when both ( , m) = (1, −1) and ( , m) = (2, −2) modes are present.In Sect.3.3, we construct synthetic frequency patterns affected by magnetic fields and test simple methods of detecting magnetic signatures.

Variations in the magnetic frequency shift with the spin parameter
We successively studied the terms F, T B , and I B 2 r φ,r involved in the expression Eq. ( 21) of the frequency shift.To do so, we used both numerical and analytical calculations, the latter being valid in the s 1 limit.This allowed us to derive and test approximate analytical formulas of the magnetic frequency shifts.

The F factor
We recall that for a given ( , m) the F factor defined by Eq. ( 22) only depends on s.The numerical determination of F requires us to compute Λ(s) and H r (θ, s) from Laplace's tidal equation, to derive H θ and H φ from (Λ(s), H r (θ, s)), and to compute the integrals involved in the F factor.These calculations are detailed in Appendix B. Approximate analytical expressions of F can also be obtained from analytic expressions of Λ and the Hough functions valid in the high-s limit (Townsend 2003).These calculations are detailed in Appendices C.1 and C.2.

The T B factor
The factor, T B , defined by Eq. ( 23) is a weighted latitudinal integral of B 2 r φ,r (θ) normalised by B 2 r .Here, we analyse its variations with the spin parameter and the magnetic field topology for the same four modes.For a given mode, the weight function, K θ (θ), only depends on the spin parameter.Figure 2 shows the latitudinal profile of the weight function of the four modes computed at s = [0.1, 1, 8, 15].We observe that all the weight functions concentrate towards the equator as s increases.This is a consequence of the well-known tendency of the Hough functions to equatorial concentration (Lee & Saio 1997;Townsend 2003).We also observe that at a fixed s this effect is more or less important, depending on the mode.The axisymmetric dipolar mode is the most concentrated, followed by the ( , m) = (2, −2) mode, the ( , m) = (1, −1) mode, and finally the r mode, which is not yet equatorially concentrated at s = 8.This behaviour is in qualitative agreement with the high-s asymptotic form of the Hough functions (Townsend 2003) that indicate that the halfwidth of H θ are 1/s for the ( , m) = (1, 0) mode, √ 3/(−ms) for the ( , m) = (2, −2) and ( , m) = (1, −1) modes, and √ 3/(s − 1) for the r-mode.
As a consequence of the equatorial concentration, the K θ average of B 2 r is increasingly weighted towards the equator so that for very large spin parameters we expect T B ≈ B 2 r φ,r (θ = π/2)/ B 2 r .In this limit, T B is thus constant with s and identical for all the modes.From a Taylor expansion of B 2 r at the equator, we can determine how T B (s) tends towards For the two sectoral prograde modes and the r mode, we obtain with α = 4 for ( = 1, m = −1), α = 8 for ( = 2, m = −2), and α = 4/3 for (k = −2, m = 1), while, for the dipolar mode axisymmetric, This approximation requires that near the equator B 2 r φ,r (θ) varies on scales larger than the width of the weight function, K θ .To the first two orders in the 1/s expansion, T B depends on B 2 r φ,r and its second derivative at the equator.The next order would involve the fourth derivative of B 2 r φ,r at the equator.These expressions also indicate that the convergence towards B 2 r φ,r (θ = π/2)/ B 2 r is faster (∝ 1/s 2 ) for the dipolar axisymmetric mode than for the other three modes (∝ 1/s).
To test these expressions, we computed T B (s) for an oblique dipolar field.This is a favourable case since the field varies on large scales around the equator.If β denotes the inclination angle of the dipole with respect to the rotation axis, the radial component of the oblique dipole reads B r (r, θ, φ) = B 0 b(r)(cos θ cos β + sin θ cos φ sin β) (Prat et al. 2020).A straightforward calculation (see Appendix D) then yields Figure 3 shows both the numerical and analytical forms of T B (s) for the four modes and for three inclination angles, β = 0 • , 47 • , 90 • , of the oblique dipole.We find that the high-s asymptotic expressions provide good approximations above some spin  parameters, say above s 5 for ( , m) = (1, −1), s 2.5 for ( , m) = (2, −2) and ( , m) = (1, 0), and s 20 for (k, m) = (−2, 1).As for the F factor, the convergence is more or less rapid depending on the mode tendency to equatorial concentration.We also observe that, as expected, T B (s) converges towards the same value, B 2 r φ,r (θ = π/2)/ B 2 r = 3 2 sin 2 β, for all the modes.
3.1.3.The I B 2 r φ,r term The magnetic frequency shift can also depend on s through the radii of the mode cavity.This dependency is contained in the product of the two terms I and B 2 r φ,r : The inner and outer cavity radii being determined by the conditions ω 0 = Min(N, S ), where S = c s r √ Λ is the Lamb frequency A2, page 5 of 18 and c s the sound speed, they are likely to depend on the mode frequency, ω 0 , and thus on s.The inner radius, r i , remains very close to the outer radius of the convective core for highorder g modes though.This is due to the sharp increase in the Brunt-Väisälä frequency at the bottom of the radiative zone.Conversely, the outer radius, r o , that is located in the radiative envelope of γ Dors does depend on the mode frequency.The vertical lines in Fig. 4 illustrate, in the case of the ( , m) = (1, −1) mode, the variations in r o as s takes three different values [2,5,15].
The variation in the outer radius of the mode cavity does not affect significantly r o r i N/r dr, the denominator of I B 2 r φ,r , because N/r peaks strongly near the bottom of the radiative zone.But it can affect the numerator, that is, the integral, 3 is dominant and increases rapidly in the star envelope up to the surface convective zone where N vanishes.Assuming a uniform magnetic field, Fig. 5 shows that this induces a sharp increase in I B 2 r φ,r between s ∼ 2 and s ∼ 6 for ( , m) = (1, −1) modes.For higher spin parameters, r o remains fixed at the bottom of the surface convective zone so that I B 2 r φ,r no longer varies with the spin parameter.Considering a dipole-like ∝ 1/r 3 decrease in B r within the envelope instead of a uniform field would not significantly modify this picture.Conversely, if the magnetic field decreases by various orders of magnitude between r i and r o , the variations in r o will not affect the weighted integral of B 2 r φ .In this case, I B 2 r φ,r depends neither on s nor on ( , m).

Approximate analytical formulas
Using the analytical forms of F and TB, we obtain the following approximate expressions for the magnetic frequency shifts of the four mode considered: ) where we introduced B 2 eq = B 2 r φ,r (θ = π/2) and D eq = dθ 2 (θ = π/2)/B 2 eq .These expressions are valid in the highs limit.They require that B 2 r φ,r varies on a length scale larger than the half-width of the mode weight function, K θ .They imply that in the high-s limit these modes probe the magnetic fields only in the equatorial region.
The analytical expressions allow us to conclude on the s dependence of the frequency shift in the case of a magnetic field buried at the bottom of the radiative zone for which the product, IB 2 eq , has a negligible variation with the spin parameter.
A2, page 6 of 18  Then, if B 2 eq 0, the magnetic frequency shift varies as 1/ω 3 0 for the two prograde sectoral modes and the r mode, and as 1/ω 5 0 for the dipolar axisymmetric modes.If B 2 eq = 0, the frequency shift rather varies as 1/ω 2 0 for the two prograde sectoral modes and the r mode, and as 1/ω 3 0 for the dipolar axisymmetric modes.We should nevertheless keep in mind that the conditions B φ /B r (N/ω 0 ) and B θ /B r (N/2Ω) 1/2 (N/ω 0 ) 1/2 are not verified at the equator in this case.
The magnetic shift scaled by ) 3 , that is, equal to s 3 F(s)T B , is shown in Fig. 6 for ( , m) = (1, −1) modes in an oblique dipolar field.The numerical computation converges towards the analytical form for high spin parameters (the formula for oblique dipolar fields is derived in Appendix D).We see that the magnetic shift is larger for lower frequencies and strong equatorial radial fields.While the magnetic shift decreases at low s, it does not vanish in practice because s has a lower limit associated with the upper limit of the g mode frequencies.
We recall that if the envelope magnetic field is dominant, the variation of I B 2 r with s (as shown in Fig. 5 for ( , m) = (1, −1)) must also be taken into account in the s dependence of the frequency shift.

Comparing the magnetic shifts of different ( , m) modes
We then compared the magnetic shifts of different ( , m) modes.This allowed us to determine the mode that produces the highest shifts but also to find relations between the different magnetic shifts that can be used to detect or confirm magnetic signatures.
Among We first compute and analyse the ratio of the magnetic shifts (Sect.3.2.1)and then discuss how they can be used to detect magnetic signatures in frequency spectra (Sect.3.2.2).

The ratio of the magnetic shifts of different ( , m) modes
We considered modes that had the same radial order, n, but a different horizontal quantum number, ( , m) or (k, m).From the dispersion relation Eq. ( 7), it follows that the product s √ Λ is identical for all these modes: Thus, knowing the spin parameter of, say, an ( = 1, m = −1) mode, s 1,−1 , we can determine the spin parameters of the other modes by solving the equation s This can be done numerically or using asymptotic forms of Λ (see Eq. (E.1)).It is not necessary to specify the common radial order to determine the relations between the spin parameters, although n can be determined by Eq. ( 34) once Ω, Π 0 , and g are specified.
A2, page 7 of 18 For a given star, the ratio of the magnetic frequency shifts is equal to the ratio of s 3 F(s)T B (s).The I B 2 r term cancels because the outer radii of the cavity is determined by s c s and is thus identical for the different modes.The ratio between the magnetic shifts of the ( = 2, m = −2) and ( = 1, m = −1) modes, denoted ω 1 (2, −2, n)/ω 1 (1, −1, n), is shown in Fig. 7.It has been computed for oblique dipolar fields of three inclination angles and for the observed interval of prograde dipolar mode spin parameters, that is, s 1,−1 ∈ [0, 20], which translates into s 2,−2 ∈ [0, 10] for constant s √ Λ.The ratio of the s 3 F(s) terms, which does not depend on the field topology, is also displayed.
We observe that the magnetic shift ratio, ω 1 (2, −2, n)/ω 1 (1, −1, n), is always smaller than one and that, for all inclination angles, it converges rapidly towards 1/2 as s increases.The asymptotic value of 1/2 is predicted by the asymptotic form of the magnetic shift of the two modes (see Eqs. ( 30) and ( 31)) and the asymptotic relation between the two spin parameters: s 2,−2 ≈ s 1,−1 /2.The rapid convergence of ω 1 (2, −2, n)/ω 1 (1, −1, n) can be understood by looking at the weight functions of the two prograde sectoral modes computed, respectively, at s 1,−1 and s 2,−2 .Figure 8 indeed shows that both functions become rapidly similar as s increases above one.This comes from the fact that at high s the prograde sectoral horizontal Hough functions behave as Gaussian functions whose widths are proportional to 1/(−ms) 1/2 .Since s 2,−2 ≈ s 1,−1 /2, the Gaussian widths of the two ( = 2, m = −2, n) and ( = 1, m = −1, n) Hough functions are very close, so that we expect and verify that the weight functions K 1,−1 θ (s 1,−1 ) and K 2,−2 θ (s 2,−2 ), and therefore the T B factors of the two modes, are also very close.Crucially, this implies that the property, ω 1 (2, −2, n) ≈ ω 1 (1, −1, n)/2, holds for all field topologies but the rare ones that probe the small differences between the two latitudinal weight functions.
The ratio between the magnetic shifts of the r mode (k = −2, m = 1) mode and the prograde dipolar mode,  The ratio of the product, F(s)s 3 , which does not depend on the field topology, is also displayed (black).interval translates into s −2,1 ∈ [6, 20] for constant s √ Λ.As shown in Fig. 10, the weight function, K θ , of the two modes are strongly concentrated towards the equator at high s.The ratio of the magnetic frequency shifts tends towards three in accordance with the analytical forms, Eqs. ( 30) and (32), and the asymptotic relation, s −2,1 ≈ 3s −1,1 .We also observe that ω 1 (−2, 1, n)/ω 1 (1, −1, n) seems to diverge as s −2,1 approaches six.Indeed, at this point, s 3 1,−1 vanishes more rapidly than F −2,2 , so that the ratio of s 3 F(s) terms goes to infinity.This s 1,−1 → 0 limit is however not relevant in practice because both magnetic shifts vanish and will thus be undetectable.It remains that, for s 1,−1 = 2, the ratio of the s 3 F(s) terms is equal to ∼7.For these low s 1,−1 , the weight function profiles (see the blue curve in Fig. 10) also shows that the r mode is more sensitive to fields concentrated at the pole than the dipolar sectoral mode.This A2, page 8 of 18 Fig.10.Comparison of the latitudinal weight functions, K θ , of the = 1, m = −1 (continuous lines) and k = −2, m = 1 (dashed lines) modes of same radial order, n.The weight functions were computed for three (s 1,−1 , s −2,1 ) couples: (2.2, 8.9; blue), (9.5, 29.5; green), and (18.6, 56.8; red).Fig. 11.Ratio of the magnetic shifts produced by the zonal dipolar mode, = 1, m = 0, and the prograde dipolar mode, = 1, m = −1, of the same radial order, n, as a function of their spin parameters, s 1,−1 (lower x axis) and s 1,0 (upper x axis).The ratio was computed for three different angles of an inclined dipolar field, β = 0 • (blue), β = 47 • (green), and β = 90 • (red).The ratio of the product, F(s)s 3 , which does not depend on the field topology, is also displayed (black).
provides an opportunity to probe magnetic fields in both polar and equatorial regions.

√
Λ.The decrease in the ratio towards high s is close to 3/(2s 1,0 ), in agreement with the analytical forms, Eqs. ( 30) and (33), and the asymptotic relation, s 1,0 = √ s 1,−1 .For 1 s 1,0 2.5, the range of observed s 1,0 in Li et al. (2020), the magnetic shift ratio is higher than one if β is low enough and is smaller than one for higher β.Consequently, depending on the field the rotational splitting, ω 1 (1, −1, n) − ω 1 (1, 0, n), can be either larger or smaller than the unmagnetised TAR predictions.Such rotational splittings were determined by Li et al. (2020) for 11 stars and it was found that they can be either higher or smaller than the TAR prediction.Magnetic fields are thus a We now propose a method of searching for and measuring magnetic fields in stars where both ( = 2, m = −2) and ( = 1, m = −1) frequency patterns are present.If frequencies with the same radial order have been identified in these two frequency patterns, it is possible to compute the difference between the ( = 1, m = −1, n) frequency and half the ( = 2, m = −2, n) frequency, denoted δK (as it involves two Kelvin modes).If magnetic effects can be treated perturbatively, it reads The first approximated equation, Eq. ( 36), comes from the relation, ω , can be determined.This can be done in principle using a seismic determination of Ω and an unmagnetised oscillation model.But, we also expect δK 0 to be small and, if it is negligible, compared to 3 4 ν 1 (1, −1, n), then Eq. (37) holds, giving direct access to the magnetic frequency shift from the two observable frequencies, ν in (1, −1, n) and ν in (2, −2, n).
The difference, δK 0 , should indeed be small because it vanishes for high-s, high-order TAR modes (as s 2,−2 = s 1,−1 /2 in this limit).If we still assume high n TAR modes but now consider finite s, δK 0 is 9.6 nHz at s 1,−1 = 5 (for a rotation period of one day) and decreases for higher s 1,−1 .We also determined A2, page 9 of 18 δK 0 with more realistic frequencies computed from a stellar model and using the TOP code (Reese et al. 2006(Reese et al. , 2009) ) to relax the TAR of the Coriolis force.The structure of this model has been computed with the CESAM evolution code (Morel 1997;Morel & Lebreton 2008).It has a mass, M = 1.4 M , a radius, R = 1.38 R , and a central hydrogen abundance, X c = 0.68.The oscillation spectrum has been computed for a rotation, ν rot = 1.28 d −1 .In Fig. 13, we show that δK 0 decreases from ∼20 nHz to ∼7 nHz between s −1,1 = 5 and s −1,1 = 7.This remains small although it is a bit higher than the same difference computed for high n TAR modes.
While δK 0 tends to vanish as n increases, the magnetic shift, ν 1 (1, −1, n), increases with n (or s), as described by Eq. ( 30).Thus, observing that δK in ≈ δK 0 + 3ν 1 (1, −1, n)/4 increases with n is already a clue of the presence of a magnetic field.Fitting δK in with Eq. ( 30) then gives access to B 2 eq and D eq .To test further the method proposed here, the properties of the difference δK 0 should be studied for various stellar models, including the effects of steep chemical gradients above the convective core, differential rotation, and centrifugal deformation.These effects are briefly discussed in Sect. 4.

Probing magnetic signatures in frequency spectra
In this section, we construct frequency spectra of prograde sectoral modes, ( = 1, m = −1) and ( = 2, m = −2), affected by a radial magnetic field in order to test simple methods of revealing magnetic signatures.
This specific period, Π 0 = 4122 s, has been calculated from a mid-main-sequence stellar model of a γ Dor star (M = 1.6 M , R = 2.13 R , X c = 0.35), which we use in this section.It has been computed with the CESAM evolution code, by including turbulent diffusion of chemicals, D = 400 cm 2 s −1 .It possesses a convective core with a radius r cc = 0.069R.
The magnetic frequency shifts, ν 1,−1 1 and ν 2,−2 1 , were determined by Eqs. ( 30) and ( 31), assuming D eq = 0 for simplicity.This stellar model and three different radial field profiles were considered to compute IB 2 eq and then the magnetic frequency shifts.The first case (profile U) is a uniform field, B r = 20 kG, between r cc , the radius of the convective core, and r = 0.5R.For the two other profiles, we consider that B r is 24.3 kG at r = r cc and decreases exponentially by one order of magnitude either at r = 0.2R (profile B+) or at r = r cc + 0.02R (profile B-).The magnetic frequency shifts, ν 1,−1 1 , computed for the frequency corresponding to s 1,−1 = 5 are 110, 37, and 9.1 nHz, for these three configurations, respectively.The magnetic frequency shift, ν 2,−2 1 , is simply ν 1,−1 1 /2 at s 2,−2 = s 1,−1 /2.These values are larger or just above the frequency resolution of Kepler for classical pulsators (1/T = 7.9 nHz for T = 4 yr).In practice, models are not able to fit observed frequencies at this level of accuracy so that the magnetic frequency shifts should rather be compared with other source of frequency deviations such as the differential rotation or the steep gradient of molecular weight overlying the convective core.
Figure 14 shows the observed period spacings, ∆P, of the ( = 1, m = −1) modes as a function of the period, P, for the different magnetic cases.The ∆P − P diagrams are frequently used to determine the period, Π 0 , and the rotation frequency, ν rot , using the TAR (e.g.Van Reeth et al. 2015;Li et al. 2019).The magnetic field does not produce noticeable signatures on the shape of the profile of ∆P: it slightly changes its mean value and its slope.Such changes can be mimicked by a modification of Π 0 and/or ν rot and may bias the determination of these parameters.To show this, we analysed the magnetically shifted frequency patterns with methods that do not take magnetic effects into account.We thus applied the method developed by Christophe et al. (2018) to the frequency sets of prograde dipolar and quadripolar modes generated with the magnetic field profiles U, B+, and B-.This method returns the period spacing, Π 0 , and the rotation frequency, ν rot , that best fit a frequency set within the TAR.Once Π 0 and ν rot were determined, we identified the radial orders, n, of the modes and the offset, g (see Eq. ( 7)).Table 1 summarises the results obtained for the different cases A2, page 10 of 18  and Fig. 15 shows the difference between the simulated frequencies and the fitted frequencies.First, we notice that the fits are in good agreement and the residuals between input and fitted frequencies do not generally exceed a few nHz.It means that a model without a magnetic field can correctly reproduce data that are indeed affected by a magnetic field.However, the values of Π 0 and ν rot are biased towards larger values.The stronger the magnetic shift is, the stronger the bias is: ( = 1, m = −1) modes are more sensitive to magnetic field than ( = 2, m = −2) modes and the effect are the largest for profile U and the smallest for profile B-.The presence of a magnetic field also lowers the apparent value of the offset, g .This effect can be strong enough to not only modify g , but also to shift the apparent values of the radial order, n (it is the case for the profile U and for = 1 with the profile B+, see Table 1).
Considering a real star for which only ( = 1, m = −1) modes are observed, biases of Π 0 and ν rot can hardly be spotted since the real values of these quantities are a priori unknown.However, for strong fields (such as with profile U, or with profile B+ to a lesser extent) we may identify the existence of a magnetic effect by analysing the departure of the observed frequencies from the TAR prediction.This difference can be revealed in a stretched échelle diagram, in which mode frequencies are plotted as a function of the period in the corotating frame, stretched with the Λ function, modulo the period spacing, Π 0 (Fig. 16).Within the TAR, modes form vertical lines in stretched échelle diagrams.We see in Fig. 16 that a magnetic field distorts the ridge in a S-shape profile with an amplitude larger than 500 s (for profile U).This signature can be compared to various effects stud- ied in Christophe et al. (2018), especially in their Fig. 5. First, we notice that the differences between complete computations and the TAR are clearly smaller (<100 s).Second, buoyancy glitches generate oscillating signatures with amplitudes similar to the ones produced by a strong magnetic field, but their periodic nature helps to discriminate between both.Finally, a strong differential rotation also produces a smooth change in the ridge shape, but with a very small amplitude and mainly induces avoided crossings between modes, which create a quite different signature to the one shown in Fig. 16.
Having both ( = 1, m = −1) and ( = 2, m = −2) frequency patterns aids the detection of seismic magnetic signatures.Table 1 indeed shows clear differences in the values of Π 0 obtained for each pattern.For the U magnetic profile, the ( = 1, m = −1) pattern provides a period spacing 21% larger than the ( = 2, m = −2) pattern.This corresponds to a difference of 922 s, which is larger than most errors in the Π 0 determination quoted in Li et al. (2020).The difference in the ν rot values is smaller, although the 0.05 d −1 discrepancy is also larger than most errors in ν rot values quoted in Li et al. (2020).For the weakest magnetic shifts (profile B-), the difference in Π 0 is at the level of the best-fitting errors indicated in Li et al. (2020).
The ( = 1, m = −1) and ( = 2, m = −2) frequency patterns can also be combined to detect and measure the magnetic frequency shifts.Using the radial order identifications provided by the fits (see Table 1), the difference, δK in , introduced in Sect.3.2.2can be computed as a function, n.In a non-magnetic star, δK in = δK 0 is small (of the order of 10 nHz) and a decreasing function of n, whereas a magnetic term increasing with n adds up in a magnetic star.As the order identifications may be wrong, differences formed by shifting the radial orders of one of the pattern ν in ( = 1, m = 1, n + n s ) − ν in ( = 2, m = 2, n)/2 were computed until the smallest positive difference was found.The result of this process is shown in Fig. 17 for the three magnetic profiles considered.We observe that the smallest positive difference is an increasing function of n for the three magnetic profiles.As explained in Sect.3.2.2, this quantity is close to 3ν 1,−1 1 /4 and can be used to derive the magnetic seismic observables, B eq and D eq .Here, δK 0 = 0 because we computed the unperturbed frequencies in the asymptotic regime.In real stars, the magnetic term will need to be higher than δK 0 for the smallest positive difference, ν in ( = 1, m = 1, n + n s ) − ν in ( = 2, m = 2, n)/2, to be an increasing function of n.Considering that δK 0 is of the order A2, page 11 of 18  17.Magnetic shift computed from the difference, δK in for the magnetic profile U (orange), B+ (red), and B-(green), as a function of the apparent radial order, n, of = 2 modes.We notice that for the profile U, the apparent radial orders are shifted compared to the real ones (see Table 1 and text for a discussion).Dashed lines are obtained when the radial order of = 1 modes is incorrectly identified compared to the radial order of = 2 modes by ±1.
of 10 nHz as in the realistic model considered in Sect.3.2.2(see Fig. 13), a magnetic signature would be detected for the magnetic profiles U and B+.

Summary and discussion
We used a perturbative method to investigate the frequency shift induced by an arbitrary magnetic field on gravito-inertial modes.From the general expression of the magnetic frequency shift, Eq. ( 10), we derived approximated simplified formulas valid for short radial wavelength gravito-inertial modes described by the TAR.This allowed us to determine the frequency dependence of the magnetic shift for the modes most often observed in γ Dor stars, the ( , m) = (1, −1), (2, −2), (1, 0) g modes and the (k, m) = (−2, 1) r mode, and to define the seismic observable quantities and their relation with the internal magnetic field.We find that the magnetic frequency shift is larger for low-frequency modes and for strong equatorial radial fields.Comparing magnetic shifts of different modes with the same radial order, n, we find that the (k = −2, m = 1) r mode is the most sensitive to magnetic field followed by the ( = 1, m = −1) g mode.We then proposed simple methods to detect magnetic signatures: the characteristic form of the error when the frequency pattern is analysed with a non-magnetic model, the discrepancy between the period spacings, Π 0 , determined, respectively, with the ( , m) = (1, −1) and the ( , m) = (2, −2) frequency patterns analysed with a nonmagnetic model, and the determination of the magnetic shift from the difference, Before starting the discussion of our results, we first have to mention a clear discrepancy between our results and those presented in Prat et al. (2019Prat et al. ( , 2020)).The magnetic shifts computed in these papers result in a sawtooth pattern in the period spacing versus period relation.This feature is incompatible with the smooth variations of the magnetic shift with s (or n) that we find in both our numerical and analytical results.We notice that in Prat et al. (2019) the equations relating the horizontal Hough functions to the radial Hough function given in their Appendix A are not correct.
In the following, we first discuss the approximations made to derive the simplified formulas of the magnetic frequency shifts, then the conditions for applying the perturbative method, and finally the origin and detectability of magnetic fields in γ Dor radiative zones.

Approximations in the derivation of the magnetic shifts
In the regime of short radial wavelengths, the effect of the radial component of the magnetic field tends to dominate over the horizontal components.As already discussed in Li et al. (2022), the contribution of the horizontal components to the magnetic shift can be neglected if k r B r k θ B θ and k r B r k φ B φ , where (k r , k θ , k φ ) are the components of the mode wave vector.When applied to high-n, high-s TAR oscillations, these conditions translate into upper limits of the ratio of the horizontal to radial components, namely, B φ /B r (N/ω 0 ) and B θ /B r (N/2Ω) 1/2 (N/ω 0 ) 1/2 , (see details in Appendix A).For a ( , m) = (1, −1) mode of spin parameter s = 5 in a M = 1.6 M , R = 2.13 R , X c = 0.35, ν rot = 1 d −1 star model, these constraints are B φ /B r 125 and B θ /B r 55 near the bottom of the radiative zone.If the horizontal field contributions were not negligible, a different frequency dependence of the magnetic shift is expected and this property can be used to identify a situation where the radial field is not dominant.Indeed, as argued by Li et al. (2022) in the case of slowly rotating stars, the magnetic shift varies as ∝1/ω 0 for purely azimuthal fields, a frequency dependence that can be distinguished from the ∝1/ω 3 0 dependence produced by radial fields.The numerical results of Dhouib et al. (2022) obtained for a rapidly rotating star with azimuthal fields also indicate that the increase in the magnetic shift with decreasing ω 0 is much slower than ∼1/ω 3 0 (see their Fig.11).We also used the TAR to describe the effects of the star rotation on g modes.The TAR has been tested against two oscillation codes, TOP (Reese 2006;Reese et al. 2009) and ACOR (Ouazzani et al. 2012), which include a full treatment of the Coriolis and centrifugal accelerations.It was found the TAR generally provides a good approximation when the centrifugal deformation is not too strong (Ballot et al. 2012;Ouazzani et al. 2017) although it fails in some specific frequency intervals where the so-called Rosette modes (Ballot et al. 2012;Takata & Saio 2013) are observed or where gravito-inertial modes resonantly interact with convective core inertial modes (Ouazzani et al. 2020).
We recall that, when necessary, the approximations made to obtain simplified forms of the perturbative magnetic frequency shifts can be avoided by directly applying Eq. ( 10) to a chosen field configuration and to unperturbed modes computed for a particular stellar model with full oscillation codes like TOP or ACOR.It may then be necessary to use the full expression of the Lorentz term, ξ 0 , L L (ξ 0 ) , derived in Prat et al. (2020).

The limits of the perturbative methods
Another limit of the present study is that perturbative methods only apply when the strength of the magnetic field is low enough.The frequency deviation produced by the magnetic field must be small compared with the unperturbed frequency (i.e.ω 1 ω 0 ).For the synthetic spectra studied in Sect.3.3, the magnetic shifts indeed remained much smaller than the unperturbed frequencies (the largest ratio is ω 1 /ω 0 = 0.048 and corresponds to the strongest field and lowest frequency considered).Non-perturbative calculations are nevertheless necessary to determine the accuracy of the perturbative method.In a recent study, Rui et al. (2023) investigated the regime of A2, page 12 of 18 strong non-perturbative magnetic fields in the particular case of a dipolar field aligned with the rotation rate.They show that perturbative methods can deviate significantly from the nonperturbative results when the field strength is close to the critical radial field above which gravity waves no longer propagate radially.As mentioned in the introduction, this critical radial field has been invoked to account for the unexpectedly low-amplitude, = 1 g modes observed in about 20% of red giants (Fuller et al. 2015).In our case, a critical radial field puts constraints on the possible amplitude of the magnetic shifts.Indeed, when a mode is observed, it means the internal radial magnetic field is everywhere smaller than the critical field, which in turn limits the amplitude of the magnetic frequency shift.For a full observed spectrum, the condition must hold for the smallest frequency because B c ∝ ω 2 .To construct the ( , m) = (1, −1) frequency patterns of Sect.3.3, we did choose radial fields such that B r < B c (ν min ).The exact value of B c , which is the subject of ongoing research (Rui & Fuller 2023;Rui et al. 2023), should depend on the mode and the field geometries.Using B c = a 4πρ(r/N)ω 2 0 with a ≈ 1/2 as in Rui & Fuller (2023) for = 1, m = −1 modes, we find B c = 24.3kG at the bottom of the wave cavity near r = r cc and B c = 21.6 kG at r = 0.5R for the stellar model considered.Both constraints limit the amplitude of the magnetic shifts.Finally, we note that an upper limit of ν 1 = ν 0 /(2a 2 ) is obtained by assuming B r = B c (ν 0 ) everywhere in the wave cavity and using Eq. ( 30).But it is very unlikely to be reached as there is no physical reason for the actual stellar radial field to be equal to the field that prevents gravity wave propagation.

Origin and detectability of magnetic fields in γ Dor radiative zones
We now discuss the possible origin of magnetic fields in γ Dor radiative zones and whether they can produce detectable magnetic frequency shifts.Dynamo action generates magnetic fields in the convective core of γ Dors.Field strengths ranging between 10 and 100 kG are expected from numerical simulations of a 2 M convective core by Brun et al. (2005) and equipartition fields determined using the mixing length theory (Cantiello et al. 2016).Numerical simulations also show that the dynamo field pervades the radiative layers just above the convective core, although its strength drops off sharply (Brun et al. 2005;Augustson et al. 2016).In Brun et al. (2005), the field intensity decreases by approximately one order of magnitude between r cc and r cc +0.02R.A 100 kG radial field at r cc with such a decrease into the overlying radiative layers produces magnetic shifts of ν 1 = 154 nHz for a ( = 1, m = −1) frequency, corresponding to s = 5 and the stellar model M = 1.6 M , R = 2.13 R , X c = 0.35 previously considered.However, this field exceeds the critical field, which is 24.3 kG at r = r cc + 0.003R for the smallest frequency of a typical pattern, n ∈ [35,65].The fact that core dynamo fields can exceed the critical field for gravity wave propagation in γ Dor stars has already been put forward in Cantiello et al. (2016, see their Fig. 8 comparing B c to a core dynamo field).These authors even suggest that this modesuppression process explains why some stars in the γ Dors or SPB instability strips are not pulsating (see also Lecoanet et al. 2022).If instead we consider as in Sect.3.3 a convective core field no larger than the critical field with the same exponential decrease, the magnetic shifts do not exceed ∼10 nHz.They would be very difficult if not impossible to detect even with Kepler data.However, the assumption of a sharp exponential decrease in the field strength does not take into account the decrease of convective core mass along the main sequence that leaves behind magnetic fields in radiative layers that were previously convective (Cantiello et al. 2016, see their Fig. 6).This process should increase the amplitude of the magnetic shifts in evolved γ Dors.
Other types of magnetic fields could be present in radiative interiors although we generally lack theoretical constraints on their strength and topology.This includes fossil magnetic fields possibly inherited from earlier evolutionary stages, either from the fully convective pre-main-sequence phase or from the stellar formation (Braithwaite & Spruit 2004;Lignières et al. 2014).Such fields could interact with the convective core as investigated by Featherstone et al. (2009).Moreover, in differentially rotating radiative zones, magnetic fields could be generated by dynamo processes triggered by the Tayler instability (Petitdemange et al. 2023) or the magneto-rotational instability (Guseva et al. 2017).
Spectroscopic observations provide constraints on surface fields.The Ap/Bp stars, a class of intermediate-mass mainsequence stars, are known to harbour surface fields with a typical value of 1 kG.The radial weight function plotted in Fig. 4 shows that envelope fields of such strength would produce a frequency shift of the same order as a 100-kG field at the bottom of the radiative zone.While interesting, this situation should remain unusual among γ Dors because the incidence of Ap/Bp stars is very low (less than 1%) in the mass range 1.4 M < M < 1.8 M (Sikora et al. 2019).It is more common to find Ap/Bp stars among SPB stars, the B star HD 43317 being an example of a star with strong surface field and g-mode oscillations (Lecoanet et al. 2022).Deep spectropolarimetric surveys also revealed that much weaker magnetic fields ∼1 G might be present at the surface of many intermediate-mass stars (Lignières et al. 2009;Blazère et al. 2016).The magnetic shift produced by such envelope fields would be very small though.
In Sect.3.3, we proposed methods to detect magnetic signatures based on the analysis of magnetically shifted asymptotic spectra.While promising, these methods should be tested with more realistic unperturbed spectra, that is, non-asymptotic spectra that include the effects of chemical composition gradients, differential rotation, and centrifugal deformation.In particular, robust generic properties of the small difference, δK 0 = ν 0 (1, −1, n) − ν 0 (2, −2, n)/2, in non-magnetic stars would be useful to detect magnetic fields.We can already notice from the study by Miglio et al. (2008) that steep chemical composition gradients seem to have a very small influence on the relation between ( = 1, n) and ( = 2, n) frequencies in slowly rotating stars (see their Fig.27).It is also known that a weak radial differential rotation has the same effect on all high-n ( , m) modes (Van Reeth et al. 2018;Takata et al. 2020) so that it will not modify the difference.While latitudinal differential rotation has been taken into account in studying gravito-inertial propagating waves (e.g.Ogilvie & Lin 2004;Prat et al. 2018), the effect of a latitudinal differential rotation on the frequencies of gravito-inertial modes has not been determined yet.However, as the ( , m) = (1, −1, n) and (2, −2, n) modes have very similar latitudinal eigenfunctions (see Fig. 8), we anticipate that the effects of a weak latitudinal differential rotation will be similar for the two modes.Strong radial or latitudinal differential rotation would have stronger effects (as shown in Christophe et al. 2018 andTakata et al. 2020, a strong radial differential rotation induces avoided crossings between modes that produce large and localised deviations in the frequency pattern) but they have not yet been detected in γ Dor stars.The methods we propose offer new possibilities to search for internal magnetic fields in existing A2, page 13 of 18 where p is an integer and the horizontal Hough functions are sin θH θ = Θ = H p (u)e −u 2 /2 (C.2) (C.3)where u = Λ 1/2 s 1/2 µ, with µ = cos θ, and H p is the "physical" Hermite polynomial of order, p, and Θ and − Θ are the notations used by Townsend (2003).(C.13) The value of Λ comes from Townsend (2020), who improved the expression given in Townsend (2003).In the following we simply use Λ = 1/9 as we find it is a better approximation for the numerical solution of Λ.In the following we approximate sin 2 θ ≈ 1, as it is the basic assumption of Townsend (2003) asymptotic solutions.For each of the four modes considered, we rewrite the integrals G 1 and G 2 using u or τ variables instead of θ.We then get exact and then asymptotic s 1 expressions of G 1 and G 2 using Wolfram alpha.The resulting asymptotic expressions for F are

C.3. Asymptotic form of the T B (s) factor
For large spin parameters, s, the weight function, K θ , is focused near the equator.Asymptotically it plays the same role as a Dirac function, so that T B (s) → B 2 r φ,r / B 2 r , as s → +∞.In order to study the next order in the high-s asymptotic regime, we expanded B 2 r φ,r around the equator using the µ = cos θ variable, namely

Fig. 5 .
Fig. 5. Variation in IK r (r) with s computed for the ( , m) = (1, −1) modes in the case of a uniform radial field (which implies that IK r (r) = I) and for a M = 1.6 M , R = 2.13 R , X c = 0.35 stellar model.

Fig. 14 .
Fig. 14.Period spacing, ∆P, of ( = 1, m = −1) modes as a function of their periods for a model with a magnetic field following the radial profiles U (orange), B+ (red), and B-(green).The non-magnetic case is also shown (black).

Fig. 16 .
Fig. 16.Stretched échelle diagram from a magnetically shifted frequency spectra (profile U), obtained by applying the method of Christophe et al. (2018).Modes are plotted twice for clarity.
Fig.17.Magnetic shift computed from the difference, δK in for the magnetic profile U (orange), B+ (red), and B-(green), as a function of the apparent radial order, n, of = 2 modes.We notice that for the profile U, the apparent radial orders are shifted compared to the real ones (see Table1and text for a discussion).Dashed lines are obtained when the radial order of = 1 modes is incorrectly identified compared to the radial order of = 2 modes by ±1.

F
Fig C.1 presents a comparison of the asymptotic expression and numerical calculation of F for the ( = 1, m = 0) mode.

Table 1 .
Christophe et al. (2018)and g fitted with the method ofChristophe et al. (2018)applied to frequency sets generated in presence of magnetic fields (profiles U, B+, B-).Notes.The table gives also the radial orders identified by the method.The first line corresponds to the real values without magnetic field.