Issue 
A&A
Volume 518, JulyAugust 2010
Herschel: the first science highlights



Article Number  A30  
Number of page(s)  8  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201014426  
Published online  26 August 2010 
Gravity modes in rapidly rotating stars
Limits of perturbative methods
J. Ballot^{1}  F. Lignières^{1}  D. R. Reese^{2}  M. Rieutord^{1}
1  Laboratoire d'Astrophysique de ToulouseTarbes, Université de Toulouse,
CNRS, 14 Avenue E. Belin, 31400 Toulouse, France
2 
LESIA, UMR8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon, France
Received 13 March 2010 / Accepted 3 May 2010
Abstract
Context. CoRoT and Kepler missions are now providing
highquality asteroseismic data for a large number of stars. Among
intermediatemass and massive stars, fast rotators are common objects.
Taking the rotation effects into account is needed to correctly
understand, identify, and interpret the observed oscillation
frequencies of these stars. A classical approach is to consider
the rotation as a perturbation.
Aims. In this paper, we focus on gravity modes, such as those occurring in Doradus, slowly pulsating B (SPB), or Be stars. We aim to define the suitability of perturbative methods.
Methods. With the twodimensional oscillation program (TOP), we
performed complete computations of gravity modes  including the
Coriolis force, the centrifugal distortion, and compressible
effects  in 2D distorted polytropic models of stars. We
started with the modes ,
n=114, and
,
n=15, 1620 of a nonrotating star, and followed these modes by
increasing the rotation rate up to 70% of the breakup rotation
rate. We then derived perturbative coefficients and determined the
domains of validity of the perturbative methods.
Results. Secondorder perturbative methods are suited to computing loworder, lowdegree mode frequencies up to rotation speeds 100 km s^{1} for typical Dor stars or 150 km s^{1} for B stars. The domains of validity can be extended by a few tens of km s^{1}
thanks to the thirdorder terms. For higher order modes, the
domains of validity are noticeably reduced. Moreover, perturbative
methods are inefficient for modes with frequencies lower than the
Coriolis frequency .
We interpret this failure as a consequence of a modification in the
shape of the resonant cavity that is not taken into account in the
perturbative approach.
Key words: asteroseismology  stars: oscillations  stars: rotation  methods: numerical
1 Introduction
CoRoT (Convection, Rotation and planetary Transits, Baglin et al. 2006) and Kepler (Borucki et al. 2007) are space missions providing uninterrupted highquality photometry time series over several months or years ideally suited for asteroseismic study. Asteroseismology provides very accurate determinations of the stellar parameters (mass, radius, age, etc.) and probes stellar structure to constrain physical processes occurring in stars. The first step towards this goal requires correctly understanding the structure of the observed oscillation spectra, and especially correctly identifying the observed modes. In the case of mainsequence (e.g. Benomar et al. 2009; Michel et al. 2008; Chaplin et al. 2010) and giant (e.g. Bedding et al. 2010; Miglio et al. 2009; Hekker et al. 2009) FGK stars with solarlike oscillations, the spectrum structure is well understood, which eases interpretation.
The spectra of classical pulsators is often noticeably more complex. For instance, the highquality observations of Scuti (e.g. Poretti et al. 2009; García Hernández et al. 2009) and Doradus (e.g. Mathias et al. 2009) stars have exhibited very rich and complex spectra of acoustic (p) and gravity (g) modes, respectively, containing several hundred  or more  modes. Interpreting their spectra is very challenging today. Indeed, these stars generally spin rapidly, so the effects of rotation on the mode frequencies must be considered.
Here, we are concerned with gravity modes, i.e., lowfrequency modes driven by the buoyancy force. They are excited and observed in a broad panel of stars, for instance, in Dor, SPB and some Be stars. The Dor stars form a class of mainsequence stars with type around F0V that can sometimes rotate rapidly (e.g. De Cat et al. 2006, and reference therein), while the rotation rate of Be stars is extreme, usually very close to their breakup limit (e.g. Frémat et al. 2006).
The effects of rotation on the oscillation modes can be treated as a perturbation where the rotation rate is the small parameter. A 1storder correction has been proposed by Ledoux (1951), 2ndorder by Saio (1981), Dziembowski & Goode (1992), or Suárez et al. (2006), and 3rdorder terms have been developed by Soufi et al. (1998). While perturbative methods are expected to be accurate enough for slowly rotating stars, their true domain of validity cannot be determined in the absence of exact calculations to compare them with.
In the past few years, calculations of p modes with both the centrifugal distortion and the Coriolis force have been performed in polytropic models of stars (Lignières et al. 2006; Reese et al. 2006) and realistic 2D stellar structures (Reese et al. 2009; Lovekin & Deupree 2008). Lignières et al. (2006) and Reese et al. (2006) have shown that, above , perturbation methods fail to reproduce lowdegree and loworder pmode frequencies ( and ) with the accuracy of CoRoT long runs. The structure of the modes is also drastically modified, and this leads to deep changes in the structure of the pmode spectrum (Lignières & Georgeot 2008; Reese et al. 2008; Lignières & Georgeot 2009).
We used an oscillation code based on Reese et al. (2006) to perform gmode calculations with a complete description of the rotational effects on the modes. In this paper, we focus on the limits of validity for perturbative methods. The models and the method are described in Sect. 2. We then derive the perturbative coefficients from the complete computations (Sect. 3), and compare the results obtained with both methods to determine and discuss the domains of validity for perturbative methods (Sect. 4) before concluding in Sect. 5.
2 Models and methods
We consider fully radiative stars for this work. Since the gravity modes are driven by the buoyancy force, they cannot exist in convective regions. SPB and Dor stars have large radiative zones with a convective core, and even a thin convective envelope for the latter. The effects of convective cores are not considered here, since we are mainly interested in the general behavior of g modes under rotation effects.
2.1 2D stellar models
As in Lignières et al. (2006) and Reese et al. (2006), we approximate the equilibrium structure of rotating stars with selfgravitating uniformlyrotating polytropes. They are described in the corotating frame by the three following equations:
(1)  
(2)  
(3) 
where p_{o} is the pressure, the density, the gravitational potential, K the polytropic constant, the polytropic index, G the gravitational constant, and the effective gravity defined as
(4) 
with s the distance to the rotation axis. Due to the centrifugal distortion, the star is not spherical and a suited surfacefitting spheroidal system of coordinates based on Bonazzola et al. (1998) has been used. Hereafter, we also classically denote r the distance to the center, and z the coordinate along the rotation axis. This equation system is numerically solved with the ESTER (Evolution STEllaire en Rotation) code as described in Rieutord et al. (2005). This is a spectral code using Chebychev polynomials in the direction, and spherical harmonics with even and m=0 in the horizontal one. We computed models decomposed on spherical harmonics up to degree . This resolution is high enough to accurately model the centrifugal effects for the maximal value of that we have considered. In the pseudoradial direction, the resolution is the same as the one we use for the frequency computation (see Sect. 2.3, n_{r}=96 generally).
To approximate a fully radiative star, we chose the polytropic index . We considered models spinning with rotation frequency between 0 and , where is the Keplerian breakup rotation rate for a star of mass M and equatorial radius .
2.2 Linearized equations for the oscillations
In the corotating frame, the equations governing the temporal evolution of small adiabatic inviscid perturbations of the equilibrium structure read, in the corotating frame,
where , p, , and are the Eulerian perturbations of density, pressure, velocity, and gravitational potential the adiabatic sound speed and N_{o} the BruntVäisälä frequency, defined as
denotes the first adiabatic exponent. In Eq. (7) we have used the structure barotropicity, ensured by the uniform rotation.
Figure 1: (Top) Map in the meridional plane of the BruntVäisälä frequency N_{o} for the model with . N_{o} has been normalized by where is the polar radius. Dashed lines show the shape of the spheroidal grid we used. (Bottom) Solid black lines show the profile of N_{o} normalized by for the nonrotating star. Red lines show N_{o} for the model along the polar (solid line) and equatorial radius (dashes). A thick red tick on the xaxis indicates the polar radius . 

