<figure><img src="image_1.png" /></figure>

HAL
open science

# The traditional approximation of rotation for rapidly rotating stars and planets

H. Dhouib, V. Prat, T. van Reeth, S. Mathis

► To cite this version:

H. Dhouib, V. Prat, T. van Reeth, S. Mathis. The traditional approximation of rotation for rapidly rotating stars and planets: II. Deformation and differential rotation. Astronomy and Astrophysics - A&A, 2021, 656, pp.A122. 10.1051/0004-6361/20211152 . cea-03475470

HAL Id: cea-03475470

https://cea.hal.science/cea-03475470

Submitted on 10 Dec 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
# The traditional approximation of rotation for rapidly rotating stars and planets

## II. Deformation and differential rotation

H. Dhouib<fnref n="1" />, V. Prat<fnref n="1" />, T. Van Reeth<fnref n="2" />, and S. Mathis<fnref n="1" />

Received 22 April 2021 / Accepted 5 October 2021

## ABSTRACT

*Context.* We examine the dynamics of low-frequency gravito-inertial waves (GIWs) in differentially rotating deformed radiation zones in stars and planets by generalising the traditional approximation of rotation (TAR). The TAR treatment was built on the assumptions that the star is spherical (i.e. its centrifugal deformation is neglected) and uniformly rotating. However, it has been generalised in our previous work by including the effects of the centrifugal deformation using a non-perturbative approach. In the meantime, TAR has been generalised in spherical geometry to take the differential rotation into account.

*Aims.* We aim to carry out a new generalisation of the TAR treatment to account for the differential rotation and the strong centrifugal deformation simultaneously.

*Methods.* We generalise our previous work by taking into account the differential rotation in the derivation of our complete analytical formalism that allows the study of the dynamics of GIWs in differentially and rapidly rotating stars.

*Results.* We derived the complete set of equations that generalises the TAR, simultaneously taking the full centrifugal acceleration and the differential rotation into account. Within the validity domain of the TAR, we derived a generalised Laplace tidal equation for the horizontal eigenfunctions and asymptotic wave periods of the GIWs, which can be used to probe the structure and dynamics of differentially rotating deformed stars with asteroseismology.

*Conclusions.* A new generalisation of the TAR, which simultaneously takes into account the differential rotation and the centrifugal acceleration in a non-perturbative way, was derived. This generalisation allowed us to study the detectability and the signature of the differential rotation on GIWs in rapidly rotating deformed stars and planets. We found that the effects of the differential rotation in early-type deformed stars on GIWs is theoretically largely detectable in modern space photometry using observations from *Kepler* and *TESS*.

**Key words.** hydrodynamics – waves – stars: rotation – stars: oscillations – methods: analytical – methods: numerical

## 1. Introduction

Understanding how angular momentum and chemicals are transported in the interiors of stars (and planets) along their evolution is one of the key challenges of modern stellar (and planetary) astrophysics. Indeed, rotation modifies their structure, their chemical stratification, their internal flows and magnetism, and their mass losses and winds (e.g. Maeder 2009; Mathis et al. 2013; Aerts et al. 2019, and references therein). In this quest, asteroseismology has bought a fundamental breakthrough by demonstrating that all stars are the seat of a strong extraction of angular momentum during their evolution in comparison with the predictions by stellar models taking the rotation into account following the standard rotational transport and mixing theory (Eggenberger et al. 2012; Marques et al. 2013; Ceillier et al. 2013; Cantiello et al. 2014; Ouazzani et al. 2019). This was first obtained thanks to mixed pulsation modes splitted by rotation propagating in evolved low- and intermediate-mass stars (Beck et al. 2012, 2014, 2018; Mosser et al. 2012; Deheuvels et al. 2012, 2014, 2015, 2020; Di Mauro et al. 2016; Triana et al. 2017; Gehan et al. 2018; Tayar et al. 2019). Then, observations of oscillation

modes in F- and A-type stars (Kurtz et al. 2014; Saio et al. 2015; Bedding et al. 2015; Keen et al. 2015; Van Reeth et al. 2015, 2016, 2018; Schmid & Aerts 2016; Murphy et al. 2016; Sowicka et al. 2017; Guo et al. 2017; Saio et al. 2018, 2021; Mombarg et al. 2019; Li et al. 2019, 2020; Ouazzani et al. 2020) and in B-type stars (Pápics et al. 2015, 2017; Triana et al. 2015; Moravveji et al. 2016; Kallinger et al. 2017; Buysschaert et al. 2018; Szewczuk & Daszyńska-Daszkiewicz 2018; Pedersen et al. 2021; Szewczuk et al. 2021) provided us new Rosetta stones to constrain the transport of angular momentum in the whole Hertzsprung-Russell diagram. More particularly, this pushes gravity- and gravito-inertial mode pulsators such as γ-Doradus and SPB stars at the forefront of this research. For instance, recent theoretical developments have demonstrated how it is possible to probe stellar internal rotation in γ-Doradus stars from their surface to their convective core (Ouazzani et al. 2020; Saio et al. 2021). These stars are rapid rotators for a large proportion of them. Therefore, it is necessary to study gravito-inertial modes. These modes are gravity modes, which propagate only in stably stratified stellar radiation zones under the action of the restoring buoyancy force in the absence of rotation, which are modified by rotation. If their

<footnote n="1">¹ Département d'Astrophysique-AIM, CEA/DRF/IRFU, CNRS/INSU, Université Paris-Saclay, Université Paris-Diderot, Université de Paris, 91191 Gif-sur-Yvette, France
e-mail: hachem.dhouib@cea.fr</footnote>

<footnote n="2">² Institute of Astronomy, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium</footnote>
frequency is super-inertial (i.e. above the inertial frequency $2\Omega$, $\Omega$ being the stellar angular velocity), they are propagating in stellar radiation zones and evanescent in convective regions. If their frequency is sub-inertial (below $2\Omega$) they propagate in an equatorial belt in radiation zones and they become propagative inertial waves in convective regions (e.g. Dintrans & Rieutord 2000; Mathis et al. 2014). The challenge of studying these waves is that the equation describing their dynamics are intrinsically bidimensional and non-separable (Dintrans et al. 1999; Prat et al. 2016, 2018; Mirouh et al. 2016). This makes the development of seismic diagnosis difficult analytically (Prat et al. 2017) or expansive in computation time when using 2D oscillation and stellar structure codes (e.g. Ouazzani et al. 2017; Reese et al. 2021) in the general case.

In this framework, the traditional approximation of rotation (TAR), which has been first introduced in geophysics (Eckart 1960) to treat the propagation of gravito-inertial waves (GIWs) in the case where the Coriolis acceleration can be neglected in front of the buoyancy force in the direction of the entropy and chemical stable stratifications, is very useful. Indeed, it allows one to consider that the velocity of gravito-inertial waves are mostly horizontal and to neglect the latitudinal component of the rotation vector in the momentum equation. That leads to decoupling the vertical and horizontal directions by keeping the non-perturbative action of the Coriolis acceleration in the latitudinal and azimuthal directions and neglecting it along the vertical one. Propagation equations become separable as in the non-rotating case (e.g. Bildsten et al. 1996; Lee & Saio 1997). In addition, the derivation of powerful and flexible seismic diagnostics using the period spacing between consecutive high radial order gravito-inertial modes becomes possible in uniformly rotating spherical stars (Bouabid et al. 2013). These period spacing provide constrains on properties of the chemical stratification and the rotation rate (through their slope) near the interface between the convective core and the radiative envelope in rapidly rotating intermediate-mass $\gamma$-Doradus (e.g. Van Reeth et al. 2015, 2016, 2018; Ouazzani et al. 2017; Christophe et al. 2018; Li et al. 2019, 2020; Saio et al. 2021) and SPB stars (e.g. Pápics et al. 2015, 2017; Moravveji et al. 2016; Szewczuk & Daszyńska-Daszkiewicz 2018; Pedersen et al. 2021), thanks to the high precision of space-based photometric observations (e.g. Aerts et al. 2019; Aerts 2021, and references therein). This allows us to build intensive forward modelling of g-mode pulsating stars to constrain their internal structural and dynamical properties.

This approximation, in its standard version, is only applicable for low-frequency GIWs propagating in strongly stably stratified zones of uniformly rotating spherical stars. In this case a set of assumptions should be verified: (i) the buoyancy force is stronger than the Coriolis acceleration (i.e. $2\Omega \ll N$, where N is the Brunt-Väisälä frequency) in the direction of stable entropy or chemical stratification, (ii) the Brunt-Väisälä frequency is larger than the frequency of the waves in the rotating frame $\omega$ ($\omega \ll N$), (iii) the rotation is assumed to be uniform and (iv) the star is assumed to be spherical, in other words, the centrifugal deformation of the star is neglected (i.e. $\Omega \ll \Omega_K \equiv \sqrt{GM/R^3}$, where $\Omega_K$ is the Keplerian critical (breakup) angular velocity, and G, M, and R are the universal constant of gravity, the mass of the star, and the stellar radius, respectively).

Recently, the TAR has been examined in deformed stars. First, Mathis & Prat (2019) have included the centrifugal acceleration for slightly deformed stars using a first-order perturbative approach. Then, this perturbative framework has been optimised for practical applications to one-dimensional stellar

