Issue 
A&A
Volume 601, May 2017



Article Number  A130  
Number of page(s)  14  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201321264  
Published online  22 May 2017 
Frequency regularities of acoustic modes and multicolour mode identification in rapidly rotating stars
^{1} Institut d’Astrophysique et Géophysique de l’Université de Liège, Allée du 6 Août 17, 4000, Liège, Belgium
email: daniel.reese@obspm.fr
^{2} School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
^{3} LESIA, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, UPMC Univ. Paris 06, Univ. Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, 92195 Meudon, France
^{4} Université de Toulouse, UPSOMP, IRAP, 31028 Toulouse, France
^{5} CNRS, IRAP, 14 avenue Édouard Belin, 31400 Toulouse, France
^{6} GEPI, Observatoire de ParisMeudon, CNRS, Université Paris Diderot, 92125 Meudon Cedex, France
^{7} High Altitude Observatory, National Center for Atmospheric Research, Boulder, CO 80307, USA
Received: 8 February 2013
Accepted: 27 January 2017
Context. Mode identification has remained a major obstacle in the interpretation of pulsation spectra in rapidly rotating stars. This has motivated recent work on calculating realistic multicolour mode visibilities in this type of star.
Aims. We would like to test mode identification methods and seismic diagnostics in rapidly rotating stars, using oscillation spectra that are based on these new theoretical predictions.
Methods. We investigate the autocorrelation function and Fourier transform of theoretically calculated frequency spectra, in which modes are selected according to their visibilities. Given that intrinsic mode amplitudes are determined by nonlinear saturation and cannot currently be theoretically predicted, we experimented with various adhoc prescriptions for setting the mode amplitudes, including using random values. Furthermore, we analyse the ratios between mode amplitudes observed in different photometric bands to see up to what extent they can identify modes.
Results. When nonrandom intrinsic mode amplitudes are used, our results show that it is possible to extract a mean value for the large frequency separation or half its value and, sometimes, twice the rotation rate, from the autocorrelation of the frequency spectra. Furthermore, the Fourier transforms are mostly sensitive to the large frequency separation or half its value. The combination of the two methods may therefore measure and distinguish the two types of separations. When the intrinsic mode amplitudes include random factors, which seems more representative of real stars, the results are far less favourable. It is only when the large separation or half its value coincides with twice the rotation rate, that it might be possible to detect the signature of a frequency regularity. We also find that amplitude ratios are a good way of grouping together modes with similar characteristics. By analysing the frequencies of these groups, it is possible to constrain mode identification, as well as determine the large frequency separation and the rotation rate.
Key words: stars: oscillations / stars: rotation / stars: interiors / stars: variables:δScuti
© ESO, 2017
1. Introduction
One of the major obstacles in interpreting the acoustic frequency spectra of rapidly rotating stars is mode identification, i.e. finding the correspondence between theoretically calculated modes and observed pulsations. Several reasons make it difficult to match the two. First and foremost is the lack of simple frequency patterns, as is found in solarlike stars. Indeed, rapid rotation leads to complex spectra with overlapping classes of pulsation modes, each with an independent frequency organisation (Lignières & Georgeot 2008). Next comes the whole problem of mode amplitudes. Most rapid rotators tend to be massive or intermediate mass stars where modes are predominantly excited by the κ mechanism. This leads to nonlinear saturation and coupling between modes, making it nearly impossible to predict the amplitudes with current theory. Further difficulties include avoided crossings between modes, and all of the theoretical and numerical challenges associated with rapid rotation. From an observational point of view, the high quality data from the CoRoT (Baglin et al. 2009; Auvergne et al. 2009) and Kepler (Borucki et al. 2009) space missions have painted a new picture of δ Scuti stars through the detection of hundreds of pulsation modes (Poretti et al. 2009; Balona et al. 2012). In a similar way, the number of detected modes has also increased for stars from other classes of rapidly rotating pulsators, and along with it the complexity of the spectra (e.g. Uytterhoeven et al. 2011). Consequently, most asteroseismic analyses have focused on interpreting the general characteristics of these spectra rather than identifying individual modes.
Various strategies have been devised to identify modes. One can, for instance, search for frequency patterns appropriate for rapid rotation. The background for this search is the discovery of asymptotically uniform frequency spacings in the numerically computed spectra of uniformly rotating polytropic models (Lignières et al. 2006; Reese et al. 2008) and differentially rotating realistic selfconsistent field (SCF) models (Reese et al. 2009a). These uniform spacings have also been modelled through asymptotic semianalytical formulas (Pasek et al. 2012). In observed spectra, recurrent frequency spacings that may correspond to the large separation or half its value have been found in some stars (García Hernández et al. 2009; García Hernández et al. 2013; Paparó et al. 2016). Moreover, García Hernández et al. (2015) show that mean density estimates based on this type of spacings (obtained via a scaling relation similar to the one in Reese et al. 2008, but based on SCF models) are compatible with independent mass and radii measurements obtained for δ Scuti stars in binary systems. Nonetheless, it is expected that various effects may contribute to hide these regular frequency patterns. First, as mentioned before, the full spectrum is a superposition of subspectra corresponding to different classes of modes and some of the uniform spacings only concern one class. This complicates their detection in the full spectrum. Also, owing to their asymptotic nature, these spacings might not be relevant to analyse the low to moderate (up to radial order n ~ 10) frequency domain, typical of most rapidly rotating pulsators. A third effect that may come into play is the presence of mixed modes in evolved stars and/or sharp sound speed gradients since they can potentially modify the regular spacings. Finally, mode selection effects that are due to the nonlinearly determined intrinsic mode amplitudes could affect the detectability of the regular patterns. As a first attempt, Reese et al. (2009b) developed a strategy to find these frequency spacings but ran into difficulties when including chaotic modes, which come from another class of modes. Lignières et al. (2010) addressed the same question with encouraging results but their analysis was restricted to the asymptotic regime and relied on simplifying assumptions regarding the spectrum of chaotic modes and the mode visibilities. In this paper, our first goal is to search for regular frequency spacings in the most realistic synthetic spectra available, using relevant frequency ranges and accurate visibility calculations. While they can provide guidance to a similar search in real data, we already know that these results must be taken with caution since the intrinsic mode amplitudes used in this paper are not realistic, but based on adhoc prescriptions.
Another strategy, which avoids this difficulty, is to constrain the identification through multicolour photometric or spectroscopic observations. Indeed, one can measure the amplitudes and the phases of a given pulsation mode in different photometric bands and then compare these by calculating amplitude ratios or phase differences. The geometry of the modes will then lead to different characteristic signatures. One important advantage of this method is that these signatures are independent of the intrinsic mode amplitude. In a similar fashion, the oscillatory movements induced by a pulsation mode cause Doppler shifts that show up as variations in the shape of spectroscopic absorption lines, known as line profile variations or LPVs. These variations are then directly related to the geometry of the mode. By comparing theoretical predictions with observations, we can then constrain the mode’s identification. In the following, we focus on multicolour mode identification.
Most of the previous theoretical investigations of amplitude ratios and phase differences have been based on mode calculations that approximate the effects of rotation. For instance, DaszyńskaDaszkiewicz et al. (2002) used a perturbative approach, whereas Townsend (2003) and DaszynskaDaszkiewicz et al. (2007) applied the traditional approximation (that is typically used when calculating gravitoinertial modes). It is only recently that these predictions have started to fully take into account the effects of rotation. First, Lignières et al. (2006) and Lignières & Georgeot (2009) calculated diskintegration factors for pulsation modes calculated in fully deformed polytropic models. Given the simplified nature of these calculations, it was not possible to calculate associated amplitude ratios or phase differences. More recently, Reese et al. (2013b, hereafter Paper I) calculated mode visibilities in realistic models of rapidly rotating stars. This work relied on a grid of Kurucz atmospheres to calculate realistic emerging intensities, thereby taking into account limb and gravity darkening. Paper I took into account the Lagrangian variations of temperature and effective gravity as well as the surface distortion from the pulsation modes. The pulsation modes were calculated using a 2D approach that takes into account centrifugal distortion. The main limitation was the adiabatic approximation, which leads to an unreliable estimate of the Lagrangian variations of the effective temperature. Nonetheless, first results were obtained in that work that will hopefully provide a qualitative insight, both into mode visibilities and multicolour amplitude ratios.
In the following two sections, we succinctly recall some of the main aspects of pulsation modes in rapidly rotating stars, as well as the basic principles behind the visibility calculations described in Paper I. We then examine the autocorrelation functions of theoretically calculated spectra that overlap the p and gmode domains. This is followed by a discussion on Fourier transforms of frequency spectra and how these complement autocorrelation functions. In Sect. 6, we show how multicolour photometric mode identification can be extended to rapidly rotating stars. Finally, a discussion concludes the paper.
2. Pulsation modes
Rapidly rotating models are calculated thanks to the selfconsistent field (SCF) method (Jackson et al. 2005; MacGregor et al. 2007). The resultant models represent zero age main sequence (ZAMS) stars with a cylindrical rotation profile, and hence a barotropic stellar structure, i.e. a structure where different thermodynamic quantities, such as the density, pressure, and temperature are constant on isopotential surfaces (deduced from the sum of the gravitational and centrifugal potentials). The pulsation modes are calculated using the twodimensional oscillation program (TOP Reese et al. 2006, 2009a). As described in Paper I, an improved treatment of the models and a slightly different mechanical boundary condition were necessary to obtain eigenfunctions appropriate for visibility calculations. Compared to Paper I, the frequency range has been extended both to higher and lower frequencies, including g modes although, for the most part, we focus on p modes.
As was shown in Lignières & Georgeot (2008, 2009), acoustic modes subdivide into several classes of modes as the rotation rate increases. Each class of mode has its own typical geometry and frequency organisation, whether regular or statistical. This behaviour stems from the gradual transition of the ray dynamics system going from being integrable to chaotic – a transition that causes different regions, associated with the different classes, to appear in the Poincaré section. Of particular interest are the island modes, the rotating counterpart to low degree modes. These modes focus around a periodic trajectory and are characterised by the quantum numbers , where is the number of nodes along the trajectory, the number of nodes perpendicular to the trajectory, and m the usual azimuthal order. As illustrated in the animation in Reese (2008), it is possible to transform these quantum numbers into the usual spherical quantum numbers (n,ℓ,m) and inversely, using the following relations based on what could be described as node conservation: (1)where ε corresponds to mode parity, i.e. ε = 0 for even modes, modes that are symmetric with respect to the equator, and ε = 1 for odd modes. We note, in particular, that the parity of corresponds to that of the mode. This simply corresponds to the fact that the node along the equatorial plane is treated as a pseudoradial node in island modes. Closely related to this is the fact that one set of values corresponds to two sets of (ℓ,m) values, depending on the parity of , as illustrated in Fig. 1. At this point, it is also useful to introduce the large frequency separation, Δ and the semilarge frequency separation, Δ/2: where ω is the frequency, indexed either by the spherical or island mode quantum numbers.
Fig. 1 Meridional crosssections of two island modes with the same values but opposing parity, as determined from the parity of . The colours indicate the Eulerian pressure perturbation divided by the square root of the equilibrium density. 
Recently, Pasek et al. (2012) found asymptotic expressions for the profile of island modes in the direction perpendicular to the periodic trajectory. Apart from a scaling, modes with the same values have the same transverse profile. This applies in particular at the stellar surface. Hence, modes with the same values have similar surface profiles in one hemisphere. They are either symmetric or antisymmetric with respect to the equator, depending on the parity of . We note that using ℓ instead of also enables us to select island modes of the same equatorial parity. In addition, the weak influence of the Coriolis force on highfrequency acoustic modes (see Reese et al. 2006) implies that ± m pairs of mode are nearly identical (with neardegenerate frequencies in the corotating frame). It follows that island modes with the same (ℓ,  m ) values will have similar surface profiles. They should then have similar visibilities (as described in the next section) and this property will play an important role in the mode identification methods described in the following.
3. Mode visibilities
Mode visibilities are obtained by perturbing the expression that gives the amount of energy radiated by a star to an observing instrument: (4)where d is the distance to the star, μ the cosine of the angle between the outward normal to the stellar surface and the observer’s direction, “Vis. Surf.” the surface visible to the observer, g_{eff} and T_{eff} the effective gravity and temperature, and I(μ,g_{eff},T_{eff}) the specific radiation intensity, multiplied by the instrument’s and/or filter’s transmission curve and integrated over the wavelength spectrum. The perturbed expression is (5)where δ denotes a Lagrangian perturbation and ℜ { ... } the real part. Variations caused by fluctuations to the boundary between the visible and hidden side of the star lead to second order effects and are therefore neglected. The Lagrangian perturbation to the specific intensity, δI, is calculated as follows: (6)The quantities I, ∂lnI/∂lnT_{eff}, ∂lnI/∂lng_{eff}, and ∂I/∂μ are calculated using a grid of Kurucz atmospheres that spans the relevant effective temperature and gravity ranges (see references and more details in Paper I). This enables us to take into account both limb and gravity darkening. The quantities δT_{eff}/T_{eff}, δg_{eff}/g_{eff} and δμ are deduced from the surface profiles of the pulsation modes as described in Paper I. We note that, since the pulsation modes are calculated using the adiabatic approximation, δT_{eff}/T_{eff} is not accessible and is therefore approximated by δT/T. As has been pointed out in Dupret et al. (2002) and Dupret et al. (2003), this can lead to poor results, since δT/T is not reliable in the outer layers when calculated adiabatically. A full nonadiabatic calculation would remedy this problem but is beyond the scope of this paper. Finally, the term δ(dS), which intervenes in the second integral in Eq. (5), is also deduced from the surface profiles of the eigenmodes, as described in Paper I. Hence, the geometrical distortions of the stellar surface, induced by the pulsation modes are fully taken into account. We note that the centrifugal deformation is also taken into account, both in the pulsation and visibility calculations.
4. Autocorrelation function of the frequency spectra
4.1. General description
We first look at the autocorrelation function of the frequency spectra using mode visibilities in CoRoT’s photometric band to set the amplitudes. These spectra have been calculated in 2 M_{⊙} stellar models with rotation rates ranging from 0.0 to 0.8 Ω_{C}, where Ω_{C} is the critical rotation rate^{1}, the rotation rate at which the centrifugal force exactly compensates the equatorial gravity, g_{eq} (e.g. Jackson et al. 2005): (7)We extract the N most visible modes in the CoRoT photometric band. The selected frequencies are given the same amplitude and are convolved by a Gaussian profile with a width of 1/15 d^{1}, after which we calculate the autocorrelation function. The width of the Gaussian profile had to be carefully selected. Indeed, a smaller width leads to signatures that are less clear given the variations of the large frequency separation in the frequency range considered here, whereas a larger width leads to a loss of accuracy in the position of peaks in the autocorrelation function.
To select modes according to how visible they are, we not only need their visibilities (as computed in Sect. 3), but also their intrinsic amplitudes, since the observed amplitude is proportional to the product of the two. However, determining the intrinsic amplitude of a mode in a classical pulsator is an unsolved theoretical problem (Goupil et al. 2005, and references therein). Accordingly, we experiment with the following adhoc ways of defining the intrinsic amplitude:

inclusion of random factors: the mode visibilities from the previous case are multiplied by random numbers;
4.2. Normalisation of the maximal displacement
Figure 2 shows the autocorrelation functions of a spectrum spanning seven large frequency separations from (n,ℓ,m) = (2,1,0) to (9,1,0), corresponding to a 2 M_{⊙} stellar model rotating at 0.5 Ω_{C}. Autocorrelation functions have been computed for four different inclination angles and for various amplitude thresholds decreasing from N = 10 to 100. We note that for the poleon configuration i = 0°, only axisymmetric modes are visible. Accordingly, high values of N implicitly lead to the assumption that ℓ> 3 modes are visible, which may be somewhat optimistic and should be treated with caution (although some publications suggest that these modes can sometimes be detected, see, for example, Poretti et al. 2009). Very clear signatures of the large frequency separation, Δ, and half its value show up for i = 0° and i = 30°. These signatures are caused by the dominant presence of island modes. We also observe that the Δ/2 signature disappears at large inclination angles (see the fourth column in Fig. 2). This is due to the cancellation of antisymmetric island modes seen in near equatoron configurations.
Another important feature of the autocorrelation function is a small peak at twice the rotation rate, which is mostly seen for high values of i and using large numbers of selected modes. This is caused by the frequency difference between prograde modes with m = −1 and their retrograde counterparts, m = 1. Indeed, as already mentioned, the weak effect of the Coriolis force induces a neardegeneracy of ± m modes in the corotating frame, regardless of whether they are island or chaotic in nature. This produces pairs of frequencies separated by 2mΩ in the observed spectra of a uniformly rotating star. Their visibilities are also very similar. The 2Ω peak is then due to m = ± 1 pairs, the most visible nonaxisymmetric modes.
Fig. 2 Autocorrelation functions of acoustic spectra in a 2 M_{⊙} stellar model rotating at 0.5 Ω_{C}. The spherical radial orders of the modes, n, range from 2 to 9, therefore spanning 7 large frequency separations. The 4 columns correspond to 4 different inclinations, i = 0° being a poleon configuration. Each row corresponds to a different number of included modes. For example, in the top row, only the 10 most visible modes are included in the frequency spectra before calculating the autocorrelation function. The vertical dotted and dashed green lines give the rotation rate, Ω, and twice its value. The vertical dashed and continuous red lines indicate the large frequency separation, Δ, and half its value (see Eqs. (2) and (3)). 
Figure 3 shows what happens with a model rotating at 0.7 Ω_{C}. This time, the large frequency separation is very similar to twice the rotation rate. As a result, the autocorrelation functions show very strong peaks at both Δ ≃ 2Ω and Δ/2 ≃ Ω, regardless of the inclination angle. This is not surprising since the corresponding regularities add up to produce strong peaks. The pulsation frequencies actually tend to cluster around points separated by Δ/2 ≃ Ω. While such regularities are then easier to spot, it is difficult to disentangle between changing the pseudoradial order, , and changing the azimuthal order. A similar coincidence occurs around 0.3Ω_{C}, where 2Ω is close to Δ/2. This also leads to strong peaks in the autocorrelation functions.
Fig. 4 Same as Fig. 2, except that the model is rotating at 0.7 Ω_{C} and the visibilities have been multiplied by random numbers between 1 and 100 (which therefore affects which modes are selected in the frequency spectra). 
It is interesting to observe that the frequency spacing Ω ≃ Δ/2 shows up quite strongly in Fig. 3, even in the equatoron (i = 90°) configuration. This is somewhat surprising since antisymmetric modes cancel themselves out, leaving only even modes that are spaced by Δ for fixed values of . Furthermore, prograde and retrograde modes with m = ± 1 are spaced by 2Ω. The explanation lies in the fact that if one considers a socalled multiplet of modes with the same values, at sufficient rotation rates, modes with consecutive m values are approximately separated by Ω, at least for small values of  m . Indeed, the advection term mΩ is much stronger than the frequency deviations in the corotating frame, which behaves as according to the numerical calculations of Reese et al. (2009a) and the analytical model of Pasek et al. (2012).
4.3. Inclusion of random factors
Here, we multiplied the mode visibilities from the previous case by random numbers, before selecting the highestamplitude modes and calculating the autocorrelation functions. The random numbers are between 1 and 100 and are uniformly distributed on a logarithmic scale. Figure 4 shows the resultant autocorrelation functions for a model rotating at 0.7 Ω_{C}. As can be seen, the signature of regularities are much less evident than previously. Nevertheless, some peaks still remain, for instance the peak around Δ for i = 30° with 30 modes and the peak around Δ/2, and to a lesser extent Δ, for i = 60° with 200 or 300 modes. Similar signatures occur around 0.3Ω_{C}, where 2Ω is close to Δ/2. However, apart from the coincidences between 2Ω and Δ or Δ/2, the multiplication of the mode amplitude by such a random factor makes the characteristic frequency separations much more difficult to extract.
4.4. Other normalisations
Figure 5 shows how the autocorrelation functions are modified when using the normalisation based on the kinetic energy. Qualitatively, this remains the same as Fig. 2. Some of the peaks stand out better, notably the 2Ω peaks for few modes. In contrast, the normalisation based on the mean surface displacement (not shown) gives poor results. The explanation for this is quite simple: by normalising by the mean surface displacement, gravity (or gravitoinertial) modes are no longer penalised. Furthermore, the (ω + mΩ)^{2} factor in Eq. (10) ends up amplifying them. Hence, the spectra of selected modes based on this normalisation are dominated by gravity modes, which do not follow the same pattern, thereby drowning out the Δ, Δ/2, Ω, and 2Ω signatures in most cases.
Fig. 5 Same as Fig. 2 (i.e. with a model rotating at 0.5Ω_{C}), but where the modes are normalised by the kinetic energy, multiplied by the appropriate power of (ω + mΩ). 
5. Fourier transform of the frequency spectrum
5.1. General description
Recently, García Hernández et al. (2009), García Hernández et al. (2013) and García Hernández et al. (2015) analysed the Fourier transforms of the frequency spectra of the δ Scuti stars HD 174936 and HD 174966, observed by CoRoT. They investigated what happens when the number of selected frequencies varies from a few tens to a few hundreds. Below, we apply the same procedure, but to our numerically calculated frequency spectra. We select modes according to their visibilities in CoRoT’s photometric band, then assign the same amplitude to the selected modes before calculating the Fourier transform of the resultant spectrum. To facilitate comparisons with the autocorrelation functions, we apply this technique to the frequency spectra spanning seven radial orders that were studied in the previous section (i.e. Figs. 2 and 3).
5.2. Normalisation of the maximal displacement
Figure 6 shows the squared modulus of the Fourier transform of frequency spectra in the model rotating at 0.5 Ω_{C}, for various numbers of selected modes and for four different inclinations. Taking the Fourier transform of a function that depends on frequency yields another function, which depends on time t; it is plotted here as a function of 1 /t to facilitate the identification of regularities. At low inclinations, peaks appear at Δ/2, with their forest of harmonics at Δ/4, Δ/6, etc. This is to be expected because the frequency spectrum behaves like a Dirac comb with a Δ/2 periodicity. At higher inclinations, a peak appears close to Δ (which is an expected regularity), but is shifted; we also recover some harmonics (especially at Δ/3), but not all of them. This is probably an effect of rotation that does not necessarily add peaks, but acts as a modulation of the amplitude of the Fourier transform. Indeed, we notice that rotation does not produce peaks at 2Ω or Ω. This is because, although there are recurrent frequency separations of 2Ω (or actually slightly smaller because of the Coriolis force), such separations are formed by pairs of frequencies rather than by a Dirac comb. Nevertheless, when the frequency spectrum is dominated by two similar subspectra, with the second one being identical to the first one, but shifted by 2Ω (this is what happens when m = ± 1 modes dominate), the Fourier transform of the full spectrum will be the Fourier transform of the subspectrum multiplied by cos^{2}(Ωt). This type of modulation can make some peaks disappear or slightly shift some broad peaks. In this kind of configuration, it is impossible to unambiguously detect the correct large separation with the Fourier transform only.
Fig. 6 Fourier transform of frequency spectra in a 2 M_{⊙} model rotating at 0.5 Ω_{C}. Each row corresponds to a different number of selected modes, as indicated on the left, and each column to a different inclination, as stated above. The red, continuous and discontinuous, vertical lines correspond to the large frequency separation, Δ, and various fractions of this value: Δ/4, Δ/3, and Δ/2. The green, dotted and dashed, vertical lines correspond to Ω and 2Ω, respectively. 
Figure 7 shows what happens with the model rotating at 0.7 Ω_{C}, where Δ nearly coincides with 2Ω. In this case, the frequency spectra take on a fairly simple form in which the frequencies cluster around points separated by Δ/2 ≃ Ω, regardless of inclination. This leads to strong peaks at Δ/2 and Δ/4 in the Fourier transforms, regardless of inclination.
5.3. Other normalisations
Figure 8 shows the effects of the alternative normalisations described in Sect. 4.1 on the Fourier transform of the spectra of the Ω = 0.5Ω_{C} and Ω = 0.7Ω_{C} models. The Fourier transforms continue to detect Δ/2 and its many harmonics in many cases, even when using random factors or the normalisation based on the mean surface displacement. However, the tests involving random factors benefit from the coincidence between Δ and 2Ω. In the absence of this type of coincidence, the Δ/2 signature is far less visible, except in a few cases where it still shows up. Overall, these tests confirm the robustness of the large and semilarge frequency separations, as detected by the Fourier transform, to different nonrandom normalisations, or when examining favourable cases where Ω coincides with Δ/2 or Δ.
Fig. 8 Same as Fig. 6, but where the mode normalisation includes random factors (top panel), or is based on the mean surface displacement (bottom panel). 
Overall success rate in finding similar modes via the multicolour mode identification scheme.
Overall, the Fourier transforms complement the autocorrelation functions quite nicely. Indeed, although both detect the separation Δ/2, only the autocorrelation functions are sensitive to 2Ω. As explained above, these frequency separations are produced by pairs of modes rather than by Dirac combs, thereby escaping detection by Fourier transforms, but not by the autocorrelation functions. Hence, this provides a simple way to distinguish between the two and to get a better grasp of the regularities present in the spectrum.
6. Multicolour mode identification
6.1. Method
We now turn our attention to multicolour mode identification. As emphasised in the introduction, multicolour mode signatures, such as the ratios between mode amplitudes observed in different photometric bands, do not depend on intrinsic mode amplitudes since these factor out. In Paper I, it was shown that amplitude ratios tend to be similar for island modes with the same (ℓ,  m ) values, even in the most rapidly rotating models. This is consistent with the fact that such modes have a similar surface structure (see Sect. 2). Hence, this raises the question as to whether we can select modes with similar properties by picking a reference mode at random in an oscillation spectrum and searching for the N other modes that produce the most similar amplitude ratios. In what follows, we use the following cost function to evaluate the proximity of a mode’s amplitude ratios to that of the reference mode: (11)where i is an index on the photometric band, d the number of photometric bands, the visibilities of the reference mode and V_{i} those of the mode being evaluated. The visibilities have been normalised using to avoid favouring a particular photometric band. Nonetheless, a small value for J implies that the V_{i}s are close to the s, and therefore that the amplitude ratios are similar.
This approach is different to that typically taken in other works such as DaszyńskaDaszkiewicz et al. (2002). Indeed, most authors compare observed amplitude ratios directly with theoretical ones. The approach described here consists in comparing observed amplitude ratios between each other. It thus bypasses limitations in the theoretical predictions. In the following, we nevertheless use theoretical amplitude ratios to test its validity.
6.2. Adiabatic case
We start with the (n,ℓ,m) = (9,4,1) mode in the model at 0.6 Ω_{C} as the reference mode and search for the nine other modes with the most similar amplitude ratios. The mode visibilities are calculated using the Geneva photometric system, which contains seven photometric bands. Before normalising the visibilities, we filter out modes where the overall visibility is more than 100 times smaller than that of the reference mode, thereby excluding most gravity modes. Figure 9 shows the results for i = 60°. These results are promising. Apart from the mode at 5.124Ω_{K}, all of the modes belong to the (ℓ,  m  ) = (4,1) family of island modes. This then leads to clear peaks in the autocorrelation function at 2Ω, Δ, and 2Δ, as well as combination peaks at Δ−2Ω, 2Δ−2Ω, and Δ + 2Ω. Another example of highly successful mode selection using the same strategy can be found in Fig. 4 of Reese et al. (2013a). However, not all cases work out so well.
Fig. 9 Amplitude ratios, frequency spectrum and ranks, associated autocorrelation function, and meridional crosssections of the reference mode (n,ℓ,m) = (9,4,1) and the 9 other most similar modes, for an inclination of i = 60°. In the first row, the left panel shows the amplitude ratios, the dashed line corresponding to the reference mode (we note that this line is hardly visible because it is mostly covered up by the solid lines from the other modes), the middle panel gives the frequencies, where a higher rank indicates a higher resemblance with the reference mode, 10 corresponding to the reference mode, and the right panel displays the autocorrelation function of the subspectrum of the selected frequencies. The second and third rows show the meridional crosssections of the selected modes. Beneath each mode, the values of (n,ℓ,m) are indicated when the mode corresponds to an island mode. Otherwise, only the azimuthal order, m, is provided. 
To get an overall picture of the mode selection method, we applied it to all of the identified island modes above a given threshold frequency in the pulsation spectrum of our 2 M_{⊙}, 0.6 Ω_{C} stellar model. Then, we quantified the mode selection success rate by finding the average number of island modes (excluding the reference mode), the average number of island modes with the same (ℓ,  m ) values as the reference mode, and the average number of island modes with the same values, all of which are subsequently divided by the total number of selected modes, i.e. N = 9, to get a success rate between 0 and 1.
Results are given in Cols. 2–4 of Table 1, for a threshold frequency ω_{threshold} = 8.0 Ω_{C}. In this case, the total number of island modes divided by the total number of modes is 0.0115. Thus, the much higher success rates obtained show that selecting modes according to similar amplitude ratios considerably increases the chances of finding island modes, provided the reference mode is an island mode. When the inclination is high, the method selects island modes of the same parity. In contrast, in near poleon configurations, it will select modes of both equatorial symmetries. This property explains the difference in the (ℓ,m) versus success rates. Repeating this test for different values of the threshold frequency shows that higher success rates are achieved for higher frequency modes, most probably because they are further into the asymptotic regime.
Overall, the method appears to be a promising way of choosing modes with similar surface distributions and hence quantum numbers. Regularities of the subspectrum of selected modes may then help to determine the azimuthal order, thanks to the frequency separations 2  m  Ω, and constrain the radial order. Nonetheless, we may still wonder if this kind of strategy will continue to work when nonadiabatic effects are taken into account. Below, we examine this question by analysing mode visibilities in which nonadiabatic effects are approximated.
6.3. Pseudo nonadiabatic effects
Nonadiabatic effects strongly modify the effective temperature variations and, hence, multicolour photometric signatures of pulsation modes. One may then wonder if these effects are able to mask the structural similarities between modes with similar quantum numbers and thus alter the promising results of the mode selection method obtained in the adiabatic case. In the absence of full nonadiabatic pulsation calculations in rapidly rotating δ Scuti stars, we use nonadiabatic calculations in the nonrotating case to derive approximate nonadiabatic visibilities in the rotating case.
In nonrotating stars, the relative effective temperature fluctuations are typically proportional to the radial displacement (12)where R^{⋆} is the radius of the nonrotating model. The quantity f_{T}exp(iφ_{T}) represents the nonadiabatic effects and generally depends on the degree ℓ and frequency ω/ Ω_{K} of the mode, as well as on the effective temperature and surface gravity. However, as illustrated in Fig. 10, it actually depends little on the harmonic degree for acoustic modes (i.e. at sufficiently high frequencies). We also note that, when described as a function of ω/ Ω_{K}, f_{T}exp(iφ_{T}) is only slightly affected by the effective gravity, except at high frequencies, where the effect is somewhat stronger. This can be seen by comparing the two models in Fig. 10.
Fig. 10 Quantities f_{T} and φ_{T} in two nonrotating models with masses 1.8 M_{⊙} and 2.0 M_{⊙}. The two models have virtually the same temperature (8143.9 K and 8137.3 K, respectively), but different log (g) values (4.1249 dex and 3.9593 dex, respectively). The different line styles correspond to different ℓ values, and the overlapping vertical dotted lines (which are hard to distinguish) indicate the position of the fundamental modes of both models. 
In a rotating star, the effective temperature and the surface gravity vary from pole to equator. Hence, to estimate nonadiabatic effects, we calculated, using the MAD code (Dupret 2001), f_{T}exp(iψ_{T}) /R^{⋆} for a set of nonrotating models along four evolutionary sequences, which span the effective temperature and surface gravity ranges of our rotating models, and for a set of modes spanning the relevant frequency range. We smoothed out the small differences in degree and interpolated according to frequency, effective temperature, and effective gravity, to have a value for f_{T}exp(iφ_{T}) /R^{⋆} at each colatitude for each frequency. The function  f_{T}  /R^{⋆} is shown in Fig. 11. This was subsequently multiplied by R(θ), the radius of the rotating model, which depends on the colatitude, θ, and then applied at the corotating pulsation frequency, (ω + mΩ)/Ω_{K}, to estimate the effective temperature variations arising in our pulsating rotating models. A quick look at the function  f_{T}  /R^{⋆} shows that it cannot simply be expressed as the product of a function that depends on θ alone and another function that depends on frequency alone, i.e. it is nonseparable in terms of these two variables. Accordingly, this distorts the θ dependence of δT_{eff}/T_{eff} as the frequency changes, thereby leading to increased scatter in the photometric signatures of a given (ℓ,  m ) family of island modes.
Fig. 11 Quantity  f_{T}  /R^{⋆} as a function of colatitude and frequency, for the 1.8 M_{⊙} rotating model. 
As in the previous section, the method to find similar modes is applied to all of the island modes above a frequency threshold of 8.0 Ω_{C} and its efficiency is quantified by computing an overall success rate. In Table 1, the adiabatic results for a rotating model, M = 1.8 M_{⊙}, Ω = 0.6 Ω_{C}, (see Cols. 5–7) are compared with the pseudo nonadiabatic ones (Cols. 8–10). In the present case, the number of island modes divided by the total number of modes is 0.0330. Table 1 shows that pseudo nonadiabatic success rates are worse than adiabatic ones, but the method remains efficient in selecting island modes when island modes are used as reference modes.
There are a number of cases where pseudo nonadiabatic calculations lead to similar or even better results than the adiabatic calculations, such as the case illustrated in Fig. 12. In the adiabatic case, five of the selected modes are not island modes, even though the reference mode is an island mode. In the pseudo nonadiabatic case, only one of the selected modes is not an island mode, and there is also one island mode with a different value of  m . The remaining modes correspond to pairs of (ℓ,  m  ) = (2,1) island modes. This gives rise to strong signatures at 2Ω, Δ, and 2Δ, as well as some combination peaks in the autocorrelation function.
Hence the technique appears to remain viable in the nonadiabatic case. Furthermore, phase differences (e.g. Balona & Evers 1999; Dupret et al. 2005) could also be exploited and may lead to higher success rates by providing supplementary constraints. This can, however, only be tested when full nonadiabatic pulsation calculations are available in rapidly rotating δ Scuti stars, since it is only in such conditions that phase differences can be obtained with a reasonable accuracy.
7. Conclusion
In this paper, we investigated different ways of applying the visibility calculations of Reese et al. (2013b) to the interpretation of pulsation spectra in rapidly rotating stars. Given the lack of a comprehensive theory on nonlinear mode saturation in δ Scuti stars, we tested various ad hoc assumptions to determine intrinsic mode amplitudes. As such, these results must be taken with caution but represent a first step towards interpreting pulsation spectra in these stars using realistic mode visibilities.
Fig. 12 Same as Fig. 9, but for the reference mode (n,ℓ,m) = (4,2,1) (and for both the adiabatic and pseudo nonadiabatic case). 
We first looked at the autocorrelation functions of frequency spectra in much the same way as was done in Lignières et al. (2010), but using the newer more realistic visibilities, applied to spectra of acoustic modes calculated in a selfconsistent way (rather than having the nonaxisymmetric modes estimated from the axisymmetric ones). Our results confirm those of Lignières et al. (2010) in the sense that it is possible to observe peaks corresponding to a mean value of the large frequency separation, Δ, half its value, Δ/2, the rotation rate, Ω, and twice the rotation rate, 2Ω, when conditions are favourable. These conditions are achieved when the frequencies span a large enough range to reinforce recurrent spacings, and when not too many modes are included. The orientation of the star is important since it will favour either Δ/2 for low inclinations (i.e. close to poleon), or Δ and 2Ω for high inclinations. The Ω spacing is visible owing to the usual multiplet structure at small rotation rates, then it disappears and becomes visible again at high rotation rates when a new of type multiplet forms. Of particular interest are the situations where 2Ω coincides with either Δ/2 or Δ, which occurs for models rotating at 0.3 Ω_{C} or 0.7 Ω_{C}, respectively. These types of situations lead to a very strong signature in the autocorrelation function owing to a simplification of the spectrum in which modes tend to cluster together. Although this is ideal for detecting these frequency separations, it also makes it more difficult to disentangle the two. Finally, we experimented with multiplying the intrinsic amplitudes by random numbers between 1 and 100, uniformly spread out on a logarithmic scale. This is used as a poor substitute for the effects of nonlinear saturation and mode coupling on the amplitudes in such stars. Despite this, a weak signature of the frequency separations remained in the favourable case, where Δ coincided with 2Ω. This gives us hope that it may be possible to detect, at least in favourable cases, these characteristic separations, and might explain the recent detection of recurrent spacings in δ Scuti stars (Mantegazza et al. 2012; Suárez et al. 2014; García Hernández et al. 2009; García Hernández et al. 2013; García Hernández et al. 2015).
We also looked at Fourier transforms of the frequency spectra. Although these also detected Δ/2 as a recurrent spacing, 2Ω escaped detection. This is because the latter is not formed by a Dirac comb, but rather by isolated pairs of frequencies. This difference between the two approaches is quite useful since it helps us distinguish between the two types of separations. Hence, we hope to interpret observed spectra by combining the two approaches and correctly identifying Δ/2 and 2Ω. Further tests based on different mode normalisations confirmed the robustness of both the autocorrelation functions and the Fourier transforms at detecting these characteristic frequency separations, as long as gravitoinertial modes were excluded from the analysis.
Finally, we turned our attention to multicolour mode identification. A key advantage of multicolour mode signatures, such as amplitude ratios, is that the intrinsic amplitudes factor out, thereby leaving a signature that is only sensitive to the geometric structure of the mode. Previous investigations into the matter had concluded that, due to the dependence of amplitude ratios on inclination and azimuthal order, it would be very difficult to identify modes in rapidly rotating stars from multicolour photometry alone (e.g. DaszyńskaDaszkiewicz et al. 2002; Townsend 2003). In contrast, we present more promising results. To achieve this, we apply a strategy that is different from the nonrotating case, a strategy which involves choosing a reference mode and searching for other modes with the most similar amplitude ratios. By repeating this procedure for different reference modes, we can group modes together into families with similar amplitude ratios. It turns out that these families have similar azimuthal orders and degrees, and their frequencies follow patterns that help to constrain the identification of these modes, as well as the rotation rate and the large frequency separation. Furthermore, by comparing modes between each other rather than with theoretical predictions, we bypass the current limitations with our theory such as the lack of full 2D nonadiabatic pulsation calculations, or the difficulties in modelling the interactions between such modes and convection. Nonetheless, we investigated, in an approximate way, nonadiabatic effects on multicolour mode amplitudes. We found that nonadiabatic effects do tend to increase the scatter between the multicolour photometric signatures of similar modes owing to its distortion of δT_{eff}/T_{eff} as a function of frequency. This makes the above mode identification strategy more difficult to apply although some of the results still remain promising. Even if it turned out to be too penalising for mode identification, it would still be possible to extract recurrent spacings, such as the rotation rate, and possibly the large frequency separation.
Of course, for multicolour mode identification to work well, highquality multicolour photometric observations of stars with numerous pulsation modes, such as δ Scuti stars, are needed. In this regard, the constellation of nanosatellites, BRITE (Kuschnig et al. 2009), is a promising source of this type of data, given that it observes in red and blue photometric bands. The PLATO mission, scheduled for launch in 2025, will contain a platform of 26 telescopes, two of which will include broadband filters and the remaining 24 will operate in white light (Catala et al. 2011; Rauer et al. 2014).
The critical rotation rate described here differs from that typically used by observers. Indeed, it is calculated using g_{eq}, the equatorial gravity of the current model, whereas most observers use the equatorial gravity of the model at breakup, which tends to be smaller owing to the increased equatorial radius. As such, the values used here convert to larger values if using the observers’ convention.
Acknowledgments
The authors thank the third referee for clear recommendations and suggestions which helped improve the manuscript. D.R.R. was financially supported through a postdoctoral fellowship from the Subside fédéral pour la recherche 2012, University of Liège, and was funded by the European Community’s Seventh Framework Programme (FP7/20072013) under grant agreement No. 312844 (SPACEINN), both of which are gratefully acknowledged. This work was granted access to the HPC resources of IDRIS under the allocation 201199992 made by GENCI (Grand Équipement National de Calcul Intensif).
References
 Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Baglin, A., Auvergne, M., Barge, P., et al. 2009, IAU Symp., 253, 71 [NASA ADS] [CrossRef] [Google Scholar]
 Balona, L. A., & Evers, E. A. 1999, MNRAS, 302, 349 [NASA ADS] [CrossRef] [Google Scholar]
 Balona, L. A., Lenz, P., Antoci, V., et al. 2012, MNRAS, 419, 3028 [NASA ADS] [CrossRef] [Google Scholar]
 Borucki, W., Koch, D., Batalha, N., et al. 2009, IAU Symp., 253, 289 [Google Scholar]
 Catala, C., Appourchaux, T., & Plato Mission Consortium. 2011, J. Phys. Conf. Ser., 271, 012084 [NASA ADS] [CrossRef] [Google Scholar]
 DaszyńskaDaszkiewicz, J., Dziembowski, W. A., Pamyatnykh, A. A., & Goupil, M.J. 2002, A&A, 392, 151 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 DaszynskaDaszkiewicz, J., Dziembowski, W. A., & Pamyatnykh, A. A. 2007, Acta Astron., 57, 11 [Google Scholar]
 Dupret, M. A. 2001, A&A, 366, 166 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dupret, M., De Ridder, J., Neuforge, C., Aerts, C., & Scuflaire, R. 2002, A&A, 385, 563 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dupret, M.A., De Ridder, J., De Cat, P., et al. 2003, A&A, 398, 677 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dupret, M.A., Grigahcène, A., Garrido, R., et al. 2005, MNRAS, 361, 476 [NASA ADS] [CrossRef] [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]
 García Hernández, A., Moya, A., Michel, E., et al. 2013, A&A, 559, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 GarcíaHernández, A., MartínRuiz, S., Monteiro, M. J. P. F. G., et al. 2015, ApJ, 811, L29 [NASA ADS] [CrossRef] [Google Scholar]
 Goupil, M.J., Dupret, M. A., Samadi, R., et al. 2005, JApA, 26, 249 [Google Scholar]
 Jackson, S., MacGregor, K. B., & Skumanich, A. 2005, ApJS, 156, 245 [NASA ADS] [CrossRef] [Google Scholar]
 Kuschnig, R., Weiss, W. W., Moffat, A., & Kudelka, O. 2009, in SolarStellar Dynamos as Revealed by Helio and Asteroseismology: GONG 2008/SOHO 21, eds. M. Dikpati, T. Arentoft, I. González Hernández, C. Lindsey, & F. Hill, ASP Conf. Ser., 416, 587 [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]
 MacGregor, K. B., Jackson, S., Skumanich, A., & Metcalfe, T. S. 2007, ApJ, 663, 560 [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]
 Paparó, M., Benkő, J. M., Hareter, M., & Guzik, J. A. 2016, ApJS, 224, 41 [NASA ADS] [CrossRef] [Google Scholar]
 Pasek, M., Lignières, F., Georgeot, B., & Reese, D. R. 2012, A&A, 546, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Poretti, E., Michel, E., Garrido, R., et al. 2009, A&A, 506, 85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rauer, H., Catala, C., Aerts, C., et al. 2014, Exper. Astron., 38, 249 [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. 2009a, A&A, 506, 189 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reese, D. R., Thompson, M. J., MacGregor, K. B., et al. 2009b, A&A, 506, 183 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reese, D. R., Lignières, F., Ballot, J., et al. 2013a, in Progress in Physics of the Sun and Stars: A New Era in Helio and Asteroseismology, eds. H. Shibahashi, & A. E. LynasGray, ASP Conf. Ser., 479, 545 [Google Scholar]
 Reese, D. R., Prat, V., Barban, C., van ’t VeerMenneret, C., & MacGregor, K. B. 2013b, A&A, 550, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Suárez, J. C., García Hernández, A., Moya, A., et al. 2014, A&A, 563, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Townsend, R. H. D. 2003, MNRAS, 343, 125 [NASA ADS] [CrossRef] [Google Scholar]
 Uytterhoeven, K., Moya, A., Grigahcène, A., et al. 2011, A&A, 534, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Overall success rate in finding similar modes via the multicolour mode identification scheme.
All Figures
Fig. 1 Meridional crosssections of two island modes with the same values but opposing parity, as determined from the parity of . The colours indicate the Eulerian pressure perturbation divided by the square root of the equilibrium density. 

In the text 
Fig. 2 Autocorrelation functions of acoustic spectra in a 2 M_{⊙} stellar model rotating at 0.5 Ω_{C}. The spherical radial orders of the modes, n, range from 2 to 9, therefore spanning 7 large frequency separations. The 4 columns correspond to 4 different inclinations, i = 0° being a poleon configuration. Each row corresponds to a different number of included modes. For example, in the top row, only the 10 most visible modes are included in the frequency spectra before calculating the autocorrelation function. The vertical dotted and dashed green lines give the rotation rate, Ω, and twice its value. The vertical dashed and continuous red lines indicate the large frequency separation, Δ, and half its value (see Eqs. (2) and (3)). 

In the text 
Fig. 3 Same as Fig. 2, except that the model is rotating at 0.7 Ω_{C}. 

In the text 
Fig. 4 Same as Fig. 2, except that the model is rotating at 0.7 Ω_{C} and the visibilities have been multiplied by random numbers between 1 and 100 (which therefore affects which modes are selected in the frequency spectra). 

In the text 
Fig. 5 Same as Fig. 2 (i.e. with a model rotating at 0.5Ω_{C}), but where the modes are normalised by the kinetic energy, multiplied by the appropriate power of (ω + mΩ). 

In the text 
Fig. 6 Fourier transform of frequency spectra in a 2 M_{⊙} model rotating at 0.5 Ω_{C}. Each row corresponds to a different number of selected modes, as indicated on the left, and each column to a different inclination, as stated above. The red, continuous and discontinuous, vertical lines correspond to the large frequency separation, Δ, and various fractions of this value: Δ/4, Δ/3, and Δ/2. The green, dotted and dashed, vertical lines correspond to Ω and 2Ω, respectively. 

In the text 
Fig. 7 Same as Fig. 6, but for a model rotating at 0.7 Ω_{C}. 

In the text 
Fig. 8 Same as Fig. 6, but where the mode normalisation includes random factors (top panel), or is based on the mean surface displacement (bottom panel). 

In the text 
Fig. 9 Amplitude ratios, frequency spectrum and ranks, associated autocorrelation function, and meridional crosssections of the reference mode (n,ℓ,m) = (9,4,1) and the 9 other most similar modes, for an inclination of i = 60°. In the first row, the left panel shows the amplitude ratios, the dashed line corresponding to the reference mode (we note that this line is hardly visible because it is mostly covered up by the solid lines from the other modes), the middle panel gives the frequencies, where a higher rank indicates a higher resemblance with the reference mode, 10 corresponding to the reference mode, and the right panel displays the autocorrelation function of the subspectrum of the selected frequencies. The second and third rows show the meridional crosssections of the selected modes. Beneath each mode, the values of (n,ℓ,m) are indicated when the mode corresponds to an island mode. Otherwise, only the azimuthal order, m, is provided. 

In the text 
Fig. 10 Quantities f_{T} and φ_{T} in two nonrotating models with masses 1.8 M_{⊙} and 2.0 M_{⊙}. The two models have virtually the same temperature (8143.9 K and 8137.3 K, respectively), but different log (g) values (4.1249 dex and 3.9593 dex, respectively). The different line styles correspond to different ℓ values, and the overlapping vertical dotted lines (which are hard to distinguish) indicate the position of the fundamental modes of both models. 

In the text 
Fig. 11 Quantity  f_{T}  /R^{⋆} as a function of colatitude and frequency, for the 1.8 M_{⊙} rotating model. 

In the text 
Fig. 12 Same as Fig. 9, but for the reference mode (n,ℓ,m) = (4,2,1) (and for both the adiabatic and pseudo nonadiabatic case). 

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.