Open with DEXTER 
Figure 2: Spectrum as a function of the degree of the radial velocity of the mode at . is the absolute value of the highest Chebychev coefficient in the decomposition on the spherical harmonics of degree of the radial component of . The spectrum is normalized to its maximum value. The different lines correspond to different spatial resolutions. 

Open with DEXTER 
The 2D distribution of the BruntVäisälä frequency is shown in Fig. 1 for the most rapidly rotating model we have considered, together with the profiles of N_{o} along the polar and equatorial radii, which are compared to the N_{o} profile of the nonrotating star. Within about the inner half of the star, the deviations from sphericity induced by the centrifugal force remain limited. We also see that N_{o} diverges at the surface of the polytrope because and p_{o} vanish there.
Looking for timeharmonic solutions of the system (5)(8), we obtain an eigenvalue problem, which we then solve using the twodimensional oscillation program (TOP). The details of this oscillation code closely follow Reese et al. (2006). The equations are projected on the spherical harmonic basis . Due to the axisymmetry of the system, the projected equations are decoupled relatively to the azimuthal order m, but in contrast to the spherical nonrotating case, they are coupled for all degrees of the same parity.
2.3 Resolutions and method accuracy
The different sources of error of our numerical method have been discussed in Valdettaro et al. (2007) and tested in a context similar to the present one in Lignières et al. (2006) and Reese et al. (2006). The numerical resolution has been chosen to ensure a sufficient accuracy for the computed frequencies. In the horizontal direction, the resolution is given by the truncation of the spherical harmonics expansion. The highest degree of the expansion is and, for most of the calculations presented here, we used L = 40+m, i.e. coupled spherical harmonics. In the pseudoradial direction , the solution has been expanded over the set of Chebychev polynomials up to n_{r}=96.
Using higher resolutions (n_{r} = 96, L = 80 and n_{r} = 144, L = 80), we find that the relative agreement of the frequencies always remains better than 5 10^{8}. It also does not affect the mode significantly as illustrated in Fig. 2 where the spectral expansion of the radial velocity component of the mode at is displayed for the three different resolutions. Figure 2 also shows that a unique  or even a few  spherical harmonics would not properly describe such an eigenmode.
2.4 Following modes with rotation
We computed the frequencies of to 3 modes in a nonrotating polytrope. We recall that without rotation the system to solve becomes decoupled with respect to , hence the modes are represented with only one spherical harmonic. A reference frequency set, , was computed from a 1D polytropic model with a radial resolution n_{r}=512.
We then followed the variation in frequency of each mode of degree , azimuthal order m_{0}, and radial order n_{0} by slowly increasing the rotation rate, step by step. The ArnoldiChebychev method requires an initial guess for the frequency, and returns the solutions that are the closest to this guess. The guess we provide is extrapolated from the results at lower rotation rates: we compute from the three last computed points a quadratic extrapolation at the desired rotation rate. For the first point (), we use the frequency obtained in 1D as a guess. For the second point, we extrapolate a guess with the asymptotic relation (Ledoux 1951).
Among the solutions found around the initial guess, we select the correct one by following this strategy:
 1.
 For each calculated mode, we determine, from its spatial spectrum (like the ones shown Fig. 2), the two dominant degrees, and .
 2.
 We compare and with the degree of the mode we are following.
 3.
 We select the solutions such that ; if none of the solutions verifies this criterion, we select the solutions such that .
 4.
 If more than one solution has been selected at this point, we consider the projection of the modes on the spherical harmonic and compare it to the projection of the mode at a lower rotation rate. The solution that gives the highest correlation is finally selected.