structure models by Henneco et al. (2021). Second, in Dhouib et al. (2021, hereafter Paper I), we have used a non-perturbative approach to include the centrifugal acceleration for strongly deformed stars and planets by studying the dynamics of low-frequency GIWs in a general spheroidal coordinate system defined by Bonazzola et al. (1998) which follows the shape of a deformed star. It is important to note that in these studies the rotation was assumed to be uniform. These semi-analytical studies demonstrated that centrifugal effects' signatures are potentially detectable. Moreover, their results are in good qualitative agreement with direct computations of gravito-inertial modes in centrifugally deformed stellar models using 2D oscillation codes that take into account the full Coriolis and centrifugal accelerations such as the Toulouse Oscillation Program (TOP) (Ballot et al. 2010, 2012) and the adiabatic code of oscillation including rotation (ACOR) code (Ouazzani et al. 2012, 2015, 2017). In the case where the centrifugal acceleration is treated as a perturbation, but where the full Coriolis acceleration was taken into account, we can refer to the Tohoku oscillation code (Lee & Saio 1987; Lee & Baraffe 1995). The chosen coordinate system is used in ACOR, TOP and in the Evolution STEllaire en Rotation (ESTER) code (Espinosa Lara & Rieutord 2013) that computes the structure and the stationary internal hydrodynamics of (differentially-)rotating early-type stars such as the observed gravito-inertial modes pulsators we are studying.

However, the solid-body rotation assumption has to be potentially abandoned along the evolution of real stars (and planets) where gradients of angular velocity can develop, both in the radial and in the latitudinal directions. They can be triggered by the redistribution of angular momentum induced by stellar winds (e.g. Zahn 1992), structural adjustments (Maeder & Meynet 2000; Decressin et al. 2009), and by transport processes (Mathis et al. 2013). First, as shown in Charbonnel et al. (2013) and Gallet & Bouvier (2015), early main sequence low-mass stars are potentially subject to a strong differential rotation. Moreover, Aerts et al. (2019) and references therein (we refer the reader to those provided in the first paragraph for F, A, and B-type stars) pointed out that the radiative envelope of intermediate-mass stars are the seat of a weak differential rotation for those observed with a small difference of rotation between the stellar surface and the core. This is also predicted by numerical simulations and models of different transport mechanisms like internal gravity waves (Rogers et al. 2013; Rogers 2015) or meridional flows (Maeder & Meynet 2000; Decressin et al. 2009; Rieutord et al. 2016; Ouazzani et al. 2019). In this framework, Ogilvie & Lin (2004) and Mathis (2009) have generalised the TAR by including the effects of general differential rotation on low-frequency GIWs propagating in spherical stars (and planets). This new formalism has been successfully applied to gravito-inertial mode pulsators observed by the Kepler (Borucki et al. 2009) space telescope. Indeed, Van Reeth et al. (2018) have applied it to derive the variation of the asymptotic period spacing in the case of a weak radial differential rotation. They evaluated the sensitivity of GIWs period spacing to the effect of such a radial shear in spherical stars and they demonstrated that they can be potentially detected if observing several modes.

Since both effects of the centrifugal deformation and differential rotation are potentially observable, we generalise here these previous works, treating the case of a general differential rotation in deformed stars (and planets). We begin in Sect. 2 with the derivation of the system of linearised hydrodynamic equations for GIWs in spheroidal coordinates which follows the shape of a differentially rotating deformed body. We choose again the coordinate system that was introduced by
Bonazzola et al. (1998) which is used in 2D numerical models of rotating stars and pulsation codes. In Sect. 3, we build the generalised TAR in this set of equations by adopting the adequate assumptions. In Sect. 4, we rewrite the oscillation equations in the form of a generalised Laplace tidal equation and deduce the asymptotic frequencies of low-frequency GIWs. In Sect. 5, we carry out a numerical exploration of the eigenvalues and Hough eigenfunctions of the generalised Laplace tidal equation within the domain of validity of the TAR. We model early-type stars with 2D ESTER models (Espinosa Lara & Rieutord 2013). In Sect. 6, we quantify the differential rotation and the centrifugal deformation combined effects on the asymptotic period spacing pattern and we discuss their detectability. Finally, we discuss and summarise our work and results in Sect. 8.

## 2. Hydrodynamic equations in differentially rotating deformed stars

### 2.1. Spheroidal geometry

As in Paper I, here, we use the spheroidal coordinate system ($\zeta, \theta, \varphi$) proposed by Bonazzola et al. (1998), where $\zeta$ is the pseudo-radial coordinate, $\theta$ the colatitude and $\varphi$ the azimuth. Following Espinosa Lara & Rieutord (2013), this new coordinate system, illustrated in Fig. 1, can be linked to the usual spherical one ($r, \theta, \varphi$) via the following mapping

$$r(\zeta, \theta) = a_i \xi \Delta \eta_i + R_i(\theta) + A_i(\xi) (\Delta R_i(\theta) - a_i \Delta \eta_i), \quad \eta_i \le \zeta \le \eta_{i+1}, \quad (1)$$

in the ith subdomain $D_{i \in [0,n-1]} \in [R_i(\theta), R_{i+1}(\theta)]$ of the spheroidal domain $D$ where $R_{i \in [0,n]}(\theta)$ are series of functions, such that $R_n(\theta) = R_s(\theta)$ is the outer boundary and $R_0(\theta) = 0$ is the centre. Additionally, $\eta_i = R_i(\theta = 0)$ are the polar radii of the interfaces between the subdomains, $\Delta\eta_i = \eta_{i+1} - \eta_i$, $\Delta R_i(\theta) = R_{i+1}(\theta) - R_i(\theta)$ and $\xi = (\zeta - \eta_i)/\Delta\eta_i$. The functions $A_{i \in [1,n-1]}(\xi) = -2\xi^3 + 3\xi^2$ and $A_0(\xi) = -1.5\xi^5 + 2.5\xi^3$ and the constants $a_i = 1$ are chosen to satisfy the boundary conditions between the different subdomains. We recall also the spheroidal base ($a_\zeta, a_\theta, a_\varphi$) defined using the mapping (Eq. (1)) (Rieutord et al. 2005; Lignières et al. 2006; Reese et al. 2006)

$$\left\{ \begin{array}{l} a_{\zeta} = \frac{\zeta^2}{r^2} e_r \\ a_{\theta} = \frac{\zeta}{r^2 r_{\zeta}} (r_{\theta} e_r + r e_{\theta}), \\ a_{\varphi} = \frac{\zeta}{r r_{\zeta}} e_{\varphi} \end{array} \right. \quad (2)$$

where $r_\zeta \equiv \partial_\zeta r$, $r_\theta \equiv \partial_\theta r$, ($e_r, e_\theta, e_\varphi$) is the usual spherical base, $\epsilon = 1 - R_{pol}/R_{eq}$ is the flatness, $R_s(\theta)$, $R_{eq}$ and $R_{pol}$ are the surface, the equatorial and the polar radii, respectively.

### 2.2. Linearised hydrodynamic equations in spheroidal coordinates

To treat the wave dynamics in differentially rotating, strongly deformed stars and planets, we derive the complete adiabatic inviscid system of equations in spheroidal coordinates. First, the linearised momentum equation for an inviscid fluid is written as

$$(\partial_t + \Omega \partial_\varphi) \left[ \frac{\zeta^2 r_\zeta}{r^2} v^\zeta + \frac{\zeta r_\theta}{r^2} v^\theta \right] = 2\Omega \frac{\zeta \sin \theta}{r} v^\varphi - \frac{1}{\rho_0} \partial_\zeta P' + \frac{\rho'}{\rho_0^2} \partial_\zeta P_0 - \partial_\zeta \Phi', \quad (3)$$

<figure><img src="image_2.png" /><figcaption>Fig. 1. Spheroidal coordinate system with $n = 13$ subdomains used to compute the equilibrium model of a star rotating at 60% of its break-up velocity corresponding to the case of the 3 $M_\odot$ ESTER model (with the central fraction in hydrogen $X_c = 0.7$) studied in Sect. 5. $\zeta$ and $\theta$ are the pseudo-radius and the colatitude, respectively. $\eta_i$ and $\Delta R_i(\theta)$ are defined in our mapping (Eq. (1)).</figcaption></figure>

$$(\partial_t + \Omega \partial_\varphi) \left[ \frac{\zeta^2 r_\theta}{r^2} v^\zeta + \frac{\zeta (r^2 + r_\theta^2)}{r^2 r_\zeta} v^\theta \right] = 2\Omega \frac{\zeta (r_\theta \sin \theta + r \cos \theta)}{r r_\zeta} v^\varphi - \frac{1}{\rho_0} \partial_\theta P' + \frac{\rho'}{\rho_0^2} \partial_\theta P_0 - \partial_\theta \Phi', \quad (4)$$

$$\frac{\zeta}{r_\zeta}(\partial_t + \Omega \partial_\varphi)v^\varphi = -2\Omega \frac{\zeta^2 \sin \theta}{r} v^\zeta - 2\Omega \frac{\zeta (r_\theta \sin \theta + r \cos \theta)}{r r_\zeta} v^\theta \\ - \frac{1}{\rho_0 \sin \theta} \partial_\varphi P' - \frac{\zeta \sin \theta}{r_\zeta} (\zeta \partial_\zeta \Omega u^\zeta + \partial_\theta \Omega u^\theta) \\ - \frac{1}{\sin \theta} \partial_\varphi \Phi', \quad (5)$$

where $v^\zeta, v^\theta$, and $v^\varphi$ are the contravariant components of the velocity field $v = v^i a_i$ and $\Omega = \Omega(\zeta, \theta)(\cos \theta e_r - \sin \theta e_\theta)$ is the differential angular velocity of the star. $\rho, \Phi$ and $P$ are the fluid density, gravitational potential and pressure, respectively. Each of these scalar quantities has been expanded as

$$X(r, \theta, \varphi, t) = X_0(r, \theta) + X'(r, \theta, \varphi, t),$$

where $X_0$ is the hydrostatic component of $X$ and $X'$ the wave's associated linear fluctuation. Following Cowling (1941), we neglect the fluctuation of the gravitational potential. Next, the linearised continuity equation is obtained

