Issue 
A&A
Volume 546, October 2012



Article Number  A11  
Number of page(s)  13  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201219716  
Published online  26 September 2012 
Regular oscillation subspectrum of rapidly rotating stars
^{1} CNRS, IRAP, 14 avenue Edouard Belin, 31400 Toulouse, France
email: pasek@irsamc.upstlse.fr
^{2} Université de Toulouse, UPSOMP, IRAP, Toulouse, France
^{3} CNRS, LPT (IRSAMC), 31062 Toulouse, France
^{4} Université de Toulouse, UPS, Laboratoire de Physique Théorique (IRSAMC), 31062 Toulouse, France
^{5} Institut d’Astrophysique et Géophysique de l’Université de Liège, Allée du 6 Août 17, 4000 Liège, Belgium
Received: 30 May 2012
Accepted: 17 August 2012
Aims. We present an asymptotic theory that describes regular frequency spacings of pressure modes in rapidly rotating stars.
Methods. We use an asymptotic method based on an approximate solution of the pressure wave equation constructed from a stable periodic solution of the ray limit. The approximate solution has a Gaussian envelope around the stable ray, and its quantization yields the frequency spectrum.
Results. We construct semianalytical formulas for regular frequency spacings and mode spatial distributions of a subclass of pressure modes in rapidly rotating stars. The results of these formulas are in good agreement with numerical data for oscillations in polytropic stellar models. The regular frequency spacings depend explicitly on internal properties of the star, and their computation for different rotation rates gives new insights on the evolution of mode frequencies with rotation.
Key words: asteroseismology / chaos / methods: analytical / stars: oscillations / stars: rotation / waves
© ESO, 2012
1. Introduction
The field of asteroseismology has now reached its age of maturity with the exploitation of space missions CoRoT (Baglin et al. 2006) and Kepler (Koch et al. 2010) that are gathering stellar light curves with high accuracy. However, there are still unresolved issues that hinder the successful pairing of light curve frequencies with pulsation modes, which is crucial to obtain detailed information on the inner structure of observed stars. One of these issues is the rapid rotation of a star around its axis, since the exact nature of rotational effects on pulsation modes is not known. In particular, the centrifugal flattening (e.g. Monnier et al. 2007) affects the spectrum of pressure modes (pmodes) in a complex way (Lignières & Georgeot 2009). This difficulty mainly concerns nonevolved massive and intermediatemass pulsating stars which are typically rapid rotators (Royer 2009). Recently though, hints of regular frequency spacings have been found in the spectrum of rapidly rotating δ Scuti stars observed with CoRoT (García Hernández et al. 2009; Mantegazza et al. 2012), and this could ease future mode identification.
The recent development of accurate numerical models has enabled progress in the comprehension of pulsation modes in rapidly rotating stars. It has been found in particular (Lignières et al. 2006; Reese et al. 2008, 2009) that in the rapidly rotating regime a subset of pmodes shows approximate regular frequency spacings in the form: (1)where frequencies ω_{n,ℓ,m} are given in the corotating frame. Quantum numbers n, ℓ and m correspond to node numbers of the mode amplitude distributions, Δ_{n}, Δ_{ℓ} and Δ_{m} are frequency regularities, and α is a constant term. The approximate formula in Eq. (1) shows a better agreement with numerical results towards highfrequencies, thus suggesting that this relation is of an asymptotic nature. It should also be noted that, from computations of diskaveraging factors, the pmodes following Eq. (1) are expected to be among the most visible ones (Lignières & Georgeot 2009). An example of such a mode can be seen in Fig. 1.
Fig. 1 Pressure amplitude on a meridian plane for a polytropic stellar model, with d the distance to the rotation axis and ρ_{0} the equilibrium density. The mode shown corresponds to n = 50, ℓ = 1 and m = 1 at a rotation rate of Ω/Ω_{K} = 0.300, where is the limiting rotation rate for which the centrifugal acceleration equals the gravity at the equator, M being the stellar mass and R_{eq} the equatorial radius. Colors/grayness (for respectively the colour/black and white version of the figure) denote pressure amplitude, from red/gray (maximum positive value) to blue/black (minimum negative value) through white (null value). The thick black line is the ray γ located in the center of the main stable island. (This figure is available in color in the electronic form.) 
The frequency spacings of Eq. (1) are notably similar to the regularities described by Tassoul’s asymptotic formula (Tassoul 1980) for low degree pmodes in nonrotating stars. Tassoul’s formula at leading order is (2)with the large frequency separation (3)where c(r) is the radially inhomogeneous sound speed, and R the stellar radius. The integer n_{s} is the node number of the radial component of the mode, while ℓ_{s} is the degree of the associated spherical harmonics, and α_{s} depends on surface properties. Tassoul’s theory has proved to be very useful for interpretating solarlike oscillations in slowly rotating stars. Indeed, the formula relates observable quantities, such as the regular frequency spacing Δ, to physical properties of stellar interiors. For rapidly rotating stars, it would be clearly desirable to gain insights on the underlying physics of the potentially observable regular spacings Δ_{n}, Δ_{ℓ} and Δ_{m} by a similar asymptotic analysis. In this paper we derive a formula for these regular frequency spacings in the asymptotic regime.
The generalization of the pmode asymptotic theory to rapidly rotating stars is not trivial. Tassoul’s theory requires separation of variables, which is no longer possible when the star is flattened by rotation. For nonseparable wave systems, a wellknown technique to obtain eigenmodes is to study the shortwavelength limit of the propagating waves. This limit gives an equation for the propagation of rays that is similar to the geometrical optics limit of electromagnetism, or the classical limit in quantum mechanics. Then, by imposing quantization conditions on the phase of waves propagating on these rays, one obtains the eigenmodes of the wave system. This technique was first developed in the context of quantum physics, and is often called semiclassical quantization.
For spherical stars, the ray limit of pressure waves has been previously used to recover the Tassoul asymptotic formula from the EinsteinBrillouinKeller (EBK) quantization of ray dynamics (Gough 1993). This analytical approach is possible only when the ray system is integrable. A dynamical system is said to be integrable when it has as many conserved quantities (energy, angular momentum, etc.) as degrees of freedom (Ott 2002). In rapidly rotating stars, there are not enough conserved quantities to ensure integrability of the ray dynamics. Indeed, in Lignières & Georgeot (2008, 2009), it has been found that acoustic rays in rotating stars have a very different dynamical behavior depending on their initial conditions in positionmomentum space (the socalled phase space). For a polytropic stellar model, the numerical integration of the equations for acoustic rays displayed various types of solutions. Indeed, one can obtain either stable rays staying on torusshaped surfaces in phase space which form structures such as stable islands, or chaotic rays that are dense and ergodic on a phase space volume.
A similar behavior has been found in many systems studied in the field of theoretical physics known as quantum chaos or wave chaos (Gutzwiller 1990). This field has among its objectives to analyze quantum (wave) systems whose classical (shortwavelength) limit is partly or fully chaotic. In this framework, one can predict the existence of some eigenfunctions (mode profiles) and energies (frequencies) of the quantum (wave) system from the different structures that are present in the classical system phase space (Percival 1973; Berry & Robnik 1984). In the stellar pulsation setting, Lignières & Georgeot (2008, 2009) found that the mixed (i.e. regular and chaotic) character of the acoustic ray dynamics in rapidly rotating stars results in a classification of pmodes in two broad families: regular modes either associated with stable islands or whispering gallery zones, and chaotic modes associated with ergodic regions in phase space. For the regular modes associated with stable islands, the socalled island modes, it is known to be possible to obtain approximate analytical solutions by solving the wave equation in the vicinity of a periodic stable ray (Babich & Buldyrev 1991). A simple application of such a method is found in modes of optical resonators, where the periodic stable light ray is a straight line between two reflecting mirrors (Kogelnik & Li 1966). These methods have been previously employed to obtain modes of more complex lasing (Tureci et al. 2002) and electronic (Zalipaev et al. 2008) cavities as well as quantum chaos systems (Vagov et al. 2009). In this paper, we apply this approach to rapidly rotating stars.
In the present analysis, we thus construct an asymptotic formula for regularities in the pmode spectrum of rapidly rotating stars. Part of the results were already presented in the short communication of Pasek et al. (2011). In the present paper we give a detailed derivation of these results, specify their domain of validity, extend them with a study of rotational splittings, and explore their astrophysical applications.
The paper is organized as follows. In Sect. 2 we present the wave equation for pmodes in rotating stars and its asymptotic limit leading to an equation for acoustic rays. In Sect. 3 we use a stable periodic solution of the ray dynamics to obtain a semianalytical formula for the associated pmodes, and to derive a formula for the associated regular frequency spacings. We then compare the results obtained from the derived formulas for mode frequencies and spatial distributions with numerical results (Sect. 4). Finally, we suggest directions on how these results could be used for the asteroseismic diagnosis of rapidly rotating stars by discussing the phenomenological implications of the theory in Sect. 5.
2. Pmodes in rotating stars and their asymptotic limit
In Sect. 2.1 we introduce the wave equation for pmodes in rotating stars. We then present the asymptotic limit of this equation in order to obtain an equation for the dynamics of acoustic rays (Sect. 2.2).
2.1. Pressure modes in rotating stars
We start with the equation for small adiabatic timeharmonic perturbations of the pressure field in a selfgravitating gas. Since we are interested in obtaining an asymptotic theory for pmodes in the highfrequency regime, we use the Cowling approximation (i.e. we neglect the perturbations of the gravitational potential), an approximation known to be valid for highfrequency perturbations in nonrotating stars (Aerts et al. 2010). We also neglect the Coriolis force. Indeed, in the highfrequency regime, the time scale associated with this force is much longer than the mode period, and thus the influence of the Coriolis force on pulsation frequencies is weak. This has been numerically checked in Lignières et al. (2006), Reese et al. (2006, 2008). In the asymptotic regime of pmodes, the oscillation frequencies are far greater than the BruntVäisälä frequency and thus we can discard the terms corresponding to gravity waves. With these assumptions, the equation for pressure perturbations is a Helmholtz equation such that (4)where is the complex amplitude associated with the pressure perturbation , f is a function of the background model, ω_{c} is the cutoff frequency of the model and c_{s} its inhomogeneous sound velocity (for a detailed derivation of this equation see Lignières & Georgeot 2009). The stellar model is not spherically symmetric due to the centrifugal distortion, but is however cylindrically symmetric with respect to the rotation axis. Therefore, we can write the pressure field as Ψ = Ψ_{m}exp(imφ) where m is an integer and φ is the azimuth angle of spherical coordinates. By inserting this expression in Eq. (4) we obtain (cf. Appendix A) (5)where d is the distance to the rotation axis. The new mode amplitude Φ_{m} is such that . We introduce a renormalized sound velocity: (6)We notice that besides its spatial dependence, also depends on ω and m, and that m is taken as a parameter for the twodimensional wave equation (Eq. (5)).
2.2. Ray limit of pmodes
In nonrotating stars, the asymptotic theory of high frequency pmodes has first been derived by Vandakurov (1967), and Tassoul (1980). The method was to use the spherical symmetry of the star to reduce the problem to a onedimensional equation in order to obtain the mode frequencies. This method is not applicable when the centrifugal force breaks the spherical symmetry of the star. In this case though, one can study the shortwavelength limit (ω → ∞) of the wave equation (Eq. (4)) (as detailed in Lignières & Georgeot 2009). This provides a Hamiltonian system describing the propagation of acoustic rays. The Hamiltonian has been derived in Lignières & Georgeot (2009) as (7)where the frequencyscaled wavevector is the projection of onto the meridional plane of the star. We notice that this expression has been derived from the shortwavelength limit of the three dimensional wave equation (Eq. (4)) and then, projected onto the corotating meridian plane. An alternative derivation would be to start from the twodimensional wave equation (Eq. (5)). In this case, the ray limit yields the same expression with the addition of the 1/4 factor of Eq. (5) that accounts for the impossibility of acoustic rays to go through the rotation axis (i.e. d = 0). Throughout the paper we use Eq. (7) as the Hamiltonian for acoustic rays.
To probe the integrability property of a dynamical system, it is convenient to use the Poincaré surface of section (PSS), a standard tool in dynamical systems theory (Gutzwiller 1990; Ott 2002) to visualize the structures in phase space. A PSS is a lower dimensional slice of phase space. The acoustic ray dynamical system in the meridional plane has two degrees of freedom, which gives a phase space of dimension four (two for positions, and two for momenta). There is one conserved quantity in the form of the acoustic wave frequency, so the dynamics belongs to a threedimensional manifold in phase space. By fixing an additional position or momentum coordinate, we obtain a twodimensional PSS which can be easily visualized. An example of such a section for our system is shown in Fig. 2. Different choices of PSS variables are possible, some of which are presented in Lignières & Georgeot (2009). We have here chosen to fix the colatitude θ = π/2, so that the PSS corresponds to the crossing of rays with the equatorial halfplane. We thus display a section in coordinates (r/R_{eq},k_{r}/ω) where k_{r} is the norm of the radial wavevector, ω the mode frequency, and R_{eq} the equatorial radius (that may be greater than the polar radius since the star is flattened by rotation).
Fig. 2 Poincaré Surface of Section (PSS) at the rotation rate Ω/Ω_{K} = 0.589 for quantum number m = 0. Each dot corresponds to the crossing of an acoustic ray with the equatorial halfplane in the (r/R_{eq},k_{r}/ω) phase space, r being the radial coordinate and k_{r} the associated momentum. R_{eq} is the equatorial radius and ω the mode frequency. Red/dark gray denotes a chaotic ray, green/light gray a whispering gallery ray, blue/black a stable island ray (see text). Upper inset is a closeup of the main stable island. (This figure is available in color in the electronic form.) 
In such plots, each dot corresponds to the crossing of an acoustic ray with the PSS. Successive dots from a single ray will form lines in integrable zones, or fill surfaces densely in chaotic zones. We see in Fig. 2 that when the rotation rate Ω/Ω_{K} (where is the limiting rotation rate) is large, different structures coexist in the system phase space: stable islands correspond to concentric circles around a stable ray, whispering gallery rays to lines near the surface, and chaotic zones to densely filled areas (for more details, see Lignières & Georgeot 2009). In this paper, we will focus on the 2periodic stable island which is the main stable island (shown in the inset of Fig. 2). The rays’ dynamics is very sensitive to the rotation rate, so the PSS will be different for each rotational velocity. Indeed, the locus in phase space of the main stable island changes as rotation increases. The major change happens when the central ray of the 2periodic stable island undergoes a bifurcation at Ω/Ω_{K} ≃ 0.26. For slow rotation rates, the central ray of the island is located on the polar axis, and through this bifurcation it transforms into two stable rays surrounding one unstable ray on the polar axis. Then, as rotation increases, the stable island will coast away from the polar axis. This bifurcation will be of some importance in the following, when we will show that one can construct approximate eigenmodes of the wave system from this 2periodic stable island. In Fig. 1, one can see an example of an island mode obtained from a fullnumerical computation, together with the central ray of the 2periodic stable island for the same rotation rate and quantum number m.
3. Semianalytical method for island modes
In this section we construct an asymptotic approximation of a subset of regular pmodes associated with a stable periodic ray. The method is based on the works of Babich and coworkers (see Babich & Buldyrev 1991 and references therein for the general formalism; and e.g. Zalipaev et al. 2008; Vagov et al. 2009, for some applications). It consists in deriving an approximation of the wave equation in the vicinity of the ray (Sect. 3.1), finding Gaussian wavepacket solutions related to the stability properties of the ray (Sect. 3.2), and then deriving the asymptotic formula for the frequencies from a quantization condition (Sect. 3.3).
3.1. Approximate wave equation in the vicinity of a stable ray
For a given rotation rate and quantum number m, we start with the central periodic ray of the main stable island (see Sect. 2). This ray must be computed by numerically evaluating the Hamiltonian equations derived from Eq. (7). In the following we will call this ray γ. The first step is to write the wave equation (Eq. (5)) in the vicinity of γ. For this, we use a local orthonormal coordinate system (s,ξ) defined as r = sT + ξN where s is the arc length along the ray, T the unit tangent vector, ξ the transverse coordinate and N the unit vector normal to T. The two basis vectors are related by the curvature κ(s) of the ray as follows: (8)In this coordinate system, the wave equation (Eq. (5)) reads (9)where the scale factors h_{s} and h_{ξ} (Arfken & Weber 2005) are (10)and (11)This yields (12)In the vicinity of the ray γ, that is for small ξ, the terms of Eq. (9) are simplified using: and (16)We then express the function Φ_{m}(s,ξ) in terms of a WentzelKramersBrillouin (WKB) ansatz as (17)where τ is an unknown function of position. The fundamental assumption underlying the theory of Babich is that, as ω → + ∞, the mode is localized on the acoustic ray and that its transverse extent scales as 1/. Such a solution can be found by assuming that the transverse variable ξ scales as (18)Then, from an expansion of Eq. (9) in ω, one obtains at the dominant order that the WKB phase in Eq. (17) depends only on s as . At the next order in ω one finds a parabolic equation for the function V_{m}: (19)with (20)where we introduced the scaled coordinate and .
3.2. Solutions of the parabolic wave equation
To find a solution to Eq. (19), we first find a solution at a fixed arc length s, and then study how this solution must evolve with s. At fixed s, the first two terms of Eq. (19) correspond to the equation for a quantum harmonic oscillator in the direction N transverse to γ. Thus, as we know from quantum mechanics (CohenTannoudji et al. 1973), a solution of this equation is a Gaussian wavepacket, transverse to the ray, that we write as (21)with Γ an unknown complexvalued function. To find the variation of this Gaussian wavepacket along the ray γ, we introduce a solution of this form in the parabolic equation (Eq. (19)) and obtain a Riccati equation for Γ (22)simple equation for the factor A(23)In the following, we show that the equation for Γ is related to the ray properties in the vicinity of γ and can thus be solved from ray dynamics computations. First, using the variables (z(s),p(s)) defined as and , Eq. (22) is transformed into the Hamiltonian system: where the (timedependent) Hamiltonian function is (26)τ being the time coordinate. From Eqs. ((24), (25)), one can derive an equation for z only to obtain (27)Equations ((24), (25)) have two independent solutions (z,p) and that are conjugate to each other. Of these two solutions, only one is physically relevant, i.e. corresponds to a localized wavepacket. According to Eq. (21) this happens if Im(Γ) > 0 for all s, and as shown in Appendix B, the sign of Im(Γ) stays constant along γ since (28)The variation along γ of the Gaussian wavepacket can now be linked to the dynamics of the acoustic rays nearby γ. Indeed, Im(Γ) has a simple expression in terms of the complex variable z, as shown in Eq. (28). This variable, on the other hand, obeys Eq. (27) which can be shown to be the same as the equation describing the deviation from γ of a ray nearby γ (see derivation in Appendix C). The Hamiltonian in Eq. (26) is thus a local integrable approximation, also known as a normal form approximation (Arnol’d 1989), to the full Hamiltonian for acoustic rays written in Eq. (7).
Now, our task is to find the two linearly independent complex conjugate solutions of Eqs. ((24), (25)). The terms in these equations depend only on quantities that are evaluated on the periodic ray γ. Therefore, these equations are periodic in s, or equivalently in τ. Equations ((24), (25)) can thus be written as: (29)where the matrix Σ (30)verifies Σ(τ + T_{γ}) = Σ(τ), being the acoustic travel time along γ. Then if (z(τ),p(τ)) is a solution of Eq. (29), so is (z(τ + T_{γ}),p(τ + T_{γ})) and the two solutions are related by the following linear map (31)where M is called the monodromy matrix (Cvitanović et al. 2010, and references therein). As (z,p) describe ray deviations from γ, the matrix M characterizes the stability of γ. As γ is stable, we know that  Tr(M)  < 2 and that the eigenvalues are of modulus one and complex conjugates of each other i.e. Λ^{ ± } = exp( ± iα) with α ∈ ] 0, π [ (cf. Appendix D), where α is called a Floquet phase or stability angle. Hence the two linearly independent solutions of Eq. (29) can be written in the form: (32)where the functions u_{ ± }(τ) are periodic with period T_{γ} and v^{ ± } are independent eigenvectors of the monodromy matrix M.
An expression of the monodromy matrix in terms of second derivatives of the action function S can be derived (Bogomolny 2006). The action function S is defined by a trajectory from the position q_{i} to q_{f} for a given energy or frequency ω (Gutzwiller 1990). For our purposes, the action is written as (33)where σ is the arclength along a ray nearby γ. If we write the monodromy matrix as (34)we can express its components from the second derivatives of the action function S with the following formulas (35)where the positions q_{i} and q_{f} are written as (s_{i},z_{i}) and (s_{f},z_{f}), z_{i} and z_{f} being respectively the initial and final transverse positions of the neighboring ray after one period, and the derivatives are evaluated on the periodic ray. From these expressions, and the simple formula giving the roots of a second degree polynomial (cf. Appendix D), we can obtain an expression for the stability angle α as (36)where (37)We thus have obtained a solution of the approximate wave equation Eq. (19) in the form of a Gaussian wavepacket (Eq. (21)), whose evolution along the ray γ is given by Eqs. ((29), (30)). It is possible to find other solutions of Eq. (19) that have a finite number of nodes in the direction transverse to γ. As for the quantum harmonic oscillator (CohenTannoudji et al. 1973), these solutions can be obtained from Eq. (21) using the annihilation â and creation operator that are defined as (Babich & Buldyrev 1991) (38)These operators have the commutation rule , where the commutator is defined as . Indeed, from the commutation relation of with the operator of Eq. (19), the functions (39)are also solutions of Eq. (19) (Babich & Buldyrev 1991). The calculation of these higherℓ solutions then leads to (40)with H_{ℓ} the Hermite polynomials of order ℓ. Finally, we can write the solutions of Eq. (5) as (41)
3.3. Quantization condition and regular frequency spacings
The quantization condition is based on the singlevaluedness of the solution presented in Eq. (41), and thus asserts that the phase accumulated by the function over one period must be a multiple of 2π. In the following, we assume without loss of generality that the eigenvalue which corresponds to the wavepacket localization is exp( + iα). The contribution of function to the dynamical phase of is obtained from Eqs. (32) and (40). We obtain that the phase accumulated over one period is (42)where α is the Floquet phase that is defined modulo 2π. For our purposes, we must also take into account the number N_{r} of multiples of 2π acquired by the phase. N_{r} can be computed by following the evolution of the eigenvector v over one period, and we verified numerically that, alternatively, N_{r} is also the winding number around γ of a ray nearby γ during one period. The last term in Eq. (42) is the Maslov phase (Gutzwiller 1990) that comes from the reflection of the wave on the boundaries. The formula for the frequencies, ω_{n,l,m}, of island modes is thus (43)or, in a form that makes the regularities more visible, (44)with the frequency regularities (45)and the constant term (46)The quantum numbers n and ℓ correspond to node numbers of the pmode, respectively in the longitudinal and transverse directions of the central ray γ as illustrated in Fig. 3. This is to be contrasted with the case of a spherical mode where the most natural labeling are the quantum numbers of spherical harmonics. Equation (43) is a semianalytical formula since the quantities T_{γ}, N_{r} and α must be computed numerically from the RungeKutta integration of the Hamiltonian equations for acoustic rays. The acoustic time, T_{γ}, is directly computed from γ itself. From the intersections of a ray nearby γ with the PSS (as can be seen in Fig. 2), we compute the monodromy matrix M that maps one intersection with the PSS to the next one. Then, by diagonalizing this matrix one obtains the Floquet phase α from its eigenvalues, and the functions z and Γ from its eigenvectors. It is important to note that for m even, only modes symmetric with respect to the rotation axis exist. Since the preceding theory does not take this phenomenon into account, the theoretical value of δ_{ℓ} is multiplied by two when the ray γ coincides with the rotation axis, i.e. for rotation rates less than the bifurcation point. Finally, it can be noted that a formula similar to Eq. (43) can also be obtained through the formalism of the Gutzwiller trace formula following the method of Miller (1975).
Fig. 3 Left: schematic representation of the mode labeling/quantum numbers used for island modes. Right: illustration of the relation between the quantum numbers of spherical and island modes. The quantum number n is the number of nodes in the longitudinal direction of the island mode (along the ray γ), ℓ is the number of nodes in the transverse direction of the mode (transverse to γ) and m is the number of azimuthal nodes. n_{s}, ℓ_{s} and m_{s} are the quantum numbers of spherical modes. In the right figure, the multiplets of island modes correspond to diagonal colored bands, whereas the multiplets of spherical modes would have the form of vertical bands. The plotted mode on the left corresponds to n = 46, ℓ = 1 and m = 0. (This figure is available in color in the electronic form.) 
4. Comparison with numerical results
As the present asymptotic theory relies on various assumptions, its relevance for stellar seismology is not guaranteed and needs to be assessed through a comparison with exact calculations of realistic stellar models. In this section, the comparison is done with numerically computed modes in uniformly rotating polytropic models of stars. The hypotheses of the asymptotic theory are the following: first, it is valid in the asymptotic regime, that is, for high enough frequencies. Second, the island modes are constructed from a stable island of acoustic ray phase space. At null rotation such a structure does not exist, so we expect that the theory fails to describe spherical mode amplitudes and frequencies. A stable island immediately appears at nonzero rotation, but its phase space volume must be high enough for an island mode to exist. This volume increases with frequency and rotation (see Lignières & Georgeot 2008, 2009, for details). Thus, for a given frequency range, the number of island modes starts from zero at small rotation rates and progressively increases as the rotation and thus the phase space volume of the stable island grows. Actually, low degree spherical modes become progressively island modes as rotation increases. Another assumption used in finding a solution to the wave equation is that the mode decays as ∝ 1/ in the direction transverse to the periodic ray. Finally, the theory also neglects the Coriolis force and the perturbations of the gravitational potential.
In the following, the asymptotic theory is compared with highly accurate computations of high frequency adiabatic modes in uniformly rotating polytropic stellar models with index N = 3, the Coriolis effect and perturbations of the gravitational potential being taken into account. The accuracy of these calculations, described in detail in Reese et al. (2006), is very high (the relative precision on the frequencies is 10^{7}) and thus does not interfere with the present comparison. A large number of modes were followed from Ω/Ω_{K} = 0 to Ω/Ω_{K} = 0.896. At zero rotation these modes are low degree ℓ_{s} ∈ { 0,1,2,3 } , high order n_{s} ∈ [21,25] modes. At higher rotation rates, they become island modes and can thus be labeled with n and ℓ, the number of nodes along and transverse to γ, respectively, as illustrated in Fig. 3. The relation between the quantum numbers at zero and high rotation rates is (Reese 2008): We remind the reader that rotational multiplets are defined, in the nonrotating case, as a set of frequencies with identical n_{s}, ℓ_{s} quantum numbers but different values of m_{s} for m_{s} ∈ [ − ℓ_{s},ℓ_{s}] . For rotating stars, we can define multiplets as frequencies with identical n and ℓ but different m ∈ Z, i.e. without any limiting value for m. The relation between the two sets of quantum numbers and the two types of multiplets is illustrated in Fig. 3. In this figure the multiplets of island modes correspond to diagonal colored bands, whereas the multiplets at zero rotation would have the form of vertical bands. We restricted ourselves to numerical modes with so, in terms of island mode quantum numbers, the range of numerically computed modes is the one given in Table 1.
Island mode quantum numbers n, ℓ, m of numerically computed modes, corresponding to n_{s} ∈ [21,25] , ℓ_{s} ∈ {0,1,2,3}, m_{s} ∈ [−ℓ_{s},ℓ_{s}] in terms of spherical mode quantum numbers.
The associated numerical frequency spacings are defined as: (50)and (51)The semianalytical asymptotic theory also requires determining the α term in Eq. (45) numerically. To test the robustness of this calculation, we checked that the frequency spacing δ_{ℓ} only weakly depends on the choice of the ray nearby γ that is used to compute α. Also, the spacings δ_{n} and δ_{ℓ} neither depend on the resolution of the background model nor on the integration parameters of the RungeKutta method.
4.1. Regular frequency spacings
According to Eq. (44), the structure of the spectrum is characterized by the two spacings δ_{n}(m) and δ_{ℓ}(m). In Fig. 4, their semianalytical and numerical values computed for m = 0 and  m  = 1 are compared as a function of the rotation rate.
Fig. 4 Comparison of frequency spacings δ_{n,ℓ} between island modes, computed from numerical simulations and semianalytical formulas for different values of Ω/Ω_{K} (the frequency spacings are normalized by with R_{p} the polar radius). Circles: δ_{n}, triangles: δ_{ℓ}, orange/gray: semianalytical results, blue/black: numerical results. Numerical results correspond to different sets of and values, with n ∈ [42,51] , ℓ ∈ { 0,1 } for m = 0; and n ∈ [42,51] , ℓ = 0, n ∈ { 42,44,46,48,50 } , ℓ = 1 for m ∈ { − 1,1 } . Upper panel: m = 0. Lower panel: m = ± 1. (This figure is available in color in the electronic form.) 
One can see that the semianalytical regularities δ_{n}, δ_{ℓ}, and the full computations of highfrequency pmodes are in good agreement for almost all rotation rates. For m = 0, around Ω/Ω_{K} ≃ 0.26, the agreement degrades significantly. In this rotation range, the ray γ in the center of the main stable island undergoes a bifurcation from one stable ray on the polar axis to two stable rays surrounding one unstable ray. When such a bifurcation occurs, the eigenvalues of the monodromy matrix become Λ^{ ± } = 1, corresponding to a Floquet phase α = 0 mod 2π (see e.g. Brack 2001). This is indeed what happens at Ω/Ω_{K} ≃ 0.26, as δ_{ℓ} ∝ α goes to zero. Such a behavior conveys the nonvalidity of the present normal form approximation for rays undergoing a bifurcation. One possibility would be to use other local approximations of the ray dynamics called uniform approximations (Schomerus & Sieber 1997). The discrepancy coming from the bifurcation is not to be found for m ≠ 0, since the stable ray stays away from the polar axis and does not undergo a bifurcation as rotation increases.
Although, as mentioned before, the theoretical and numerical frequency spacings are not expected to match for slow rotation rates, the discrepancies remain small in this rotation range. This is due to the fact that, as Ω/Ω_{K} approaches zero, the stable ray is along the polar axis and, according to the expression of δ_{n}, this implies that δ_{n} = Δ/2, i.e. half the large separation defined in Eq. (3). Now, using the first order of Tassoul’s formula and the quantum numbers conversion rules Eqs. ((47)–(49)), it is easy to see that, at zero rotation, is expected to be close to half the large separation. Note also that the doubling of the numerical values observed in Fig. 4 at small rotation rates is due to the small separation that appears at the next order of Tassoul’s theory. Concerning δ_{ℓ}, the semianalytical calculations indicates that δ_{ℓ} goes to 2δ_{n}, that is Δ, for slow rotation rates. Again, Tassoul’s formula applied to the n, ℓ quantum numbers shows that is close to Δ. For these reasons, the frequency spacings δ_{n} and δ_{ℓ} converge to the results of the first order of Tassoul’s formula, though their derivation is formally not possible for nonrotating stars.
In order to investigate the drift between the spectra of different m as rotation increases, we consider the frequency spacing: (52)Figure 5 displays a comparison between the numerical and semianalytical values of δ_{m} for ℓ = 0 and  m  = 1. As expected, the agreement is not good at small rotation rates. Using the first order of Tassoul’s formula to approximate the numerical results at zero rotation, is found to be close to  m  Δ/2 when Ω/Ω_{K} = 0. This is not compatible with the asymptotic theory of the island mode since it predicts that δ_{m} goes to zero when ω goes to infinity. Indeed, δ_{n}(m) depends on m/ω because and the ray path γ, given by the Hamiltonian Eq. (7), both depend on m/ω. An alternative explanation is to consider the spatial distribution of island modes of fixed m and ℓ: one finds that increasing ω produces both larger derivatives along the stable ray associated with a higher node number n, and larger transverse derivatives because the transverse extent scales as 1/. Thus, the contribution of the azimuthal derivatives becomes negligible in the wave equation (Eq. (4)). We also verified that δ_{m} displayed in Fig. 5 diminishes when n is increased. Thus, for rotation rates such that the numerical modes are not fully island modes, they behave more like spherical modes and shows clear discrepancies with the asymptotic results.
Fig. 5 Comparison of frequency spacing δ_{m} between island modes, computed from numerical simulations and semianalytical formula Eq. (53) for different values of Ω/Ω_{K} (the frequency spacings are normalized by with R_{p} the polar radius). Orange/gray: semianalytical results of Eq. (53) for n = 42 and m = 1, blue/black: numerical results. The numerical results correspond to different sets of δ_{m} for n ∈ [42,51] , ℓ = 0 and m ∈ { − 1,1 } . (This figure is available in color in the electronic form.) 
By contrast, at high rotation rates, an approximate analytical formula for δ_{m} is derived in the following and shown to closely reproduce the numerical results. Starting from (53)we first assume that n is large enough to neglect β(m) − β(0). From Eq. (45), δ_{n}(m) is equal to 2π/T_{γ}(m) where (54)The dependence of T_{γ} on m is explicit in the integrand but is implicit in the integration path γ. In the following, the variation of the location of γ with m/ω is assumed to be negligible. Then, an expansion in 1/ω of the integrand in Eq. (54) leads to: (55)Hence we obtain an approximate expression for δ_{n}(m) of the form: (56)If we insert the previous expression for δ_{n}(m) in Eq. (53) and neglect β(m) − β(0), we have (57)Finally, normalizing by and replacing ω by yields (58)where we have used the fact that stays constant in the frequency range considered here, and is known to be close to . The previous expression will be made more precise by renormalizing the value of c_{s} by to take into account that ω_{c} is not negligible, and indeed is of the order of ω, close to the stellar surface. In Fig. 6, the numerical values of δ_{m}(ℓ = 0) as well as results for the last term in Eq. (58) are plotted as a function of m/ for Ω/Ω_{K} = 0.419, showing a good agreement. This behavior is valid for rotation rates higher than Ω/Ω_{K} ≃ 0.4. It must also be noted that in the numerical calculations by Reese et al. (2009), using more realistic stellar models, the asymptotic dependency was also found empirically.
Fig. 6 Frequency spacings δ_{m}(ℓ = 0) = ω_{n,ℓ = 0,m} − ω_{n,ℓ = 0,m = 0}, normalized by ω_{p}, as a function of m/. The integers n, ℓ and m are the quantum numbers of island modes. The rotation rate is Ω/Ω_{K} = 0.419. Blue/black dots: numerical results. Orange/gray dashed line: semianalytical results for the last term in Eq. (58). Numerical modes included are for quantum numbers n ∈ {42,44,46,48,50}, ℓ = 0 and m ∈ [−3,3] . (This figure is available in color in the electronic form.) 
4.2. Pressure amplitudes of island modes
In this section, we compare the results obtained from the semianalytical formula for mode spatial distributions Eq. (41) with results from full numerical computations. Equatorial cuts of the semianalytical modes can be expressed as a function of , where r_{0} is the radial position of the ray γ, while the value of Γ is obtained from the eigenvectors of the monodromy matrix M. In Fig. 7, the equatorial cuts of semianalytical and numerical modes are plotted for different rotational velocities and quantum numbers ℓ and m. The chosen modes are representative of the different behaviors observed.
Fig. 7 Examples of normalized amplitudes distributions (real part of ) on the equator as a function of position r/R_{eq} (where R_{eq} is the equatorial radius of the stellar model) for semianalytical and numerical modes. Modes displayed are for Ω/Ω_{K} = 0.300 (left column), and Ω/Ω_{K} = 0.707 (right column). Quantum numbers (n,ℓ,m) are as follows. Left column (50,0,0), (50,1,0), (50,0,1), (50,1,1). Right column (50,0,0), (50,1,0), (50,2,0), (50,3,0). Blue/black continuous line: numerical results; orange/gray dashed line: semianalytical results. (This figure is available in color in the electronic form.) 
Discrepancies between semianalytical and numerical results are mainly due to edge effects. This occurs when the transverse extent of the mode (which scales as 1/) reaches either the polar axis for small rotations, or the surface near the equator for high rotations. Finally, avoided crossings can also contravene an accurate prediction for mode amplitudes since the amplitudes of crossing modes will be linear combinations of all the modes contributing to the crossing. Hence, modes undergoing an avoided crossing can differ significantly from Eq. (41) (cf. third panel in Fig. 7). Overall, there is nevertheless a good agreement between the semianalytical and numerical results for mode spatial distributions, showing the validity of Eq. (41).
5. Phenomenology and observables for asteroseismology
In this section, we show that the asymptotic theory provides a simple understanding of the evolution of the island mode spectrum with rotation. Then, the physical content of the potentially observable frequency spacings δ_{n}, δ_{ℓ} and δ_{m} is discussed.
Figure 8 displays the global evolution of all the numerical frequencies considered in the observer’s frame, whose island mode quantum numbers can be found in Table 1 (or equivalently n_{s} ∈ [21,25] , ℓ_{s} ∈ [0,3], m_{s} ∈ [−ℓ_{s},ℓ_{s}] in spherical modes quantum numbers). The first phenomenon that can be noticed is a global decrease of frequencies with rotation. This effect is simply due to the increasing volume of the star when it is spinning rapidly.
Fig. 8 Normalized frequencies ω/ω_{p} as a function of normalized rotation velocity Ω/Ω_{K} in the observer’s frame. Different colors correspond to different values of  m  . Blue: m = 0, green:  m  = 1, red:  m  = 2, cyan:  m  = 3. Modes included are for the quantum numbers n_{s} ∈ [21,23] , ℓ_{s} ∈ [0,3] , m_{s} ∈ [ − ℓ_{s},ℓ_{s}] (see Table 1). (This figure is available in color in the electronic form.) 
Besides this global effect, the evolution of the spectrum’s organization can be inferred from the evolution of frequency spacings δ_{n}, δ_{ℓ} and δ_{m}. The spacing δ_{n} stays almost constant from null up to high rotations, its value remaining close to half the large frequency separation of the spherical model. If a large number of island modes are detected in an observed spectrum, δ_{n} should be easily extracted from the data. By contrast, the rapid evolution of δ_{ℓ} with rotation will strongly modify the spectrum’s organization. This is shown in Fig. 9 where for clarity only a few m = 0 modes have been displayed: the ℓ = 0,n ∈ [43,46] and ℓ = 1,n ∈ [42,44] modes (or equivalently the (n_{s},ℓ_{s}) ∈ { (21,1),(21,2),(21,3),(22,0),(22,1),(22,2),(23,0) } modes).
Fig. 9 Normalized frequencies ω/ω_{p} as a function of normalized rotation velocity Ω/Ω_{K}. The quantum numbers of the displayed modes are: ℓ = 0, n ∈ [43,46] and ℓ = 1, n ∈ [42,44] with m = 0. Blue: ℓ = 0, green: ℓ = 1. For clarity, the corresponding degrees of the spherical harmonics are also written. We have outlined the large frequency separation Δ, the small frequency separation δ = ω_{ns,ℓs} − ω_{ns − 1,ℓs + 2} and the δ_{n}, δ_{ℓ} spacings with arrows. (This figure is available in color in the electronic form.) 
Starting from the usual structure at zero rotation involving the large and small separations of Tassoul’s theory, the spectrum reorganization induced by the decrease of δ_{ℓ} can also be viewed as an increase of the small separation δ. Then, above Ω/Ω_{K} ≃ 0.45, the structure of the m = 0 spectrum remains practically unchanged.
Now, to illustrate the evolution of the spectra of different m, Fig. 10 displays the n = 44, ℓ = 0, m ∈ { − 2, − 1,1,2 } mode frequencies as a function of the rotation rate together with the n ∈ [43,46] , ℓ = 0, m = 0 frequencies. The main feature of this evolution is the decrease of δ_{m} from ≃ at zero rotation to very small values at high rotations. When multiplets of island modes are defined as in Sect. 4 and Fig. 3, they show no regularity at small rotation rates. Note, however, that the splitting ω_{m} − ω_{ − m} is always very close to − 2mΩ because the effects of the Coriolis force are negligible. By contrast, at high rotation rates, as δ_{m} vanishes above Ω/Ω_{K} ≃ 0.45, the m ∈ { − 2, − 1,0,1,2 } modes clearly form a regular multiplet, as can be seen in Fig. 10, where deviations from strict Ω spacings are due to the (m/ term. Since for such rotation rates the structure of the m = 0 spectrum remains unchanged, the evolution of the whole spectrum in the observer’s frame is dominated by the advection term mΩ.
Fig. 10 Normalized frequencies ω/ω_{p} as a function of normalized rotation velocity Ω/Ω_{K}. The quantum numbers of the displayed modes are: n = 44, ℓ = 0, m ∈ { − 2, − 1,1,2 } with the n ∈ [43,46] , ℓ = 0, m = 0 modes. Blue: m = 0, red:  m  ∈ { 1,2 } . For clarity, the corresponding degrees of the spherical harmonics are also written. Continuous line: frequencies in the observer frame, dashed line: in the rotating frame. We have outlined the large frequency separation Δ and the δ_{m} spacing with arrows. (This figure is available in color in the electronic form.) 
The global evolution of mode frequencies in the observer’s frame shown in Fig. 8 also presents some particular events: a first clustering of mode frequencies occurs around Ω/Ω_{K} ≃ 0.25 and then a second one around Ω/Ω_{K} ≃ 0.56. Both phenomena can be understood from the asymptotic theory. According to the asymptotic formulas Eqs. ((43)–(46)), crossings of mode frequencies will happen when δ_{ℓ}/δ_{n}, or equivalently α/π, has a rational value. Though the asymptotic theory predicts true eigenvalue crossings, it is known that these crossings will be avoided if the two modes are of the same symmetry class (Landau & Lifshitz 1977). As can be seen in Fig. 4, δ_{ℓ} becomes equal to δ_{n} at some rotation rate around Ω/Ω_{K} ≃ 0.25 where the spectrum for a given m simplifies to ω_{n,m} = δ_{n}(m)n + β. The degeneracy occurs between modes of a different symmetry class, and the rotation rate at which it occurs depends only weakly on the m values of the modes, if m is small. This property translates itself into a clustering of the full spectrum in the observer’s frame because it turns out that this rotation rate is close to δ_{n}/2; and δ_{m}, that decreases from an initial value of mδ_{n} to zero at high rotation, is around mδ_{n}/2 at this intermediate rotation. The second frequency clustering close to Ω/Ω_{K} ≃ 0.56 is related to the fact that δ_{m} vanishes at high rotation. In this regime, the different m spectra are expected to collapse onto a single spectrum in the rotating frame but not in the observer’s frame. However, when Ω is equal to δ_{n}, the near degeneracy of the m spectra produces the frequency clustering observed at Ω/Ω_{K} ≃ 0.56.
One of the interests of asymptotic theories in asteroseismology is to gain physical insights into seismic observables such as δ_{n}, δ_{l}, and δ_{m}. In the following, we briefly discuss this point with emphasis on the differences and similarities with the physical content of the large and small separation from Tassoul’s theory. The spacing δ_{n} depends only on the acoustic travel time along the acoustic ray γ. We expect T_{γ} to be dominated by the time spent in the subsurface region where the sound speed is much smaller than in the interior. While the path of the ray varies with rotation, δ_{n} remains approximately proportional to the mean density as shown by Reese et al. (2008). On the other hand, the δ_{ℓ} spacing depends also on the second derivatives of the sound speed transverse to the ray γ, an information integrated all along the ray. As long as the path of the stable ray goes through the central region of the star, the island mode frequencies should be sensitive to the chemical stratification and thus the age of the star. However, after the bifurcation at Ω/Ω_{K} ≃ 0.26, the ray path progressively avoids the central region and the island modes do not contain this information anymore. Another interesting property of δ_{ℓ} (or δ_{ℓ}/δ_{n}) is that it is very sensitive to rotation as long as Ω/Ω_{K} ≤ 0.35. Finally, for high rotation rates (Ω/Ω_{K} ≥ 0.40), the value of δ_{m}, that can be detected through the irregularity of multiplets, also gives an information on rotation since it is proportional to (Eq. (58)) where the distance of the ray to the rotation axis d strongly depends on the rotation rate.
6. Conclusions
In this paper, we derived an asymptotic formula for frequencies that predicts and describes regular spacings in the pmode spectrum of rapidly rotating stars. The derivation relied on finding a stable periodic solution of the acoustic ray dynamics, and obtaining an expression for the modes that are localized around this ray, the socalled island modes. The method thus provides a formula for the island modes frequencies, as well as a formula for the mode spatial distributions. We compared these semianalytical formulas with results from numerical computations of highfrequency oscillations in rotating polytropic stellar models. The frequency spectrum is characterized by the three spacings δ_{n}, δ_{ℓ} and δ_{m}. The agreement was shown to be good for δ_{n} and δ_{ℓ} at almost all rotation rates, while δ_{m} shows significant discrepancies at low rotation rates. The spacing δ_{n} stays almost constant at all rotation rates with a value that is close to half the large frequency separation of the nonrotating model. On the other hand, the rapid decrease of δ_{ℓ} strongly modifies the spectrum’s organization up to Ω/Ω_{K} ≃ 0.4, while above that rotation rate, δ_{ℓ} remains approximately constant. For such high rotation rates, the spacing δ_{m} nearly vanishes, thus in the observer’s frame the evolution of the whole spectrum is dominated by the advection term mΩ. We have also seen that the combined evolution of these frequency spacings with rotation leads to particular events such as true or near degeneracies, that can significantly simplify the spectrum. In addition to these new insights on the evolution of the island mode spectrum with rotation, the asymptotic theory provides semianalytical formulas for the regular spacings, in particular simple formulas for δ_{n} and δ_{ℓ}.
The present asymptotic theory should be useful for different aspects of stellar seismology in the presence of rapid stellar rotation. The regular frequency spacings are potentially observable and our results provide guidance to look for them in data. While investigations dedicated to the search for regularities are necessary (e.g. Lignières et al. 2010), we expect that the easiest quantities to detect in an island mode spectrum are δ_{n} at any rotation rate, 2mΩ at small rotation rates, and Ω at high rotation. For modeling pulsations, the asymptotic theory provides a new approach, complementary to numerical computations. One of its advantages is to give a quick estimate of frequency spacings for a given stellar model, which in turn can be used to search for patterns in numerically computed spectra. In the same spirit, the semianalytical amplitude distributions might provide a useful approximation for calculating mode visibilities and spectral signatures.
The asymptotic theory in itself can be improved and extended in various ways. We have already mentioned that the method needs to be refined at the rotation rate where the bifurcation of the stable ray occurs, using uniform approximations of the ray dynamics. It would also be interesting to predict analytically the rotation rate of the bifurcation for a given sequence of stellar models. The present method also assumes that the modes are governed by local dynamics around the stable ray. This assumption can be tested with a numerical EBK method applied to the tori of the stable island (Bohigas et al. 1993), although this method is complicated to implement in pratice. In this paper, we left aside the determination of the actual number of modes that are described by the asymptotic theory. An estimate of such a number can be obtained by computing systematically the phase space volume of stable islands for different rotations (e.g. Lignières & Georgeot 2009). Then, knowing the value of δ_{ℓ} that gives the mean distance between island mode frequencies, or the mean density of these modes, one could compute the ℓ^{max} of the modes that satisfy our formulas. Another aspect that we have not modeled is avoided crossing in spite of the fact that it will induce important deviations, especially at low frequencies. Strong gradients of the sound speed will also produce deviations from the asymptotic theory. A technique called raysplitting, that has already been used successfully in quantum chaos (Blümel et al. 1996), could account for this effect. Finally, a similar technique could be applied to asymptotic gravity modes, that were shown recently to have connections with ray theory (Ballot et al. 2011).
Acknowledgments
We thank J. Ballot for his help at various stages of this work. We also thank the ANR project SIROCO for funding and CALMIP (“CALcul en MIdiPyrénées”) for the use of their supercomputer. M.P., F.L. and D.R.R. acknowledge the KITP staff of UCSB for their warm hospitality during the research program “Asteroseismology in the Space Age”. This research was supported in part by the National Science Foundation under Grant No. PHY1125915. D.R.R. acknowledges financial support through a postdoctoral fellowship from the “Subside fédéral pour la recherche 2011”, University of Liège.
References
 Aerts, C., ChristensenDalsgaard, J., & Kurtz, D. 2010, Asteroseismology, Astronomy and Astrophysics Library (Springer) [Google Scholar]
 Arfken, G., & Weber, H. 2005, Mathematical Methods For Physicists (Elsevier) [Google Scholar]
 Arnol’d, V. 1989, Mathematical methods of classical mechanics, Graduate texts in mathematics (SpringerVerlag) [Google Scholar]
 Babich, V., & Buldyrev, V. 1991, Shortwavelength diffraction theory: asymptotic methods, Springer series on wave phenomena (SpringerVerlag) [Google Scholar]
 Baglin, A., Auvergne, M., Barge, P., et al. 2006, in ESA Spec. Publ. 1306, eds. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, 33 [Google Scholar]
 Ballot, J., Lignières, F., Prat, V., Reese, D. R., & Rieutord, M. 2011 [arXiv:1109.6856] [Google Scholar]
 Berry, M. V., & Robnik, M. 1984, J. Phys. A Math. Gen., 17, 2413 [NASA ADS] [CrossRef] [Google Scholar]
 Blümel, R., Antonsen, Jr., T. M., Georgeot, B., Ott, E., & Prange, R. E. 1996, Phys. Rev. Lett., 76, 2476 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Bogomolny, E. 2006, in Frontiers in Number Theory, Physics, and Geometry I, eds. P. Cartier, B. Julia, P. Moussa, & P. Vanhove (Springer Berlin Heidelberg), 3 [Google Scholar]
 Bohigas, O., Tomsovic, S., & Ullmo, D. 1993, Phys. Rep., 223, 43 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Brack, M. 2001, in Atomic Clusters and Nanoparticles, eds. C. Guet, P. Hobza, & F. Spiegelman, 161 [Google Scholar]
 CohenTannoudji, C., Diu, B., & Laloë, F. 1973, Mécanique quantique I, Collection enseignement des sciences (Hermann) [Google Scholar]
 Cvitanović, P., Artuso, R., Mainieri, R., Tanner, G., & Vattay, G. 2010, Chaos: Classical and Quantum (Copenhagen: Niels Bohr Institute) [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]
 Gough, D. O. 1993, in Astrophysical Fluid Dynamics – Les Houches 1987, eds. J.P. Zahn, & J. ZinnJustin, 399 [Google Scholar]
 Gutzwiller, M. 1990, Chaos in classical and quantum mechanics, Interdisciplinary applied mathematics (SpringerVerlag) [Google Scholar]
 Koch, D. G., Borucki, W. J., Basri, G., et al. 2010, ApJ, 713, L79 [NASA ADS] [CrossRef] [Google Scholar]
 Kogelnik, H., & Li, T. 1966, Appl. Opt., 5, 1550 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Landau, L., & Lifshitz, E. 1977, Quantum mechanics: nonrelativistic theory (ButterworthHeinemann) [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]
 Lignières, F., Georgeot, B., & Ballot, J. 2010, Astron. Nachr., 331, 1053 [NASA ADS] [CrossRef] [Google Scholar]
 Mantegazza, L., Poretti, E., Michel, E., et al. 2012, A&A, 542, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Miller, W. H. 1975, J. Chem. Phys., 63, 996 [NASA ADS] [CrossRef] [Google Scholar]
 Monnier, J. D., Zhao, M., Pedretti, E., et al. 2007, Science, 317, 342 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Ott, E. 2002, Chaos in Dynamical Systems (Cambridge University Press) [Google Scholar]
 Pasek, M., Georgeot, B., Lignières, F., & Reese, D. R. 2011, Phys. Rev. Lett., 107, 121101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Percival, I. C. 1973, J. Phys. B At. Mol. Phys., 6, L229 [NASA ADS] [CrossRef] [Google Scholar]
 Reese, D. 2008, J. Phys. Conf. Ser., 118, 012023 [NASA ADS] [CrossRef] [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, 506, 189 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Royer, F. 2009, in Lecture Notes in Physics The Rotation of Sun and Stars (Berlin Springer Verlag), 765, 207 [Google Scholar]
 Schomerus, H., & Sieber, M. 1997, J. Phys. A Math. Gen., 30, 4537 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Tassoul, M. 1980, ApJS, 43, 469 [NASA ADS] [CrossRef] [Google Scholar]
 Tureci, H. E., Schwefel, H. G. L., Stone, A. D., & Narimanov, E. E. 2002, Opt. Exp., 10, 752 [NASA ADS] [CrossRef] [Google Scholar]
 Vagov, A., Schomerus, H., & Zalipaev, V. V. 2009, Phys. Rev. E, 80, 056202 [NASA ADS] [CrossRef] [Google Scholar]
 Vandakurov, Y. V. 1967, AZh, 44, 786 [NASA ADS] [Google Scholar]
 Zalipaev, V. V., Kusmartsev, F. V., & Popov, M. M. 2008, J. Phys. A Math. Gen., 41, 065101 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Wave equation in the meridional plane
In this section, we derive the twodimensional wave equation in the meridional plane of the star. We start from the threedimensional wave equation (Eq. (4)) in spherical coordinates (r,θ,φ), and from the expression of the mode amplitude as Ψ = Ψ_{m}exp(imφ) we obtain the following wave equation: (A.1)where θ is the colatitude. We want to cancel out of this equation the terms multiplied by first order derivatives of Ψ_{m} in order to obtain a twodimensional Helmholtzlike equation. We thus introduce the ansatz Ψ_{m}(r,θ) = β(r,θ)Φ_{m}(r,θ) in Eq. (A.1) to obtain that β(r,θ) must satisfy (A.2)to cancel out terms in , and (A.3)to cancel out terms in . This leads to the solution (A.4)and we choose B = 1 to yield Eq. (5).
Appendix B: Expression of the imaginary part of Γ
In this section, we derive the expression for the imaginary part of Γ using the Wronskian of Eq. (27). First, we write Eq. (27) as (B.1)Then, we shall use the Abel’s differential equation identity that says that for an equation of the following type (B.2)with the Wronskian defined as (B.3)and y_{1}, y_{2} the two independent solutions of Eq. (B.2), there is the identity (B.4)With the change of variable we can transform Eq. (B.1) into a form similar to Eq. (B.2) with P(u) = 0. It follows that the Wronskian stays constant for all u: (B.5)and then, coming back to the s variable, we obtain that (B.6)or equivalently (B.7)stays constant for all s such that . From Eqs. (B.6), (B.5) and the fact that Γ = p/z we thus obtain that (B.8)yielding (B.9)and therefore (B.10)Hence this shows that Im [Γ] keeps a constant sign along the ray when s varies. Thus if Γ is chosen such that it localizes the function at s_{0}, then the function will stay localized around the ray for all s. Since the localization is obtained for Im [Γ] > 0 or equivalently i, we choose the Wronskian such that (B.11)
Appendix C: Analogy between the function z and the deviation of nearby rays
In this section, we show that the equation satisfied by the function z, Eq. (27), is the same as the equation describing the deviation of two nearby rays. It will thus be possible to compute the evolution of z from the ray dynamics. To be an acoustic ray of the system, the ray γ must correspond to an extremum of the action S defined as (C.1)The length element dσ can be expressed in the coordinate system centered on the ray γ defined in Sect. 3.1 as (C.2)yielding (C.3)The action function S can therefore be expressed in this coordinate system as (C.4)with the corresponding Lagrangian (C.5)The Lagrangian L becomes for ξ small: (C.6)where and derivatives in ξ are evaluated at ξ = 0. As is known from classical mechanics, for the action to be extremum, the associated Lagrangian must satisfy the EulerLagrange equation (C.7)Keeping only quadratic terms in ξ from this equation yields (C.8)This equality must be valid for all ξ so we can deduce that: (C.9)and (C.10)We thus obtain the equation for the deviation between nearby rays as (C.11)and this equation is the same as Eqs. ((27), (B.1)) that are satisfied by the function z.
Appendix D: Eigenvectors and eigenvalues of the monodromy matrix M
In this section, we derive the necessary formulas to express the stability angle α and the function Γ in terms of the elements of the monodromy matrix M. The characteristic polynomial of a 2 × 2 matrix M is
(D.1)Then, it is known that a monodromy matrix is a symplectic matrix, which implies that det(M) = 1. For a stable ray, we have  Tr(M)  < 2, and thus the roots of this polynomial become (D.2)We thus obtain the expression for the stability angle α as (D.3)and we obtain Eq. (36).
The eigenvector v_{1} will then be given by the following equation: (D.4)which corresponds to Recalling that Γ(s) = p(s)/z(s) we write (D.5)and thus we obtain that (D.6)This yields (D.7)and with the same procedure we obtain the expression of Γ_{2}: (D.8)
All Tables
Island mode quantum numbers n, ℓ, m of numerically computed modes, corresponding to n_{s} ∈ [21,25] , ℓ_{s} ∈ {0,1,2,3}, m_{s} ∈ [−ℓ_{s},ℓ_{s}] in terms of spherical mode quantum numbers.
All Figures
Fig. 1 Pressure amplitude on a meridian plane for a polytropic stellar model, with d the distance to the rotation axis and ρ_{0} the equilibrium density. The mode shown corresponds to n = 50, ℓ = 1 and m = 1 at a rotation rate of Ω/Ω_{K} = 0.300, where is the limiting rotation rate for which the centrifugal acceleration equals the gravity at the equator, M being the stellar mass and R_{eq} the equatorial radius. Colors/grayness (for respectively the colour/black and white version of the figure) denote pressure amplitude, from red/gray (maximum positive value) to blue/black (minimum negative value) through white (null value). The thick black line is the ray γ located in the center of the main stable island. (This figure is available in color in the electronic form.) 

In the text 
Fig. 2 Poincaré Surface of Section (PSS) at the rotation rate Ω/Ω_{K} = 0.589 for quantum number m = 0. Each dot corresponds to the crossing of an acoustic ray with the equatorial halfplane in the (r/R_{eq},k_{r}/ω) phase space, r being the radial coordinate and k_{r} the associated momentum. R_{eq} is the equatorial radius and ω the mode frequency. Red/dark gray denotes a chaotic ray, green/light gray a whispering gallery ray, blue/black a stable island ray (see text). Upper inset is a closeup of the main stable island. (This figure is available in color in the electronic form.) 

In the text 
Fig. 3 Left: schematic representation of the mode labeling/quantum numbers used for island modes. Right: illustration of the relation between the quantum numbers of spherical and island modes. The quantum number n is the number of nodes in the longitudinal direction of the island mode (along the ray γ), ℓ is the number of nodes in the transverse direction of the mode (transverse to γ) and m is the number of azimuthal nodes. n_{s}, ℓ_{s} and m_{s} are the quantum numbers of spherical modes. In the right figure, the multiplets of island modes correspond to diagonal colored bands, whereas the multiplets of spherical modes would have the form of vertical bands. The plotted mode on the left corresponds to n = 46, ℓ = 1 and m = 0. (This figure is available in color in the electronic form.) 

In the text 
Fig. 4 Comparison of frequency spacings δ_{n,ℓ} between island modes, computed from numerical simulations and semianalytical formulas for different values of Ω/Ω_{K} (the frequency spacings are normalized by with R_{p} the polar radius). Circles: δ_{n}, triangles: δ_{ℓ}, orange/gray: semianalytical results, blue/black: numerical results. Numerical results correspond to different sets of and values, with n ∈ [42,51] , ℓ ∈ { 0,1 } for m = 0; and n ∈ [42,51] , ℓ = 0, n ∈ { 42,44,46,48,50 } , ℓ = 1 for m ∈ { − 1,1 } . Upper panel: m = 0. Lower panel: m = ± 1. (This figure is available in color in the electronic form.) 

In the text 
Fig. 5 Comparison of frequency spacing δ_{m} between island modes, computed from numerical simulations and semianalytical formula Eq. (53) for different values of Ω/Ω_{K} (the frequency spacings are normalized by with R_{p} the polar radius). Orange/gray: semianalytical results of Eq. (53) for n = 42 and m = 1, blue/black: numerical results. The numerical results correspond to different sets of δ_{m} for n ∈ [42,51] , ℓ = 0 and m ∈ { − 1,1 } . (This figure is available in color in the electronic form.) 

In the text 
Fig. 6 Frequency spacings δ_{m}(ℓ = 0) = ω_{n,ℓ = 0,m} − ω_{n,ℓ = 0,m = 0}, normalized by ω_{p}, as a function of m/. The integers n, ℓ and m are the quantum numbers of island modes. The rotation rate is Ω/Ω_{K} = 0.419. Blue/black dots: numerical results. Orange/gray dashed line: semianalytical results for the last term in Eq. (58). Numerical modes included are for quantum numbers n ∈ {42,44,46,48,50}, ℓ = 0 and m ∈ [−3,3] . (This figure is available in color in the electronic form.) 

In the text 
Fig. 7 Examples of normalized amplitudes distributions (real part of ) on the equator as a function of position r/R_{eq} (where R_{eq} is the equatorial radius of the stellar model) for semianalytical and numerical modes. Modes displayed are for Ω/Ω_{K} = 0.300 (left column), and Ω/Ω_{K} = 0.707 (right column). Quantum numbers (n,ℓ,m) are as follows. Left column (50,0,0), (50,1,0), (50,0,1), (50,1,1). Right column (50,0,0), (50,1,0), (50,2,0), (50,3,0). Blue/black continuous line: numerical results; orange/gray dashed line: semianalytical results. (This figure is available in color in the electronic form.) 

In the text 
Fig. 8 Normalized frequencies ω/ω_{p} as a function of normalized rotation velocity Ω/Ω_{K} in the observer’s frame. Different colors correspond to different values of  m  . Blue: m = 0, green:  m  = 1, red:  m  = 2, cyan:  m  = 3. Modes included are for the quantum numbers n_{s} ∈ [21,23] , ℓ_{s} ∈ [0,3] , m_{s} ∈ [ − ℓ_{s},ℓ_{s}] (see Table 1). (This figure is available in color in the electronic form.) 

In the text 
Fig. 9 Normalized frequencies ω/ω_{p} as a function of normalized rotation velocity Ω/Ω_{K}. The quantum numbers of the displayed modes are: ℓ = 0, n ∈ [43,46] and ℓ = 1, n ∈ [42,44] with m = 0. Blue: ℓ = 0, green: ℓ = 1. For clarity, the corresponding degrees of the spherical harmonics are also written. We have outlined the large frequency separation Δ, the small frequency separation δ = ω_{ns,ℓs} − ω_{ns − 1,ℓs + 2} and the δ_{n}, δ_{ℓ} spacings with arrows. (This figure is available in color in the electronic form.) 

In the text 
Fig. 10 Normalized frequencies ω/ω_{p} as a function of normalized rotation velocity Ω/Ω_{K}. The quantum numbers of the displayed modes are: n = 44, ℓ = 0, m ∈ { − 2, − 1,1,2 } with the n ∈ [43,46] , ℓ = 0, m = 0 modes. Blue: m = 0, red:  m  ∈ { 1,2 } . For clarity, the corresponding degrees of the spherical harmonics are also written. Continuous line: frequencies in the observer frame, dashed line: in the rotating frame. We have outlined the large frequency separation Δ and the δ_{m} spacing with arrows. (This figure is available in color in the electronic form.) 

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.