Figure 3: Evolution with the rotation rate of the frequency of the modes and , with a zoom on a region where an avoided crossing between these two modes occurs. At a given rotation rate, a diamond (respectively, a triangle) indicates the mode is dominated by the component (). Crosses indicate that the modes are hardly discernible: they are both dominated by and their structures are very similar. 

Open with DEXTER 
The second difficulty comes from the socalled avoided crossings. Two modes with the same m and the same parity cannot have the same frequency. This implies that the two curves associated to their evolution with cannot cross each other. Figure 3 illustrates this phenomenon with modes and : the frequencies get closer and closer, but since the curves cannot cross, the modes exchange their properties. During an avoided crossing, the two modes have the mixed properties of the two initial modes. With our modefollowing method, when the coupling is strong and the avoided crossing takes long, the method can follow the wrong branch. For instance, in the case illustrated in Fig. 3, if the program follows the mode, it continues sometimes on the branch instead of jumping to the other branch.
3 Perturbative coefficients
The approach used to determine the perturbative coefficients in this paper is very close to the one of Reese et al. (2006).
3.1 Determining perturbative coefficients
In the perturbative approach, frequencies are developed as a function of the rotation rate, .
For instance, to the 3rd order, it reads
where is the frequency for the nonrotating case, and the perturbative coefficients. The bar denotes the normalization and . We normalize the frequencies by since the polar radius is expected to be a slowly varying function of in real stars, as opposed to .
The coefficients
can be numerically calculated from the complete computations since they are directly linked to the jth derivative of the function
at .
However, to improve the accuracy, we use symmetry properties of the problem: changing
in ,
one easily shows that
(11) 
We define
and get
We note that vanish for odd j.
We compute and on a grid of k points from to and use the Eqs. (15) and (16) to calculate the terms with the (k1)thorder interpolating polynomials. The determination of the coefficients is then accurate to the (2k1)thorder in . In practice we use a typical step 10^{3} and k=4. In Eq. (14), we use , making it totally independent of the 1D solutions.
By explicitly expressing the dependence on m of the perturbative coefficients, Eq. (10) becomes
The form of the 1st order comes from Ledoux (1951), the 2nd order from Saio (1981), and the 3rd is derived from Soufi et al. (1998). We have verified that the derived coefficients fit these relations with a very good accuracy (see below) and list them in Table 1.
Table 1: Perturbative coefficients (see development Eq. (17)) for g modes with frequency , radial order , and in a polytropic stellar model with an index .
To know the coefficients for another normalization, for instance for
,
one can use the following development:
(18) 
The perturbed frequencies in this new normalization then express
From our models we have computed .
3.2 Coefficient accuracy and comparisons with previous works
The zerothorder coefficients were compared to the 1D computations and we find agreement within 10^{9}. We also compared our results to previous frequency computations of in a nonrotating polytropic model performed by ChristensenDalsgaard & Mullan (1994) with a totally different method. We renormalized their results for g modes (Table 4 of their paper) to their dynamical frequency (Eq. (3.2) of their paper). The relative differences with our results do not exceed 5 10^{8}.
The choice for the step is important for the accuracy of the terms . Ideally, we should choose as small as possible, but when it is too small, the numerical noise, produced by the uncertainties on the computed (Sect. 2.3), drastically increases. We then chose the value of to have the best tradeoff. These uncertainties on determinations were taken into account for the estimated accuracy of the coefficients , and .
The 1storder perturbative coefficients are expressed with integrals of the eigenmodes in the nonrotating model (Ledoux 1951). We then computed these terms with our 1D eigensolutions and compared them to . The results are consistent within 10^{8}.
An explicit computation of 2nd and 3rdorder coefficients requires calculating the 1st and 2ndorder corrections of the eigenfunctions, which is not so straightforward. It is the reason we performed a direct numerical computation of these coefficients. The numerical errors we estimated for and are generally around 10^{5} and always less than 10^{4}. We then checked the consistency of our computations with the 2ndorder calculations of Saio (1981) for g modes with n=1 to 3. In this work, all frequencies were normalized by the dynamical frequency of the nonrotating polytrope. By noticing that with , and using the relation (19), we were able to compare these results with ours. We get a good qualitative agreement with absolute differences better than 10^{2}, which is reasonable relative to the lower accuracy of Saio's computations. It gives an interesting consistency check for our calculations. Overall, the perturbative coefficients listed in Table 1 have been determined with high accuracy.
4 Domains of validity of perturbative approaches
From the previously computed coefficients, we calculated mode frequencies with the 1st to 3rdorder perturbative approximations for rotation rates ranging from to and compared them to complete computations. Figure 4 illustrates such a comparison by showing the evolution of the frequencies of the seven m components of an mode together with their 2ndorder perturbative approximation. We clearly observe that the agreement between both approaches at low rotation progressively disappears as the rotation increases.
Figure 4: Evolution with the rotation rate of the frequencies of the components of the (, n=16) multiplet obtained with a complete computation (solid line) and with the 2ndorder perturbative approximation (dashed lines). 