$$(\partial_t + \Omega \partial_\varphi) \rho' = - \frac{\zeta^2 \partial_\zeta \rho_0}{r^2 r_\zeta} v^\zeta - \frac{\zeta \partial_\theta \rho_0}{r^2 r_\zeta} v^\theta \\ - \frac{\zeta^2 \rho_0}{r^2 r_\zeta} \left[ \frac{\partial_\zeta (\zeta^2 v^\zeta)}{\zeta^2} + \frac{\partial_\theta (\sin \theta v^\theta)}{\zeta \sin \theta} + \frac{\partial_\varphi v^\varphi}{\zeta \sin \theta} \right]. \quad (6)$$

Then, the linearised energy equation in the adiabatic limit is derived

$$(\partial_t + \Omega \partial_\varphi) \left( \frac{1}{\Gamma_1} \frac{P'}{P_0} - \frac{\rho'}{\rho_0} \right) = \frac{N^2}{\|g_{\text{eff}}\|^2} \mathbf{v} \cdot \mathbf{g}_{\text{eff}}, \quad (7)$$
where $\Gamma_1 = (\partial \ln P_0 / \partial \ln \rho_0)_S$ ($S$ being the macroscopic entropy) is the adiabatic exponent, $g_{\text{eff}} = -\nabla\Phi_0 + \frac{1}{2}\Omega^2\nabla(r^2\sin^2\theta) = \nabla P_0 / \rho_0$ is the background effective gravity which includes the centrifugal acceleration and $N^2$ the squared Brunt–Väisälä frequency (or buoyancy frequency) given by

$$N^2(\zeta, \theta) = g_{\text{eff}} \cdot \left( -\frac{1}{\Gamma_1} \frac{\nabla P_0}{P_0} + \frac{\nabla \rho_0}{\rho_0} \right). \qquad (8)$$

Finally, we expand the wave's velocity and fluctuations on Fourier series both in time and in azimuth

$$v(\zeta, \theta, \varphi, t) \equiv \sum_{\omega^{\text{in}}, m} \{u(\zeta, \theta) \exp[i(\omega^{\text{in}}t - m\varphi)]\}, \qquad (9)$$

$$X'(\zeta, \theta, \varphi, t) \equiv \sum_{\omega^{\text{in}}, m} \{\bar{X}(\zeta, \theta) \exp[i(\omega^{\text{in}}t - m\varphi)]\}, \qquad (10)$$

where $m$ is the azimuthal order and $\omega^{\text{in}}$ is the wave frequency in an inertial reference frame. In a differentially rotating region, the waves are Doppler-shifted due to the differential rotation so we can define the wave frequency $\omega$ in the rotating reference frame as

$$\omega(\zeta, \theta) = \omega^{\text{in}} - m\Omega(\zeta, \theta). \qquad (11)$$

### 3. Generalised TAR

#### 3.1. Approximation on the stratification profile: $N^2(\zeta, \theta) \approx N^2(\zeta)$

As in Paper I, in order to obtain a separable system of equations when applying the TAR, we assume here that the Brunt–Väisälä frequency $N^2$ depends mainly on the pseudo-radius $\zeta$. This implies also that the hydrostatic pressure $P_0$ and the hydrostatic density $\rho_0$ depend mainly on $\zeta$. We present the validity domain of this approximation in Sect. 5.1 using two-dimensional ESTER stellar models.

#### 3.2. The TAR with centrifugal acceleration in differentially rotation stars

