Issue 
A&A
Volume 627, July 2019



Article Number  A64  
Number of page(s)  9  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201935462  
Published online  01 July 2019 
Period spacings of gravity modes in rapidly rotating magnetic stars
I. Axisymmetric fossil field with poloidal and toroidal components
^{1}
AIM, CEA, CNRS, Université ParisSaclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 GifsurYvette, France
email: vincent.prat@cea.fr
^{2}
LESIA, Observatoire de Paris, PSL University, CNRS, Sorbonne Université, Univ. Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, 92195 Meudon, France
^{3}
Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
^{4}
Dept. of Astrophysics, IMAPP, Radboud University Nijmegen, 6500 GL Nijmegen, The Netherlands
Received:
13
March
2019
Accepted:
14
May
2019
Context. Stellar magnetic fields are often invoked to explain the missing transport of angular momentum observed in models of stellar interiors. However, the properties of an internal magnetic field and the consequences of its presence on stellar evolution are largely unknown.
Aims. We study the effect of an axisymmetric internal magnetic field on the frequency of gravity modes in rapidly rotating stars to check whether gravity modes can be used to detect and probe such a field.
Methods. Rotation is taken into account using the traditional approximation of rotation and the effect of the magnetic field is computed using a perturbative approach. As a proof of concept, we compute frequency shifts due to a mixed (i.e. with both poloidal and toroidal components) fossil magnetic field for a representative model of a known magnetic, rapidly rotating, slowly pulsating Btype star: HD 43317.
Results. We find that frequency shifts induced by the magnetic field scale with the square of its amplitude. A magnetic field with a nearcore strength of the order of 150 kG (which is consistent with the observed surface field strength of the order of 1 kG) leads to signatures that are detectable in period spacings for highradialorder gravity modes.
Conclusions. The predicted frequency shifts can be used to constrain internal magnetic fields and offer the potential for a significant step forward in our interpretation of the observed structure of gravitymode period spacing patterns in rapidly rotating stars.
Key words: asteroseismology / waves / stars: magnetic field / stars: oscillations / stars: rotation
© V. Prat et al. 2019
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The development of helio and asteroseismology in recent decades has drastically improved our knowledge of stellar interiors (see e.g. Aerts et al. 2010, 2019; Chaplin & Miglio 2013; Hekker & ChristensenDalsgaard 2017). In particular, seismic constraints on the internal rotation of red giants and subgiants (Beck et al. 2012; Mosser et al. 2012; Deheuvels et al. 2012, 2014, 2015; Triana et al. 2017; Gehan et al. 2018) and intermediatemass mainsequence stars (Kurtz et al. 2014; Saio et al. 2015; Murphy et al. 2016; Van Reeth et al. 2016; Aerts et al. 2017; Ouazzani et al. 2017) highlight a weak coretosurface rotation contrast, which requires a much stronger transport of angular momentum than is predicted by current models (e.g. Marques et al. 2013; Ceillier et al. 2013; Cantiello et al. 2014; Spada et al. 2016; Ouazzani et al. 2019). For a few rapidly rotating, pulsating γ Doradus (γ Dor) stars, Van Reeth et al. (2018) also found very weak differential rotation, but the majority of the 37 analysed stars revealed uniform radial rotation profiles. Two physical processes are good candidates to explain the missing transport of angular momentum: internal gravity waves (Lee & Saio 1993; Zahn et al. 1997; Talon & Charbonnel 2005; Pantillon et al. 2007; Mathis 2009; Lee et al. 2014; Rogers 2015) and stable or unstable magnetic fields (Spruit 1999; Heger et al. 2005; Mathis & Zahn 2005; Fuller et al. 2019). We refer the reader to Aerts et al. (2019) for a thorough review.
The improvements in the instrument sensitivity and the increasing number of stars included in spectropolarimetric campaigns allow us to detect magnetic fields at the surface of stars using the Zeeman effect (see e.g. Donati & Landstreet 2009). For the massive O and Btype stars, highresolution spectropolarimetric surveys such as Magnetism in Massive Stars (MiMeS; Wade et al. 2016) and B fields in OB stars (BOB; Morel et al. 2014) have detected a largescale magnetic field at the surface of approximately 7% of the analysed stars. Most of the detected fields around intermediatemass and massive stars have a simple geometry, usually an inclined dipole, and a polar field strength between 300 G and a few tens of kG. These fields are probably of fossil origin, since their properties do not scale with stellar parameters or rotation (Mestel 1999; Neiner et al. 2015; Emeriau & Mathis 2015), in contrast with convective dynamo fields detected at the surface of lowmass stars. However, a convective dynamo may still be present in the core (e.g. Brun et al. 2005; Augustson et al. 2016).
Theoretical and numerical studies show that the stability of fossil magnetic fields requires that they extend deep within the stellar radiative envelope and have a mixed configuration with both poloidal and toroidal components (Tayler 1973, 1980; Markey & Tayler 1973; Braithwaite & Spruit 2004; Braithwaite 2007, 2009; Duez et al. 2010a). The internal magnetic field of a massive star with a convective core and radiative envelope, however, is still largely unknown because stellar interiors are opaque to spectropolarimetric measurements. Indirect information on the magnetic field can be obtained by constraining the transport of angular momentum and chemical species it may induce (e.g. Briquet et al. 2012). Another promising way to constrain the structure of the internal magnetic field is to look at its effect on the oscillation modes of stars, in particular lowfrequency gravity (g) modes (e.g. Buysschaert et al. 2018).
In nonrotating stars, the internal magnetic field induces the splitting of gmode frequencies of the same angular degree but different azimuthal orders (Ledoux & Simon 1957). This has been applied for a model representative of slowly pulsating Btype (SPB) stars with a purely poloidal dipolar field (e.g. Hasan et al. 2005).
In slowly rotating, weakly magnetic stars, rotation and the magnetic field can both be treated as perturbations. This has been mainly done for pressure (p) modes in the Sun and rapidly oscillating Ap (roAp) stars, with either an axisymmetric toroidal field (Gough & Taylor 1984) or an oblique (not aligned with the rotation axis) dipolar field (Dziembowski & Goode 1985; Gough & Thompson 1990; Goode & Thompson 1992; Shibahashi & Takata 1993). A similar treatment was applied to β Cephei by Shibahashi & Aerts (2000).
When the full effects of rotation need to be taken into account, that is for rapidly rotating stars (such as γ Dor, SPB, or Be stars) or very lowfrequency modes (such as Rossby modes, also called r modes), more complex formalisms have to be used. Schenk et al. (2002) proposed such a formalism to describe the coupling between modes of rotating stars due to external forces. Morsink & Rezania (2002) used this formalism to investigate the effect of a general, nonperturbative magnetic field on r modes.
The rotationally perturbed gmode pulsation frequencies in rapidly rotating stars are better described using the traditional approximation of rotation (TAR) than a perturbative approach. This approximation assumes that the radial component of the Coriolis force is negligible with respect to the buoyancy force and that radial displacements are limited by buoyancy and are small compared to horizontal displacements (Eckart 1960; Townsend 2003). Thus, the horizontal component of the rotation vector is neglected. The TAR also assumes that the star is spherically symmetric, which is reasonable only for moderately fast rotators because of centrifugal deformation. This allows for separation of variables in spherical coordinates and more efficient computation of eigenmodes than when considering the full effects of rotation. We refer the reader to Mathis et al. (2008) and Ouazzani et al. (2017), and references therein for a complete discussion of the TAR.
A standard tool used to interpret spectra of gravity modes is period spacing patterns, i.e. the morphology of period differences between modes of consecutive radial orders and the same angular degree and azimuthal order. In chemically homogeneous, nonrotating stars, we expect to have a constant asymptotic period spacing, as derived from the asymptotic relations by Tassoul (1980). However, in rotating stars, period spacing values seen in the inertial frame are not constant, and their slope as a function of the period is related to the rotation (Bouabid et al. 2013; Van Reeth et al. 2015a,b; Ouazzani et al. 2017). As an alternative approach to fitting the gradient of the period spacing pattern, Christophe et al. (2018) proposed a method of recovering uniform period spacings by stretching the pulsation periods.
Recently, Buysschaert et al. (2018) considered the magnetic, rapidly rotating SPB star HD 43317, following the earlier studies of Pápics et al. (2012), Briquet et al. (2013), and Buysschaert et al. (2017). They computed magnetic splittings due to an axisymmetric, purely poloidal dipolar field in the nonrotating case (following the perturbative formalism of Hasan et al. 2005) and found that they were negligible with respect to rotational splittings, thus justifying the perturbative approach. The authors combined photometric and spectroscopic time series to perform mode identification and modelled the star using a grid of nonrotating, nonmagnetic equilibrium stellar structure models. They found tentative evidence for a low amount of convective core overshooting, which was interpreted as being caused by the largescale magnetic field in HD 43317 (see also Briquet et al. 2012). However, as already mentioned above, such a purely poloidal magnetic field would be unstable. In addition, rapid rotation most likely influences the effect of the magnetic field on stellar oscillations. The goal of the present study is to develop a perturbative description of gravitoinertial modes (i.e. gravity modes affected by rotation) in the presence of a stable, mixed (i.e. with both poloidal and toroidal components), axisymmetric magnetic field within the TAR.
We present how the perturbation theory is modified in the presence of rotation in Sect. 2. We show how we derived the expression of the frequency shifts induced by an axisymmetric magnetic field in Sect. 3. Then, we apply our new formalism to a representative stellar structure model of HD 43317 found by Buysschaert et al. (2018) in Sect. 4. We discuss the results and prospects in Sect. 5.
2. Firstorder perturbation theory of rotating stars
In the uniformly rotating case, assuming a Lagrangian displacement of the form , where ω is the angular frequency in the corotating frame and t the time, the linearised momentum equation can be written as
where ℬ(ξ) = 2Ω ∧ ξ is the Coriolis operator, Ω = Ωe_{z} is the rotation vector, and 𝒞(ξ) describes forces that do not depend on the frequency. In the nonmagnetic case it reduces to pressure and buoyancy forces, 𝒞_{0}(ξ) = − ∇δP/ρ + δρ∇P/ρ^{2}, where δP and δρ are respectively the pressure and density Eulerian perturbations around the equilibrium values P and ρ. In the TAR, the horizontal component of the rotation vector is neglected, so that the Coriolis operator reduces to ℬ(ξ) = 2Ωcos θe_{r} ∧ ξ, where θ is the colatitude and e_{r} is the radial unit vector.
The firstorder perturbation by an additional force that does not depend on ω is obtained by writing ω = ω_{0} + εω_{1}, ξ = ξ_{0} + εξ_{1}, and 𝒞 = 𝒞_{0} + ε𝒞_{1}, where ε ≪ 1. At the zeroth order in ε, Eq. (1) leads to the unperturbed equation
and the firstorder equation reads
Given that iℬ and 𝒞_{0} are Hermitian, Eq. (2) allows us to simplify the scalar product of the previous equation with ξ_{0}, which leads to
where the scalar product is defined by
and the asterisk (^{*}) denotes the complex conjugate.
3. Magnetic frequency shifts
In the TAR, unperturbed eigenmodes for gravitoinertial waves are given in spherical coordinates (r, θ, φ) by
where m is the azimuthal order, and H_{r}, H_{θ}, and H_{φ} are radial, latitudinal, and azimuthal Hough functions, respectively (Hough 1898; Lee & Saio 1997; Townsend 2003). Their definitions are given in Appendix A.
In the present work we consider a magnetic field that is weak enough such that the effect of the unperturbed Lorentz force on the equilibrium state is negligible (e.g. Duez et al. 2010b). In addition, for simplicity we assume for our analytical calculations that the field is axisymmetric and dipolar:
For our numerical calculations (see Sect. 4), we use a field defined by
with ψ being the stream function, which satisfies
where e_{φ} is the unit vector in the azimuthal direction, R is the stellar radius, λ ≃ 35.89 is the smallest nonnegative constant such that B vanishes at the stellar surface, μ_{0} is the vacuum permeability, α is an amplitude scaling factor, and j_{1} (respectively y_{1}) is the firstorder spherical Bessel function of the first (respectively second) kind. Defining A such that ψ = A(r) sin ^{2}θ, we can write b_{r} = 2A/r^{2}, b_{θ} = −A′/r, and b_{φ} = λA/(rR), where the prime symbol (′) denotes a total radial derivative. One of the properties of this field is that the specific perturbed Lorentz force δF_{L}/ρ (see Eq. (11)) remains finite at the surface even though ρ becomes very small, so the perturbative approach can be used a priori.
This field aims to represent a fossil field in the stellar radiative zone, although its origin is not explicitly used further. Since purely poloidal and toroidal fields are known to be unstable (Tayler 1973; Markey & Tayler 1973), we use the mixed configuration proposed by Duez et al. (2010a), which is based on energy^{1} and stability arguments. This configuration has been demonstrated to be stable, using threedimensional simulations (Duez et al. 2010a), and to have properties similar to those obtained using numerical simulations of the formation of fossil fields (Braithwaite & Spruit 2004; Braithwaite & Nordlund 2006; Braithwaite 2008). Although this configuration has been derived in the nonrotating case, numerical simulations and theoretical calculations show that similar magnetic equilibria can be obtained in the rotating case (Duez 2011; Emeriau & Mathis 2015) even if their formation time can be modified (Braithwaite & Cantiello 2013). Generally, these equilibria are not axisymmetric, but rather oblique dipoles. However, as a first step, we restrict our study to an axisymmetric configuration, which is mathematically simpler. Here the field extends down to the centre, but in reality the field would be more complex due to its interaction with the dynamo field likely present in the convective core (Featherstone et al. 2009). For magnetic F and Atype stars, the fossil field might also interact with the dynamo in the subsurface convective envelope (Augustson et al. 2013).
The induction equation in the ideal magnetohydrodynamical approximation implies that the Eulerian perturbation to the magnetic field due to the oscillation displacement is
The perturbed Lorentz force reads
Equation (4) implies that magnetic frequency shifts δω = εω_{1} are proportional to ⟨ξ_{0}, δF_{L}/ρ⟩ = ⟨ξ_{0}, ε𝒞_{1}(ξ_{0})⟩. Therefore, the frequency shifts are proportional to the square of the magnetic field amplitude.
The computation of the magnetic frequency shifts involves a large number of terms. A significant fraction of them are zero for symmetry reasons. Hasan et al. (2005) considered only one term, which they claimed was dominant for highradialorder g modes, i.e. for modes with an absolute value of the radial order^{2} n much larger than the angular degree ℓ. Since the effect of rotation may change this, we choose to keep all nonzero terms. They are listed in Appendix B.
The first term of the denominator in Eq. (4) involves the scalar product
For highradialorder g modes, the first (radial) term is much smaller than the second (horizontal) one, and can thus be neglected, as done in Hasan et al. (2005). For lowradialorder modes (i.e. n < 5), though, this introduces significant errors. Note that in the model for HD 43317 (Buysschaert et al. 2018), some modes had such low radial orders.
The second term of the denominator in Eq. (4) yields
Since this term scales with the spin factor 2Ω/ω, it is negligible for lowradialorder modes, but not for highradialorder ones. When it is negligible, Eq. (4) implies
as in the nonrotating case, where Hasan et al. (2005) found that
This defines the splitting coefficient S_{c}, which is proportional to in their study, with
where ρ_{c} is the central density of the star. It follows that the effect of the magnetic field is stronger at lower frequencies, i.e. longer periods or higher radial orders. This is consistent with the fact that those frequencies are closer to the Alfvén frequency, which characterises the propagation of magnetic waves.
Highradialorder modes have a large radial wavenumber, which allows us to perform a Jeffreys–Wentzel–Kramers–Brillouin smallwavelength analysis (see e.g. Unno et al. 1989, in the context of stellar oscillations). When the poloidal component of the magnetic field is much larger than the toroidal component or of the same order of magnitude, which is the case here, the dominant term of the numerator in Eq. (4) is proportional to
This is consistent with the result obtained in the nonrotating case by Hasan et al. (2005). However, for subinertial (ω < 2Ω) gravitoinertial waves, which are often excited in rapidly rotating stars (e.g. Neiner et al. 2012; Moravveji et al. 2016; Saio et al. 2018), Hough functions differ significantly from spherical harmonics and our more general formalism is needed. These waves are trapped in an equatorial belt. In contrast, when the toroidal component is much larger than the poloidal component, the analysis predicts that eight other terms of the numerator in Eq. (4) could have a significant impact on the perturbation (see Appendix B.2).
4. Application to HD 43317: proof of concept
As a proof of concept, we now apply our new formalism to a representative stellar model of HD 43317, based on the analysis by Buysschaert et al. (2018), which is a rapidly rotating, magnetic Btype star exhibiting g modes. The photometric and spectroscopic analysis of this star by Pápics et al. (2012) indicated a solarlike metallicity Z = 0.014, an effective temperature T_{eff} = 17 350 ± 750 K, a surface gravity log g = 4.0 ± 0.1 dex, and a rotation period P_{rot} = 0.897673(4) d (about 60% of the Roche critical rotation rate). The spectropolarimetric analysis (Briquet et al. 2013; Buysschaert et al. 2017) showed the presence of a dipolar surface magnetic field of 1312 ± 332 G, and Buysschaert et al. (2018) estimated the obliquity angle as β = 81 ± 6°. The best model computed by Buysschaert et al. (2018) has a stellar mass M_{⋆} = 5.8 M_{⊙}, a central hydrogen mass fraction X_{c} = 0.54, and an exponential convective core overshooting parameter f_{ov} = 0.004. The stellar model was evaluated using a gridbased approach with the Modules for Experiments in Stellar Astrophysics (MESA) onedimensional stellar structure and evolution code (Paxton et al. 2018), and the eigenmodes were computed with the GYRE oscillation code (Townsend & Teitler 2013; Townsend et al. 2018).
The magnetic field obtained from Eq. (8) with the density profile of the MESA model is illustrated in Fig. 1. We note that the obliquity of the field with respect to the rotation axis is not accounted for in the present study. In the following, unless mentioned otherwise, we fix the maximum value for the internal magnetic field to B_{max} = 1.5 × 10^{5} G, which is approximately 100 times the value of the surface magnetic field detected by spectropolarimetry. This ratio is consistent with numerical simulations of magnetic equilibria (Braithwaite 2008), but our value of B_{0} (≃B_{max}) is significantly larger than the values (26.1 kG and 82.4 kG) used to compute magnetic splittings in Buysschaert et al. (2018).
Fig. 1. Representation of the magnetic field used in the numerical computations. The stream lines represent the poloidal component of the field (their line width is proportional to the amplitude). The background coloured contours represent the toroidal component. In this plot, the field is normalised so that the maximum of its norm is unity. Both poloidal and toroidal components have a maximum amplitude of order unity. This field is scaled using the amplitude parameter α mentioned in Eq. (9). 
We numerically computed the frequency perturbations of dipole (ℓ = 1) and quadrupole (ℓ = 2) modes. For nonaxisymmetric modes, it is important to distinguish between the frequency in the corotating frame ω and the frequency in the inertial frame ω_{i}, where
Buysschaert et al. (2018) incorrectly computed the magnetic frequency shifts of the identified g modes of HD 43317. A confusion concerning the numbering of GYRE output files led the authors to use incorrect eigenmodes in their computations. Our Table 1 gives a corrected version of their Table 3, which compares rotational splittings with magnetic splittings in the nonrotating case. We also add the magnetic frequency shifts in the TAR computed in the present paper. Although the magnetic frequency shifts in the TAR sometimes show significant differences from those computed in the nonrotating case, most of the time they are of the same order of magnitude, and remain small compared to the rotational shifts, which is consistent with our perturbative approach.
Comparison of rotational frequency shifts, magnetic frequency shifts in the nonrotating case, and magnetic frequency shifts in presence of rotation.
More generally, to check whether the perturbative approach is relevant, we compare the Alfvén frequency and the pulsation frequency ω. The Alfvén frequency is given by
where k is the wave vector, and must be much smaller than the pulsation frequency in the corotating frame. For highradialorder modes, B ⋅ k ≃ B_{r}k_{r} = b_{r}(r)k_{r} cos θ. Therefore, the magnitude of the Alfvén frequency is highest along the rotation axis, but highradialorder modes are subinertial, and thus trapped in the region where cos θ < ω/(2Ω) (see e.g. Prat et al. 2016). The criterion reads
To go further, we use the rough estimate k_{r} ∼ n/R. The quantity is highest at the centre of the star, but gravity waves are mostly sensitive to the nearcore region. In this region, , which finally leads to
This criterion is verified for all the modes computed in this work, which justifies the perturbative treatment a priori.
In the following sections we study more deeply the magnetic frequency shifts of zonal (Sect. 4.1), prograde (Sect. 4.2), and retrograde modes (Sect. 4.3). We investigate the effect of rotation in Sect. 4.4.
4.1. Zonal modes
Figure 2 represents for ℓ = 1 and m = 0 the period spacing ΔP between modes of consecutive radial orders as a function of the period P = 2π/ω_{i}. In a nonrotating model, the period spacings should be mostly constant with dips related to the chemical stratification near the core (Miglio et al. 2008). In rotating stars, the slope of the gmode period spacing pattern is related to the nearcore rotation rate (Bouabid et al. 2013; Van Reeth et al. 2016; Ouazzani et al. 2017). As can be seen in Fig. 2, the slope is not affected by the magnetic field, but the depth and the spacing of the dips show a significant modification at longer periods (i.e. low frequencies), where the effects of a magnetic field are stronger, as expected from Eq. (14). In addition, the magnetic field induces the presence of peaks close to the dips, thus forming a sawtooth pattern. A similar behaviour is observed for ℓ = 2. Figure 2 also illustrates the fact the magnetic shifts scale with the square of the field strength.
Fig. 2. Period spacings of g modes with radial orders from −1 (left) to −74 (right) as a function of the period in the inertial frame for ℓ = 1, m = 0, and different magnetic field strengths. The vertical black bar on the left represents a typical observational error bar of 250 s (Van Reeth et al. 2015b). 
According to Van Reeth et al. (2018), the effect of differential rotation on gravitoinertial modes is similar to a modification of the Brunt–Väisälä profile, which creates additional dips in the period spacings. Therefore, magnetic signatures are different from the signatures of radial differential rotation or stratification. This is a key property that may help to constrain the structure of internal stellar magnetic fields.
4.2. Prograde modes
Figure 3 represents the gmode period spacing of modes with ℓ = 1 and m = 1 and shows that lowfrequency modes are strongly affected by the magnetic field, similarly to zonal modes. The main difference is that period spacings of prograde modes tend to zero for high radial orders while their period remains finite, which makes them potentially harder to distinguish from each other.
4.3. Retrograde modes
For retrograde modes, Eq. (18) may lead to negative frequencies. Observationally, negative frequencies cannot be distinguished from positive ones. Therefore, retrograde modes are split into two series of modes. The corresponding period spacings are plotted in Fig. 4. Lowradialorder modes (which have short periods in the corotating frame) are only slightly affected by the magnetic field. In contrast, highradialorder modes (which have long periods in the corotating frame, but moderate ones in the inertial frame) show clear signatures of the magnetic field.
Fig. 4. Period spacings of g modes with radial orders from −1 (bottom left) to −26 (top right) and −44 (right) to −74 (bottom) as a function of the period in the inertial frame for ℓ = 1, m = −1, and B_{0} = 1.5 × 10^{5} G. We note that modes from n = −1 to n = −11 were observed for HD 43317. The vertical black bar (very small here) represents a typical observational error bar of 250 s. 
4.4. Effect of rotation
In this section, we present the computed modes and the associated magnetic frequency shifts with 50% and 150% of the measured rotation rate of HD 43317 (around 30% and 90% of the critical rotation rate, respectively). The obtained period spacings are plotted in Fig. 5. It is clear from this figure that increasing the rotation rate significantly decreases the amplitude of the expected magnetic signatures. This can be easily explained: when the Brunt–Väisälä frequency is much higher than the Coriolis frequency 2Ω, which is usually the case, the lower bound for the frequency of gravitoinertial waves ω_{−} is close to 2Ωcos θ (see e.g. Prat et al. 2016). The Alfvén frequency is usually lower than ω_{−}. When increasing the rotation rate, ω_{−} also increases and moves further away from the Alfvén frequency. Hence, gravitoinertial waves become less sensitive to the presence of the magnetic field with increasing rotation velocity. This is also consistent with the criterion derived in Eq. (21). It is therefore best to search for the signatures of a magnetic field in period spacing patterns of slow rotators.
Fig. 5. Same as Fig. 2, for B_{0} = 1.5 × 10^{5} G and two different rotation rates: Ω = 0.5Ω_{⋆} and Ω = 1.5Ω_{⋆} (around 30% and 90% of the critical rotation rate, respectively). 
For the slower rotation rate considered here, Fig. 5 displays a negative period spacing value, with a higherradialorder mode having a shorter period than a lowerradialorder mode. In reality, this could lead to an avoided crossing between consecutive modes (e.g. Lignières et al. 2006), which would require a nonperturbative treatment of the magnetic field. This would also make the detection of period spacing patterns significantly more difficult.
5. Discussion
In the present work we investigated the effect of a mixed, axisymmetric, internal largescale magnetic field, which presumably (but not necessarily) is of fossil origin, on the oscillation frequencies of gravitoinertial modes in the traditional approximation of rotation. The numerical application to a model of a SPB star showed that highradialorder modes exhibit significantly distinct gmode period spacing patterns because of the magnetic field compared to nonmagnetic stars. In particular, the magnetic field reveals itself by a sawtooth pattern in the period spacing morphology, rather than the typical dips that occur due to mode trapping caused by a μgradient left behind by the shrinking convective core (Miglio et al. 2008). Thus, the search for such predicted patterns is a way to discover internal magnetic fields from gravity modes. In addition, computations at different rotation rates highlighted that magnetic signatures decrease with rotation. As a consequence, it is crucial to take rotation into account when computing magnetic frequency shifts. Otherwise, the strength of the magnetic field needed to explain observed signatures would be drastically underestimated. It is noteworthy that the two Kepler SPB stars that were modelled in detail are both slow to moderate rotators. Their modelling revealed shortcomings in their frequency fits for the modes with higher periods (hence higher orders) (Moravveji et al. 2015, 2016). Our work provides a good reason to revisit the modelling of these stars, assuming the presence of a magnetic field.
The predicted magnetic signatures seem to be related (at least partly) to the dips in the period spacings created by chemical gradients. These gradients are a consequence of stellar evolution and are thus not present in stars on the zeroage main sequence (ZAMS). Recently, Mombarg et al. (2019) determined the asteroseismic ages of 37 γ Dor stars, finding several of those to be near the ZAMS, even though they show significant dips in their gmode period spacing patterns. Therefore, it would be interesting to determine whether the signatures of a largescale magnetic field are measurable for such stars. More generally, we intend to perform a systematic exploration of magnetic effects as a function of stellar parameters, such as mass, core mass, metallicity, and magnetic field strength in the nearcore region, where g modes are most sensitive.
In order to verify that the TAR is valid in the considered regime and that the magnetic field is weak enough to use the perturbation theory, it would be interesting to compare the obtained frequencies to twodimensional computations of modes, for example with the Twodimensional Oscillation Program (TOP; Reese et al. 2006) or the Adiabatic Code of Oscillation including Rotation (ACOR; Ouazzani et al. 2012). However, neither of these oscillation codes accounts for the magnetic field at the moment. Another application of these codes would be to apply the perturbative theory to gravitoinertial modes computed in centrifugally deformed stars.
In the case of HD 43317, magnetic signatures could not be extracted because too few mode frequencies were identified, especially for highradialorder modes, where the effect of the magnetic field is strongest. We note that this star was observed with CoRoT. In principle, stars observed with Kepler should allow many more modes to be detected. Another possible reason for the low number of identified frequencies is that only a fraction of computed modes may actually be excited. In addition, some modes may not be visible because of surface cancellation effects and the inclination of the star. However, some stars do exhibit gravity modes with radial orders down to −50 (see Van Reeth et al. 2015b for F stars and Pápics et al. 2017 for B stars). When a series of these modes is observed, our results show that magnetic signatures could be extracted for a field of a reasonable strength. Van Reeth et al. (2015b) give estimates of the error on the period spacings for γ Dor stars that range from roughly 20 to 1000 s (an average of 250 s was used in Figs. 2–5), which is in some cases smaller than the expected magnetic signatures.
Although our mixed (poloidal and toroidal) magnetic configuration is a significant improvement compared to previous studies, it is still relatively simple (confined axisymmetric dipole). First, a more realistic work would require considering a magnetic field that is not confined within the star. The formalism derived in this work is still valid a priori, but since the density becomes very small at the surface, surface effects are to be expected (see e.g. Bigot et al. 2000). Modelling these effects would require a nonperturbative treatment of the magnetic field (e.g. Morsink & Rezania 2002). A similar treatment would also be needed to investigate the effect of a very strong internal magnetic field. Second, many stars, including HD 43317, have a dipolar magnetic field inclined with respect to the rotation axis. Our formalism needs to be generalised to the case of a nonaxisymmetric field to predict the influence of the obliquity angle on the magnetic signatures, but this is beyond the scope of the present paper and is the subject of future work. Finally, general threedimensional nonaxisymmetric configurations should be studied. They will be of natural interest for fields generated by dynamos (e.g. Brun et al. 2005; Augustson et al. 2016; Brun & Browning 2017, and references therein), complex fossil fields (Donati et al. 2006; Braithwaite 2008; Kochukhov et al. 2011), and unstable fields (e.g. Braithwaite 2006, 2007; Brun 2007; Zahn et al. 2007). Seismic magnetic signatures could also be used to help detect and constrain magnetic fields of Vegalike stars (Lignières et al. 2009).
Another possible followup of this work is to generalise it to the case of differential rotation. A full treatment of differential rotation would probably make the present formalism unusable. However, nonaxisymmetric magnetic fields^{3} are known to inhibit differential rotation (Moss 1992; Spruit 1999). Thus, a perturbative treatment of differential rotation might be sufficient.
Broderick & Narayan (2008) and Duez & Mathis (2010) have demonstrated that dipolar fields should be favoured for fossil fields, as shown by the observations, because they are the equilibrium states with the lowest energies.
Acknowledgments
The research leading to these results received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreements No. 647383: SPIRE with PI S.M. and No. 670519: MAMSIE with PI C.A.). V.P. and S.M. acknowledge support from the CNES PLATO grant at CEA/DAp. V.P. acknowledges the International Space Science Institute (ISSI) for supporting the SoFAR international team. The authors thank the referee for the useful comments.
References
 Aerts, C., ChristensenDalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology, Astronomy and Astrophysics Library (Heidelberg: SpringerVerlag) [Google Scholar]
 Aerts, C., Van Reeth, T., & Tkachenko, A. 2017, ApJ, 847, L7 [NASA ADS] [CrossRef] [Google Scholar]
 Aerts, C., Mathis, S., & Rogers, T. 2019, ARA&A, 57 [Google Scholar]
 Augustson, K. C., Brun, A. S., & Toomre, J. 2013, ApJ, 777, 153 [Google Scholar]
 Augustson, K. C., Brun, A. S., & Toomre, J. 2016, ApJ, 829, 92 [NASA ADS] [CrossRef] [Google Scholar]
 Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Bigot, L., Provost, J., Berthomieu, G., Dziembowski, W. A., & Goode, P. R. 2000, A&A, 356, 218 [NASA ADS] [Google Scholar]
 Bouabid, M.P., Dupret, M.A., Salmon, S., et al. 2013, MNRAS, 429, 2500 [NASA ADS] [CrossRef] [Google Scholar]
 Braithwaite, J. 2006, A&A, 453, 687 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Braithwaite, J. 2007, A&A, 469, 275 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Braithwaite, J. 2008, MNRAS, 386, 1947 [NASA ADS] [CrossRef] [Google Scholar]
 Braithwaite, J. 2009, MNRAS, 397, 763 [NASA ADS] [CrossRef] [Google Scholar]
 Braithwaite, J., & Cantiello, M. 2013, MNRAS, 428, 2789 [NASA ADS] [CrossRef] [Google Scholar]
 Braithwaite, J., & Nordlund, Å. 2006, A&A, 450, 1077 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Braithwaite, J., & Spruit, H. C. 2004, Nature, 431, 819 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Briquet, M., Neiner, C., Aerts, C., et al. 2012, MNRAS, 427, 483 [NASA ADS] [CrossRef] [Google Scholar]
 Briquet, M., Neiner, C., Leroy, B., & Pápics, P. I., & MiMeS Collaboration 2013, A&A, 557, L16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Broderick, A. E., & Narayan, R. 2008, MNRAS, 383, 943 [Google Scholar]
 Brun, A. S. 2007, Astron. Nachr., 328, 1137 [NASA ADS] [CrossRef] [Google Scholar]
 Brun, A. S., & Browning, M. K. 2017, Liv. Rev. Sol. Phys., 14, 4 [Google Scholar]
 Brun, A. S., Browning, M. K., & Toomre, J. 2005, ApJ, 629, 461 [NASA ADS] [CrossRef] [Google Scholar]
 Buysschaert, B., Neiner, C., Briquet, M., & Aerts, C. 2017, A&A, 605, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Buysschaert, B., Aerts, C., Bowman, D. M., et al. 2018, A&A, 616, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cantiello, M., Mankovich, C., Bildsten, L., ChristensenDalsgaard, J., & Paxton, B. 2014, ApJ, 788, 93 [Google Scholar]
 Ceillier, T., Eggenberger, P., García, R. A., & Mathis, S. 2013, A&A, 555, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chaplin, W. J., & Miglio, A. 2013, ARA&A, 51, 353 [NASA ADS] [CrossRef] [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]
 Deheuvels, S., García, R. A., Chaplin, W. J., et al. 2012, ApJ, 756, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Deheuvels, S., Doğan, G., Goupil, M. J., et al. 2014, A&A, 564, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Ballot, J., Beck, P. G., et al. 2015, A&A, 580, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Donati, J.F., & Landstreet, J. D. 2009, ARA&A, 47, 333 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Donati, J.F., Howarth, I. D., Jardine, M. M., et al. 2006, MNRAS, 370, 629 [NASA ADS] [CrossRef] [Google Scholar]
 Duez, V. 2011, Astron. Nachr., 332, 983 [NASA ADS] [CrossRef] [Google Scholar]
 Duez, V., & Mathis, S. 2010, A&A, 517, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Duez, V., Braithwaite, J., & Mathis, S. 2010a, ApJ, 724, L34 [Google Scholar]
 Duez, V., Mathis, S., & TurckChièze, S. 2010b, MNRAS, 402, 271 [NASA ADS] [CrossRef] [Google Scholar]
 Dziembowski, W., & Goode, P. R. 1985, ApJ, 296, L27 [NASA ADS] [CrossRef] [Google Scholar]
 Eckart, C. 1960, Hydrodynamics of Oceans and Atmospheres (New York: Pergamon Press) [Google Scholar]
 Emeriau, C., & Mathis, S. 2015, in New Windows on Massive Stars, eds. G. Meynet, C. Georgy, J. Groh, & P. Stee, IAU Symp., 307, 373 [NASA ADS] [Google Scholar]
 Featherstone, N. A., Browning, M. K., Brun, A. S., & Toomre, J. 2009, ApJ, 705, 1000 [NASA ADS] [CrossRef] [Google Scholar]
 Fuller, J., Piro, A. L., & Jermyn, A. S. 2019, MNRAS, 485, 3661 [NASA ADS] [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., & Taylor, P. P. 1984, Mem. Soc. Astron. It., 55, 215 [NASA ADS] [Google Scholar]
 Gough, D. O., & Thompson, M. J. 1990, MNRAS, 242, 25 [Google Scholar]
 Hasan, S. S., Zahn, J.P., & ChristensenDalsgaard, J. 2005, A&A, 444, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 [NASA ADS] [CrossRef] [Google Scholar]
 Hekker, S., & ChristensenDalsgaard, J. 2017, A&ARv, 25, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hough, S. S. 1898, Phil. T. Roy. Soc. London A, 191, 139 [Google Scholar]
 Kochukhov, O., Lundin, A., Romanyuk, I., & Kudryavtsev, D. 2011, ApJ, 726, 24 [Google Scholar]
 Kurtz, D. W., Saio, H., Takata, M., et al. 2014, MNRAS, 444, 102 [NASA ADS] [CrossRef] [Google Scholar]
 Ledoux, P., & Simon, R. 1957, Ann. Astrophys., 20, 185 [NASA ADS] [Google Scholar]
 Lee, U., & Saio, H. 1993, MNRAS, 261, 415 [Google Scholar]
 Lee, U., & Saio, H. 1997, ApJ, 491, 839 [NASA ADS] [CrossRef] [Google Scholar]
 Lee, U., Neiner, C., & Mathis, S. 2014, MNRAS, 443, 1515 [NASA ADS] [CrossRef] [Google Scholar]
 Lignières, F., Rieutord, M., & Reese, D. 2006, A&A, 455, 607 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lignières, F., Petit, P., Böhm, T., & Aurière, M. 2009, A&A, 500, L41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Markey, P., & Tayler, R. J. 1973, MNRAS, 163, 77 [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., & Zahn, J.P. 2005, A&A, 440, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., Talon, S., Pantillon, F.P., & Zahn, J.P. 2008, Sol. Phys., 251, 101 [NASA ADS] [CrossRef] [Google Scholar]
 Mestel, L. 1999, in Stellar Magnetism, (Oxford: Oxford University Press), Int. Ser. Monographs Phys., 99 [Google Scholar]
 Miglio, A., Montalbán, J., Noels, A., & Eggenberger, P. 2008, MNRAS, 386, 1487 [NASA ADS] [CrossRef] [Google Scholar]
 Mombarg, J. S. G., Van Reeth, T., Pedersen, M. G., et al. 2019, MNRAS, 485, 3248 [NASA ADS] [CrossRef] [Google Scholar]
 Moravveji, E., Aerts, C., Pápics, P. I., Triana, S. A., & Vandoren, B. 2015, A&A, 580, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Moravveji, E., Townsend, R. H. D., Aerts, C., & Mathis, S. 2016, ApJ, 823, 130 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Morel, T., Castro, N., Fossati, L., et al. 2014, The Messenger, 157, 27 [NASA ADS] [Google Scholar]
 Morsink, S. M., & Rezania, V. 2002, ApJ, 574, 908 [NASA ADS] [CrossRef] [Google Scholar]
 Moss, D. 1992, MNRAS, 257, 593 [NASA ADS] [Google Scholar]
 Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012, A&A, 548, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Murphy, S. J., Fossati, L., Bedding, T. R., et al. 2016, MNRAS, 459, 1201 [NASA ADS] [CrossRef] [Google Scholar]
 Neiner, C., Floquet, M., Samadi, R., et al. 2012, A&A, 546, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Neiner, C., Mathis, S., Alecian, E., et al. 2015, in Polarimetry, eds. K. N. Nagendra, S. Bagnulo, R. Centeno, & M. Jesús Martínez González, IAU Symp., 305, 61 [Google Scholar]
 Ouazzani, R.M., Dupret, M.A., & Reese, D. R. 2012, A&A, 547, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ouazzani, R.M., Salmon, S. J. A. J., Antoci, V., et al. 2017, MNRAS, 465, 2294 [NASA ADS] [CrossRef] [Google Scholar]
 Ouazzani, R.M., Marques, J. P., Goupil, M., et al. 2019, A&A, 626, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pantillon, F. P., Talon, S., & Charbonnel, C. 2007, A&A, 474, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pápics, P. I., Briquet, M., Baglin, A., et al. 2012, A&A, 542, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pápics, P. I., Tkachenko, A., Van Reeth, T., et al. 2017, A&A, 598, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Paxton, B., Schwab, J., Bauer, E. B., et al. 2018, ApJS, 234, 34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Prat, V., Lignières, F., & Ballot, J. 2016, A&A, 587, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reese, D., Lignières, F., & Rieutord, M. 2006, A&A, 455, 621 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rogers, T. M. 2015, ApJ, 815, L30 [NASA ADS] [CrossRef] [Google Scholar]
 Saio, H., Kurtz, D. W., Takata, M., et al. 2015, MNRAS, 447, 3264 [NASA ADS] [CrossRef] [Google Scholar]
 Saio, H., Kurtz, D. W., Murphy, S. J., Antoci, V. L., & Lee, U. 2018, MNRAS, 474, 2774 [Google Scholar]
 Schenk, A. K., Arras, P., Flanagan, É. É., Teukolsky, S. A., & Wasserman, I. 2002, Phys. Rev. D, 65, 024001 [NASA ADS] [CrossRef] [Google Scholar]
 Shibahashi, H., & Aerts, C. 2000, ApJ, 531, L143 [NASA ADS] [CrossRef] [Google Scholar]
 Shibahashi, H., & Takata, M. 1993, PASJ, 45, 617 [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. 1999, A&A, 349, 189 [NASA ADS] [Google Scholar]
 Talon, S., & Charbonnel, C. 2005, A&A, 440, 981 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tassoul, M. 1980, ApJS, 43, 469 [NASA ADS] [CrossRef] [Google Scholar]
 Tayler, R. J. 1973, MNRAS, 161, 365 [NASA ADS] [CrossRef] [Google Scholar]
 Tayler, R. J. 1980, MNRAS, 191, 151 [NASA ADS] [Google Scholar]
 Townsend, R. H. D. 2003, MNRAS, 340, 1020 [NASA ADS] [CrossRef] [Google Scholar]
 Townsend, R. H. D., & Teitler, S. A. 2013, MNRAS, 435, 3406 [NASA ADS] [CrossRef] [Google Scholar]
 Townsend, R. H. D., Goldstein, J., & Zweibel, E. G. 2018, MNRAS, 475, 879 [Google Scholar]
 Triana, S. A., Corsaro, E., De Ridder, J., et al. 2017, A&A, 602, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial Oscillations of Stars (Tokyo: University of Tokyo Press) [Google Scholar]
 Van Reeth, T., Tkachenko, A., Aerts, C., et al. 2015a, A&A, 574, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Van Reeth, T., Tkachenko, A., Aerts, C., et al. 2015b, ApJS, 218, 27 [NASA ADS] [CrossRef] [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]
 Wade, G. A., Neiner, C., Alecian, E., et al. 2016, MNRAS, 456, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, H., Boyd, J. P., & Akmaev, R. A. 2016, Geosci. Model Dev., 9, 1477 [NASA ADS] [CrossRef] [Google Scholar]
 Zahn, J.P., Talon, S., & Matias, J. 1997, A&A, 322, 320 [Google Scholar]
 Zahn, J.P., Brun, A. S., & Mathis, S. 2007, A&A, 474, 145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Hough functions
The radial Hough function (Hough 1898) is defined by H_{r}(θ) = f(μ = cos θ), where f is the solution of the eigenvalue equation
and ν = 2Ω/ω is the spin factor. In the nonrotating case, the eigenvalue λ reduces to −ℓ(ℓ+1), where ℓ is the angular degree of the mode, and H_{r} reduces to the classical associated Legendre polynomial .
The latitudinal and azimuthal Hough functions are respectively given by
and
where ′ denotes here a total latitudinal derivative. It is obvious from these definitions that H_{φ} has the same parity as H_{r} with respect to θ = π/2, whereas H_{θ} has the opposite parity.
In the present work, we computed Hough functions using the implementation based on Chebyshev polynomials proposed by Wang et al. (2016).
Appendix B: Nonzeroaverage terms of the Lorentz work
In this section, we list all the nonzeroaverage terms of the work of the Lorentz force δF_{L} ⋅ ξ^{*} (for convenience, we take here the product with μ_{0}). Interestingly, all these terms involve either only poloidal components of the magnetic field, or only the toroidal component, and we group them accordingly. In addition, all terms have a purely radial part multiplied by a purely latitudinal part, and the prime symbol (′) denotes a total derivative, either radial or latitudinal depending on the context.
B.1. Poloidal terms
Defining A = [(rb_{θ})′ + b_{r}], the terms involving poloidal components are
B.2. Toroidal terms
The terms that involve the toroidal component are
The first eight terms in the previous equation are those predicted to have a significant impact on frequency shifts for a strong toroidal field.
All Tables
Comparison of rotational frequency shifts, magnetic frequency shifts in the nonrotating case, and magnetic frequency shifts in presence of rotation.
All Figures
Fig. 1. Representation of the magnetic field used in the numerical computations. The stream lines represent the poloidal component of the field (their line width is proportional to the amplitude). The background coloured contours represent the toroidal component. In this plot, the field is normalised so that the maximum of its norm is unity. Both poloidal and toroidal components have a maximum amplitude of order unity. This field is scaled using the amplitude parameter α mentioned in Eq. (9). 

In the text 
Fig. 2. Period spacings of g modes with radial orders from −1 (left) to −74 (right) as a function of the period in the inertial frame for ℓ = 1, m = 0, and different magnetic field strengths. The vertical black bar on the left represents a typical observational error bar of 250 s (Van Reeth et al. 2015b). 

In the text 
Fig. 3. Same as Fig. 2, but for ℓ = 1, m = 1, and B_{0} = 1.5 × 10^{5} G. 

In the text 
Fig. 4. Period spacings of g modes with radial orders from −1 (bottom left) to −26 (top right) and −44 (right) to −74 (bottom) as a function of the period in the inertial frame for ℓ = 1, m = −1, and B_{0} = 1.5 × 10^{5} G. We note that modes from n = −1 to n = −11 were observed for HD 43317. The vertical black bar (very small here) represents a typical observational error bar of 250 s. 

In the text 
Fig. 5. Same as Fig. 2, for B_{0} = 1.5 × 10^{5} G and two different rotation rates: Ω = 0.5Ω_{⋆} and Ω = 1.5Ω_{⋆} (around 30% and 90% of the critical rotation rate, respectively). 

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.