Open with DEXTER 
To define the domains of validity of perturbative approaches, we fix the maximal departure
allowed between the perturbed frequencies
and the ``exact'' ones
.
For each mode and each approximation order, we define the domain of validity
,
such that
.
The precision of the observed frequencies
can be related to the normalized error
through
From this expression, we see that, for a fixed precision , the normalized error depends on the dynamical frequency of the star considered. We thus display the domains of validity of the perturbative approximations for two types of stars with different dynamical frequencies, a typical Dor star, and a typical B star. The Dor star is such that , , i.e. , while the B star has , , i.e. . For the frequency precision, we chose , which corresponds to the resolution of an oscillation spectrum after a hundred days. This is the typical accuracy for a CoRoT long run. Accordingly, the normalized error is equal to 1.6 10^{3} for the Dor star and 9.3 10^{3} for the B star. It must be noted that the domains of validity will not be affected by the numerical errors of (cf. Sect. 3.2) and (cf. Sect. 2.3) as the values of remain at least an order of magnitude higher.
Figure 5: Evolution of the frequencies of modes (top to bottom). Frequencies are computed in the corotating frame. Perturbative approximations have been tested for a typical Dor (left panels) and for a B star (right panels). Green/red/blue parts of curves indicate that 1st/2nd/3rd order is sufficient to reproduce complete calculations within an error . Error bars on the righthand side of each panel show and 10 . Magenta lines indicate . For each plot, the bottom xaxis and left yaxis show dimensionless units, whereas the top xaxis and right yaxis show physical units. 