By assuming the following frequencies hierarchy imposed by the TAR: $\omega \ll N$ and $2\Omega \ll N$ (we verify this hierarchy in Sect. 7) which leads to a mostly horizontal velocity equation because of the strong buoyancy restoring force action in the vertical direction ($|u^\zeta| \ll |u^\theta|$, $|u^\varphi|$), the approximation $N^2(\zeta, \theta) \approx N^2(\zeta)$, the Cowling approximation ($\Phi' = 0$), and the anelastic approximation (acoustic waves are filtered out), we can rewrite accordingly the energy equation (Eq. (7)) and the three components of the momentum equation (Eqs. (3)–(5)). First, we rewrite the energy equation as

$$\frac{\tilde{\rho}}{\rho_0} = \frac{1}{i\omega} \frac{N^2}{g_{\text{eff}}} \frac{\zeta^2}{r^2 r_{\zeta}} u^{\zeta}, \qquad (12)$$

where we neglect the term $(1/\Gamma_1)P'/P_0$ in Eq. (7) using the anelastic approximation and we simplify the expression of the squared Brunt–Väisälä frequency (Eq. (8)) using the approximation (3.1). Then, we simplify the radial momentum equation (Eq. (3)) to

$$i\omega\partial_{\zeta}\tilde{W} + N^2\zeta^2\mathcal{A}u^{\zeta} = 0, \qquad (13)$$

where we neglect the term $(P'/\rho_0^2)\partial_{\zeta}\rho_0$ thanks to the anelastic approximation and the vertical component of the acceleration

and of the Coriolis acceleration, which are dominated by the buoyancy term since we assume that $\omega \ll N$ and $2\Omega \ll N$ within the TAR. Subsequently, the latitudinal component of the momentum equation (Eq. (4)) reduces to

$$i\omega\zeta\mathcal{B}u^{\theta} - 2\Omega\zeta C u^{\varphi} + \partial_{\theta}\tilde{W} = 0, \qquad (14)$$

where we neglect the terms involving the vertical wave velocity since $|u^\zeta| \ll \{|u^\theta|, |u^\varphi|\}$ because of the strong stable stratification. Finally, the azimuthal component of the momentum equation (Eq. (5)) simplifies, for the same reasons, into

$$i\omega\zeta\mathcal{D}u^{\varphi} + \left(2\Omega\zeta C + \frac{\sin\theta}{r_{\zeta}}\partial_{\theta}\Omega\right)u^{\theta} - \frac{im}{\sin\theta}\tilde{W} = 0. \qquad (15)$$

The coefficients $\mathcal{A}$, $\mathcal{B}$, $\mathcal{C}$, and $\mathcal{D}$ are defined in Table 1 and $\tilde{W} = \tilde{P}/\rho_0$ is the normalised pressure. We thus obtain a system of equations which has the same mathematical form as in the case of: (i) spherically symmetric uniformly rotating (Bildsten et al. 1996; Lee & Saio 1997) and differentially rotating (Mathis 2009) stars, (ii) weakly deformed uniformly rotating (Mathis & Prat 2019) stars and (iii) strongly deformed uniformly rotating (Paper I) stars. We thus still manage to partially decouple the vertical and horizontal components of the velocity. By solving the system formed by Eqs. (14) and (15) we can express $u^\theta$ and $u^\varphi$ as a function of $\tilde{W}$ as follows

$$\begin{aligned} u^{\theta}(\zeta, x) &= \mathcal{L}_{\omega^{\text{in}}, m}^{\theta} \left[ \tilde{W}(\zeta, x) \right] \\ &= -i \frac{1}{\zeta} \frac{1}{\omega(\zeta, x)} \frac{1}{\mathcal{B}(\zeta, x)} \frac{1}{\sqrt{1-x^2}} \\ &\quad \times \left[ (1-x^2) \left( 1 + \frac{\nu^2(\zeta, x)C^2(\zeta, x)/\mathcal{B}(\zeta, x) + \mathcal{F}(\zeta, x)}{\mathcal{E}(\zeta, x)\mathcal{B}(\zeta, x) - \mathcal{F}(\zeta, x)} \right) \partial_x \right. \\ &\quad \left. + m \frac{\nu(\zeta, x)C(\zeta, x)}{\mathcal{E}(\zeta, x) - \mathcal{F}(\zeta, x)} \right] \tilde{W}, \end{aligned} \qquad (16)$$

$$\begin{aligned} u^{\varphi}(\zeta, x) &= \mathcal{L}_{\omega^{\text{in}}, m}^{\varphi} \left[ \tilde{W}(\zeta, x) \right] \\ &= \frac{1}{\zeta} \frac{1}{\omega(\zeta, x)} \frac{1}{\mathcal{E}(\zeta, x) - \mathcal{F}(\zeta, x)} \frac{1}{\sqrt{1-x^2}} \\ &\quad \times \left[ (1-x^2) \left( \frac{\nu(\zeta, x)C(\zeta, x)}{\mathcal{B}(\zeta, x)} + \frac{\mathcal{F}(\zeta, x)}{\nu(\zeta, x)C(\zeta, x)} \right) \partial_x + m \right] \tilde{W}, \end{aligned} \qquad (17)$$

where $\mathcal{E}$ and $\mathcal{F}$ are defined in Table 1, $x = \cos\theta$ is the reduced latitudinal coordinate and

$$\nu(\zeta, \theta) = \frac{2\Omega(\zeta, \theta)}{\omega(\zeta, \theta)}, \qquad (18)$$

is the spin parameter.

The structure of the new equations that we obtain in the spheroidal differentially rotating case is similar to the one with uniform rotation (Paper I). But we can point out two major differences. First, the wave frequency $\omega$ and the spin parameter $\nu$ are no longer uniform but they now depend on the pseudo-radius $\zeta$ and the reduced latitudinal coordinate $x$. That is why it is more convenient to index the operators and the variables by $\omega^{\text{in}}$ (independent of $\zeta$ and $\theta$) instead of using the spin parameter $\nu$ as in Paper I. Second, a new coefficient $\mathcal{F}$ is involved in the derivation of the dynamics of GWs that models the action of differential rotation.
<table><caption>**Table 1.** Terms involved in the derivation of the generalised Laplace tidal equation in the general case of spheroidal geometry and in the particular case of spherical geometry with differential rotation.</caption><thead><tr><th>Terms</th><th>Spheroidal</th><th>Spherical</th></tr></thead><tbody><tr><td>A</td><td><sup>1</sup>&frasl;<sub>r<sup>2</sup>r<sub>ζ</sub></sub></td><td><sup>1</sup>&frasl;<sub>ζ<sup>2</sup></sub></td></tr><tr><td>B</td><td><sup>r<sup>2</sup> + (1 - x<sup>2</sup>)r<sub>x</sub><sup>2</sup></sup>&frasl;<sub>r<sup>2</sup>r<sub>ζ</sub></sub></td><td>1</td></tr><tr><td>C</td><td><sup>-(1 - x<sup>2</sup>)r<sub>x</sub> + rx</sup>&frasl;<sub>rr<sub>ζ</sub></sub></td><td>x</td></tr><tr><td>D</td><td><sup>1</sup>&frasl;<sub>r<sub>ζ</sub></sub></td><td>1</td></tr><tr><td>E</td><td><sup>D - <sup>v<sup>2</sup>C<sup>2</sup></sup>&frasl;<sub>B</sub></sup>&frasl;<sub>v</sub></td><td>1 - v<sup>2</sup>x<sup>2</sup></td></tr><tr><td>F</td><td>-<sup>v</sup>&frasl;<sub>ω</sub><sup>C D</sup><sub>B</sub>(1 - x<sup>2</sup>)∂<sub>x</sub>Ω</td><td>-<sup>v x</sup>&frasl;<sub>ω</sub>(1 - x<sup>2</sup>)∂<sub>x</sub>Ω</td></tr></tbody></table>

## 4. Dynamics of low-frequency gravito-inertial waves

We now derive the generalised Laplace tidal equation for the normalised pressure $\tilde{W}$, which allows us to compute the frequencies and periods of low-frequency GWs and to build the corresponding seismic diagnostics. Applying the anelastic approximation and the approximation $N^2(\zeta, \theta) \approx N^2(\zeta)$ in the continuity equation (Eq. (6)) simplifies it into

$$ \zeta \partial_{\zeta} \rho_0 u^{\zeta} + \rho_0 \left[ \frac{\partial_{\zeta} (\zeta^2 u^{\zeta})}{\zeta} + \frac{1}{\sin \theta} \partial_{\theta} (\sin \theta u^{\theta}) - \frac{im}{\sin \theta} u^{\varphi} \right] = 0. \quad (19) $$

### 4.1. JWKB approximation

Under the assumption that $\omega \ll N$, each scalar field and each component of $u$ can be expanded using the two-dimensional Jeffreys-Wentzel-Kramers-Brillouin JWKB approximation (Fröman & Fröman 1965). In this case, the spatial structure of the waves can be described by the product of a rapidly oscillating plane-like wave function in the pseudo-radial direction multiplied by a slowly varying envelope:

$$ \tilde{W}(\zeta, \theta) = \sum_k \left\{ w_{\omega^{\text{in}}km}(\zeta, \theta) \frac{A_{\omega^{\text{in}}km}}{k^{1/2}} \exp \left[ i \int_{V;\omega^{\text{in}}km} d\zeta \right] \right\}, \quad (20) $$

$$ u^j(\zeta, \theta) = \sum_k \left\{ \hat{u}_{\omega^{\text{in}}km}^j(\zeta, \theta) \frac{A_{\omega^{\text{in}}km}}{k^{1/2}} \exp \left[ i \int_{V;\omega^{\text{in}}km} d\zeta \right] \right\}, \quad (21) $$

with $j \equiv \{\zeta, \theta, \varphi\}$, $k$ is the index of a latitudinal eigenmode (cf. Sect. 4.4) and $A_{\omega^{\text{in}}km}$ is the amplitude of the wave. Substituting the expansion given in Eqs. (20) and (21) into Eqs. (13), (16) and (17), the final pseudo-radial, latitudinal and azimuthal components of the velocity are obtained:

$$ \hat{u}_{\omega^{\text{in}}km}^{\zeta}(\zeta, x) = \frac{k_{V;\omega^{\text{in}}km}(\zeta)}{N^2(\zeta)} \frac{\omega_{km}(\zeta, x)}{\zeta^2 \mathcal{A}(\zeta, x)} w_{\omega^{\text{in}}km}(\zeta, x), \quad (22) $$

$$ \hat{u}_{\omega^{\text{in}}km}^{\theta}(\zeta, x) = \mathcal{L}_{\omega^{\text{in}}m}^{\theta}[w_{\omega^{\text{in}}km}(\zeta, x)], \quad (23) $$

$$ \hat{u}_{\omega^{\text{in}}km}^{\varphi}(\zeta, x) = \mathcal{L}_{\omega^{\text{in}}m}^{\varphi}[w_{\omega^{\text{in}}km}(\zeta, x)]. \quad (24) $$

### 4.2. Approximation: $\mathcal{A}(\zeta, \theta) \approx \mathcal{A}(\zeta)$

As in Paper I, to be able to introduce the eigenvalues $\Lambda_{\omega^{\text{in}}km}$ which depends only on $\zeta$ and derive the generalised Laplace tidal equation, we have to do a partial separation between the pseudo-radial and latitudinal variables in the pseudo-radial component of the velocity (Eq. (22)). So, we have to assume that the coefficient $\mathcal{A}$ depends mainly on the pseudo-radius $\zeta$. The validity domain of this approximation is presented in Sect. 5.1 using two-dimensional ESTER stellar models. Unlike the uniformly rotating case (Paper I), this assumption alone is insufficient to partially separate the variables in Eq. (22) in the differentially rotating case, so it is mandatory to assume a third approximation.

### 4.3. Approximation: $\Omega(\zeta, \theta) \approx \Omega(\zeta)$

The wave frequency in the rotating frame (Eq. (11)) depends on $\zeta$ and $\theta$. This dependency comes from the differential rotation $\Omega(\zeta, \theta)$. Therefore, in order to preform the partial separation of variables in Eq. (22), the angular velocity $\Omega$ should have a dependency only on $\zeta$. We examine in detail the validity of this approximation for two-dimensional ESTER stellar models in Sect. 5.1.1.

### 4.4. Generalised Laplace tidal equation

Substituting the expansion given in Eqs. (20) and (21) into Eqs. (13), (16) and (17), then replacing them into the continuity equation (Eq. (19)) we obtain the generalised Laplace tidal operator (GLTO)

$$ \mathcal{L}_{\omega^{\text{in}}m} = \omega \partial_x \left[ \frac{1}{\omega \mathcal{B}(\mathcal{E}-\mathcal{F})} \left( \mathcal{E} + \frac{\nu^2 C^2}{\mathcal{B}} \right) (1-x^2) \partial_x \right] \\ - \frac{m\mathcal{F}}{\nu\mathcal{C}(\mathcal{E}-\mathcal{F})} \partial_x + m\omega \partial_x \left( \frac{\nu\mathcal{C}}{\omega\mathcal{B}(\mathcal{E}-\mathcal{F})} \right) - \frac{m^2}{(\mathcal{E}-\mathcal{F})(1-x^2)}, \quad (25) $$

thus the generalised Laplace tidal equation (GLTE) for the normalised pressure $w_{\omega^{\text{in}}km}$ is written as

$$ \mathcal{L}_{\omega^{\text{in}}m} [w_{\omega^{\text{in}}km}] = -\Lambda_{\omega^{\text{in}}km}(\zeta) w_{\omega^{\text{in}}km}, \quad (26) $$

where $\Lambda_{\omega^{\text{in}}km}$ are the eigenvalues deduced from the following dispersion relation:

$$ k_{V;\omega^{\text{in}}km}^2(\zeta) = \frac{N^2(\zeta)\mathcal{A}(\zeta)}{\omega_{km}^2(\zeta)} \Lambda_{\omega^{\text{in}}km}(\zeta), \quad (27) $$

and $w_{\omega^{\text{in}}km}$ the generalised Hough functions (eigenfunctions) of the GLTE. We choose to define our latitudinal ordering number $k$ to enumerate, for each $(\omega^{\text{in}}, m)$, the infinite set of solutions as in Lee & Saio (1997), Mathis (2009) and Mathis & Prat (2019).

The GLTO (Eq. (25)) reduces to the Laplace tidal operator in the case of a uniformly rotating deformed star (Paper I). In this case, $\omega$ and $\nu$ are constants so we can take them out of the derivatives and $\mathcal{F} = 0$; thus we obtain

$$ \mathcal{L}_{vm} = \partial_x \left[ \frac{1}{\mathcal{B}} \left( 1 + \frac{\nu^2 C^2}{\mathcal{E}\mathcal{B}} \right) (1 - x^2) \partial_x \right] \\ + \left( m\nu \partial_x \left( \frac{C}{\mathcal{E}\mathcal{B}} \right) - \frac{m^2}{\mathcal{E}(1 - x^2)} \right). \quad (28) $$

Furthermore, the GLTO (Eq. (25)) reduces to the Laplace tidal operator derived by Mathis (2009) in the case of a spherical differentially rotating star. In this case, we can replace coefficients
A, B, C, D, and E by their analytical expressions in the particular case of a spherical geometry presented in Table 1 of Paper I, so the new coefficient F defined here can be written in this particular case as

$$F(\zeta, x) = -\frac{vx}{\omega}(1 - x^2)\partial_x\Omega, \qquad (29)$$

thus we obtain

$$L_{\omega^{in}m} = \omega \partial_x \left[ \frac{1-x^2}{\omega G} \partial_x \right] - \frac{m \partial_x \Omega}{\omega G} (1-x^2) \partial_x + \left( m \omega \partial_x \left( \frac{vx}{\omega G} \right) - \frac{m^2}{G(1-x^2)} \right), \qquad (30)$$

with

$$G(\zeta, x) = \mathcal{E}(\zeta, x) - \mathcal{F}(\zeta, x) = 1 - v^2 x^2 + \frac{vx}{\omega}(1 - x^2)\partial_x \Omega. \qquad (31)$$

### 4.5. Asymptotic frequency and period spacing of low-frequency GIWs

By substituting the dispersion relation (Eq. (27)) into the following quantisation condition in the vertical direction defined by Unno et al. (1989), Gough (1993), and Christensen-Dalsgaard (1997) and used in Paper I

$$\int_{\zeta_1}^{\zeta_2} k_{V;\omega^{in}nkm} d\zeta = (n + 1/2)\pi, \qquad (32)$$

where $\zeta_1$ and $\zeta_2$ are the turning points of the Brunt-Väisälä frequency N and n is the radial order, we compute numerically the asymptotic frequencies of low-frequency GIWs in deformed differentially rotating stars in the inertial frame. We describe in detail the used numerical method in Sect. 6. In the case of a weak differential rotation, we can extract, analytically, the expression of the asymptotic frequencies in the rotating frame

$$\omega_{nkm} = \frac{\int_{\zeta_1}^{\zeta_2} N(\zeta) \sqrt{\mathcal{A}(\zeta) \Lambda_{\omega^{in}km}(\zeta)} d\zeta}{(n + 1/2)\pi}, \qquad (33)$$

with

$$\omega_n^{in} = \omega_{nkm}^{in} = \omega_{nkm} + m\Omega_{av}, \qquad (34)$$

where $\Omega_{av}$ is here the averaged rotation. Then by applying a first-order Taylor development following Bouabid et al. (2013), we can generalise his expression of the corresponding period spacing obtained in the spherical case

$$\Delta P_{km} = P_{nkm} - P_{n-1km} \approx \frac{2\pi^2}{\int_{\zeta_1}^{\zeta_2} N \sqrt{\mathcal{A} \Lambda_{\omega^{in}km} d\zeta (1 + \alpha)}}, \qquad (35)$$

with

$$\alpha = \frac{1}{2} \frac{\int_{\zeta_1}^{\zeta_2} N \sqrt{\mathcal{A} \Lambda_{\omega^{in}km} \frac{d \ln \Lambda_{\omega^{in}km}}{d \ln \omega^{in}}} d\zeta}{\int_{\zeta_1}^{\zeta_2} N \sqrt{\mathcal{A} \Lambda_{\omega^{in}km} d\zeta}}. \qquad (36)$$

### 5. Application to rapidly and differentially rotating early-type stars

As in Paper I, we use here also ESTER models (Espinosa Lara & Rieutord 2013) to implement our equations. We use 3 $M_\odot$ stellar models with a hydrogen mass fraction in the core $X_c = 0.7$ rotating at the equator at [0%, 90%] of the Keplerian break-up rotation rate.

#### 5.1. Domain of validity of the generalised TAR

We study, within the framework of ESTER models, the validity domain of the three approximations (3.1), (4.2) and (4.3) that are necessary to build the TAR in the case of rapidly differentially rotating deformed bodies:

$$N^2(\zeta, \theta) \approx N^2(\zeta), \qquad (37)$$

$$\mathcal{A}(\zeta, \theta) \approx \mathcal{A}(\zeta), \qquad (38)$$

$$\Omega(\zeta, \theta) \approx \Omega(\zeta), \qquad (39)$$

The validity of the approximations (37) and (38) is well discussed in Paper I (the $N^2$ and $\mathcal{A}$ profiles are the same in the uniformly and differentially rotating cases). Here, we focus on determining the validity domain of the approximation (39) by applying the same method used to evaluate the other two assumptions that we recall here.

#### 5.1.1. Validity of $\Omega(\zeta, \theta) \approx \Omega(\zeta)$

To visualise the problem, first we illustrate in Fig. 2 the function $\Omega(\zeta, \theta)$ computed with an ESTER model (3 $M_\odot$, $X_c = 0.7$) rotating at the equator at $\Omega/\Omega_K = 20\%$. We can see that the angular velocity is equal to 14.66 µHz at the core, and at the surface we get 11.41 µHz at the pole and 12.22 µHz at the equator. Then, to evaluate the committed error by adopting the approximation

$$\Omega(\zeta, \theta) \approx \Omega(\zeta), \qquad (40)$$

we represent in Fig. 3 the angular velocity profile as a function of the pseudo-radius $\zeta$ for different values of the colatitude $\theta_f$, at a fixed rotation rate at the equator $\Omega = 0.2\Omega_K$. Visually, we can notice that the latitudinal gradient of the angular velocity is less important than the radial gradient. In order to confirm this observation, we show, in the bottom panels, the corresponding relative error $\Delta\Omega_{\theta_f}(\zeta)$ between the exact value $\Omega(\zeta, \theta_f)$ given by the model and the approximated model $\Omega_{approx}(\zeta)$ (the weighted average of the exact value over the colatitude $\theta$, $\bar{\Omega}(\zeta)$). We define these quantities as follows:

$$\delta\Omega_{\theta_f}(\zeta) = \frac{\Omega_{approx}(\zeta) - \Omega(\zeta, \theta_f)}{\Omega(\zeta, \theta_f)}, \qquad (41)$$

with

$$\Omega_{approx}(\zeta) = \bar{\Omega}(\zeta) = \frac{\int_{0}^{\pi/2} \Omega(\zeta, \theta) \sin \theta d\theta}{\int_{0}^{\pi/2} \sin \theta d\theta} = \int_{0}^{\pi/2} \Omega(\zeta, \theta) \sin \theta d\theta, \qquad (42)$$

where $\theta_f$ is a fixed value of the colatitude $\theta$ and we used the symmetry propriety of $\Omega$ with respect to the equator. We can conclude that this approximation is valid with a level of confidence of 90% for all pseudo-radii $\zeta$ and for all rotation rates $\Omega/\Omega_K \in [0\%, 90\%]$ since the error rate is always lower than the maximal error rate fixed to 10%. We also recall here that within the physics treated in the ESTER code the differential rotation results from the baroclinic meridional circulation (Zahn 1992; Maeder & Zahn 1998; Mathis & Zahn 2004; Rieutord 2006) and an isotropic viscous transport. A weaker dependence of $\Omega$ on $\theta$ can be enforced by horizontal turbulent transport (Zahn 1992; Maeder 2003; Mathis et al. 2004, 2018; Park et al. 2021) or magnetic fields (Mestel & Weiss 1987).
<figure><img src="image_3.png" /><figcaption>Fig. 2. Angular velocity profile $\Omega(\zeta, \theta)$ of an ESTER model ($3 M_\odot$, $X_c = 0.7$) rotating at the equator at $\Omega/\Omega_K = 20\%$ (the isolines are represented in black).</figcaption></figure>

### 5.1.2. Validity domain of the TAR

We determine the validity domain of the three approximations (37)–(39) as a function of the rotation rate and the pseudo-radius by varying systematically the normalised rotation rate of the star $\Omega/\Omega_K$ and by calculating for each case the associated maximum relative error. Afterwards, we fix a maximum error rate equal to 10% and we deduce the pseudo-radius limit $\zeta_{limit}$ where the maximum relative error exceeds this threshold and we decide that the approximations become invalid. Physically, it means that the variation of the structural properties and rotation profile are weak in the latitudinal direction and vary mostly with the pseudo-radius that allows us to build the TAR. The best method in the future to improve the control of the adopted assumptions would be to compute asymptotic frequencies using the latitudinally averaged quantities within the TAR and to compare them to frequencies computed using 2D oscillation codes applied to deformed models (e.g. Ouazzani et al. 2015, 2017; Reese et al. 2021). In Fig. 4, we display the pseudo-radius limit $\zeta_{limit}$ as a function of the rotation rate $\Omega/\Omega_K$ for 3 solar masses ESTER models with $X_c = 0.7$ and $X_c = 0.2$. This curve defines the validity domains of the approximations (37)–(39) which consequently define the validity domain of the TAR. We note that the influence of the hydrogen mass fraction in the core $X_c$ on the validity domain of the TAR is very weak. It is therefore possible to apply the generalised TAR to differentially rotating early-type stars rotating up to 20% of the Keplerian critical angular velocity. This limit is the same as the one found in Paper I, so the differential rotation have no influence on the validity domain of the TAR. Chemical gradients created near the core along the evolution are not treated in 2D ESTER models yet. Their treatment can be foreseen by adding them ‘by hand’ to explore their effects

<figure><img src="image_4.png" /><figcaption>Fig. 3. Profile of the angular velocity $\Omega$ and the relative error $\delta\Omega_{\theta t}$ of the approximation $\Omega(\zeta, \theta) \approx \Omega(\zeta)$ as a function of $\zeta$ at different colatitudes $\theta_t$ using an ESTER model rotating at the equator at 20% (above) and 90% (below) of the Keplerian break-up rotation rate (the light orange area indicate the margin of error which we allow).</figcaption></figure>

as this has been done for instance for acoustic glitches by Reese et al. (2014, 2021).

### 5.2. Eigenvalues and Hough functions

We solve the GLTE for different pseudo-radii, wave frequency in the inertial frame, and rotation rates within the scope of the defined validity domain using an implementation based on Chebyshev polynomials (Wang et al. 2016) as in Paper I. Here, we are no longer able to calculate the spectrum of the GLTE as a function of the spin parameter $v$ because in the differentially rotating case, $v$ is no longer a constant value but it becomes a function of $\zeta$ and $\theta$. So to make the comparison between the new results and our previous results clear, we define the weighted average of the spin parameter

$$ \bar{v} = \int_{0}^{1} \int_{0}^{\pi/2} v(\zeta, \theta) \sin \theta d\theta d\zeta. \quad (43) $$

Figure 5 shows the eigenvalues $\Lambda_{\omega^{in}km}$ as a function of $\omega^{in}$ for $m = 1$ and $\zeta = \zeta_{limit} = 1$ at $\Omega/\Omega_K = 0.2$. Since $\Omega$ points in the direction of $\theta = 0$ and the oscillations are proportional to $e^{i(\omega t - m\phi)}$, a prograde (retrograde) oscillation correspond to a
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Validity domain of the TAR within the framework of 3 $M_\odot$, $X_c = 0.7$ (above) and 3 $M_\odot$, $X_c = 0.2$ (below) ESTER models as a function of the rotation rate $\Omega/\Omega_K$ and the pseudo-radius $\zeta$ (with 90% degree of confidence).</figcaption></figure>

positive (negative) value of the product $mv$ which depends on the spatial coordinates. This figure reveals that the differential rotation of the star causes at $\zeta = 1$ a modest gradual shift in the eigenvalues. We recall here the two main families of eigenvalues. The first one being gravity-like solutions ($\Lambda_{\omega^{\text{in}}km} \ge 0$). They exist for any value of $\omega^{\text{in}}$ and we attach positive $k$'s to them. They correspond to gravity waves ($g$ modes) modified by rotation. The second one is Rossby-like solutions. They exist only for specific values of wave frequency in the inertial reference frame ($\omega^{\text{in}} \in [\omega^{\text{in}}, \omega^{\text{in}}]$) and we attach negative $k$'s to them. They appear only in rotating stars and they correspond to Rossby modes ($r$ modes) if they are retrograde and have positive eigenvalues ($\Lambda_{\omega^{\text{in}}km} > 0$) and to overstable convective modes if they are prograde and have negative eigenvalues ($\Lambda_{\omega^{\text{in}}km} < 0$).

Contrary to the uniformly rotating case we obtain here eigenvalues which depend considerably on the pseudo-radius $\zeta$ as shown in Fig. 6 where we represent the spectrum of the GLTE as a function of the pseudo-radius at $\Omega/\Omega_K = 0.2$, $m = 1$ and $\omega^{\text{in}} = 1.12 \times 10^{-4} \text{ rad s}^{-1}$. More precisely we can see that the differential rotation influence especially the region close to the convective core ($0.153 \le \zeta \le 0.35$) where the rotation is maximal and the radial gradient of the angular velocity is the highest (cf. Fig. 2). We can see that this dependency becomes very low

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Spectrum of the GLTE as a function of the wave frequency in the inertial frame $\omega^{\text{in}}$ at $\zeta = 1$ and $m = 1$ for $\Omega = 0.2\Omega_K$ in the case of a uniformly rotating star (above) and of a differentially rotating star (below). Blue (respectively, orange) dots correspond to even (respectively, odd) eigenfunctions. The vertical lines correspond to the mean values of the spin parameter $\bar{\nu}$.</figcaption></figure>

in the external region $0.5 \le \zeta \le 1$ where we obtain eigenvalues very close to the ones in the uniformly rotating case.

We focus now on the influence of the differential rotation on the generalised Hough functions $w_{km}$ which varies with the pseudo-radius $\zeta$ and the horizontal coordinate $x$. This dependence is illustrated in Fig. 7 at $\Omega/\Omega_K = 20\%$ for prograde dipole $\{k = 0, m = 1\}$ modes with $\omega^{\text{in}} = 10^{-4} \text{ rad s}^{-1}$ and retrograde Rossby $\{k = -2, m = -1\}$ modes with $\omega^{\text{in}} = -5.5 \times 10^{-5} \text{ rad s}^{-1}$ for four different cases: Fig. 7a represents our starting point where the rotation is uniform (Paper I); in Fig. 7b, we take into account the pseudo-radial differential rotation; in Fig. 7c, we take into account the latitudinal differential rotation and Fig. 7d represents our final point where we take into account the full differential rotation profile. As the pseudo-radius decreases from the surface ($\zeta = 1$) to the edge of the radiative zone ($\zeta = 0.153$), Rossby-like solutions are slightly modified, whereas gravity-like solutions considerably change. We can also see that the major modification comes from the differential rotation gradient along the pseudo-radius $\zeta$ (Fig. 7b) while the differential rotation according to the colatitude $\theta$ modifies marginally the gravity-like solutions and has almost no impact on the Rossby-like solution (Fig. 7c). On top of that, we can see that the Hough functions at the surface are not too affected in all the cases and that the strongest variation occur in the interior of the star where the pseudo-radial gradient of the rotation is high.

Overall, we can see clearly in Fig. 7 that the differential rotation strongly influences and modifies the eigenfunctions of the GLTE by introducing a new dependency on the pseudo-radial coordinate $\zeta$, whereas in the spherically symmetric uniformly rotating case they were only dependant on the latitudinal coordinate $\theta$ while they were also slightly dependent on $\zeta$ in the uniformly rotating deformed case. In fact, the gravity-like and the Rossby-like solutions shift as the distance from the border of
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Spectrum of the GLTE as a function of the pseudo-radius for $m = 1$ and $\omega^{\text{in}} = 1.12 \times 10^{-4} \text{ rad s}^{-1}$ at $\Omega/\Omega_{\text{K}} = 0.2$ in the case of a uniformly rotating star (above) and of a differentially rotating star (below). Blue (respectively, orange) dots correspond to even (respectively, odd) Hough functions.</figcaption></figure>

the radiative zone of the star ($\zeta = 0.153$) to its surface ($\zeta = 1$) increases. More specifically, we observe that the gravity-like solutions migrate onwards, away from the equator ($x = 0$), causing a broadening of its shape. This corresponds to the modification of the region of propagation of gravito-inertial waves by the differential rotation (Mathis 2009; Mirouh et al. 2016; Prat et al. 2018). Indeed, since the angular velocity $\Omega(\zeta)$ decreases with the distance to the centre, the equatorial trapping of sub-inertial ($\bar{v} > 1$) gravito-inertial modes become less and less important.

# 6. Asymptotic seismic diagnosis

## 6.1. Asymptotic period spacing pattern

In order to compute the period spacing patterns, we adapt the method developed in Henneco et al. (2021) and used in Paper I. First, we calculate for each radial order $n$ of a given mode ($k, m$), the corresponding asymptotic frequencies in the inertial frame $\omega_{nkm}^{\text{in}}$ using Eqs. (11) and (33). Then, we can calculate the asymptotic periods in the inertial frame ($P_{nkm}^{\text{in}} = 2\pi/\omega_{nkm}^{\text{in}}$) and the corresponding period spacing which is shown in Fig. 8 for $\{k = 0, m = 1\}$, $\{k = 1, m = 0\}$, and $\{k = 0, m = -1\}$ modes. The periods represented here are calculated for radial orders between $n = 5$ and $n = 50$ but the Cowling approximation that we adopted here to develop our formalism is valid only for high radial orders (Bouabid et al. 2013; Ouazzani et al. 2017). That is why we hatched the region with low radial order ($n < 20$), so we can identify where the Cowling approximation is potentially invalid. We find a net decrease in the period spacing $\Delta P_{nkm}^{\text{in}}$ caused by the differential rotation. The centrifugal and the differential rotation effects are moderately significant when assessed within the validity domain of the generalised TAR nonetheless the global

characteristics of the period spacing pattern are conserved. This is in good agreement with the results obtained by Ouazzani et al. (2017). They computed gravito-inertial modes and their period spacing (represented in Fig. 6 of their article) using the ACOR 2D oscillation code with deformed stellar models (see Ouazzani et al. 2015 for details on their method), which take the centrifugal acceleration into account following the method of deformation of acoustic models by Roxburgh (2006). Ballot et al. (2012) also studied the effect of the centrifugal acceleration on the period spacing. They found, for a sectoral mode ($k = 0$), a discrepancy between the period spacing (represented in Fig. 2 of their article) computed using the spherical TAR and the complete computations using TOP with spheroidal models. They demonstrate that this discrepancy originates in the centrifugal distortion of the 2D models. But for other modes, they do not see a considerable effect of the centrifugal deformation on the period spacing pattern in agreement with our results.

Even though the centrifugal deformation effects predicted by the TAR are weak in intermediate-mass stars, they are theoretically detectable for some radial orders of several modes using observations from TESS (Transiting Exoplanet Survey Satellite; Ricker et al. 2014) and Kepler (Paper I). This is why we do not neglect them and we study the differential rotation and the deformation effects simultaneously. To quantify the induced variation and evaluate its detectability, we compute the frequency differences.

## 6.2. Detectability of the differential rotation effect

To evaluate the detectability of the differential rotation effect with space-based photometric observations, we compute first the frequency differences $\Delta f_{\text{centrifugal}}$ between asymptotic frequencies calculated in the standard TAR ($\text{TAR}_s$) and those calculated in the generalised TAR with a uniform rotation ($\text{TAR}_{g_u}$) and a differential rotation ($\text{TAR}_{g_d}$) through Eq. (33):

$$ \Delta f(n) = |f_{\text{TAR}_x}(n) - f_{\text{TAR}_s}(n)|, \quad (44) $$

with $x$ and $y$ two parameters that allow us to choose the effect to be evaluated:

$$ (x, y) = \begin{cases} (g_u, s) & \text{for the effect of the centrifugal acceleration;} \\ (g_u, g_d) & \text{for the effect of the differential rotation;} \\ (g_d, s) & \text{for the two effects simultaneously.} \end{cases} $$

Then, by comparing the obtained frequency differences with the frequency resolutions ($f_{\text{res}} = 1/T_{\text{obs}}$) of Kepler and TESS light curves covering quasi-continuously observation times of $T_{\text{obs}} = 4$ years and $T_{\text{obs}} = 351$ days, respectively, we can deduce the radial orders $n$'s for which the frequency differences are expected to be theoretically detectable. In Fig. 9, we display these results for $\{k = 0, m = 1\}$, $\{k = 1, m = 0\}$, and $\{k = 0, m = -1\}$ modes rotating at $0.2 \Omega_{\text{K}}$. In this figure, we show the detectability of the centrifugal effect, the differential rotation effect and of the two effects taken simultaneously into account. Unlike the centrifugal effect, we can see clearly that the effect of the differential rotation for these modes is, theoretically, largely detectable for all radial order higher than 15 using TESS and Kepler observations (not represented here because its observation time is larger than the one of TESS so the detectability using Kepler is guaranteed). The detectability of the centrifugal effect is discussed in details in Paper I.

As we pointed out in Paper I, this is a formal evaluation due to the presence of large correlations between the different parameters of our stellar models. This will mask the effect of the
<figure><img src="image_8.png" /><figcaption>Fig. 7. Generalised Hough functions (normalised) as a function of the horizontal coordinate x at different pseudo-radii from $\zeta = 0.153$ (blue) to $\zeta = 1$ (yellow) at $\Omega/\Omega_K = 0.2$ for gravity-like solutions $\{k = 0, m = 1\}$ with $\omega^{\text{in}} = 10^{-4} \text{ rad s}^{-1}$ corresponding to $\bar{v} = 10.62$ (the solid lines) and Rossby-like solutions $\{k = -2, m = -1\}$ with $\omega^{\text{in}} = -5.5 \times 10^{-5} \text{ rad s}^{-1}$ corresponding to $\bar{v} = 6.48$ (the dotted lines) in the case of a star: (a) uniformly rotating, (b) differentially rotating according to the pseudo-radius, (c) differentially rotating according to the colatitude and (d) differentially rotating according to the pseudo-radius and the colatitude.</figcaption></figure>

centrifugal acceleration and the differential rotation in forward asteroseismic modelling analyses, even when it is theoretically detectable:

$$ \Delta f_i > f_{\text{res}}, i \equiv \{\text{centrifugal, differential}\}. \quad (45) $$

## 7. Evaluation of the terms hierarchy imposed by the TAR within differentially rotating deformed stars

The hierarchy of terms imposed by the TAR can be summarised by the following frequency hierarchy:

$$ 2\Omega \ll N, \quad (46) $$

$$ \omega \ll N, \quad (47) $$

which ensures the other hierarchies. Using the Brunt–Väisälä frequency and rotation profiles from ESTER models and using the asymptotic frequencies calculated in Sect. 6, we compare these frequencies and we discuss whether the TAR is still valid in differentially rotating deformed stars as in Paper I.

### 7.1. The strong stratification assumption: $2\Omega \ll N$

We evaluate the term $N/2\Omega$ using $3 M_\odot, X_c = 0.7$ ESTER models for rotation rates $\Omega/\Omega_K \in [0.1, 0.9]$. As shown in Fig. 10, the strong stratification assumption is valid only in the radiative zone away from the border between the convective core and the radiative envelope ($\zeta \ge 0.2$) at $\Omega/\Omega_K \le 0.2$. Beyond this critical value, the Brunt–Väisälä frequency and the frequency of rotation have a close order of magnitude so the strong stratification approximation is no longer valid.

### 7.2. The low frequency assumption: $\omega \ll N$

We evaluate the term $N/\omega$ using the asymptotic frequencies for $\{k = 0, m = 1\}, \{k = 1, m = 0\}, \text{ and } \{k = 0, m = -1\}$ modes calculated in Sect. 6. As shown in Fig. 11, the low frequency assumption is valid in all the space domain for the $\{k = 0, m = 1\}$ mode only for high radial orders ($n > 20$). Below this critical value, the Brunt–Väisälä frequency and the wave frequency can

have a very close order of magnitude near the transition layer between the convective core and the radiative envelope. But far from this interface, the low frequency approximation is valid for all the modes.

## 8. Discussion and conclusions

This work is a continuation of our previous study where we generalised the TAR in the case of strongly deformed, rapidly and uniformly rotating stars and planets. In this perspective, we study the possibility of carrying out a new generalisation of the TAR that abandons the assumption of uniform rotation and takes into account the radial and latitudinal differential rotation. We approached this exploration by deriving the generalised Laplace Tidal equation in a spheroidal coordinate system. We relied mainly on two assumptions that will define the validity domain of the generalised TAR (Paper I). The equation that we derive has a similar form to the one that is obtained when the TAR is applied to weakly rotating spherical stars (Lee & Saio 1997), differentially rotating spherical stars (Mathis 2009; Van Reeth et al. 2018), moderately rotating weakly deformed stars (Mathis & Prat 2019; Henneco et al. 2021) and uniformly rotating strongly deformed stars (Paper I). So with this new formalism, we can study GIWs in the radiative region of all types of stars and planets.

We apply this general formalism to rapidly and differentially rotating early-type stars using 2D ESTER models ($M = 3 M_\odot, X_c = 0.7$) where we found that the signature of the differential rotation effect in the period spacing patterns is stronger than the signature of the centrifugal effect for the $\{k = 0, m = 1\}$ mode at $\Omega/\Omega_K = 20\%$, typically up to a factor twenty. This is caused essentially by the presence of a strong pseudo-radial $\Omega$-gradients in stellar interior of early-type stars. This new generalisation of the TAR can be used to study the dissipation of stellar and planetary tides induced by low-frequency GIWs in rapidly and differentially rotating deformed stars and planets and the angular momentum transport with a formalism that can be directly implemented in ESTER models.

The next step will be the inclusion of the magnetic field in a non-perturbative way within the generalised TAR formalism.
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

**Fig. 8.** Period spacing pattern in the inertial frame computed for {$k = 0, m = 1$} (top), {$k = 1, m = 0$} (middle), and {$k = 0, m = -1$} (bottom) modes at $\Omega/\Omega_K = 20\%$ in a spherical (red) and a deformed (blue) star uniformly rotating (at the weighted mean of the rotation rate $\bar{\Omega} = 12.86$ µHz) and in a deformed star differentially rotating (yellow). The hatched pink area indicates where the Cowling approximation is potentially invalid. (The fluctuations in the period spacing pattern are caused by the numerical noise which is introduced by the numerical derivatives of the mapping with respect to $\zeta$ and $\theta$ used in the resolution of the GLTE.)

So far, Prat et al. (2019, 2020) and Van Beeck et al. (2020) have focused on the case where magnetic fields are weak enough to be treated within a perturbative treatment to study the effects of a magnetic field on the seismic parameters of *g* modes which become magneto-gravito inertial modes (Mathis & de Brye 2011, 2012). So a possible follow-up of this work (Paper III) is to take into account the stellar magnetic fields and generalise the TAR framework to the case of differentially rotating strongly deformed magnetic stars.

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

**Fig. 9.** Detectable radial orders *n* for {$k = 0, m = 1$} (top), {$k = 1, m = 0$} (middle), {$k = 0, m = -1$} (bottom) modes at a rotation rate $\Omega/\Omega_K = 20\%$ based on the frequency resolution of TESS. Blue, yellow and red dots represent respectively the centrifugal effect, the differential rotation effect and the sum of the two effects (the hatched pink area indicates where the Cowling approximation is potentially invalid).

<figure><img src="image_15.png" /></figure>

**Fig. 10.** $N/2\Omega$ term as a function of the pseudo-radius $\zeta$ at different rotation rates $\Omega/\Omega_K$ (the hatched area represents the convective region of the star).
<figure><img src="image_16.png" /><figcaption>**Fig. 11.** $N/\omega$ term as a function of the pseudo-radius $\zeta$ and the radial order $n$ for the {$k=0, m=1$} (top), {$k=1, m=0$} (middle) and {$k=0, m=-1$} (bottom) modes at $\Omega/\Omega_K = 20\%$.</figcaption></figure>

**Acknowledgements.** We thank the referee for their constructive comments and suggestions that allow us to improve our article. We thank the ESTER code developers for their efforts and for making their codes publicly available. We are also grateful to Pr. C. Aerts who kindly commented on the manuscript and suggested improvements. H.D., V.P. and S.M. acknowledge support from the European Research Council through ERC grant SPIRE 647383. H.D. and S.M. acknowledge support from the CNES PLATO grant at CEA/DAP. TVR gratefully acknowledges support from the Research Foundation Flanders (FWO) under grant agreement No. 12ZB620N.

## References

Aerts, C. 2021, Rev. Mod. Phys., 93, 015001

Aerts, C., Mathis, S., & Rogers, T. M. 2019, ARA&A, 57, 35

Ballot, J., Lignières, F., Reese, D. R., & Rieutord, M. 2010, A&A, 518, A30

Ballot, J., Lignières, F., Prat, V., Reese, D. R., & Rieutord, M. 2012, in Progress in Solar/Stellar Physics with Helio- and Asteroseismology, eds. H. Shibahashi, M. Takata, & A. E. Lynas-Gray, ASP Conf. Ser., 462, 389

Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55

Beck, P. G., Hambleton, K., Vos, J., et al. 2014, A&A, 564, A36

Beck, P. G., Kallinger, T., Pavlovski, K., et al. 2018, A&A, 612, A22

Bedding, T. R., Murphy, S. J., Colman, I. L., & Kurtz, D. W. 2015, Eur. Phys. J. Web Conf., 101, 01005

Bildsten, L., Ushomirsky, G., & Cutler, C. 1996, ApJ, 460, 827

Bonazzola, S., Gourgoulhon, E., & Marck, J.-A. 1998, Phys. Rev. D, 58, 104020

Borucki, W., Koch, D., Batalha, N., et al. 2009, in Transiting Planets, eds. F. Pont, D. Sasselov, & M. J. Holman, IAU Symp., 253, 289

Bouabid, M. P., Dupret, M. A., Salmon, S., et al. 2013, MNRAS, 429, 2500

Buysschaert, B., Aerts, C., Bowman, D. M., et al. 2018, A&A, 616, A148

Cantiello, M., Mankovich, C., Bildsten, L., Christensen-Dalsgaard, J., & Paxton, B. 2014, ApJ, 788, 93

Ceillier, T., Eggenberger, P., García, R. A., & Mathis, S. 2013, A&A, 555, A54

Charbonnel, C., Decressin, T., Amard, L., Palacios, A., & Talon, S. 2013, A&A, 554, A40

Christensen-Dalsgaard, J. 1997, Stellar Oscillations (CiteSeerX)

Christophe, S., Ouazzani, R. M., Ballot, J., et al. 2018, in SF2A-2018: Proceedings of the Annual Meeting of the French Society of Astronomy and Astrophysics, eds. P. Di Matteo, F. Billebaud, F. Herpin, et al., 129

Cowling, T. G. 1941, MNRAS, 101, 367

Decressin, T., Mathis, S., Palacios, A., et al. 2009, A&A, 495, 271

Deheuvels, S., García, R. A., Chaplin, W. J., et al. 2012, ApJ, 756, 19

Deheuvels, S., Doğan, G., Goupil, M. J., et al. 2014, A&A, 564, A27

Deheuvels, S., Ballot, J., Beck, P. G., et al. 2015, A&A, 580, A96

Deheuvels, S., Ballot, J., Eggenberger, P., et al. 2020, A&A, 641, A117

Dhouib, H., Prat, V., Van Reeth, T., & Mathis, S. 2021, A&A, 652, A154

Di Mauro, M. P., Ventura, R., Cardini, D., et al. 2016, ApJ, 817, 65

Dintrans, B., & Rieutord, M. 2000, A&A, 354, 86

Dintrans, B., Rieutord, M., & Valdettaro, L. 1999, J. Fluid Mech., 398, 271

Eckart, C. 1960, Hydrodynamics of Oceans and Atmospheres (Oxford: Pergamon Press)

Eggenberger, P., Montalbán, J., & Miglio, A. 2012, A&A, 544, L4

Espinosa Lara, F., & Rieutord, M. 2013, A&A, 552, A35

Fröman, N., & Fröman, P. O. 1965, JWKB Approximation: Contributions to the Theory (Amsterdam: North-Holland Publishing Company)

Gallet, F., & Bouvier, J. 2015, A&A, 577, A98

Gehan, C., Mosser, B., Michel, E., Samadi, R., & Kallinger, T. 2018, A&A, 616, A24

Gough, D. 1993, in Les Houches Session XLVII, eds. J. P. Zahn, & J. Zinn-Justin (Amsterdam: Elsevier), 399

Guo, Z., Gies, D. R., & Matson, R. A. 2017, ApJ, 851, 39

Henneco, J., Van Reeth, T., Prat, V., et al. 2021, A&A, 648, A97

Kallinger, T., Weiss, W. W., Beck, P. G., et al. 2017, A&A, 603, A13

Keen, M. A., Bedding, T. R., Murphy, S. J., et al. 2015, MNRAS, 454, 1792

Kurtz, D. W., Saio, H., Takata, M., et al. 2014, MNRAS, 444, 102

Lee, U., & Baraffe, I. 1995, A&A, 301, 419

Lee, U., & Saio, H. 1987, MNRAS, 225, 643

Lee, U., & Saio, H. 1997, ApJ, 491, 839

Li, G., Van Reeth, T., Bedding, T. R., Murphy, S. J., & Antoci, V. 2019, MNRAS, 487, 782

Li, G., Van Reeth, T., Bedding, T. R., et al. 2020, MNRAS, 491, 3586

Lignières, F., Rieutord, M., & Reese, D. 2006, A&A, 455, 607

Maeder, A. 2003, A&A, 399, 263

Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (Springer), 827

Maeder, A., & Meynet, G. 2000, A&A, 361, 159

Maeder, A., & Zahn, J.-P. 1998, A&A, 334, 1000

Marques, J. P., Goupil, M. J., Lebreton, Y., et al. 2013, A&A, 549, A74

Mathis, S. 2009, A&A, 506, 811

Mathis, S. 2013, in Transport Processes in Stellar Interiors, eds. M. Goupil, K. Belkacem, C. Neiner, F. Lignières, & J. J. Green, 865, 23

Mathis, S., & de Brye, N. 2011, A&A, 526, A65

Mathis, S., & de Brye, N. 2012, A&A, 540, A37

Mathis, S., & Prat, V. 2019, A&A, 631, A26

Mathis, S., & Zahn, J. P. 2004, A&A, 425, 229

Mathis, S., Palacios, A., & Zahn, J. P. 2004, A&A, 425, 243

Mathis, S., Neiner, C., & Tran Minh, N. 2014, A&A, 565, A47

Mathis, S., Prat, V., Amard, L., et al. 2018, A&A, 620, A22

Mestel, L., & Weiss, N. O. 1987, MNRAS, 226, 123

Mirouh, G. M., Baruteau, C., Rieutord, M., & Ballot, J. 2016, J. Fluid Mech., 800, 213

Mombarg, J. S. G., Van Reeth, T., Pedersen, M. G., et al. 2019, MNRAS, 485, 3248

Moravveji, E., Townsend, R. H. D., Aerts, C., & Mathis, S. 2016, ApJ, 823, 130

Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012, A&A, 548, A10

Murphy, S. J., Fossati, L., Bedding, T. R., et al. 2016, MNRAS, 459, 1201

Ogilvie, G. I., & Lin, D. N. C. 2004, ApJ, 610, 477

Ouazzani, R. M., Dupret, M. A., & Reese, D. R. 2012, A&A, 547, A75

Ouazzani, R. M., Roxburgh, I. W., & Dupret, M. A. 2015, A&A, 579, A116

Ouazzani, R.-M., Salmon, S. J. A. J., Antoci, V., et al. 2017, MNRAS, 465, 2294
H. Dhouib et al.: The traditional approximation of rotation for rapidly rotating stars and planets. II.

Ouazzani, R. M., Marques, J. P., Goupil, M. J., et al. 2019, A&A, 626, A121

Ouazzani, R. M., Lignières, F., Dupret, M. A., et al. 2020, A&A, 640, A49

Pápics, P. I., Tkachenko, A., Aerts, C., et al. 2015, ApJ, 803, L25

Pápics, P. I., Tkachenko, A., Van Reeth, T., et al. 2017, A&A, 598, A74

Park, J., Prat, V., Mathis, S., & Bugnet, L. 2021, A&A, 646, A64

Pedersen, M. G., Aerts, C., Pápics, P. I., et al. 2021, Nat. Astron., 5, 715

Prat, V., Lignières, F., & Ballot, J. 2016, A&A, 587, A110

Prat, V., Mathis, S., Lignières, F., Ballot, J., & Culpin, P. M. 2017, A&A, 598, A105

Prat, V., Mathis, S., Augustson, K., et al. 2018, A&A, 615, A106

Prat, V., Mathis, S., Buysschaert, B., et al. 2019, A&A, 627, A64

Prat, V., Mathis, S., Neiner, C., et al. 2020, A&A, 636, A100

Reese, D., Lignières, F., & Rieutord, M. 2006, A&A, 455, 621

Reese, D. R., Lara, F. E., & Rieutord, M. 2014, in Precision Asteroseismology, eds. J. A. Guzik, W. J. Chaplin, G. Handler, & A. Pigulski, 301, 169

Reese, D. R., Mirouh, G. M., Espinosa Lara, F., Rieutord, M., & Putigny, B. 2021, A&A, 645, A46

Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2014, in Space Telescopes and Instrumentation 2014: Optical, Infrared, and Millimeter Wave, eds. J. M. Oschmann, Jr., M. Clampin, G. G. Fazio, & H. A. MacEwen, SPIE Conf. Ser., 9143, 914320

Rieutord, M. 2006, A&A, 451, 1025

Rieutord, M., Dintrans, B., Lignières, F., Corbard, T., & Pichon, B. 2005, in SF2A-2005: Semaine de l'Astrophysique Française, eds. F. Casoli, T. Contini, J. M. Hameury, & L. Pagani, 759

Rieutord, M., Espinosa Lara, F., & Putigny, B. 2016, J. Comput. Phys., 318, 277

Rogers, T. M. 2015, ApJ, 815, L30

Rogers, T. M., Lin, D. N. C., McElwaine, J. N., & Lau, H. H. B. 2013, ApJ, 772, 21

Roxburgh, I. W. 2006, A&A, 454, 883

Saio, H., Kurtz, D. W., Takata, M., et al. 2015, MNRAS, 447, 3264

Saio, H., Bedding, T. R., Kurtz, D. W., et al. 2018, MNRAS, 477, 2183

Saio, H., Takata, M., Lee, U., Li, G., & Van Reeth, T. 2021, MNRAS, 502, 5856

Schmid, V. S., & Aerts, C. 2016, A&A, 592, A116

Sowicka, P., Handler, G., Dębski, B., et al. 2017, MNRAS, 467, 4663

Szewczuk, W., & Daszyńska-Daszkiewicz, J. 2018, MNRAS, 478, 2243

Szewczuk, W., Walczak, P., & Daszyńska-Daszkiewicz, J. 2021, MNRAS, 503, 5894

Tayar, J., Beck, P. G., Pinsonneault, M. H., García, R. A., & Mathur, S. 2019, ApJ, 887, 203

Triana, S. A., Moravveji, E., Pápics, P. I., et al. 2015, ApJ, 810, 16

Triana, S. A., Corsaro, E., De Ridder, J., et al. 2017, A&A, 602, A62

Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial Oscillations of Stars (Tokyo: University of Tokyo Press)

Van Beeck, J., Prat, V., Van Reeth, T., et al. 2020, A&A, 638, A149

Van Reeth, T., Tkachenko, A., Aerts, C., et al. 2015, A&A, 574, A17

Van Reeth, T., Tkachenko, A., & Aerts, C. 2016, A&A, 593, A120

Van Reeth, T., Mombarg, J. S. G., Mathis, S., et al. 2018, A&A, 618, A24

Wang, H., Boyd, J. P., & Akmaev, R. A. 2016, Geosci. Model Dev., 9, 1477

Zahn, J. P. 1992, A&A, 265, 115