Issue 
A&A
Volume 631, November 2019



Article Number  A26  
Number of page(s)  12  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201935639  
Published online  15 October 2019 
The traditional approximation of rotation, including the centrifugal acceleration for slightly deformed stars
AIM, CEA, CNRS, Université ParisSaclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 GifsurYvette Cedex, France
email: stephane.mathis@cea.fr
Received:
8
April
2019
Accepted:
15
August
2019
Context. The traditional approximation of rotation (TAR) is a treatment of the dynamical equations of rotating and stably stratified fluids in which the action of the Coriolis acceleration along the direction of the entropy (and chemicals) stratification is neglected, while assuming that the fluid motions are mostly horizontal because of their inhibition in the vertical direction by the buoyancy force. This leads to the neglect of the horizontal projection of the rotation vector in the equations for the dynamics of gravitoinertial waves (GIWs) that become separable, such as in the nonrotating case, while they are not separable in the case in which the full Coriolis acceleration is taken into account. This approximation, first introduced in geophysical fluid dynamics for thin atmospheres and oceans, has been broadly applied in stellar (and planetary) astrophysics to study lowfrequency GIWs that have short vertical wavelengths. The appoximation is now being tested thanks to direct 2D oscillation codes, which constrain its domain of validity. The mathematical flexibility of this treatment allows us to explore broad parameter spaces and to perform detailed seismic modelling of stars.
Aims. The TAR treatment is built on the assumptions that the star is spherical (i.e. its centrifugal deformation is neglected) and uniformly rotating while an adiabatic treatment of the dynamics of the waves is adopted. In addition, their induced gravitational potential fluctuations is neglected. However, it has been recently generalised with including the effects of a differential rotation. We aim to carry out a new generalisation that takes into account the centrifugal acceleration in the case of deformed stars that are moderately and uniformly rotating.
Methods. We construct an analytical expansion of the equations for the dynamics of GIWs in a spheroidal coordinates system by assuming the hierarchies of frequencies and amplitudes of the velocity components adopted within TAR in the spherical case.
Results. We derive the complete set of equations that generalises TAR by taking the centrifugal acceleration into account. As in the case of a differentially rotating spherical star, the problem becomes 2D but can be treated analytically if we assume the anelastic and JWKB approximations, which are relevant for lowfrequency GIWs. This allows us to derive a generalised Laplace tidal equation for the horizontal eigenfunctions and asymptotic wave periods, which can be used to probe the structure and dynamics of rotating deformed stars thanks to asteroseismology. A first numerical exploration of its eigenvalues and horizontal eigenfunctions shows their variation as a function of the pseudoradius for different rotation rates and frequencies and the development of avoided crossings.
Key words: hydrodynamics / methods: analytical / stars: oscillations / stars: rotation / stars: interiors
© S. Mathis and V. Prat 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 traditional approximation of rotation (TAR) was first introduced to study the dynamics of the shallow Earth atmosphere and oceans (e.g. Eckart 1960; Gerkema et al. 2008; Zeitlin 2018). In particular, it is broadly used for the understanding of the propagation of waves in stably stratified (in entropy and chemical composition) and rotating atmospheric and oceanic layers. There, inertiagravity waves, which are often called gravitoinertial waves (GIWs) in stellar physics (e.g. Dintrans et al. 1999), propagate under the combined action of buoyancy force and Coriolis acceleration. Assuming that the buoyancy force is stronger than the Coriolis acceleration (i.e. 2Ω ≪ N, where Ω is the angular velocity and N the BruntVaïsälä frequency) in the direction of stable entropy or chemical stratification, this leads to waves whose horizontal velocities are higher than their vertical velocities. This allows us to neglect the terms involving the latitudinal component of the rotation vector Ω_{H} = Ωsinθ, where θ is the colatitude, in the linearised hydrodynamical equations. In this framework, the (Poincaré) wave equation, which is 2D and nonseparable in the general case, becomes separable (e.g. Gerkema & Shrira 2005). Therefore, scalar quantities (e.g. the fluctuations of density, pressure, temperature, and entropy of the waves) and the velocity components can be expressed as products of radial functions; special latitudinal functions, the socalled Hough functions that reduce to Legendre polynomials in the nonrotating case (Hough 1898; LonguetHiggins 1968); and Fourier series in time and azimuth.
In stellar physics, TAR and its flexibility has been extensively used to study lowfrequency GIWs, including Rossby waves, both in single and double stars (e.g. Berthomieu et al. 1978; Lee & Saio 1987, 1989, 1997; Townsend 2003; Mathis 2009). This treatment allows stellar physicists to derive powerful seismic diagnoses for the period spacing between consecutive high radial order gravitoinertial modes in uniformly and differentially rotating spherical stars (Bouabid et al. 2013; Ouazzani et al. 2017; Van Reeth et al. 2018). With the advent of space asteroseismology using high precision photometry (Aerts et al. 2010), the period spacing derived within the TAR framework gives access to properties of chemical stratification and to the rotation rate near the convective core of rapidly rotating intermediatemass γ Doradus stars for more than 60 stars (Van Reeth et al. 2015a,b, 2016, 2018; Aerts et al. 2017; Christophe et al. 2018; Li et al. 2019). These quantities constitute key ingredients to improve our knowledge of stellar structure, evolution, internal angular momentum transport, and to calibrate stellar models (e.g. Pedersen et al. 2018; Aerts et al. 2018, 2019; Ouazzani et al. 2019, and references therein). However, in addition to assuming that 2Ω ≪ N, other hypotheses are carried out to apply TAR in stellar interiors (e.g. Lee & Saio 1997; Townsend 2003). First, the rotation is assumed to be uniform. Next, we assume that the centrifugal distortion of the star can be neglected, i.e. , where G, M, and R are the universal constant of gravity, mass of the star, and stellar radius, respectively, and Ω_{K} is the Keplerian critical angular velocity. The fluctuation of the gravitational potential of the waves is neglected following Cowling (1941). Finally, the motions of the waves are assumed to be adiabatic.
Different approaches can be considered to go beyond these assumptions. We can choose to use 2D oscillation codes (Reese et al. 2006; Ballot et al. 2010; Ouazzani et al. 2012). However, they are not available yet to the whole community while their needed computation time and resources can prevent detailed seismic modelling at this stage. We can alternatively build a 2D asymptotic theory for GIWs that go beyond the TAR (Prat et al. 2016, 2018). However, the associated derivation of the needed asymptotic seismic diagnosis is still in its infancy (Prat et al. 2017) and should be developed. Finally, in the case of strongly stratified stellar radiation zones for which 2Ω ≪ N, we can try to improve TAR by using the hierarchy between frequencies and the corresponding properties of motions. This has been done successfully to include the effects of general differential rotation on lowfrequency GIWs (Ogilvie & Lin 2004; Mathis 2009). In this case, the problem becomes 2D and nonseparable as in the general case in which the full Coriolis acceleration is taken into account even in the case of a “shellular” rotation that depends on the radius alone. However, the problem can be treated analytically in the case of lowfrequency GIWs following Ogilvie & Lin (2004) and Mathis (2009) who considered rapidly oscillating waves in the vertical direction by assuming the anelastic approximation where acoustic waves are filtered out. These authors introduced generalised 2D Hough functions that depend on the latitude and radius, the latter acting only as a parameter. The results of this method have been successfully applied in Van Reeth et al. (2018) to derive the variation of the asymptotic period spacing in the case of a weak radial differential rotation as observed in intermediatemass stars using asteroseismology (Kurtz et al. 2014; Saio et al. 2015; Murphy et al. 2016; Aerts et al. 2017).
In this theoretical work, we consider the case of “moderately” rapidly rotating stars (or planets). Their shape becomes a slightly deformed spheroid because of the action of the centrifugal acceleration. This case has been studied using 2D oscillation modes numerical computation in Ouazzani et al. (2017). These authors demonstrated that the results obtained using TAR and the related assumptions (i.e. studying uniformly rotating spherical stars) are in qualitative agreement with the complete treatment of the Coriolis acceleration when using 1D spherical structure models (see their Figs. 1 and 3). In addition, they showed that this latter treatment is also in qualitative agreement with a 2D treatment that takes the full Coriolis and centrifugal accelerations into account^{1} (see their Fig. 6), but with some quantitative differences that should be explained. Indeed, while the period spacings of individual modes is different, the global properties such as the mean value of the period spacing, the number of modes, the extend and slope of the pattern for each group of fixed azimuthal order are similar. This could be understood by invoking that studied lowfrequency gravitoinertial modes are propagating in (and sounding) deep stellar layers, which are less influenced by the centrifugal acceleration than the surface (e.g. Ballot et al. 2010). In this framework, it becomes important to study if the TAR could be generalised to take into account the effects of the centrifugal acceleration in the case of slightly deformed moderately rotating stars (this work) before considering the more extreme cases of stars rotating close to their breakup velocity. This could have several key applications such as new seismic diagnosis, the study of the transport of angular momentum by GIWs (e.g. Lee & Saio 1993; Mathis et al. 2008; Mathis 2009; Lee et al. 2014), and the evaluation of tidal dissipation (Ogilvie & Lin 2004, 2007; Braviner & Ogilvie 2014) in deformed stars (and planets).
In this work, we focus on the first goal. First, we introduce in Sect. 2 the formalism already introduced in the literature to study oscillation modes in moderately deformed rotating stars (or planets; e.g. Smeyers & Denis 1971; Saio 1981; Lee 1993; Lee & Baraffe 1995). We consider the simplest case of a uniform rotation to disentangle the effects of the deformation from those of differential rotation. In Sect. 3, we show how to generalise TAR in this configuration and in Sect. 4 we study the dynamics of corresponding lowfrequency GIWs. In Sect. 5, we derive the periods of GIWs, which can be used to probe the chemical composition of stars and their internal rotation. In Sect. 6, we use the new formalism to determine how the solutions of the Laplace tidal equation are affected by the centrifugal deformation. Finally, we give in Sect. 7 the conclusions of this theoretical work and we discuss its perspectives and future applications.
2. Dynamical equations in moderately deformed stars
We follow the formalism presented by Lee & Baraffe (1995) (see also Smeyers & Denis 1971; Saio 1981; Lee 1993) to describe the dynamics of stellar oscillation modes in slightly deformed stars. We work in a spheroidal system of coordinates (a,θ,φ) to take into account the centrifugal deformation of the star. The origin of this coordinate system (a = 0) is the centre of the deformed star, θ is the colatitude with θ = 0 on the rotation axis, and φ is the azimuth. The spheroidal coordinates system is related to the spherical coordinates (r,θ,φ) by a mapping
where ε is a function describing the centrifugal perturbation of the hydrostatic balance. Appendix A provides its derivation. This perturbation scales as Ω^{2}, where Ω is the angular velocity of the star, which is assumed to be uniform. We define the unitvector basis attached to the spheroidal coordinates system as follows:
The dynamical equations for waves are derived in these spheroidal coordinates. As in studies of TAR in spherical stars, we focus in this work on adiabatic oscillations.
First, the linearised NavierStockes equation is written as
where
and
We introduce the Lagrangian displacement ξ of the oscillation, gravific potential Φ, pressure P, and density ρ. The scalar quantities (Φ, P, ρ) are expended as the sum of their hydrostatic components () and their wave fluctuations (). Finally, the Coriolis acceleration operator is given by
where
We adopt the Cowling approximation in which the perturbations of the gravific potential induced by pulsations are neglected (Cowling 1941).
Next, the linearised continuity equation is obtained
The linearised energy equation in the adiabatic limit is derived
where (S being the macroscopic entropy) is the adiabatic exponent. It allows us to identify the squared BruntVaïsälä frequency
Finally, the wave’s displacement and fluctuations () are expanded on Fourier series both in time and in azimuth
where ω is the frequency in the corotating frame and m the azimutal degree.
3. The TAR with centrifugal acceleration
We consider each component of the momentum equation to identify the hierarchy of the different terms and the corresponding simplifications when using the TAR in deformed stars.
First, the spheroidal radial component can be written as
where and we explain the frequency ratios ω/N and 2Ω/N. When using TAR, we focus on lowfrequency waves for which ω ≪ N. Their propagation can be studied within the anelastic approximation in which acoustic waves are filtered out. Equations (13) and (9) can be simplified accordingly by neglecting the terms and , respectively (more details are provided in Sects. 2 and 3 in Mathis 2009). In addition, TAR can be applied only in the case of “strong” stratification for which 2Ω ≪ N. In this case, the buoyancy force dominates the radial components of the Coriolis acceleration and of the acceleration of the wave. Therefore, the radial momentum equation can be simplified to
as in the case of spherical stars and for the same reasons.
Next, we examine the latitudinal component of the momentum equation, which we write
In the case of lowfrequency GIWs, the wave displacement is mostly horizontal, i.e. ξ_{a} ≪ {ξ_{θ}, ξ_{φ}}, because of the strong stable stratification. This allows us to neglect the term (ξ_{a}/ξ_{θ})∂_{θ}ε that couples horizontal and vertical directions. The previous equation simplifies onto
where
In the case of a uniform rotation ε(a,θ) = ε_{0}(a) + ε_{2}(a)P_{2}(cosθ) (see Eq. (A.24)). Therefore, ∂_{θ}ε ∝ cosθsinθ and the term tanθ∂_{θ}ε is regular.
For the same reasons, the azimuthal component of the momentum equation
reduces to
As in the spherical case, we thus obtain decoupled equations for the vertical and horizontal components of the displacement. We can thus solve the system formed by Eqs. (16) and (20) and express and as a function of the normalised pressure W′ as follows:
with
We can identify that the structure of the equations in the spheroidal case is very similar to those in the usual spherical case (e.g. Lee & Saio 1997; Townsend 2003). The most important difference is that the coefficients 𝒜, ℬ, 𝒞, and 𝒟 are function of a and θ through ε(a,θ), while they reduce to functions that only depend on θ in the spherical case with 𝒜 = ℬ = 𝒞 = 1 and 𝒟 = 1 − ν^{2} cos^{2}θ. This situation is similar to the case in which differential rotation is taken into account (Mathis 2009; Van Reeth et al. 2018) and we see in the next section that it is possible to solve the problem by adopting the same method. We identify, as in the uniformly rotating spherical case, the socalled spin parameter ν = 2Ω/ω. The regime ν > 1 (ν < 1) corresponds to sub (super) inertial waves in which ω < 2Ω (ω > 2Ω).
Following Lee & Saio (1997) and Mathis (2009), we introduce the reduced latitudinal coordinate x = cosθ; Eqs. (16) and (20) transform into
4. Dynamics of lowfrequency gravitoinertial waves
As in Mathis (2009) and Van Reeth et al. (2018), we focus from now on lowfrequency GIWs. Our goal is to derive the socalled Poincaré equation for the normalised pressure (W′), which allows us to compute their frequencies and periods and to built the corresponding seismic diagnosis. Using the anelastic approximation again, the continuity equation becomes , where we recall the relation u = ∂_{t} ξ = iω ξ between the velocity (u) and the Lagrangian displacement (ξ). The continuity equation (Eq. (8)) simplifies into
where
Lowfrequency GIWs are rapidly oscillating with short wavelengths along the (vertical) direction, which are very small compared to the characteristic lengths of variation of the background quantities. This allows us, following Mathis (2009), to use the JeffreysWentzelKramersBrillouin (JWKB) approximation (Fröman & Fröman 2005) along the vertical and to expand the pressure fluctuation and the components of the displacement as
where j ≡ {r,θ,φ}, k is the index of a latitudinal eigenmode (see Eq. (35)), and A_{νkm} is the amplitude of the wave. Using Eqs. (14), (27), and (28), we can obtain the following polarisation relations:
Using the vertical momentum equation Eqs. (14) and (29) becomes
where the JWKB approximation allows us to neglect in front of the dominant term .
Using Eqs. (27) and (28) allows us to obtain the equation for w_{νkm}
where we identify the dispersion relation for lowfrequency GIWs within TAR
This equation is in fact the Poincaré partial differential equation for GIWs. The combined use of the TAR and JWKB approximations, allows us to transform it in a linear secondorder ordinary differential equation on x with only a parametric dependence on a. This result is very similar to those obtained in the case of differential rotation in Ogilvie & Lin (2004), Mathis (2009), and Van Reeth et al. (2018). Therefore, Eq. (35) can be seen as a generalised Laplace tidal equation (and operator) for generalised Hough functions when taking the centrifugal acceleration into account. In Table 1, we recall the expressions of all the involved coefficients. Their firstorder linearisation in ε is derived in Appendix B.1. Since the eigenvalues and eigenfunctions of the generalised Laplace tidal equation vary with the pseudoradius, we choose to define our latitudinal ordering number k as in Lee & Saio (1997) by considering the eigenvalues and eigenfunctions at the centre where they are not affected by the centrifugal acceleration since the mapping we choose in this work (see Eq. (A.26)) is such that ε → 0 for a → 0. We recall that Lee & Saio (1997) ordered the eigenvalues of the nondeformed Hough functions as follows: for the eigenvalues that exist for any value of ν, they attached positive k including zero with Λ_{ν = 0, km}(a = 0) = (m+k)(m+k+1) which corresponds to the l(l+1) eigenvalues of the spherical harmonics (with l = m+k), which are the horizontal eigenfunctions for nonradial pulsations in the nonrotating case (i.e. ν = 0); for the eigenvalues that only exist when ν> 1, they used negative k in such a way that Λ_{ν, −1, m}(a = 0) > Λ_{ν, −2, m}(a = 0) > ⋯. As we see in the next section, avoided crossings can appear when looking at the variation of the eigenvalues along the pseudoradius (see Fig. 9). The corresponding k as defined by Lee & Saio (1997) then changes.
Centrifugal terms involved in the generalised Laplace tidal equation.
As in the uniformly and differentially rotating spherical cases, two classes of waves are identified: those for which 𝒟 > 0 in the whole spheroidal shell and that propagate at all colatitudes and those for which 𝒟 vanishes in the spheroidal cavity at a critical colatitude θ_{c} such that 𝒟(a,θ_{c}) = 0, which depends on a. For this second class, waves propagate only within an equatorial belt where θ > θ_{c}(a) (some examples of these waves can be found in Ballot et al. 2010; Prat et al. 2016, 2018, while an approximate value for θ_{c}(a) is derived in Eq. (B.35)). These first and second classes of waves correspond to the superinertial (ω > 2Ω) and subinertial waves (ω < 2Ω), respectively, in the case of the uniformly rotating spherical case.
This description within the TAR of lowfrequency GIWs propagating in deformed bodies can be applied to the study of the transport of angular momentum they induce (Mathis 2009), tidal dissipation in stably stratified stellar and planetary layers (Ogilvie & Lin 2004, 2007; Braviner & Ogilvie 2014; Fuller et al. 2016), and the seismology of rapidly rotating stars (e.g. Van Reeth et al. 2018).
In this framework, it is finally interesting to derive the asymptotic frequencies of lowfrequency GIWs and the corresponding periods as in Van Reeth et al. (2018). Indeed, the latter allow asteroseismologist to probe the internal rotation of regions where these waves propagate, for instance the radiative layers close to the convective core of intermediatemass stars (Van Reeth et al. 2016, 2018; Ouazzani et al. 2017).
5. Asymptotic seismic diagnosis
Following the same method that is used in the case of spherical uniformly and differentially rotating stars, we can derive the eigenfrequencies of lowfrequency GIWs by doing a vertical (i.e. in this equation along a) quantisation as follows:
where we introduce n the vertical order, while a_{t1} and a_{t2} are the inner and outer turning point for which the BruntVäisälä frequency (N) vanishes (Berthomieu et al. 1978; Tassoul 1980; Bouabid et al. 2013). Using the previously derived dispersion relation (Eq. (36)), we get the asymptotic expression for the frequencies of lowfrequency GIWs
and the corresponding period
As in the spherical case, we can thus compute the period spacing ΔP = P_{nkm} − P_{n − 1km}, which allows us to probe the internal rotation of stars but taking into account the flattening of stars by the centrifugal acceleration. The slight differences are the use of spheroidal coordinates and the variation with a of the eigenvalues Λ_{νkm} of the generalised Laplace tidal operator.
The interest of the equations derived for the deformation of the star (Eq. (A.26)), the generalised Hough functions (Eq. (35)), and the eigenfrequencies and periods (Eqs. (38) and (39)) is that they can be easily integrated for a large number of stars. This is a great asset when it is necessary to compute complete grids of stellar models to perform detailed seismic modellings (e.g. Pedersen et al. 2018, and Fig. 1).
Fig. 1. Methodology for seismic modelling of rotating deformed stars using the TAR. Red and blue boxes are related to the structure of stars and lowfrequency GIWs, respectively. 
6. Application
In this section, we numerically solve the linearised and generalised Laplace tidal equation at fixed frequency, as presented in Appendix B.1. We consider a ZAMS 1.5 M_{⊙} stellar model with a solar metallicity computed with the MESA 1D stellar evolution code (Paxton et al. 2018). First, we compute ε from the perturbation of the gravitational potential ϕ′, following the procedure described in Appendix A. Equation (A.17) implies that ϕ′ is proportional to Ω^{2} and can thus be normalised by R^{2}Ω^{2}, as represented in Fig. 2. The deformation function ε(a, θ) is illustrated in Fig. 3 for Ω/Ω_{K} = 0.2. For this rotation rate, its maximum absolute value is of the order of 2%. In addition, we observe that ε < 0. This can be easily understood since by definition (see Eq. (1)) r = a(1+ε), where r is the usual spherical radius and a the pseudoradius with a > r because of the action of the centrifugal acceleration. Since ε scales with the square of the rotation rate, this value is slightly lower than 10% for Ω/Ω_{K} = 0.4. A priori, the perturbative approach is thus valid for these rotation rates.
Fig. 2. Normalised pertubation of the gravitational potential ϕ′(a, θ)/(R^{2}Ω^{2}). 
Fig. 3. Deformation function ε(a, θ) for Ω/Ω_{K} = 0.2. 
We then solve the generalised Laplace tidal equation for different pseudoradii, spin factors, and rotation rates using an implementation based on Chebyshev polynomials similar to that by Wang et al. (2016). At the centre (a = 0), this is equivalent to solving the unperturbed classical Laplace tidal equation. The corresponding spectrum as a function of the spin factor is shown for m = −2 in Fig. 4. As expected, this is consistent with Fig. 1 of Lee & Saio (1997). This spectrum features Rossbylike solutions which have mostly negative eigenvalues and are nonexistent for ν < 1, and gravitylike solutions which have positive eigenvalues. A notable property of the Rossbylike part of the spectrum for ν > 1 is that every odd solution can be associated with an even solution that has a very close eigenvalue.
Fig. 4. Spectrum of the generalised Laplace tidal equation as a function of the spin factor ν at a = 0 for Ω = 0.2Ω_{K} and m = −2. Blue (respectively, orange) dots correspond to even (respectively, odd) eigenfunctions. 
Spectra at different pseudoradii (illustrated in Fig. 5) show significant effects of the centrifugal deformation. First, Rossbylike solutions for ν < −1 are now also grouped by two (one odd solution and one even solution) with very close eigenvalues. Second, two previously grouped Rossbylike solutions for ν > 1 now have very different behaviours and probably cause avoided crossings with other Rossbylike solutions. Finally, the behaviour of gravitylike solutions near ν = 1 becomes less regular. As the rotation rate increases, the last two effects become stronger. This can be seen in Fig. 6 for Ω/Ω_{K} = 0.4. In particular, the two previously mentioned Rossbylike solutions almost behave as gravitylike solutions for ν > 1.
Fig. 6. Same as Fig. 4, but for Ω = 0.4Ω_{K} and at different pseudoradii: a = 0.5R (left) and a = R (right). 
We focus now on how the modified Hough function w_{νkm} varies with the pseudoradius a and horizontal coordinate x. This dependence is illustrated for m = −2 and ν = 2 in Fig. 7 at two different rotation rates. Again, the solution at a = 0 is not perturbed by the centrifugal acceleration, and is thus the same for all rotation rates. As the pseudoradius increases, Rossbylike solutions are slightly modified, whereas gravitylike solutions drastically change. However, both types of solutions seem to experience an avoided crossing, which is visible through the change in the number of nodes. This behaviour is qualitatively the same for the two rotation rates studied in this work. Quantitatively, the dispersion of the eigenfunctions as a function of the pseudoradius is larger for larger rotation rates. The same eigenfunctions are plotted in Fig. 8 for m = −2 and ν = 10. Again, gravitylike solutions clearly experience an avoided crossing, but this is no longer visible for Rossbylike solutions.
Fig. 7. Solutions of the generalised Laplace tidal equation at different pseudoradii from a = 0 (dark blue) to a = R (yellow) for m = −2, ν = 2 and Ω/Ω_{K} = 0.2 (left), and Ω/Ω_{K} = 0.4 (right). Dashed lines correspond to gravitylike solutions with k = 0 (and an avoided crossing with the k = 2 mode), whereas solid lines correspond to Rossbylike solutions with k = −4 (and an avoided crossing with the k = −2 mode). 
The avoided crossings can be highlighted by plotting the eigenvalues as a function of the pseudoradius (see Fig. 9). For m = −2 and ν = 2, the avoided crossing of Rossbylike solutions is clearly visible. In contrast, that of gravitylike solutions is not because it occurs at very small pseudoradii. For m = −2 and ν = 10, it occurs at larger pseudoradii and is thus visible.
Fig. 9. Spectrum of the generalised Laplace tidal equation as a function of the pseudoradius at Ω = 0.2Ω_{K} and m = −2 for ν = 2 (left) and ν = 10 (right). Blue (respectively, orange) dots correspond to even (respectively, odd) eigenfunctions. 
We finally compute the integral
at a fixed spin factor for different rotation rates to investigate how mode frequencies could be affected by the centrifugal acceleration. The obtained values are shown in Table 2. The value of the integral slightly increases with the rotation rate, which has an effect on the quantisation condition Eq. (37) and thus on the mode frequencies.
In addition to these frequency shifts, the changes in the wave structure induced by the centrifugal deformation observed in this section may also have an impact on the transport of angular momentum and on the tidal dissipation they induce (e.g. Mathis 2009; Braviner & Ogilvie 2015).
7. Conclusions
In this theoretical article, we generalise the TAR to the case in which the deformation of a star (or planet) by the centrifugal acceleration is taken into account. We identify that the mathematical complexity introduced by the centrifugal acceleration is very similar to that appearing when applying the TAR to differentially rotating spherical stars (Ogilvie & Lin 2004; Mathis 2009; Van Reeth et al. 2018). Combining the TAR with the anelastic and JWKB approximations, we derive a generalised tidal Laplace equation, which is a secondorder linear ordinary differential equation in x = cosθ (θ being the colatitude) only with a parametric dependence on a of its coefficients. The problem thus reduces to a classical SturmLiouville problem as in the case of uniformly rotating spherical stars. It allows us to derive the asymptotic frequencies of lowfrequency GIWs and the corresponding periods and period spacings. These can be used as a seismic probe of stellar interiors and rotation in moderately rapidly rotating deformed stars. In addition, the derived formalism can be used to study the angular momentum transport and tidal dissipation induced by lowfrequency GIWs in stars and planets. We carried out a first numerical exploration of the eigenvalues and horizontal eigenfunctions of the generalised Laplace tidal equation following the methodology presented in Fig. 1. We find that both gravity and Rossbylike wave eigenfunctions are affected by the centrifugal acceleration and vary with the pseudoradii with a stronger deformation of the gravitylike solutions when compared to the spherical case. In this context, we see that both types of solutions are affected by avoided crossing phenomena. The next step will be to implement our equations in stellar evolution and oscillation codes (as we did in the differentially rotating case in Van Reeth et al. 2018) by comparing the obtained results with direct computations with 2D oscillation codes (Reese et al. 2006; Ballot et al. 2010; Ouazzani et al. 2012, 2017) and to examine if the TAR can be generalised to the case of strongly deformed stars.
The centrifugal acceleration is taken into account in the 2D oscillation code ACOR (Ouazzani et al. 2012), while the 2D rotating stellar model is built by deforming a 1D initial spherical model following the iterative method proposed in Roxburgh (2006).
Acknowledgments
We thank the referee for her/his constructive report which has allowed us to improve our work. S. M. and V. P. acknowledge support by ERC SPIRE 647383 and by CNES PLATO & GOLF grants at CEASaclay. The latter authors thank C. Aerts and C. Neiner for fruitful exchanges and K. Augustson for providing the 1.5 M_{⊙} stellar model used in this work. S. M. dedicates this work to the memory of his fatherinlaw and friend J. Neiner.
References
 Aerts, C., ChristensenDalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology (Springer Science+Business Media B.V.) [Google Scholar]
 Aerts, C., Van Reeth, T., & Tkachenko, A. 2017, ApJ, 847, L7 [NASA ADS] [CrossRef] [Google Scholar]
 Aerts, C., Molenberghs, G., Michielsen, M., et al. 2018, ApJS, 237, 15 [NASA ADS] [CrossRef] [Google Scholar]
 Aerts, C., Mathis, S., & Rogers, T. 2019, ARA&A, 57, 35 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Alvan, L., Strugarek, A., Brun, A. S., Mathis, S., & Garcia, R. A. 2015, A&A, 581, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ballot, J., Lignières, F., Reese, D. R., & Rieutord, M. 2010, A&A, 518, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Berthomieu, G., Gonczi, G., Graff, P., Provost, J., & Rocca, A. 1978, A&A, 70, 597 [Google Scholar]
 Bouabid, M.P., Dupret, M.A., Salmon, S., et al. 2013, MNRAS, 429, 2500 [NASA ADS] [CrossRef] [Google Scholar]
 Braviner, H. J., & Ogilvie, G. I. 2014, MNRAS, 441, 2321 [NASA ADS] [CrossRef] [Google Scholar]
 Braviner, H. J., & Ogilvie, G. I. 2015, MNRAS, 447, 1141 [NASA ADS] [CrossRef] [Google Scholar]
 Chandrasekhar, S. 1933, MNRAS, 93, 390 [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]
 CohenTannoudji, C., Diu, B., & Laloe, F. 1986, Quantum Mech., 2, 626 [NASA ADS] [Google Scholar]
 Cowling, T. G. 1941, MNRAS, 101, 367 [NASA ADS] [Google Scholar]
 Dintrans, B., Rieutord, M., & Valdettaro, L. 1999, J. Fluid Mech., 398, 271 [NASA ADS] [CrossRef] [Google Scholar]
 Eckart, C. 1960, Hydrodynamics of Oceans and Atmospheres (Pergamon Press) [Google Scholar]
 Fröman, N., & Fröman, P. O. 2005, Physical Problems Solved by the PhaseIntegral Method (Cambridge, UK: Cambridge University Press), 228 [Google Scholar]
 Fuller, J., Luan, J., & Quataert, E. 2016, MNRAS, 458, 3867 [NASA ADS] [CrossRef] [Google Scholar]
 Gerkema, T., & Shrira, V. I. 2005, J. Fluid Mech., 529, 195 [NASA ADS] [CrossRef] [Google Scholar]
 Gerkema, T., Zimmerman, J. T. F., Maas, L. R. M., & van Haren, H. 2008, Rev. Geophys., 46, RG2004 [NASA ADS] [CrossRef] [Google Scholar]
 Hough, S. S. 1898, Philos. Trans. R. Soc. London Ser. A, 191, 139 [Google Scholar]
 Kurtz, D. W., Saio, H., Takata, M., et al. 2014, MNRAS, 444, 102 [NASA ADS] [CrossRef] [Google Scholar]
 Lee, U. 1993, ApJ, 405, 359 [NASA ADS] [CrossRef] [Google Scholar]
 Lee, U., & Baraffe, I. 1995, A&A, 301, 419 [NASA ADS] [Google Scholar]
 Lee, U., & Saio, H. 1987, MNRAS, 224, 513 [NASA ADS] [CrossRef] [Google Scholar]
 Lee, U., & Saio, H. 1989, MNRAS, 237, 875 [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]
 Li, G., Bedding, T. R., Murphy, S. J., et al. 2019, MNRAS, 482, 1757 [NASA ADS] [CrossRef] [Google Scholar]
 LonguetHiggins, M. S. 1968, Philos. Trans. R. Soc. London Ser. A, 262, 511 [Google Scholar]
 Mathis, S. 2009, A&A, 506, 811 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathis, S., & Zahn, J.P. 2004, A&A, 425, 229 [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]
 Murphy, S. J., Fossati, L., Bedding, T. R., et al. 2016, MNRAS, 459, 1201 [NASA ADS] [CrossRef] [Google Scholar]
 Ogilvie, G. I., & Lin, D. N. C. 2004, ApJ, 610, 477 [NASA ADS] [CrossRef] [Google Scholar]
 Ogilvie, G. I., & Lin, D. N. C. 2007, ApJ, 661, 1180 [NASA ADS] [CrossRef] [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. J., et al. 2019, A&A, 626, A121 [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]
 Pedersen, M. G., Aerts, C., Pápics, P. I., & Rogers, T. M. 2018, A&A, 614, A128 [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]
 Prat, V., Mathis, S., Lignières, F., Ballot, J., & Culpin, P.M. 2017, A&A, 598, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Prat, V., Mathis, S., Augustson, K., et al. 2018, A&A, 615, A106 [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]
 Roxburgh, I. W. 2006, A&A, 454, 883 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Saio, H. 1981, ApJ, 244, 299 [NASA ADS] [CrossRef] [Google Scholar]
 Saio, H., Kurtz, D. W., Takata, M., et al. 2015, MNRAS, 447, 3264 [NASA ADS] [CrossRef] [Google Scholar]
 Smeyers, P., & Denis, J. 1971, A&A, 14, 311 [NASA ADS] [Google Scholar]
 Sweet, P. A. 1950, MNRAS, 110, 548 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Tassoul, J.L. 1978, Theory of Rotating Stars (Princeton: University Press) [Google Scholar]
 Tassoul, M. 1980, ApJS, 43, 469 [NASA ADS] [CrossRef] [Google Scholar]
 Townsend, R. H. D. 2003, MNRAS, 340, 1020 [NASA ADS] [CrossRef] [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]
 Wang, H., Boyd, J. P., & Akmaev, R. A. 2016, Geosci. Model Dev., 9, 1477 [NASA ADS] [CrossRef] [Google Scholar]
 Zahn, J. P. 1966, Ann. Astrophys., 29, 313 [NASA ADS] [Google Scholar]
 Zahn, J.P. 1992, A&A, 265, 115 [NASA ADS] [Google Scholar]
 Zeitlin, V. 2018, Geophysical Fluid Dynamics: Understanding (almost) Everything with Rotating Shallow Water Models (OUP Oxford) [Google Scholar]
Appendix A: Deformation of a moderately rotating star
The objective of this appendix is to determine the deformation of an isobar in the case of a moderately and uniformly rotating star in which the centrifugal acceleration is a linear perturbation of the order of ; Ω_{K} is the Keplerian critical angular velocity, M and R are the stellar mass and radius, respectively, and G is the universal constant of gravity.
The first step is to calculate the perturbation of the gravitational potential ϕ on the sphere of radius r. We expand ϕ on the orthogonal basis formed by the Legendre polynomials
where ϕ_{0} = −GM_{r}(r)/r is the gravitational potential of the nonrotating star (M_{r} being the mass inside the sphere of radius r) and ϕ_{1} its perturbation induced by the centrifugal acceleration. We follow the method of linearisation of the hydrostatic balance developed in Sweet (1950), Zahn (1966, 1992), and Mathis & Zahn (2004), and expand the pressure and the density around the sphere in the same way that ϕ
Then, we take the hydrostatic equation
is the centrifugal acceleration, which derives from a potential
in the case considered of a uniform rotation. The hydrostatic balance thus becomes
which we expand to the first order as
Taking the curl of Eq. (A.6), we also have
The equipotential for (ϕ+U), the isodensity and isobar thus coincide. As a consequence P can be written as a function of (ϕ+U) as
When linearised to the firstorder, we get
since dℱ/dϕ_{0} = dP_{0}/dϕ_{0} = −ρ_{0}. This leads to
which provides us the perturbation of density
where g_{0} = GM_{r}(r)/r^{2}. Next, we insert the modal expansion of ρ_{1} and those of the centrifugal potential, i.e.
and
This yields the modal amplitude of the density fluctuation over the sphere
We insert this expression in the perturbed Poisson equation ∇^{2}ϕ_{l} = 4πGρ_{l} and we retrieve the Sweet (1950) and Zahn (1966) result
where l = {0,2}. The applied boundary conditions are written as
R being the surface radius of the star (or the planet).
Taking the latitudinal component of the hydrostatic balance (Eq. (A.10)) finally provides us the radial functions of the pressure fluctuation expansion
We introduce the radial coordinate of the isobar
Taking the Taylor expansion of P to first order, we have
By definition the pressure is constant on the isobar. We conclude that
where we used Eq. (A.19) and the zerothorder hydrostatic balance dP_{0}/dr = −ρ_{0}g_{0}.
We finally introduce the pseudoradial coordinate a defined in Eq. (1), i.e.
where ε is also expanded on Legendre polynomials as
In contrast to Lee (1993), we do not use the ChandrasekharMilne expansion (e.g. Chandrasekhar 1933; Tassoul 1978) to compute a because it leads to infinite ε at the centre. Instead, we use a simple mapping such that r is equal to the deformed stellar radius r_{s}(θ) when a = R and that r ≃ a near the centre. The simplest mapping verifying these conditions is
At first order, this leads to
Appendix B: Perturbative analytical solutions
B.1. Propagative waves with fixed frequencies
We consider any propagative waves with nonquantised frequencies (ω) and thus spin parameter (ν = 2Ω/ω). This is for instance the case of tidally excited waves (e.g. Ogilvie & Lin 2004, 2007; Braviner & Ogilvie 2014), where the tidal frequency is fixed by the difference between the angular velocities of the primary and the orbit of the companion, and progressive waves (e.g. Alvan et al. 2015).
As described above in Appendix A, the structure of a moderately rotating body where the centifugal acceleration can be treated as a linear perturbation is the linear combination of the nonrotating structure and a perturbation (ε) of the order of . Therefore, we can make a corresponding linear expansion in ε of the generalised Laplace tidal operator (ℒ_{νkm}) derived in Eq. (35) and of its eigenvalue (Λ_{νkm}) using the linear perturbation theory as in quantum mechanics (CohenTannoudji et al. 1986). We obtain
where
with
is the usual Laplace tidal operator in the spherical case with its Hough eigenfunctions (Θ_{νkm}) and eigenvalues (), and is its firstorder centrifugal correction
with
and
These coefficients are obtained using the firstorder expansion of 𝒞 and 𝒟 given in Table 1 that become
while 𝒜, ℬ, and ℰ are unchanged.
The corresponding eigenvalues and eigenfunctions are derived
with their respective centrifugal perturbations
Using the polarisation relations (Eq. (33)), the same expansion is done for the JWKB amplitudes of the Lagrangian displacement as follows:
where j ≡ {r,θ,φ}.
We get in the vertical direction
where
We obtain for the latitudinal and azimuthal directions
and
where
with
This analytical solution using firstorder perturbative method can be of great interest to study propagative lowfrequency GIWs in moderately rotating stars (their excitation, their propagation, their damping, and the potential angular momentum transport they induce; e.g. Mathis 2009). It can also be used to compute tidal dissipation in slightly deformed stably stratified stellar and planetary layers (Braviner & Ogilvie 2014) such as in the case of Saturn (Fuller et al. 2016).
Finally, the firstorder linear perturbation theory at fixed given frequency can be used to compute the modification of the critical colatitude θ_{c} for which
Using the linearisation of 𝒞 given in Eq. (B.8), this becomes
Posing , where is the classical critical colatitude (e.g. Lee & Saio 1997) and scales with ε, we obtain at the first order in ε
which leads to
In the case of uniform rotation
and thus
Using the expression for as a function of ν, the critical colatitude finally reads
If we make the rough assumption that ε_{2} ≈ U_{2}/(ag_{0}) ≈ Ω^{2}a^{2}/3(GM_{a}/a) (using Eq. (A.26) with neglecting ϕ_{2}), it leads to a slight broadening of the equatorial belt where subinertial GIWs are propagating towards the surface.
B.2. Oscillation eigenmodes
Studying oscillation eigenmodes induces the use of the quantisation as derived in Sect. 5. As discussed in the previous section, it would be relevant to expand the eigenfrequencies and the corresponding spin parameters as a combination of their values in the spherical case and a centrifugal correction as follows:
where and . Such an expansion should then be introduced in the expression of the vertical wave number (k_{V; νkm}; Eq. (36)) and of (Eq. (B.4)), which would lead to complex implicit equations to solve. Therefore, we advocate to solve directly the generalised Laplace tidal equation (Eqs. (35) and (38)) when studying the case of oscillation eigenmodes.
All Tables
All Figures
Fig. 1. Methodology for seismic modelling of rotating deformed stars using the TAR. Red and blue boxes are related to the structure of stars and lowfrequency GIWs, respectively. 

In the text 
Fig. 2. Normalised pertubation of the gravitational potential ϕ′(a, θ)/(R^{2}Ω^{2}). 

In the text 
Fig. 3. Deformation function ε(a, θ) for Ω/Ω_{K} = 0.2. 

In the text 
Fig. 4. Spectrum of the generalised Laplace tidal equation as a function of the spin factor ν at a = 0 for Ω = 0.2Ω_{K} and m = −2. Blue (respectively, orange) dots correspond to even (respectively, odd) eigenfunctions. 

In the text 
Fig. 5. Same as Fig. 4, but at different pseudoradii: a = 0.5R (left) and a = R (right). 

In the text 
Fig. 6. Same as Fig. 4, but for Ω = 0.4Ω_{K} and at different pseudoradii: a = 0.5R (left) and a = R (right). 

In the text 
Fig. 7. Solutions of the generalised Laplace tidal equation at different pseudoradii from a = 0 (dark blue) to a = R (yellow) for m = −2, ν = 2 and Ω/Ω_{K} = 0.2 (left), and Ω/Ω_{K} = 0.4 (right). Dashed lines correspond to gravitylike solutions with k = 0 (and an avoided crossing with the k = 2 mode), whereas solid lines correspond to Rossbylike solutions with k = −4 (and an avoided crossing with the k = −2 mode). 

In the text 
Fig. 8. Same as Fig. 7, but for ν = 10. 

In the text 
Fig. 9. Spectrum of the generalised Laplace tidal equation as a function of the pseudoradius at Ω = 0.2Ω_{K} and m = −2 for ν = 2 (left) and ν = 10 (right). Blue (respectively, orange) dots correspond to even (respectively, odd) eigenfunctions. 

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.