Open with DEXTER 
We have determined the domains of validity of 1st, 2nd, and 3rdorder methods for lowdegree modes. Specifically, we considered modes with n=1 to 14, and and 3 loworder (n=1 to 5), and highorder (n=16 to 20) modes. The domains of validity are shown in Fig. 5 for both types of stars. Overall, the domains of validity extend to higher rotation rates for B stars than for Dor stars. This is simply due to the increase in the normalized tolerance . Besides, we observe distinct behaviors in the high and lowfrequency ranges.
In the highfrequency range, 2ndorder perturbative methods give satisfactory results up to 100 km s^{1} for Dor stars and up to 150 km s^{1} for B stars. The 3rdorder terms improve the results and increase the domains of validity by a few tens of km s^{1}. These results are to be contrasted with those found for p modes where the domains of validity are restricted to lower rotation rates. For Scuti stars, which have similar stellar parameters to Dor, Reese et al. (2006) find 5070 km s^{1} as a limit for perturbative methods. In addition, the 3rdorder terms do not improve the perturbative approximation in this case, as p modes are weakly sensitive to the Coriolis force. The rather good performance of perturbative methods at describing highfrequency g modes indicates in particular that the 2ndorder term gives a reasonable description of the centrifugal distortion. This might be surprising considering the significant distortion of the stellar surface ( at ). Actually, the energy of g modes is concentrated in the inner part of the star where the deviations from sphericity remain small (as shown in Fig. 1top). As a result, g modes ``detect'' a much weaker distortion that is then amenable to a perturbative description. A particular feature that induces a strong deviation from the perturbative method concerns mixed pressuregravity modes that arise as a consequence of the centrifugal modification of the stellar structure. For example, we found that, above a certain rotation rate, the mode becomes a mixed mode with a pmode character in the outer lowlatitude region associated with a drop in the BruntVäisälä frequency N_{o} (see Fig. 1).
The domains of validity of perturbative methods are strongly reduced in the lowfrequency range. For Dor stars, 2ndorder perturbative methods are only valid below 50 km s^{1}. Indeed, a striking feature of Fig. 5 is that perturbative methods fail to recover the correct frequencies in the inertial regime (delimited by a magenta curve). In particular, we observe that, although increasing the tolerance between the left ( Dor) and right (B star) panels subtantially extends the domains of validity in the regime, very little improvement is observed in the regime.
In the following, we argue that the failure of the perturbative method in the subinertial regime is related to changes in the mode cavity that are not taken into account by the perturbative method. Indeed, we observed that modes in the inertial regime do not explore the polar region and that the angular size of this forbidden region increases with . This is illustrated in Fig. 6 for a particular mode. Such a drastic change in the shape of the resonant cavity has a direct impact on the associated mode frequency. As perturbative methods totally ignore this effect, they cannot provide an accurate approximation of the frequencies in this regime.
This interpretation is supported by the analytical expression of the forbidden region determined by Dintrans & Rieutord (2000) for gravitoinertial modes. Indeed for frequencies
,
modes are mixed gravityinertial modes, since the Coriolis force
becomes a restoring force. In the context of their spherical model, and
within the anelastic approximation and the Cowling approximation, they
have shown that gravitoinertial waves with a frequency
only propagate in the region where
This implies that, when , a critical latitude appears above which waves cannot propagate. Even though this expression does not strictly apply to our nonspherical geometry, we have overplotted the critical surfaces with the energy distributions of our eigenmodes (see Fig. 6 for an illustration). For the nonrotating case, there is only a small circle close to the center, corresponding to the classical turning point . For the mode with , the polar forbidden region delineated by agrees pretty well with the energy distribution of our complete computation.
5 Conclusion
In the present work, we have computed accurate frequencies for g modes in polytropic models of uniformly spinning stars. We started from high and lowfrequency, lowdegree ( ) g modes of a nonrotating star and followed them up to . This allowed us to provide a table of numericallycomputed perturbative coefficients up to the 3rd order for a polytropic stellar structure (with index ). This table can serve as a reference for testing the implementation of perturbative methods. We then determined the domains of validity of perturbative approximations. For the highfrequency (loworder) modes, 2ndorder perturbative methods correctly describe modes up to for Dor stars and up to for B stars. The domains of validity can be extended by a few tens of km s^{1} with 3rdorder terms. However, the domains of validity shrink at low frequency. In particular, perturbative methods fail in the inertial domain because of a modification in the shape of the resonant cavity.
Figure 6: (Left) Meridional distribution of kinetic energy of the g mode in a nonrotating star. To enhance the contrast, it is scaled by r^{2}. (Right) The same for . Magenta lines indicate the critical surface . 

Open with DEXTER 
In a next step, we plan to compare our complete computations with the socalled traditional approximation, which is also extensively used to determine gmode frequencies (e.g Lee & Saio 1997; Berthomieu et al. 1978). We will also analyze how rotation affects the regularities of the spectrum  such as the period spacing  and compare it to the predictions of the perturbative and traditional methods. In the present study, we have focused on lowdegree modes, but a more complete exploration clearly needs to be performed. In particular, we might look for the singular modes predicted by Dintrans & Rieutord (2000). It requires to take care of dissipative processes, which play an important role in this case.
AcknowledgementsThe authors acknowledges support through the ANR project Siroco. Many of the numerical calculations were carried out on the supercomputing facilities of CALMIP (``CALcul en MIdiPyrénées''), which is gratefully acknowledged. The authors also warmly thank Boris Dintrans for discussions and useful comments on this work. D.R.R. gratefully acknowledges support from the CNES (``Centre National d'Études Spatiales'') through a postdoctoral fellowship.
References
 Baglin, A., Michel, E., Auvergne, M., & The COROT Team 2006, in Proc. of SOHO 18/GONG 2006/HELAS I, Beyond the spherical Sun, ESA SP, 624 [Google Scholar]
 Bedding, T. R., Huber, D., Stello, D., et al. 2010, ApJ, 713, 176 [Google Scholar]
 Benomar, O., Baudin, F., Campante, T. L., et al. 2009, A&A, 507, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Berthomieu, G., Gonczi, G., Graff, P., Provost, J., & Rocca, A. 1978, A&A, 70, 597 [Google Scholar]
 Bonazzola, S., Gourgoulhon, E., & Marck, J. 1998, Phys. Rev. D, 58, 104020 [NASA ADS] [CrossRef] [Google Scholar]
 Borucki, W. J., Koch, D. G., Lissauer, J., et al. 2007, in Transiting Extrapolar Planets Workshop, ed. C. Afonso, D. Weldrake, & T. Henning, ASP Conf. Ser., 366, 309 [Google Scholar]
 Chaplin, W. J., Appourchaux, T., Elsworth, Y., et al. 2010, ApJ, 713, 169 [Google Scholar]
 ChristensenDalsgaard, J., & Mullan, D. J. 1994, MNRAS, 270, 921 [NASA ADS] [Google Scholar]
 De Cat, P., Eyer, L., Cuypers, J., et al. 2006, A&A, 449, 281 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dintrans, B., & Rieutord, M. 2000, A&A, 354, 86 [NASA ADS] [Google Scholar]
 Dziembowski, W. A., & Goode, P. R. 1992, ApJ, 394, 670 [NASA ADS] [CrossRef] [Google Scholar]
 Frémat, Y., Neiner, C., Hubert, A., et al. 2006, A&A, 451, 1053 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 García Hernández, A., Moya, A., Michel, E., et al. 2009, A&A, 506, 79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hekker, S., Kallinger, T., Baudin, F., et al. 2009, A&A, 506, 465 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ledoux, P. 1951, ApJ, 114, 373 [NASA ADS] [CrossRef] [Google Scholar]
 Lee, U., & Saio, H. 1997, ApJ, 491, 839 [NASA ADS] [CrossRef] [Google Scholar]
 Lignières, F., & Georgeot, B. 2008, Phys. Rev. E, 78, 016215 [NASA ADS] [CrossRef] [Google Scholar]
 Lignières, F., & Georgeot, B. 2009, A&A, 500, 1173 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lignières, F., Rieutord, M., & Reese, D. 2006, A&A, 455, 607 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lovekin, C. C., & Deupree, R. G. 2008, ApJ, 679, 1499 [NASA ADS] [CrossRef] [Google Scholar]
 Mathias, P., Chapellier, E., Bouabid, M., et al. 2009, in Am. Inst. Phys. Conf. Ser., ed. J. A. Guzik, & P. A. Bradley, 1170, 486 [Google Scholar]
 Michel, E., Baglin, A., Auvergne, M., et al. 2008, Science, 322, 558 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Miglio, A., Montalbán, J., Baudin, F., et al. 2009, A&A, 503, L21 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
 Poretti, E., Michel, E., Garrido, R., et al. 2009, A&A, 506, 85 [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]
 Reese, D., Lignières, F., & Rieutord, M. 2008, A&A, 481, 449 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reese, D. R., MacGregor, K. B., Jackson, S., Skumanich, A., & Metcalfe, T. S. 2009, A&A, 509, 183 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rieutord, M., Corbard, T., Pichon, B., Dintrans, B., & Lignières, F. 2005, in SF2A2005: Semaine de l'Astrophysique Française, ed. F. Casoli, T. Contini, J. M. Hameury, & L. Pagani, 759 [Google Scholar]
 Saio, H. 1981, ApJ, 244, 299 [NASA ADS] [CrossRef] [Google Scholar]
 Soufi, F., Goupil, M. J., & Dziembowski, W. A. 1998, A&A, 334, 911 [NASA ADS] [Google Scholar]
 Suárez, J. C., Goupil, M. J., & Morel, P. 2006, A&A, 449, 673 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tassoul, M. 1980, ApJS, 43, 469 [NASA ADS] [CrossRef] [Google Scholar]
 Valdettaro, L., Rieutord, M., Braconnier, T., & Fraysse, V. 2007, J. Comput. Appl. Math., 205, 382 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Table 1: Perturbative coefficients (see development Eq. (17)) for g modes with frequency , radial order , and in a polytropic stellar model with an index .
All Figures
Figure 1: (Top) Map in the meridional plane of the BruntVäisälä frequency N_{o} for the model with . N_{o} has been normalized by where is the polar radius. Dashed lines show the shape of the spheroidal grid we used. (Bottom) Solid black lines show the profile of N_{o} normalized by for the nonrotating star. Red lines show N_{o} for the model along the polar (solid line) and equatorial radius (dashes). A thick red tick on the xaxis indicates the polar radius . 

Open with DEXTER  
In the text 
Figure 2: Spectrum as a function of the degree of the radial velocity of the mode at . is the absolute value of the highest Chebychev coefficient in the decomposition on the spherical harmonics of degree of the radial component of . The spectrum is normalized to its maximum value. The different lines correspond to different spatial resolutions. 

Open with DEXTER  
In the text 
Figure 3: Evolution with the rotation rate of the frequency of the modes and , with a zoom on a region where an avoided crossing between these two modes occurs. At a given rotation rate, a diamond (respectively, a triangle) indicates the mode is dominated by the component (). Crosses indicate that the modes are hardly discernible: they are both dominated by and their structures are very similar. 

Open with DEXTER  
In the text 
Figure 4: Evolution with the rotation rate of the frequencies of the components of the (, n=16) multiplet obtained with a complete computation (solid line) and with the 2ndorder perturbative approximation (dashed lines). 

Open with DEXTER  
In the text 
Figure 5: Evolution of the frequencies of modes (top to bottom). Frequencies are computed in the corotating frame. Perturbative approximations have been tested for a typical Dor (left panels) and for a B star (right panels). Green/red/blue parts of curves indicate that 1st/2nd/3rd order is sufficient to reproduce complete calculations within an error . Error bars on the righthand side of each panel show and 10 . Magenta lines indicate . For each plot, the bottom xaxis and left yaxis show dimensionless units, whereas the top xaxis and right yaxis show physical units. 

Open with DEXTER  
In the text 
Figure 6: (Left) Meridional distribution of kinetic energy of the g mode in a nonrotating star. To enhance the contrast, it is scaled by r^{2}. (Right) The same for . Magenta lines indicate the critical surface . 

Open with DEXTER  
In the text 
Copyright ESO 2010
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.