# A diagnostic diagram for γ Doradus variables and slowly pulsating B-type stars

M. Takata, R.-M. Ouazzani, H. Saio, S. Christophe, J. Ballot, V. Antoci, S. J.

A. J. Salmon, K. Hijikawa

► To cite this version:

M. Takata, R.-M. Ouazzani, H. Saio, S. Christophe, J. Ballot, et al. A diagnostic diagram for γ Doradus variables and slowly pulsating B-type stars. Astronomy and Astrophysics - A&A, 2020, 635, pp.A106. 10.1051/0004-6361/201936297. hal-02553345

HAL Id: hal-02553345

https://hal.sorbonne-universite.fr/hal-02553345

Submitted on 24 Apr 2020

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

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

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

Distributed under a Creative Commons Attribution 4.0 International License
# A diagnostic diagram for γ Doradus variables and slowly pulsating B-type stars

M. Takata<fnref n="1" />,<fnref n="2" />, R.-M. Ouazzani<fnref n="1" />, H. Saio<fnref n="3" />, S. Christophe<fnref n="1" />, J. Ballot<fnref n="4" />, V. Antoci<fnref n="5" />, S. J. A. J. Salmon<fnref n="6" />, and K. Hijikawa<fnref n="2" />

Received 12 July 2019 / Accepted 3 February 2020

## ABSTRACT

Pulsating variables of γ Doradus type (γ Dor) and slowly pulsating B-type (SPB) stars are found on and near the main sequence with typical periods varying between one and several days, making them rather hard to detect from the ground. It is only with space missions such as CoRoT and *Kepler* that we became truly capable of determining their oscillation frequencies with enough precision to perform in-depth analyses. Here we present an efficient and easy-to-implement seismic tool, in which the frequency (*v*) and the square root of the frequency difference (√Δ*v*) are plotted against each other as the abscissa and the ordinate, respectively. This allows us to immediately (1) perform mode identification; (2) estimate the average rotation rate and the characteristic period of gravity modes; and (3) recognise certain physical effects, including buoyancy glitches and avoided crossings. This diagnostic tool can only be applied to prograde sectoral *g* modes. To validate the tool presented here, we used stellar models and also applied it to three γ Dor (KIC 12066947, KIC 5608334 and KIC 4846809) and one SPB star (KIC 3459297), all observed with *Kepler*. Furthermore, we show that the rotation rates determined using this new tool are consistent with the results of previous studies.

**Key words.** asteroseismology – stars: oscillations – stars: rotation – methods: data analysis

## 1. Introduction

Gamma Doradus (γ Dor) stars are a class of pulsating variable stars located on and near the main sequence with spectral types from late-A to early-F. The corresponding mass range is between 1.4 and 2.0 $M_⊙^1$. They show multi-periodic variability with amplitudes of ≤0.1 mag and periods of about one day. The average measured projected rotational velocity, $v \sin i$, is ~70–100 km s⁻<fnref n="1" />, which indicates that the surface rotation rate is of the same order as the pulsation periods. This means that rotation significantly affects the low-frequency oscillation modes in these stars, making them perfect objects for the study of how rotation affects the internal structure and evolution of these stars.

There has been substantial progress in the understanding of γ Dor stars, thanks to space missions such as CoRoT (Baglin et al. 2006) and *Kepler* (Borucki et al. 2010), which delivered long and continuous observations with unprecedented precision of a vast number of stars. Contrary to ground-based data, space observations do not suffer from daily aliasing issues, which pose a serious problem when studying pulsating stars with periods close to one day. It is only thanks to multiple years of data from *Kepler* that oscillation periods in γ Dor stars have finally been

resolved, leading to a true revolution in the study of these stars (e.g. Van Reeth et al. 2015). While *Kepler* data have revealed a significant number of oscillation modes in many cases, their exact identification is absolutely crucial for in-depth analyses.

Oscillation frequencies in rapid rotators are fundamentally different from those of slow rotators. For the latter, the rotation can be regarded as a small perturbation and so the average internal rotation rate can be estimated by measuring the rotational splittings in the frequency spectrum (Cowling & Newing 1949; Ledoux 1951). In the case of rapid rotators, the period versus period-spacing (*P*–Δ*P*) diagram is known to provide useful information about the mode identification and the internal rotation rate (e.g. Bouabid et al. 2013; Ouazzani et al. 2017). The asymptotic formula of the period spacing (Ballot et al. 2012) serves as a useful theoretical basis for extracting information from the frequency spectrum. Nevertheless, the estimation of the rotation rates of these stars is not a simple task. Van Reeth et al. (2016), for example, use a large grid of evolutionary models with an assumed range of rotation rates to find the best-fitting parameters for their sample of stars, including the rotation rates. On the other hand, Christophe et al. (2018) propose a model-independent method based on the discrete Fourier transform. The basic idea is that the pulsation periods will be evenly spaced even when the star is rapidly rotating if the period coordinate is stretched appropriately. Li et al. (2019) developed a different model-independent method in the same theoretical

<footnote n="1">¹ 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</footnote>

<footnote n="2">² Department of Astronomy, School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
e-mail: takata@astron.s.u-tokyo.ac.jp</footnote>

<footnote n="3">³ Astronomical Institute, Graduate School of Science, Tohoku University, Sendai 980-8578, Japan</footnote>

<footnote n="4">⁴ IRAP, Université de Toulouse, CNRS, UPS, CNES, 14 Avenue Édouard Belin, 31400 Toulouse, France</footnote>

<footnote n="5">⁵ Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark</footnote>

<footnote n="6">⁶ STAR Institute, Université de Liège, Allée du 6 Août 19, 4000 Liège, Belgium</footnote>

<footnote n="1">¹ For more details about the history of γ Dor stars we refer the reader to Saio et al. (2018a). Although the prototype was discovered in 1963 (Cousins & Warren 1963), the class of γ Dor pulsators was established only three decades later (Balona et al. 1994).</footnote>
framework as Van Reeth et al. (2016) and Christophe et al. (2018), which is described in detail in Sect. 2.1. Li et al. (2020) then applied the method to as many as 611 $\gamma$ Dor stars in the Kepler field. The central problem of all the methods mentioned above is the parameter sweep in a multidimensional space, which is generally time-consuming and computationally intense.

In this paper, we show that the internal rotation rate of $\gamma$ Dor and slowly pulsating B-type (SPB) stars can be estimated in a model-independent way by using a simple linear relation between the oscillation frequency and the square-root of the frequency difference. This method can only be applied to prograde sectoral $g$ modes, which is fortunate as observations (e.g. Saio et al. 2018a; Li et al. 2019) show that $\gamma$ Dor stars predominantly pulsate in these modes. We confine ourselves to single stars in this paper because tidal effects in binaries (or multiple systems) could complicate the picture significantly. Prograde sectoral $g$ modes consist of the waves that propagate in the same longitudinal direction as the rotation, with large amplitudes confined to the equatorial region. Their main restoring force is buoyancy, while the Coriolis force plays an essential role in the equatorial confinement. These modes, which are found in $\gamma$ Dor stars, are also excited in SPB stars (Waelkens 1991; Salmon et al. 2014). These stars are known as gravity-mode pulsators on the main-sequence with the mass between 3–9 $M_{\odot}$. The second most dominant component in the frequency spectrum of $\gamma$ Dor stars is that of $r$ modes (Papaloizou & Pringle 1978; Van Reeth et al. 2016; Saio et al. 2018b), which appear at lower frequencies than prograde sectoral $g$ modes. Rossby ($r$) modes are mostly composed of toroidal waves (Rossby waves) in each spherical (or equipotential) layer, which are restored by the Coriolis force. The (phase) propagation of the modes is opposite to the rotation in the longitudinal direction (in the co-rotating frame), while the wave motions in adjacent layers are weakly coupled with each other through the influence of the buoyancy. Observations suggest that these modes are detected not only in $\gamma$ Dor stars, but also in other various types of stars, including chemically peculiar stars and accreting white dwarfs (Saio 2018, 2019).

The structure of this paper is as follows: the method is described in detail in Sect. 2. It is then applied to (a) synthetic model frequencies and (b) observed frequencies of four stars in Sect. 3. Section 4 is devoted to discussions and conclusion.

## 2. Method

### 2.1. Theoretical framework

The present study focuses on eigenmodes of stellar oscillations, specifically, prograde sectoral $g$ modes that have been detected in $\gamma$ Dor and SPB stars. These modes have significantly longer periods than the dynamical time scale of the star and they are strongly affected by the Coriolis force. This is because their rotation periods are of the same order as the oscillation periods. These properties imply that the modes are composed of gravito-inertial waves, which are waves that are subject to the influence of both buoyancy and Coriolis forces.

The equations that govern the low-frequency eigenmodes of linear adiabatic oscillation in rapidly and uniformly<fnref n="2" /> rotating stars can be formulated to a good level of precision in the framework of the traditional approximation of rotation (Eckart 1960). We neglect (1) the perturbation to the gravitational potential (the Cowling approximation after Cowling 1941), (2) the

centrifugal deformation (the assumption of spherical symmetry of the equilibrium structure), and (3) the horizontal component of the rotation vector. The derived equations are separable in the spherical coordinates ($r, \theta, \phi$), with the direction of $\theta = 0$ aligned with the rotation axis. While the $\phi$ parts of the eigenfunctions are described by the sinusoidal function, $e^{im\phi}$, where $m$ is the integral index called the azimuthal order, the $\theta$ parts come from the solution of the Laplace tidal equation (the Hough function). On the other hand, the radial parts of the governing equations are of the same form as those of the spheroidal modes of non-rotating stars under the Cowling approximation. The only difference from the non-rotating case is that the term $\ell(\ell + 1)$, where $\ell$ is the spherical degree, is replaced with the eigenvalue $\lambda_{\ell,m}$ of the Laplace tidal equation. Assuming further that the radial wavelength of the oscillation is much shorter than the local scale height of the equilibrium structure, we may apply the asymptotic analysis (e.g. Vandakurov 1967) to derive the condition for the eigenfrequency of the low-frequency modes as

$$ \int_{r_1}^{r_2} k_r dr = (n + \alpha)\pi, \quad (1) $$

where the meanings of the symbols are given as follows: $k_r$ is the radial wave number of the oscillation; $r_1$ and $r_2$ represent the inner and outer boundaries of the region where the constituent waves are propagative; $n$ is the radial (integral) order of modes; $\alpha$ is the phase shift introduced at the inner and outer boundaries of the propagative region. The expression for the radial wave number $k_r$ is given by

$$ k_r = \frac{\sqrt{\lambda_{\ell,m}(s) N(r)}}{2\pi v_{\text{co}} r}. \quad (2) $$

Here $N(r)$ is the Brunt–Väisälä frequency, which is a function of the radius $r$, and $v_{\text{co}}$ is the cyclic frequency of the oscillation in the co-rotating frame. The eigenvalue $\lambda_{\ell,m}(s)$ depends on not only $\ell$ and $m$, but also the spin parameter $s$, which is defined by

$$ s = \frac{2v_{\text{rot}}}{v_{\text{co}}}. \quad (3) $$

The index $\ell$ (spherical degree) of $\lambda_{\ell,m}(s)$ means that $\lambda_{\ell,m}(s) \to \ell(\ell + 1)$ as $s \to 0$, while the angular part of the eigenfunctions cannot be described by any single spherical harmonic for $s \neq 0$. In other words, the $\ell$ of eigenmodes in rotating stars indicates the degree of spherical harmonics that specify the angular dependence only in the limit of no rotation. We note that $v_{\text{co}}$ is related to the oscillation frequency in the inertial frame, $v$, by

$$ v_{\text{co}} = v - m v_{\text{rot}}. \quad (4) $$

In this paper, we adopt the convention that positive (negative) values of $m$ correspond to prograde (retrograde) modes. Equation (1) implies that the oscillation period in the co-rotating frame is given by

$$ v_{\text{co}}^{-1} = \frac{(n + \alpha) P_0}{\sqrt{\lambda_{\ell,m}(s)}}, \quad (5) $$

in which the parameter $P_0$ is defined by

$$ P_0 = 2\pi^2 \left( \int_{r_1}^{r_2} \frac{N}{r} dr \right)^{-1}. \quad (6) $$

Equation (5) is in general an implicit expression because $s$ depends on $v_{\text{co}}$. In the absence of rotation, the series of high-order gravity modes with a given spherical degree $\ell$ are equidistantly spaced by a constant period spacing of $P_0 / \sqrt{\ell(\ell + 1)}$.

<footnote n="2">² While it is assumed that the rotation is uniform, we can account for the radial dependence of the rotation rate without changing the essential part of the analysis (see Sect. 4 and Appendix E).</footnote>
<figure><img src="image_2.png" /><figcaption>Fig. 1. Relative differences in the square roots of the eigenvalues of the Laplace tidal equation $\sqrt{\lambda_{m,m}}$ of the prograde sectoral g modes from the asymptotic values of m for m = 1–5 as functions of the spin parameter s. The vertical solid line indicates s = 1, the limit of the sub-inertial regime, s > 1, in which the effect of the Coriolis force is important, and where the approximation of $\sqrt{\lambda_{m,m}} \approx m$ is valid.</figcaption></figure>

provided that the composition profiles have no steep gradient, and that the dependence of $r_1$ and $r_2$ on the frequency is negligible. Although the spectrum of the low-frequency oscillations in rapidly-rotating stars does not show constant period spacings, the parameter $P_0$ still characterises the structure of the spectrum through Eq. (5). In this sense, $P_0$ can be regarded as a characteristic period of the gravity-mode oscillations<fnref n="3" />.

## 2.2. Prograde sectoral g modes and the $v - \sqrt{\Delta v}$ diagram

We derive the relation between the frequency and the square root of the frequency difference for prograde sectoral g modes with a given m. Supposing that there is a list of (observed or synthetic) frequencies sorted in the ascending order with the radial orders $n_k$ ($k = 1, 2, ...$) and the same m, we can take the difference of Eq. (5) multiplied by $\sqrt{\lambda_{l,m}}$ between the (k + 1)th and kth frequencies to obtain

$$ \Delta_k \left( \sqrt{\lambda_{l,m}} \nu_{co}^{-1} \right) = (-\Delta_k n) P_0, \quad (7) $$

in which $\Delta_k q$ generally means the difference in q between the (k+1)th and kth frequencies in the list, except for $\Delta_k n = -(n_{k+1} - n_k)(\ge 1)$. When we derive Eq. (7), we assume that $r_1, r_2$ and $\alpha$ are independent of the frequency. For the prograde sectoral g modes, which have $\ell = m > 0$, there is an approximate relation for large s,

$$ \sqrt{\lambda_{m,m}(s)} \approx m \quad (8) $$

(e.g. Bildsten et al. 1996; Townsend 2003). The accuracy of Eq. (8) is depicted in Fig. 1. We may further introduce the linearisation relation,

$$ \Delta_k (\nu_{co}^{-1}) \approx -\Delta_k \nu (\nu_{k+\frac{1}{2}} - m \nu_{\text{rot}})^{-2}, \quad (9) $$

under the assumption of $\Delta_k \nu \ll \nu_{k+\frac{1}{2}}$, where $\Delta_k \nu = \nu_{k+1} - \nu_k$ and $\nu_{k+\frac{1}{2}} = (\nu_k + \nu_{k+1})/2$, with $\nu_k$ being the kth frequency in the list.

With the help of Eqs. (8) and (9), we can reduce Eq. (7) to

$$ \sqrt{\frac{\Delta_k \nu}{m \Delta_k n}} \approx \sqrt{P_0} \left( \frac{\nu_{k+\frac{1}{2}}}{m} - \nu_{\text{rot}} \right). \quad (10) $$

In Eq. (10), $\Delta_k \nu$ (the difference between the two adjacent frequencies) and $\nu_{k+\frac{1}{2}}$ (the average frequency of the two) are both observables that are computed from the same pair of mode frequencies, while m, $\Delta_k n$, $P_0$ and $\nu_{\text{rot}}$ are not known in advance. With estimates of m and $\Delta_k n$, we can construct a diagram, based on only the two observables, in which the abscissa and the ordinate are given by $\nu_{k+\frac{1}{2}}/m$ and $\sqrt{\Delta_k \nu/(m\Delta_k n)}$, respectively. We call this the diagram of $\nu$ versus $\sqrt{\Delta_k \nu}$, or simply the $v-\sqrt{\Delta_k \nu}$ diagram. Although mode identification (i.e., the knowledge of $\Delta_k n$) is not always straightforward, this tool can be used to constrain $\Delta_k n$ of these prograde modes as described in Sect. 3. According to Eq. (10), each pair of the two adjacent frequencies should correspond to a point of a common straight line in the diagram. The slope and the abscissa intercept (the intersection point with the horizontal line of $\sqrt{\Delta_k \nu/(m\Delta_k n)} = 0$) of the linear fit are given by $\sqrt{P_0}$ and $\nu_{\text{rot}}$, respectively. Using this result inversely, we can estimate $\nu_{\text{rot}}$ and $P_0$ by plotting the modes in the diagram and fitting them with a straight line. The simplest procedure is to perform linear least-squares fitting under the assumption that the errors of each point are uncorrelated with each other. This assumption is acceptable for rough estimates. However, as shown in the Appendix A, $\Delta_k \nu$ and $\Delta_{k+1} \nu$ are correlated with each other since they both depend on $\nu_{k+1}$.

## 2.3. Revised method of estimating $\nu_{\text{rot}}$ and $P_0$

The formula given by Eq. (10) can be improved by discarding the two approximations that are given by Eqs. (8) and (9). The revised formula can directly be obtained from Eq. (7) as

$$ f_k (\nu_{\text{rot}}) \sqrt{\frac{\Delta_k \nu}{m \Delta_k n}} = \sqrt{P_0} \left( \frac{\nu_{k+\frac{1}{2}}}{m} - \nu_{\text{rot}} \right), \quad (11) $$

in which

$$ f_k (\nu_{\text{rot}}) = \left[ \frac{-1}{m \Delta_k \nu} \Delta_k \left( \frac{\sqrt{\lambda_{m,m}(s)}}{\nu_{\text{co}}} \right)^{\frac{1}{2}} \right] (\nu_{k+\frac{1}{2}} - m \nu_{\text{rot}}). \quad (12) $$

The correction factor $f_k$, which is a function of $\nu_{\text{rot}}$, is practically close to 1 because Eqs. (8) and (9) are generally good approximations. This property can be utilised to construct the following iterative procedure to obtain better estimates of $\nu_{\text{rot}}$ and $P_0$:

1. set $f_k$ to $f_k^{(1)} = 1$ and perform the least-squares fitting based on Eq. (11) to estimate $\nu_{\text{rot}} = \nu_{\text{rot}}^{(1)}$ and $P_0 = P_0^{(1)}$ (cf. Appendix A); this step is called iteration 1.

2. calculate $f_k$ by substituting $\nu_{\text{rot}} = \nu_{\text{rot}}^{(1)}$ into Eq. (12), and perform the least-squares fitting again to estimate $\nu_{\text{rot}} = \nu_{\text{rot}}^{(2)}$ and $P_0 = P_0^{(2)}$ based on Eq. (11); this step is called iteration 2.

3. repeat the same procedure as step 2 to obtain $\nu_{\text{rot}}^{(i+1)}$ and $P_0^{(i+1)}$ from $\nu_{\text{rot}}^{(i)}$ (iteration $i+1$) for $i = 2, 3, ...$ until both of $|\nu_{\text{rot}}^{(i)} - \nu_{\text{rot}}^{(i+1)}|$ and $|P_0^{(i)} - P_0^{(i+1)}|$ are small enough. We use the numerical values of $\lambda_{m,m}(s)$ to evaluate Eq. (12).

## 2.4. Error estimates

The errors in $\nu_{\text{rot}}$ and $P_0$ obtained by the procedure in Sect. 2.3 can originate from (1) the observational errors in the frequency

<footnote n="3">³ Christophe et al. (2018) and Ouazzani et al. (2019) introduce $P_0$ as the buoyancy radius.</footnote>
measurement and (2) the systematic errors in Eq. (11), which arise from differential rotation, steep chemical composition gradient (discussed in Sect. 3.1), and deviation from the traditional approximation of rotation in the asymptotic limit. Because the relative errors in the frequencies of $\gamma$ Dor and SPB stars in the *Kepler* data are typically of the order of $10^{-5}$ or below, we regard that the dominant uncertainties come from (2). In order to estimate the systematic errors accurately, we need a detailed comparison between the frequencies computed under the traditional approximation of rotation in the asymptotic limit and those determined by two-dimensional computations that take full account of the effects of the rapid rotation. We leave this expensive comparison to future work, however, for a first approximation, we determine uncertainties in $v_{\text{rot}}$ and $P_0$ from the residual sum of squares of the least-squares fitting. These estimates are based on the assumptions that (a) the uncertainties in the frequencies are not correlated and that (b) they follow a Gaussian distribution with the same standard deviation. For the sake of simplicity, we assume that there is no error correlation, but we do point out that the systematic errors in Eq. (11) are likely to be highly correlated. We set the default confidence level to 99% for the uncertainties in the estimates of $v_{\text{rot}}$ and $P_0$. More details are given in Appendix A.

## 3. Results

As described in Sect. 1, the $v-\sqrt{\Delta\nu}$ diagram method was developed based on *Kepler* observations. It indicates that most of the observed frequencies in $\gamma$ Dor stars can be identified as prograde sectoral *g* modes and *r* modes. Therefore, we first address these two types of modes in this section, while the other types are discussed in Sect. 3.3.2.

Since analysing $\gamma$ Dor and SPB stars requires very similar techniques, in this section we choose to concentrate on the former, while the latter is discussed in Appendix B.

### 3.1. Tests based on models

We construct the $v-\sqrt{\Delta\nu}$ diagrams and estimate $v_{\text{rot}}$ and $P_0$ for the following three models of $\gamma$ Dor stars (all having the same mass of $1.86 M_{\odot}$), which are described in detail in Christophe et al. (2018): (1) Model A (with smooth profiles of chemical composition); (2) Model B (with steep gradients of the chemical composition profiles just above the convective core); (3) Model A with differential rotation (DR) ($v_{\text{rot}} = 15$ µHz and 7 µHz in the convective core and the envelope, respectively). For each model, we used the same set of eigenfrequencies as Christophe et al. (2018), that is, those of prograde sectoral *g* modes with *m* = 1 and the *g* mode radial orders between −50 and −20. All mode frequencies were computed with the ACOR oscillation code (Ouazzani et al. 2012, 2015). For the sake of simplicity, we took only the spherically-symmetric part of the centrifugal deformation into account. In Figs. 2–4 we show the $v-\sqrt{\Delta\nu}$ diagrams for these models. The corresponding estimates of $v_{\text{rot}}$ and $P_0$ are summarised in Table 1.

The $v-\sqrt{\Delta\nu}$ diagrams in Figs. 2–4 clearly show that all the open circle symbols are well aligned and fit a straight line. This confirms that Eq. (10) is a good approximation in all of the models. In models (1) and (2), the spin parameter of the eigenmodes ranges between 1 and 3, for which Eq. (8) is accurate to only about 10–20% (see Fig. 1). The inaccuracy of Eq. (8) results in the deviation of $f_k$ in Eq. (12) from 1. This is the main reason for the differences between the dashed and solid lines in Figs. 2

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** *v* versus $\sqrt{\Delta\nu}$ diagram for model (1) (Model A). Open circles describe the modes with an azimuthal order *m* = 1 and $\Delta_k n = 1$, which implies consecutive radial orders. The blue dashed and black solid lines describe the output of iterations 1 and 5, respectively. The open square symbol shows the abscissa intercept with the solid line, which gives the converged estimate of the rotation frequency.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Same as Fig. 2, but for model (2) (Model B).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Same as Fig. 2, but for model (3). This model includes differential rotation. A detailed explanation of the symbols and the dotted line labelled $\times\sqrt{2}$ respectively, is given in the main text.</figcaption></figure>
<table><caption>**Table 1.** Estimates of $v_{\text{rot}}$ and $P_0$ for the three theoretical models from the synthetic eigenfrequencies of the prograde sectoral g modes with the azimuthal order of $m = 1$.</caption><thead><tr><th rowspan="2"></th><th colspan="2">(1) Model A</th><th colspan="2">(2) Model B</th><th colspan="2">(3) Model A with DR</th></tr><tr><th><i>v</i><sub>rot</sub></th><th><i>P</i><sub>0</sub> [10<sup>3</sup> s]</th><th><i>v</i><sub>rot</sub></th><th><i>P</i><sub>0</sub> [10<sup>3</sup> s]</th><th><i>v</i><sub>rot</sub></th><th><i>P</i><sub>0</sub> [10<sup>3</sup> s]</th></tr></thead><tbody><tr><td>Iteration 1</td><td>6.64 &mu;Hz<br>(0.574 d<sup>-1</sup>)</td><td>4.03</td><td>6.4 &mu;Hz<br>(0.55 d<sup>-1</sup>)</td><td>3.8</td><td>11.76 &mu;Hz<br>(1.016 d<sup>-1</sup>)</td><td>5.60</td></tr><tr><td>Converged value</td><td>6.944 &plusmn; 0.008 &mu;Hz<br>(0.6000 &plusmn; 0.0007 d<sup>-1</sup>)</td><td>4.49 &plusmn; 0.01</td><td>6.7 &plusmn; 0.1 &mu;Hz<br>(0.58 &plusmn; 0.01 d<sup>-1</sup>)</td><td>4.3 &plusmn; 0.1</td><td>11.80 &plusmn; 0.04 &mu;Hz<br>(1.020 &plusmn; 0.004 d<sup>-1</sup>)</td><td>5.70 &plusmn; 0.06</td></tr><tr><td>C18</td><td>6.95 &mu;Hz<br>(0.600 d<sup>-1</sup>)</td><td>4.495</td><td>6.67 &mu;Hz<br>(0.576 d<sup>-1</sup>)</td><td>4.233</td><td>11.74 &mu;Hz<br>(1.014 d<sup>-1</sup>)</td><td>5.592</td></tr><tr><td>True value</td><td>7 &mu;Hz<br>(0.605 d<sup>-1</sup>)</td><td>4.579</td><td>7 &mu;Hz<br>(0.605 d<sup>-1</sup>)</td><td>4.453</td><td>7-15 &mu;Hz<br>(0.605-1.296 d<sup>-1</sup>)</td><td>4.579</td></tr></tbody></table>

**Notes.** The process to estimate $v_{\text{rot}}$ and $P_0$ converges after 5, 5, and 4 iterations for models (1), (2), and (3), respectively, with the requirement that the relative differences in $v_{\text{rot}}$ and $P_0$ are both less than $10^{-4}$ between two successive iterations. The errors that are given for the converged values correspond to the 99% confidence intervals. The results of Christophe et al. (2018) based on the discrete Fourier transform are also given and described in the table as C18. The true value of $v_{\text{rot}}$ for Model A with differential rotation (DR) indicates the range of the rotation rate in the model.

and 3 or, equivalently, the differences between the outputs of iteration 1 and the converged values for these two models (see Table 1).

The converged values for the rotation rate based on the iterative linear least-squares fitting described in Sect. 2.3 are generally consistent with the results based on the discrete Fourier transform (DFT) map developed by Christophe et al. (2018), as shown in Table 1. Although the differences between the converged estimates and the true values are generally small, the true values of $v_{\text{rot}}$ and $P_0$ are found outside the given confidence intervals for all cases, which indicates that the unaccounted systematic errors of Eq. (11) can be significantly larger than the uncertainties derived here.

As an example of the unaccounted error sources, we may list the steep gradients in the chemical composition profiles just above the convective core in model (2). This cannot be treated by the standard asymptotic analysis because it assumes that the wavelength of the oscillation is much shorter than the scale height of the equilibrium structure. The corresponding variation in the Brunt-Väisälä frequency (buoyancy glitch) is known to cause a wavy structure in the period spacing (e.g. Miglio et al. 2008; Bouabid et al. 2013; Kurtz et al. 2014). Such an oscillatory behaviour in $\Delta P$ caused by a glitch can also be found in $\sqrt{\Delta v}$ as shown in Fig. 3.

The case of model (3) (see Fig. 4) can be used to check what happens if the method is applied to the case of differential rotation, although the uniform rotation is assumed in the (standard) traditional approximation of rotation. Broadly speaking, we find no critical difference from Fig. 2 because most of the points are aligned. This is, at first glance, consistent with the argument in Appendix E, which claims that the effect of weak differential rotation (in the radial direction) can be taken into account only by replacing the uniform rotation frequency by its average weighted by $N/r$ over the propagation region. However, the weak differential rotation cannot be assumed for model (3) since the variation between 7 and 15 µHz is large. This means that for a strongly differentially rotating star, we cannot simply interpret the derived $v_{\text{rot}}$, even though the value of 11.80 ± 0.04 µHz for model (3) is clearly within the range of values in the model.

It is remarkable that our method allows us to highlight the two avoided crossings in Fig. 4 at a frequency of ~15.5 µHz and ~17.5 µHz, respectively. It should be noted that methods using

the traditional approximation framework cannot account for this phenomenon. This is shown by the filled circles deviating from the solid line (base line). Since these points cannot be interpreted in the framework developed in Sect. 2, we simply discard them in the fitting to estimate $v_{\text{rot}}$ and $P_0$. Lee & Saio (1989) demonstrated that an avoided crossing can be caused by the horizontal component of the rotation vector, which is neglected in the traditional approximation of rotation. However, there may be other reasons too. In principle, it can occur due to any neglected terms that make the governing equations variable-inseparable into the radial and angular parts. In the case of differential rotation, those that include the gradient of the rotation frequency can also induce the avoided crossing.

Another important difference in model (3) is that some modes are missing, meaning that not all modes are of consecutive radial order. More specifically, the open diamonds aligned with the blue dotted line labelled "$\times\sqrt{2}$" in Fig. 4 correspond to the jumps in the radial order $n$ with $\Delta_k n = 2$. This demonstrates that this method allows us to not only determine $v_{\text{rot}}$ and $P_0$ but also $\Delta_k n$. When a point is found on the line that has the same abscissa intercept as the base line but a slope steeper than the base line by factor of $\sqrt{j}$ (with $j = 2, 3, ...$), we may interpret this as a missing $j-1$ radial order(s) with $\Delta_k n = j$. With the value of $\Delta_k n$, we can calculate the left-hand side of Eq. (10), which can be used in the fitting. This correction factor of $1/\sqrt{\Delta_k n}$ is indicated by each of the vertical cyan lines in Fig. 4. We only use the open circles in Fig. 4 to estimate $v_{\text{rot}}$ and $P_0$.

We apply the same method and construct the $v-\sqrt{\Delta v}$ diagram for $r$ modes of model (1) (see Fig. 5). In sharp contrast to Fig. 2, these modes are not distributed along a straight line with a positive slope. This very clear difference can be used to distinguish $r$ modes from prograde sectoral g modes. Thus, the diagram is also useful for mode identification in this respect.

## 3.2. Application to real stars

### 3.2.1. KIC 12066947

The first $\gamma$ Dor star analysed here is KIC 12066947. The observed frequencies and the Lomb-Scargle periodogram are provided in Table A.1 and Fig. 6 of Christophe et al. (2018), respectively. All peaks were extracted using the Aarhus
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Diagram of $v$ versus $\sqrt{\Delta v}$ for $r$ modes of model (1) with $k = -2$ and $m = -1$ (cf. Lee & Saio 1997).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 6.** $v$ versus $\sqrt{\Delta v}$ diagram for KIC 12066947. Here we used the frequencies between 27 µHz and 35 µHz which were already identified to be prograde sectoral $g$ modes with $m = 1$.</figcaption></figure>

University ECHO (Extraction of Coherent Oscillations) pipeline described in Antoci et al. (2019). Analyses show that the extracted frequencies are separated in two clusters, Cluster A between 20 and 24 µHz and Cluster B between 27 and 35 µHz.

The $v-\sqrt{\Delta v}$ diagram in Fig. 6 is constructed using 14 frequencies in Cluster B, which could be identified as $m = 1$ modes. It can be seen that except for one peak, all modes follow the three lines that are displayed. The solid line is the base line, while the two dotted lines labelled “×√3” and “×√6” indicate slopes that are larger than the base line by factors √3 and √6, respectively.

This is consistent with the interpretation that these frequencies are those of the prograde sectoral $g$ modes. The outlier at $(v_{k+\frac{1}{2}}, \sqrt{\Delta_k v}) = (28.4, 1.08)$ is discarded when we perform the least-squares fitting, which converges after three iterations. Table 2 provides the converged estimates of $v_{\text{rot}}$ and $P_0$ with their 99% and $1\sigma$ confidence intervals. These results are consistent with those of Van Reeth et al. (2016) and Christophe et al. (2018) within $1\sigma$, while those of Li et al. (2019) are within 99% confidence intervals. We note that the uncertainties of Li et al. (2019) are smaller than the others because they include not only prograde sectoral $g$ modes but also $r$ modes in their analysis.

<table><caption>**Table 2.** Estimates of $v_{\text{rot}}$ and $P_0$ for KIC 12066947.</caption><thead><tr><th></th><th><i>v</i><sub>rot</sub> [µHz]</th><th><i>v</i><sub>rot</sub> [d<sup>-1</sup>]</th><th><i>P</i><sub>0</sub> [10<sup>3</sup> s]</th></tr></thead><tbody><tr><td>This work</td><td></td><td></td><td></td></tr><tr><td>99% CI</td><td>24.9 ± 0.1</td><td>2.150 ± 0.009</td><td>4.1 ± 0.1</td></tr><tr><td>1σ CI</td><td>24.88 ± 0.04</td><td>2.150 ± 0.003</td><td>4.08 ± 0.04</td></tr><tr><td>R16</td><td>25.00<sup>+0.09</sup><sub>-0.10</sub></td><td>2.160 ± 0.008</td><td>4.171 ± 0.099</td></tr><tr><td>C18</td><td>24.95 ± 0.05</td><td>2.156 ± 0.004</td><td>4.181 ± 0.063</td></tr><tr><td>L19</td><td>24.99 ± 0.02</td><td>2.159 ± 0.002</td><td>4.17 ± 0.03</td></tr></tbody></table>

**Notes.** The results of the fitting to Eq. (11) are provided for two different confidence intervals (CI) of 99% and $1\sigma$.

**References.** R16: Van Reeth et al. (2016); C18: Christophe et al. (2018); L19: Li et al. (2019).

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** $v$ versus $\sqrt{\Delta v}$ for the observed frequencies of KIC 12066947. Here we used the peaks between 20 µHz and 24 µHz.</figcaption></figure>

The $v-\sqrt{\Delta v}$ diagram constructed using the frequencies in Cluster A is shown in Fig. 7. It is clear that we fail to observe the modes showing the expected pattern for prograde sectoral $g$ modes. We therefore conclude that these peaks are those of other modes, which is consistent with Van Reeth et al. (2016) who identify these as $r$ modes (cf. Fig. 5).

### 3.2.2. KIC 5608334

Saio et al. (2018a) study KIC 5608334, which shows remarkable frequency groups (fg) in the amplitude spectrum (cf. Fig. 3 in the paper). We adopted their mode identification of the four frequency groups, fg1–fg4, as the prograde sectoral $g$ modes with $m = 1–4$, respectively. From 192 frequencies that are determined with the ECHO pipeline of Antoci et al. (2019) and listed in Table A.1 of Saio et al. (2018a), we first excluded those with amplitudes smaller than 5 ppm (all of which belong to fg3 and fg4). We chose to do so because the small peaks show significant outliers in the $v-\sqrt{\Delta v}$ diagram, which may imply that many of them do not correspond to prograde sectoral $g$ modes or additional effects such as glitches are present. The physical origin of them should be examined in a separate study. In addition, if multiple frequencies are found in a very narrow frequency range for all of fg1–fg4 (cf. Fig. 4 of Saio et al. 2018a), we chose only the one with the highest amplitude and discarded the rest. The number of the remaining frequencies is 79 in total, and is equal to
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** v versus $\sqrt{\Delta v}$ diagram for KIC5608334 constructed from the peaks in fg1 (m = 1).</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 10.** Same as Fig. 8, but for fg3 (m = 3).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 9.** Same as Fig. 8, but for fg2 (m = 2).</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Same as Fig. 8, but for fg4 (m = 4).</figcaption></figure>

<table><caption>**Table 3.** Estimates of $v_{\text{rot}}$ and $P_0$ for KIC 5608334.</caption><thead><tr><th></th><th><i>v</i><sub>rot</sub> [μHz]</th><th><i>v</i><sub>rot</sub> [d<sup>-1</sup>]</th><th><i>P</i><sub>0</sub> [10<sup>3</sup> s]</th></tr></thead><tbody><tr><td>fg1 (m = 1)</td><td>26.0 ± 0.1</td><td>2.25 ± 0.01</td><td>4.4 ± 0.1</td></tr><tr><td>fg2 (m = 2)</td><td>25.9 ± 0.2</td><td>2.24 ± 0.02</td><td>4.2 ± 0.3</td></tr><tr><td>fg3 (m = 3)</td><td>25.8 ± 0.3</td><td>2.23 ± 0.02</td><td>4.1 ± 0.3</td></tr><tr><td>fg4 (m = 4)</td><td>25.9 ± 0.7</td><td>2.23 ± 0.06</td><td>4.2 ± 0.8</td></tr><tr><td>fg1-fg4</td><td>25.86 ± 0.09</td><td>2.234 ± 0.008</td><td>4.2 ± 0.1</td></tr></tbody></table>

24, 30, 14 and 11 for fg1–fg4, respectively. We then constructed the $v-\sqrt{\Delta v}$ diagrams for each frequency group, which are shown in Figs. 8–11. Table 3 provides the corresponding results of the least-squares fitting, where only the open circles in Figs. 8–11 were used.

There are a few possible physical reasons for the discarded points, which are indicated by filled circles in Figs. 8 and 9: they might originate from an avoided crossing (as in Fig. 4) or from a contamination of other types of frequencies (e.g. those of tesseral modes or combination frequencies) than those of the prograde sectoral g modes. These outliers are caused by the frequencies with small amplitudes: they disappear if we exclude those with amplitudes smaller than 100 and 20 ppm in fg1 and fg2, respectively. We observe a small deviation of some of the open circles from the base lines in Figs. 8–11. It is possible that they also originate from an avoided crossing or the influence of the nonlinear mode interaction, as discussed in Saio et al. (2018a). Although the estimates of $v_{\text{rot}}$ and $P_0$ for fg1–fg4 are completely consistent with each other within the errors, we cannot conclude uniform rotation only from these results because the theoretical model of Saio et al. (2018a) indicates that there is little difference in the propagation cavity of the modes for the different

groups, so the estimates for fg1–fg4 correspond to essentially the same physical quantities. Based on these considerations, the last row of Table 3 gives the estimates in the case where all the points with open circles in Figs. 8–11 are simultaneously taken into account in the fitting.

Saio et al. (2018a) estimate the internal rotation rate of the star as $2.20 d^{-1}$ (25.5 μHz) by comparing the observed period spacings with theoretical models by eye. Considering the roughness of their estimate, we regard the small deviation in our results outlined in Table 3 as non-significant. Saio et al. (2018a) also note an isolated peak at $2.2397 d^{-1}$ (25.922 μHz) in the frequency spectrum, which could be due to rotational modulation.
If this is indeed the case, the derived surface rotation rate agrees well with the average internal rotation rate in Table 3 (cf. Appendix E). Li et al. (2020) also analyse this star, but based only on peaks belonging to fg1. They obtain $v_{\text{rot}} = 26.0 \pm 0.1$ µHz ($2.25 \pm 0.01$ d$^{-1}$) and $P_0 = (4.4 \pm 0.2) \times 10^3$ s. These are completely consistent with our results in Table 3.

### 3.3. Applicability of the method

#### 3.3.1. Slow rotation

While the fundamental approximation of the present analysis, $\sqrt{\lambda} \approx m$ (cf. Eq. (8)), can be justified only when the spin parameter satisfies $s \gg 1$, we show in Appendix C that in principle the $v-\sqrt{\Delta v}$ diagram method may also work for $s \lesssim 1$. However, we should stress that it is essential not to use the approximation of $\sqrt{\lambda} \approx m$ in the case of slow rotation, but to utilise the realistic values shown in Fig. 1.

#### 3.3.2. Mode identification

As described above, the $v-\sqrt{\Delta v}$ diagram allows us to distinguish prograde sectoral *g* modes from *r* modes (cf. Sect. 3.1), however other types of modes may complicate the mode identification process. In this section, we discuss possible issues that may influence mode identification.

Retrograde modes with $\ell = 1$ and $m = -1$ behave differently from the prograde modes when plotted in both, the $v-\sqrt{\Delta v}$ as well as the $P-\Delta P$ diagrams (cf. Bouabid et al. 2013). More specifically, the sign of the slope is different: negative for retrograde and positive for prograde modes in the $v-\sqrt{\Delta v}$ diagram.

One relevant question may refer to how we discriminate prograde sectoral *g* modes with different azimuthal order *m*. We believe that this is not a real problem. As in the case of KIC 5608334 (Saio et al. 2018a), these modes with different *m* usually appear in different groups in the frequency spectrum, with smaller amplitudes for the group with higher frequencies. Since the effect of the geometrical cancellation is more significant for the modes with larger *m*, we may identify those groups as the modes with *m* = 1, 2, ... in the order of frequency. In particular, if there is only one group observed, and the $v-\sqrt{\Delta v}$ diagram confirms that it consists of the prograde sectoral *g* modes, it is quite likely that the group is that of *m* = 1 modes.

Some possible confusion can be caused by the zonal modes ($m=0$), which are claimed to be observed in rare cases (cf. Van Reeth et al. 2016; Li et al. 2020). In Fig. 12 we show the $\ell=1$ zonal modes of model (1) in the $v-\sqrt{\Delta v}$ diagram. We clearly observe very similar alignment of the points on a straight line to that in Fig. 2. If we erroneously identify these as prograde sectoral *g* modes with $m=1$, we obtain the rotation rate of $v_{\text{rot}} = 3.6$ µHz (0.31 d$^{-1}$), which is quite different from the true value of 7 µHz (0.60 d$^{-1}$). This example shows that it is dangerous to apply the method blindly. We therefore suggest for the mode identification to use the information about the amplitude in addition to the $v-\sqrt{\Delta v}$ diagram. This point is demonstrated in a particular case. Van Reeth et al. (2015) detect two series of period spacings in KIC 4846809 (see Fig. 15 of the paper), in which the shorter-period (higher-frequency) series have significantly larger amplitudes. Both series cannot be prograde sectoral *g* modes, because frequencies in the higher-frequency group (with higher azimuthal order, *m*) are more likely to have smaller amplitudes due to the more severe effect of the geometrical cancellation. We therefore argue that the higher-frequency (shorter-

<figure><img src="image_13.png" /><figcaption>Fig. 12. Diagram of $v$ versus $\sqrt{\Delta v}$ for zonal modes of model (1). These modes come from the $\ell=1$ modes in the limit of no rotation.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>Fig. 13. Diagram of $v$ versus $\sqrt{\Delta v}$ for KIC 4846809. The frequencies between 17 µHz and 27 µHz are determined with the ECHO pipeline of Antoci et al. (2019). The converged estimates for the parameters are $v_{\text{rot}} = 15.7 \pm 0.1$ µHz (1.36 ± 0.01 d$^{-1}$) and $P_0 = (5.8 \pm 0.3) \times 10^3$ s.</figcaption></figure>

period) series of KIC 4846809 corresponds to the prograde sectoral *g* modes with $m=1$, which can be confirmed in the $v-\sqrt{\Delta v}$ diagram shown in Fig. 13. This identification is consistent with Van Reeth et al. (2016). The lower frequencies (longer periods) series cannot be identified by our method completely, while Van Reeth et al. (2016) regard them as zonal modes.

Other types of modes, such as $\ell=2$ and $m=\pm 1$, can also introduce confusion similar to that caused by the zonal modes. They appear as a series with a negative slope in the $P-\Delta P$ diagram, which translates into one with a positive slope in the $v-\sqrt{\Delta v}$ diagram. However, they are unlikely to have much larger amplitude than prograde sectoral *g* modes when the stellar rotation is fairly high and are, apparently, identified in only few stars (Li et al. 2020). Such exceptional cases should carefully be investigated in a separate study.

Having considered all of the above, we may propose a general procedure for mode identification as follows:

1. The first step is to identify frequency groups, each comprising a cluster in the frequency spectrum of a given star. If this
is not possible, it may be the case that the star rotates slowly and that frequency splittings may be identified.

2. If one or more groups are present, mode identification is required. From a practical point of view, it seems to be a good initial guess that the group with the largest amplitude corresponds to the prograde sectoral *g* modes of *m* = 1. For these modes we expect to observe larger amplitudes because the geometrical cancellation effect is significantly smaller in comparison to other types of modes. We may find an exception for *r* modes but only if the star rotates rapidly, or more precisely if the spin parameter satisfies *s* ≥ 1 (Saio et al. 2018a).

3. If the modes of the first group do not align in the $v-\sqrt{\Delta v}$ diagram and show a positive slope consistent with prograde sectoral *g* modes, then the frequency group with the second-largest amplitudes should be examined.

4. This process should be repeated until prograde sectoral *g* modes with *m* = 1 are found (or all of the frequency groups have been checked). If there are two (or more) frequency groups with similar amplitudes at each step, we may assume that the one with the highest frequency originates from the prograde sectoral *g* modes with *m* = 1.

5. Once the identification of the first group is successful (i.e. *m* = 1), one may look for *m* = 2, 3, ... modes that should be located around the frequencies of the first group multiplied by *m*.

We regard this procedure as a general guideline that is applicable to the majority of the stars. The further details of the mode identification should be examined carefully based on a large sample of stars.

## 4. Discussion and conclusion

We propose the $v-\sqrt{\Delta v}$ diagram as an asteroseismic tool to estimate the average rotation rate ($v_{\text{rot}}$) and the characteristic period of the gravity-mode oscillations ($P_0$) of $\gamma$ Dor and SPB stars from the frequencies of the prograde sectoral *g* modes. Our method is endorsed by the fact that the prograde sectoral *g* modes are most commonly observed in (rapidly rotating) $\gamma$ Dor stars and SPB stars (e.g. Van Reeth et al. 2015; Pápics et al. 2017). Although the method is based on the same asymptotic formula of the low-frequency oscillation under the traditional approximation of rotation as in previous works (Van Reeth et al. 2016; Christophe et al. 2018; Li et al. 2019, 2020), the biggest advantage here is in its efficiency and simplicity. It is easy to get rough estimates of $v_{\text{rot}}$ and $P_0$ by plotting the diagram and applying a linear fit, which the majority of the modes follow ideally. The estimates of $v_{\text{rot}}$ and $P_0$ can be derived from the abscissa intercept and (the square of) the slope of the line, respectively. If no clear alignment is visible, it is unlikely that the observed frequencies correspond to prograde sectoral *g* modes. The diagram thus also allows us to perform mode identification. In addition, this tool can also be used to discuss the physics of oscillation modes including the buoyancy glitch (Fig. 3) and the avoided crossing (Fig. 4). The first rough estimates of $v_{\text{rot}}$ and $P_0$ can be improved by just a linear least-squares fitting with a few iterations. It should be stressed that the method is free from any time-intensive parameter search in two- (or more) dimensional space. There is no need to constrain the ranges of any parameters in advance or give any initial guesses. These properties are important, especially when analysing a large sample of stars.

The method proposed in this paper depends on the eigenvalue of the Laplace tidal equation of the prograde sectoral *g* modes being approximately equal to the square of the azimuthal

order, $m^2$, if the spin parameter *s* is much larger than 1 (cf. Fig. 1). This implies that the mode frequencies are nearly independent of $v_{\text{rot}}$ in the co-rotating frame (cf. Eq. (5)). This is because the dispersion relation of the constituent waves (internal equatorial Kelvin waves) is independent of $v_{\text{rot}}$ when *s* ≫ 1 and also because they propagate only along the equator in the same limit. More details about the physical properties of the modes, which may be called Kelvin *g* modes, are discussed in Appendix D. While the mode frequencies are insensitive to $v_{\text{rot}}$ in the co-rotating frame, we may infer $v_{\text{rot}}$ from those in the inertial frame through the advection effect induced by rotation (Doppler shift), which is represented by the second term on the right-hand side of Eq. (4). The $v-\sqrt{\Delta v}$ diagram allows us to represent this effect graphically.

While it is standard to assume that the rotation rate is constant in the traditional approximation of rotation, this is not necessary in the asymptotic limit. In the case where the rotation rate weakly depends on the radius, the obtained estimate of $v_{\text{rot}}$ should be interpreted as the average of the internal rotation rate weighted by $N/r$ (the Brunt-Väisälä frequency divided by the radius), as is proposed by Ouazzani et al. (2019) and justified in Appendix E.

We confirm that the estimate of $v_{\text{rot}}$ is highly correlated with that of $P_0$ (cf. Van Reeth et al. 2016; Christophe et al. 2018): the correlation coefficient *a*, which is defined in Eq. (A.29), satisfies *a* ≥ 0.98 in most cases presented in this paper. Such high values can be understood based on the $v-\sqrt{\Delta v}$ diagram. Supposing we have obtained the best-fit line on the diagram, if we slightly increase (or decrease) the abscissa intercept ($v_{\text{rot}}$) from the best-fit value and perform a linear fit under this constraint, the slope of the fitted line ($\sqrt{P_0}$) is necessarily steeper (or shallower) than the best-fit line because all the points in the diagram are located in the upper right side of the abscissa intercept (the domain of $v/m > v_{\text{rot}}$ and $\sqrt{\Delta v} > 0$). The reason for the high correlation coefficient is that we take only prograde sectoral *g* modes into account. We can, in principle, reduce the correlation by including other types of modes, if available (cf. Li et al. 2019). In this context, a natural extension of the present method would be to include the other types of modes by deriving similar formulae to Eq. (11).

We adopt conservative estimates of the formal uncertainties in the present analysis, corresponding to 99% confidence intervals calculated based on the goodness of the fit. Nevertheless, Table 1 shows evidence that these values can significantly underestimate the true uncertainties. This is probably a common problem in any method based on the asymptotic frequency formula under the traditional approximation of rotation and it should be examined in detail in a future work. In order to evaluate the total uncertainty, it is necessary to improve Eq. (11), which is the fundamental relation of the present work. This can be done by taking into account additional physical effects, including centrifugal deformation, correction to the traditional approximation of rotation, higher-order terms in the asymptotic frequency formula, and the possible steep variation of the chemical composition profiles just outside the convective core.

**Acknowledgements.** We thank NASA and the *Kepler* team for their revolutionary data. M. T. thanks B. Mosser and his colleagues for their support and hospitality at the Paris Observatory, where this project started. The anonymous referee is acknowledged for valuable remarks, which helped improving this paper. M. T. is partially supported by JSPS KAKENHI Grant Number JP18K03695. Funding for the Stellar Astrophysics Centre is provided by The Danish National Research Foundation (Grant agreement no.: DNRF106). S. J. A. J. S. is funded by the Wallonia-Brussels Federation ARC grant for Concerted Research Actions.
References

Antoci, V., Cunha, M. S., Bowman, D. M., et al. 2019, MNRAS, 490, 4040

Baglin, A., Auvergne, M., Boisnard, L., et al. 2006, in 36th COSPAR Scientific Assembly, COSPAR Meeting, 36

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

Balona, L. A., Krisciunas, K., & Cousins, A. W. J. 1994, MNRAS, 270, 905

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

Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977

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

Chandrasekhar, S., & Lebovitz, N. R. 1964, ApJ, 140, 1517

Christophe, S., Ballot, J., Ouazzani, R.-M., Antoci, V., & Salmon, S. J. A. J. 2018, A&A, 618, A47

Cousins, A. W. J., & Warren, P. R. 1963, Mon. Notes Astron. Soc. S. Afr., 22, 65

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

Cowling, T. G., & Newing, R. A. 1949, ApJ, 109, 149

Cox, J. P. 1980, Theory of Stellar Pulsation, Princeton Series in Astrophysics (Princeton: Princeton University Press)

Dziembowski, W. A., Daszyńska-Daszkiewicz, J., & Pamyatnykh, A. A. 2007, MNRAS, 374, 248

Eadie, W. T., Drijard, D., James, F. E., Roos, M., & Sadoulet, B. 1971, Statistical Methods in Experimental Physics (Amsterdam: North Holland Publishing Co.)

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

Gill, A. 1982, Atmosphere-ocean Dynamics (New York: Academic Press)

Gough, D. O. 1993, in Astrophysical Fluid Dynamics – Les Houches 1987, eds. J. P. Zahn, & J. Zinn-Justin (Amsterdam: North Holland Publishing Co.), 399

Gough, D. O., & Sekii, T. 2002, MNRAS, 335, 170

Goupil, M. J., Mosser, B., Marques, J. P., et al. 2013, A&A, 549, A75

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

Ledoux, P. 1951, ApJ, 114, 373

Lee, U., & Saio, H. 1989, MNRAS, 237, 875

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

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

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

Longuet-Higgins, M. S. 1968, Philos. Trans. R. Soc. London A, 262, 511

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

Miglio, A., Montalbán, J., Noels, A., & Eggenberger, P. 2008, MNRAS, 386, 1487

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

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

Ouazzani, R.-M., Salmon, S. J. A. J., Antoci, V., et al. 2017, MNRAS, 465, 2294

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

Papaloizou, J., & Pringle, J. E. 1978, MNRAS, 182, 423

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

Saio, H. 2018, ArXiv e-prints [arXiv:1812.01253]

Saio, H. 2019, MNRAS, 487, 2177

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

Saio, H., Kurtz, D. W., Murphy, S. J., Antoci, V. L., & Lee, U. 2018b, MNRAS, 474, 2774

Salmon, S. J. A. J., Montalbán, J., Reese, D. R., Dupret, M.-A., & Eggenberger, P. 2014, A&A, 569, A18

Scuflaire, R., Théado, S., Montalbán, J., et al. 2008, Ap&SS, 316, 83

Thomson, W. 1863, Philos. Trans. R. Soc. London, 153, 583

Thomson, W. 1880, Proc. R. Soc. Edinb., 10, 92

Townsend, R. H. D. 2003, MNRAS, 340, 1020

Van Reeth, T., Tkachenko, A., Aerts, C., et al. 2015, ApJS, 218, 27

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

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

Vandakurov, Y. V. 1967, AZh, 44, 786

Waelkens, C. 1991, A&A, 246, 453
# Appendix A: Linear least-squares fitting with correlated errors

## A.1. Simplifying assumptions

Here we describe the method of estimating $v_{\text{rot}}$ and $P_0$ with their uncertainties based on Eq. (11) in detail. We paid special attention to the treatment of the error correlation. The following three assumptions are made for the sake of simplicity: (1) the errors $\{\epsilon_k\}$ of the mode frequencies $\{v_k\}$ follow a Gaussian distribution with the same (but unknown) standard deviation, $\sigma$, and those of different modes are not correlated with each other; supposing that the number of the mode frequencies is $K+1$, and that $\epsilon$ indicates the $(K+1)$-element vector whose $k$th element is equal to $\epsilon_k$, this assumption implies

$$E(\epsilon\epsilon^T) = \sigma^2 I_{K+1}, \qquad (A.1)$$

in which $v^T$, $E(a)$ and $I_k$ generally mean the transpose of $v$, the expected value of $a$ and the identity matrix of size $k$, respectively; (2) because the relation $v_{k+1/2} \gg |\Delta v_k|$ holds for any data sets we use in this paper, the frequency errors affect only the vertical position of the modes (symbols) in the $v-\sqrt{\Delta v}$ diagram; (3) since $f_k$, which is defined by Eq. (12), is generally close to 1, the frequency errors have negligible impact on $f_k$.

## A.2. Formulation in the vector form

We assume for the moment that the $K+1$ mode frequencies have consecutive radial orders and the same azimuthal order $m$, so that there are $K$ relations in the form of Eq. (11). This assumption is not essential (cf. the second paragraph of Appendix A.3). Neglecting the second- and higher-order terms in $\{\epsilon_k\}$, we may rewrite those relations as

$$y = X\theta + e, \qquad (A.2)$$

in which $e$ represents the error term caused by $\{\epsilon_k\}$. The definitions of the symbols in Eq. (A.2) are given as follows: the $k$th element of the $K$-element vector $y$ is given by

$$y_k = f_k(v_{\text{rot}}) \sqrt{\frac{\Delta_k v}{m \Delta_k n}} \quad (k = 1, 2, \dots, K); \qquad (A.3)$$

the $K \times 2$ matrix $X$ is defined by

$$X = \begin{pmatrix} 1 & v_{1+1/2}/m \\ \vdots & \vdots \\ 1 & v_{K+1/2}/m \end{pmatrix}; \qquad (A.4)$$

the 2-element vector $\theta$ is provided by

$$\theta = \left( \frac{-v_{\text{rot}} \sqrt{P_0}{\sqrt{P_0}}}{\sqrt{P_0}} \right); \qquad (A.5)$$

the $k$th element of the $K$-element vector $e$ is given by

$$e_k = h_k (\epsilon_k - \epsilon_{k+1}) \quad (k = 1, 2, \dots, K), \qquad (A.6)$$

in which

$$h_k = -\frac{y_k}{2\Delta_k v} \quad (k = 1, 2, \dots, K). \qquad (A.7)$$

We note that $e$ can be expressed by

$$e = HQ\epsilon, \qquad (A.8)$$

where $K \times K$ diagonal matrix $H$ and $K \times (K+1)$ matrix $Q$ are defined by

$$H = \begin{pmatrix} h_1 & & \\ & \ddots & \\ & & h_K \end{pmatrix} \qquad (A.9)$$

and

$$Q = \begin{pmatrix} 1 & -1 & & \\ & \ddots & \ddots & \\ & & 1 & -1 \end{pmatrix}, \qquad (A.10)$$

respectively.

## A.3. Error correlation

The variance-covariance matrix of the error vector $e$ can be calculated from Eqs. (A.1) and (A.8) as

$$C = E(ee^T) = \sigma^2 C_0, \qquad (A.11)$$

in which $C_0 = HQQ^T H$. We separate $C_0$ from $C$ because $C_0$ can be calculated from the mode frequencies, whereas $\sigma$ is unknown. The fact that $C$ is a tridiagonal matrix means that any two successive elements of $e$ are correlated with each other. The inverse of $C$ can be expressed by

$$C^{-1} = \sigma^{-2} C_0^{-1}, \qquad (A.12)$$

in which the $(i, j)$ element of $C_0^{-1}$ is given by

$$(C_0^{-1})_{i,j} = \frac{ij}{h_i h_j} \left[ \frac{1}{\max(i, j)} - \frac{1}{K+1} \right] \quad (i, j = 1, 2, \dots, K). \qquad (A.13)$$

It has so far been assumed that the $K$ relations in the form of Eq. (11) are all included in the analysis. We may have neglected some of them because they produced outliers in the $v-\sqrt{\Delta v}$ diagram. The possible reasons include avoided crossings, which cannot be described by the traditional approximation of rotation, and contamination of other types of modes than the prograde sectoral $g$ modes. In this case, the expression of $C$ (and its inverse) should be modified accordingly. For example, if we discard the $k$th element of Eq. (A.2) for $1 < k < K$, the remaining relations are separated into two groups, one with the $1, \dots, (k-1)$th elements and the other with the $(k+1), \dots, K$th elements. Since $e_{k-1}$ is not correlated with $e_{k+1}$ (cf. Eq. (A.6)), there is no error correlation between the two groups. The modified expression of $C_0^{-1}$ is composed of two diagonal blocks both in the form of Eq. (A.13) with $K$ replaced by the block size, which is equal to $k-1$ and $K-k$ for the first and second blocks, respectively. This treatment can easily be generalised as follows: supposing that the relations are separated into $J$ different groups with no error correlation, $C_0^{-1}$ should be composed of $J$ diagonal blocks in the form of Eq. (A.13). Based on this, it is possible to include the mode frequencies with different azimuthal order $m$ in the analysis because the errors of any two relations with different $m$ in the form of Eq. (11) are not correlated with each other. The different values of $m$ correspond to different diagonal blocks of $C_0^{-1}$, and each can generally consist of multiple diagonal blocks in the form of Eq. (A.13).
### A.4. Linear least-squares fitting

The problem of the linear least-squares fitting in the presence of error correlation is discussed in the literature (e.g. Eadie et al. 1971; Gough & Sekii 2002). The key idea is to multiply Eq. (A.2) by an appropriate $K \times K$ regular matrix from the left side so that any two elements of the resultant error vector are not correlated with each other. Once the error correlation disappears, we can adopt the standard procedure of the least-squares fitting. The procedure can be rewritten in terms of the original variables. In fact, we should find such $\theta$ that minimises

$$ (\mathbf{y} - \mathbf{X}\theta)^T C_0^{-1} (\mathbf{y} - \mathbf{X}\theta) \quad (\text{A.14}) $$

for given $\mathbf{y}$, $\mathbf{X}$ and $C_0^{-1}$. The solution is given by

$$ \hat{\theta} = \mathbf{A}\mathbf{X}^T C_0^{-1} \mathbf{y}, \quad (\text{A.15}) $$

where $\mathbf{A}$ is a $2 \times 2$ matrix defined by

$$ \mathbf{A} = (\mathbf{X}^T C_0^{-1} \mathbf{X})^{-1}. \quad (\text{A.16}) $$

Vector $\hat{\theta}$ is the best linear unbiased estimator of $\theta$, whose variance-covariance matrix is provided by

$$ \operatorname{cov}(\hat{\theta}) = \sigma^2 \mathbf{A}. \quad (\text{A.17}) $$

On the other hand, the residual sum of squares divided by $K-2$ yields the minimum variance unbiased estimator of $\sigma^2$,

$$ \hat{\sigma}^2 = \frac{1}{K-2} (\mathbf{y} - \mathbf{X}\hat{\theta})^T C_0^{-1} (\mathbf{y} - \mathbf{X}\hat{\theta}). \quad (\text{A.18}) $$

The corresponding estimator of $\nu_{\text{rot}}$ is naturally found to be

$$ \hat{\nu}_{\text{rot}} = -\frac{\hat{\theta}_1}{\hat{\theta}_2}. \quad (\text{A.19}) $$

Neglecting the second- and higher-order terms in $\{\epsilon_k\}$, we can show that $\hat{\nu}_{\text{rot}}$ is unbiased. Its variance is accordingly calculated as

$$ \operatorname{var}(\hat{\nu}_{\text{rot}}) = \sigma^2 \hat{\nu}_{\text{rot}}^2 \left( \frac{A_{11}}{\hat{\theta}_1^2} - \frac{2A_{12}}{\hat{\theta}_1 \hat{\theta}_2} + \frac{A_{22}}{\hat{\theta}_2^2} \right). \quad (\text{A.20}) $$

The uncertainty in $\hat{\nu}_{\text{rot}}$ can be estimated based on the fact that

$$ t = \frac{\hat{\nu}_{\text{rot}} - \nu_{\text{rot}}}{s(\hat{\nu}_{\text{rot}})} \quad (\text{A.21}) $$

follows Student's $t$-distribution with $K-2$ degrees of freedom. Here, $s(\hat{\nu}_{\text{rot}})$ stands for the standard error of $\hat{\nu}_{\text{rot}}$ defined by

$$ s(\hat{\nu}_{\text{rot}}) = \hat{\sigma} \hat{\nu}_{\text{rot}} \left( \frac{A_{11}}{\hat{\theta}_1^2} - \frac{2A_{12}}{\hat{\theta}_1 \hat{\theta}_2} + \frac{A_{22}}{\hat{\theta}_2^2} \right)^{\frac{1}{2}}. \quad (\text{A.22}) $$

The $q\%$ confidence interval of $\nu_{\text{rot}}$ is provided by

$$ \hat{\nu}_{\text{rot}} - t_0 s(\hat{\nu}_{\text{rot}}) < \nu_{\text{rot}} < \hat{\nu}_{\text{rot}} + t_0 s(\hat{\nu}_{\text{rot}}), \quad (\text{A.23}) $$

in which $t_0$ is such a value that the probability of $|t| \le t_0$ is equal to $q/100$. Alternatively, $t_0$ can be set (for $K > 4$) to

$$ t_0 = \sqrt{\frac{K-2}{K-4}}, \quad (\text{A.24}) $$

which corresponds to one standard deviation.

The way of estimating $P_0$ is analogous. All we need to do is replace the expressions for the estimator, its variance and the standard error with

$$ \hat{P}_0 = \hat{\theta}_2^2, \quad (\text{A.25}) $$

$$ \operatorname{var}(\hat{P}_0) = 4\sigma^2 \hat{P}_0 A_{22} \quad (\text{A.26}) $$

and

$$ s(\hat{P}_0) = 2\hat{\sigma} \sqrt{\hat{P}_0 A_{22}}, \quad (\text{A.27}) $$

respectively.

Finally, the covariance and the correlation coefficient between $\hat{\nu}_{\text{rot}}$ and $\hat{P}_0$ can be computed as

$$ \operatorname{cov}(\hat{\nu}_{\text{rot}}, \hat{P}_0) = 2\sigma^2 \hat{\nu}_{\text{rot}} \hat{P}_0 \left( \frac{A_{12}}{\hat{\theta}_1 \hat{\theta}_2} - \frac{A_{22}}{\hat{\theta}_2^2} \right) \quad (\text{A.28}) $$

and

$$ a = \frac{\operatorname{cov}(\hat{\nu}_{\text{rot}}, \hat{P}_0)}{\sqrt{\operatorname{var}(\hat{\nu}_{\text{rot}}) \operatorname{var}(\hat{P}_0)}}, \quad (\text{A.29}) $$

respectively.

## Appendix B: Application to SPB stars

In this section, we apply the $\nu - \sqrt{\Delta\nu}$ diagram method to SPB stars. First we validate the method based on a typical evolutionary model, and then apply it to KIC 3459297.

### B.1. Validation based on a model

As in the case of Model A for $\gamma$ Dor stars (cf. Sect. 3.1), we use the stellar evolution code CLES (Scuflaire et al. 2008) to construct a $4 M_{\odot}$ model of an SPB star. The model parameters are summarised in Table B.1. The details on how this model is calculated are the same as in Model A (cf. Christophe et al. 2018), however we note that turbulent diffusion is included to approximate the effect of rotationally-induced mixing, which smoothes the chemical composition profiles outside the convective core. The eigenfrequencies of the prograde sectoral $g$ modes with $m = 1$ are computed using the ACOR oscillation code (Ouazzani et al. 2012, 2015) with a uniform rotation frequency of $7 \mu\text{Hz}$ ($0.60 \text{ d}^{-1}$). The range of radial orders is between -50 and -18. The results are shown in Fig. B.1.

We confirm that also in the case of an SPB star all modes follow a nearly perfect linear trend, which is similar to model (1) shown in Fig. 2. Our estimates for the rotational frequency $\nu_{\text{rot}}$ and the characteristic period $P_0$ are provided in Table B.2. Strictly speaking, these estimates are slightly smaller than the true values by 0.9% and 4% for $\nu_{\text{rot}}$ and $P_0$, respectively, which is similar to the case of Model A listed in Table 1. We find that this is an intrinsic limitation of the method, which relies on the asymptotic formula under the traditional approximation of rotation. Nevertheless, since the two parameters are recovered reasonably well, we consider that the method can be applied to not only $\gamma$ Dor stars, but SPB stars as well.
<figure><img src="image_15.png" /><figcaption>Fig. B.1. Diagram of $\nu$ versus $\sqrt{\Delta\nu}$ for a typical model of SPB stars.</figcaption></figure>

<table><caption>Table B.1. Parameters of a typical evolutionary model of an SPB star.</caption><thead><tr><th>Mass</th><th>Radius</th><th>$T_{\text{eff}}$</th><th>Luminosity</th><th>Age</th></tr></thead><tbody><tr><td>4.00 $M_{\odot}$</td><td>3.19 $R_{\odot}$</td><td>1.38 × 10<sup>4</sup> K</td><td>331 $L_{\odot}$</td><td>101 Myr</td></tr><tr><th>$X_c$</th><th>$Y_{\text{ini}}$</th><th>$Z_{\text{ini}}$</th><th>$\alpha_{\text{MLT}}$</th><th>$D_t$</th></tr><tr><td>0.360</td><td>0.270</td><td>0.0142</td><td>1.8</td><td>700 cm<sup>2</sup> s<sup>-1</sup></td></tr></tbody></table>

Notes. The columns describe the following parameters: $T_{\text{eff}}$ is the effective temperature; $X_c$ is the hydrogen mass fraction at the centre; $Y_{\text{ini}}$ and $Z_{\text{ini}}$ are the initial mass fractions of helium and heavy elements, respectively; $\alpha_{\text{MLT}}$ is the mixing length parameter of convection; $D_t$ is the coefficient of turbulent diffusion.

<table><caption>Table B.2. Estimates of $v_{\text{rot}}$ and $P_0$ for the evolutionary model of an SPB star, the parameters of which are summarised in Table B.1.</caption><thead><tr><th></th><th>$v_{\text{rot}}$ [µHz]</th><th>$v_{\text{rot}}$ [d<sup>-1</sup>]</th><th>$P_0$ [10<sup>3</sup> s]</th></tr></thead><tbody><tr><td>Iteration 1</td><td>6.88</td><td>0.595</td><td>7.55</td></tr><tr><td>Iteration 4</td><td>6.94 ± 0.01</td><td>0.600 ± 0.001</td><td>7.78 ± 0.04</td></tr><tr><td>True value</td><td>7</td><td>0.605</td><td>8.08</td></tr></tbody></table>

## B.2. KIC 3459297

We use the frequencies of KIC 3459297 that are listed in Table A.2 of Christophe et al. (2018). While those frequencies distribute in three separate groups (cf. Fig. 8 of Christophe et al. 2018), we ignore the lowest-frequency group below 1.5 µHz because it probably originates from instrumental effects or the artefacts introduced during data analysis. The group between 9 and 14 µHz and that between 23 and 31 µHz are referred to as fg1 and fg2, with 24 and 8 frequencies, respectively. We identify the frequencies in fg1 and fg2 as those of prograde sectoral g modes with $m = 1$ and $m = 2$, respectively.

The $\nu-\sqrt{\Delta\nu}$ diagram is shown in Fig. B.2. As in the case of KIC 5608334, which is discussed in Sect. 3.2.2, we discard four frequencies (10.219, 10.354, 11.233 and 24.903 µHz) from the list because there is another frequency very close to each of them with higher amplitude. We also exclude the two red points in Fig. B.2 in the $v_{\text{rot}}$ and $P_0$ fitting because they do not follow the same linear trend as the other modes. These peaks might indicate the presence of a buoyancy glitch (cf. Fig. 3) or an avoided crossing (cf. Fig. 4).

<figure><img src="image_16.png" /><figcaption>Fig. B.2. Diagram of $\nu$ versus $\sqrt{\Delta\nu}$ for KIC 3459297. The red and blue circles correspond to the frequencies in the first group (with $m = 1$) and the second group (with $m = 2$), respectively. The dashed and solid lines are the fitting results of all the modes in iterations 1 and 5, respectively. The two red filled circles are excluded from the fit.</figcaption></figure>

<table><caption>Table B.3. Estimates of $v_{\text{rot}}$ and $P_0$ for KIC 3459297.</caption><thead><tr><th></th><th>$v_{\text{rot}}$ [µHz]</th><th>$v_{\text{rot}}$ [d<sup>-1</sup>]</th><th>$P_0$ [10<sup>3</sup> s]</th></tr></thead><tbody><tr><td>fg1 (m = 1)</td><td>6.80 ± 0.09</td><td>0.588 ± 0.008</td><td>6.9 ± 0.3</td></tr><tr><td>fg2 (m = 2)</td><td>7 ± 1</td><td>0.6 ± 0.1</td><td>9 ± 4</td></tr><tr><td>fg1 and fg2</td><td></td><td></td><td></td></tr><tr><td>99% CI</td><td>7.0 ± 0.1</td><td>0.61 ± 0.01</td><td>7.7 ± 0.3</td></tr><tr><td>1σ CI</td><td>7.00 ± 0.05</td><td>0.605 ± 0.004</td><td>7.7 ± 0.1</td></tr><tr><td>P17</td><td>7.3 ± 0.5</td><td>0.63 ± 0.04</td><td>8 ± 1</td></tr><tr><td>C18</td><td>6.85 ± 0.07</td><td>0.592 ± 0.006</td><td>7.0 ± 0.2</td></tr></tbody></table>

Notes. The values in the second and third row correspond to the 99% confidence intervals.

References. P17: Pápics et al. (2017); C18: Christophe et al. (2018).

The estimates for $v_{\text{rot}}$ and $P_0$ are given in Table B.3, where the second and third row are obtained by using the peaks in fg1 and fg2, respectively. Our results (fifth and sixth row of Table B.3) clearly reproduce the values published by Pápics et al. (2017) and Christophe et al. (2018). Pápics et al. (2017) identify the frequencies in fg2 as those of prograde sectoral g modes with $m = 2$, whereas Christophe et al. (2018) do not take those frequencies into account. In Fig. B.2 we show that the blue symbols corresponding to the modes in fg2 follow a similar linear trend as the red circles ($m = 1$ modes), supporting the mode identification of Pápics et al. (2017). The estimates for $v_{\text{rot}}$ are consistent within $2\sigma$ among the three independent analyses.

## Appendix C: Case of slow rotation

The method to determine $v_{\text{rot}}$ and $P_0$ based on the $\nu-\sqrt{\Delta\nu}$ diagram is inspired by the approximate expression for the eigenvalue of the Laplace tidal equation, $\lambda$, which is given by Eq. (8). Because this relation is accurate only when the spin parameter $s$ is much larger than 1, it is important to check whether the same method can be applied even for ($0 \le s \le 1$), which implies the case of slow rotation.

For this purpose, we reduce the rotation frequency of model (1) (cf. Sect. 3.1) from 7 µHz to 3 µHz, and calculate the
<figure><img src="image_17.png" /><figcaption>**Fig. C.1.** Diagram of $v$ versus $\sqrt{\Delta v}$ for Model As that rotates with the frequency of 3 µHz.</figcaption></figure>

<table><caption>**Table C.1.** Estimates of $v_{\text{rot}}$ and $P_0$ for Model As.</caption><thead><tr><th></th><th><i>v</i><sub>rot</sub> [µHz]</th><th><i>v</i><sub>rot</sub> [d<sup>-1</sup>]</th><th><i>P</i><sub>0</sub> [10<sup>3</sup> s]</th></tr></thead><tbody><tr><td>Iteration 1</td><td>2.32</td><td>0.20</td><td>3.4</td></tr><tr><td>Iteration 6</td><td>3.08 ± 0.05</td><td>0.266 ± 0.005</td><td>4.60 ± 0.04</td></tr><tr><td>True value</td><td>3</td><td>0.259</td><td>4.579</td></tr></tbody></table>

frequencies of prograde sectoral *g* modes with *m* = 1 and radial orders between -31 and -17. The corresponding range of the spin parameter is between 0.70 and 0.37.

Figure C.1 shows the $v-\sqrt{\Delta v}$ diagram for this slowly-rotating model (Model As). Although Eq. (8) is not very accurate for these modes (cf. Fig. 1), we still observe that all modes follow the linear trend in the diagram. We thus cannot find any qualitative difference from the fast-rotating case in Fig. 2.

Table C.1 provides the estimates for $v_{\text{rot}}$ and $P_0$, which are larger than the true values by only 3% and 0.5%, respectively. Comparing these results with the case of fast rotation in Table 1, we find the following: (1) larger relative uncertainties for the case of slow rotation by about an order of magnitude for $v_{\text{rot}}$ and multiple factors for $P_0$; (2) larger differences between the outputs of the initial and last iterations in spite of the similar number of iterations necessary for the convergence (5 and 6 for model (1) and Model As, respectively). The reason for point (1) is the difference in the number of frequencies that are considered in the analyses. If we extend, to include all the modes between $n = -17$ and $n = -52$, the uncertainties of $v_{\text{rot}}$ and $P_0$ are reduced to values similar to model (1). Point (2) is simply because the deviation from Eq. (8) is more significant. We may still regard that this approximation remains reasonably good even in the limit of slow rotation. In fact, the limiting value of $\sqrt{\lambda} = \sqrt{m(m+1)}$ for the non-rotating case is different from the approximation of $\sqrt{\lambda} = \sqrt{m}$ by only about 20% for $m = 1$. Therefore, the factor $f_k$ given by Eq. (12) changes during iterations by maximum this amount. This property helps the iterative procedure to be robust. However, in the case of slow rotation it is crucial to use the accurate numerical values of $\lambda$, which are shown in Fig. 1, rather than its approximate form of $m^2$ for $s \gg 1$. A failure to do so will impact the accuracy of $v_{\text{rot}}$ and $P_0$.

The method also works for a slow rotation because the asymptotic formula for the oscillation frequencies under the

traditional approximation of rotation (cf. Eq. (5)) can be justified even in the limit of $|s| \ll 1$. In fact, it can generally be shown that, for $|s| \ll 1$, $\lambda$ is given by

$$ \lambda \approx \ell(\ell+1) - ms \quad (C.1) $$

(Longuet-Higgins 1968), from which we can derive

$$ v_{n,l,m} - v_{n,l,0} \approx m \left( 1 - \frac{1}{\ell(\ell+1)} \right) v_{\text{rot}}. \quad (C.2) $$

Equation (C.2) is the standard formula of the rotational splitting for the high-order *g* modes in the case of very slow rotation (cf. Cowling & Newing 1949; Ledoux 1951). Of course, it is more practical in that case to use the rotational splittings rather than the $v-\sqrt{\Delta v}$ method to determine $v_{\text{rot}}$, while $P_0$ can be measured directly from the period spacing.

In conclusion, we find no serious problem in applying the method of the $v-\sqrt{\Delta v}$ diagram to the case of slow rotation if prograde sectoral *g* modes are properly identified.

## Appendix D: Physics of the prograde sectoral *g* modes

### D.1. Constituent waves

Prograde sectoral *g* modes are composed of internal equatorial Kelvin waves when the spin parameter *s* is much larger than 1 (cf. Townsend 2003). This can be understood from the common physical properties between the Kelvin waves and the prograde sectoral *g* modes, which are described in meteorology and oceanography literature (e.g. Gill 1982). These waves are trapped in the equatorial region, with their displacement predominantly parallel to the equator. Their horizontal propagation is in the same longitudinal direction as the rotation, while they also propagate in the radial direction because of the buoyancy force. Unlike the other gravito-inertial waves, they are restored only by the buoyancy force, whereas the Coriolis force does not influence the wave propagation except for controlling the range of the equatorial confinement by contributing to the force balance (geostrophic balance) in the latitudinal direction. This explains why the dispersion relation does not depend on the rotation rate. Having understood the preceding properties, we may call the prograde sectoral *g* modes: “Kelvin *g* modes”. We note that Kelvin modes are sometimes referred to as *f* modes, oscillation modes of non-rotating stars that consist of surface gravity waves (e.g. Chandrasekhar & Lebovitz 1964; Cox 1980). There should be no confusion between the two types of Kelvin modes of totally different physical characters. The Kelvin waves in the rotating case originate from Thomson (1880), while Thomson (1863) studies the oscillation modes of a liquid sphere (a homogeneous incompressible spherical body), which correspond to the surface gravity modes.

### D.2. Variation in $\lambda$

When there is no rotation (with the spin parameter of *s* = 0), the Laplace tidal equation is reduced to the standard associated Legendre differential equation, so that we obtain the eigenvalue of $\lambda_{m,m} = m(m+1)$ for the prograde sectoral *g* modes. As demonstrated in the previous studies (e.g. Bildsten et al. 1996; Lee & Saio 1997), $\lambda_{m,m}$ monotonically decreases as *s* becomes large, and is asymptotically equal to $m^2$ as $s \to \infty$.

We may interpret the transition from $\lambda_{m,m} = m(m+1)$ for $s = 0$ to $\lambda_{m,m} = m^2$ for $s \to \infty$ as being caused by the change
in the geometry. When $s = 0$, the system is totally spherically symmetric, so that the horizontal propagation of the waves can reach all over the (two-dimensional) sphere. From a mathematical point of view, the eigenvalue, $m(m+1)$, is that of the angular part of the Laplace operator, whereas it is physically determined by the quantisation condition that the constituent waves are in phase before and after travelling around the sphere (e.g. Gough 1993). On the other hand, the propagation region of the waves degenerates towards the equator (one-dimensional sphere) in the limit of $s \to \infty$. In addition, the geostrophic balance forces the latitudinal ($\theta$) component of the displacement to be zero, leaving only the longitudinal ($\phi$) component. Therefore, only the longitudinal waves along the equator are allowed. The eigenvalue of $\lambda_{m,m} = m^2$ can again be obtained by the quantisation condition, which requires the same phase of the waves before and after going around the equator. We can formally show that the Laplace tidal equation for $p'$ (the Eulerian perturbation to the pressure) is reduced to $-m^2p' = -\lambda p'$, which means

$$ \frac{\partial^2 p'}{\partial \phi^2} = -\lambda p', \qquad (D.1) $$

at the equator ($\theta = \pi/2$) under the condition of the geostrophic balance. Equation (D.1) is equivalent to the spatial part of the wave equation at the equator.

### D.3. Contribution of the spheroidal and toroidal components

We may discuss the decrease in $\lambda$ as $s$ increases from a different point of view than the one in Appendix D.2.

We first note that $\lambda$ can be regarded as a measure of the divergence of the horizontal components of the displacement,

$$ \nabla_h \cdot \xi_h = \frac{1}{\sin \theta} \frac{\partial}{\partial \theta} (\sin \theta \xi^{\theta}) + \frac{1}{\sin^2 \theta} \frac{\partial \xi^{\phi}}{\partial \phi}, \qquad (D.2) $$

in which $\xi^{\theta}$ and $\xi^{\phi}$ are the $\theta$ and $\phi$ components of the horizontal displacement $\xi_h$, respectively. Under the traditional approximation of rotation, the angular dependence of the eigenfunctions is described by the Hough functions, $\Theta(\cos \theta)$, which are the solutions of the Laplace tidal equation. Based on this, the following can be shown:

$$ \lambda^2 = \frac{\int_{-1}^{1} |\nabla_h \cdot \xi_h|^2 d\mu}{\Xi_h^2 \int_{-1}^{1} \Theta^2(\mu) d\mu}, \qquad (D.3) $$

where $\Xi_h$ indicates the (common) radial dependence of $\xi^{\theta}$ and $\xi^{\phi}$. The denominator on the right-hand side of Eq. (D.3) is just a normalisation factor.

Next, we may suppose to express $\xi_h$ by a sum of two parts, the spheroidal and toroidal components, which are linear combinations of the terms in the form of

$$ \nabla_h Y_{\ell,m} = \left( e_{\theta} \frac{\partial}{\partial \theta} + e_{\phi} \frac{1}{\sin \theta} \frac{\partial}{\partial \phi} \right) Y_{\ell,m} \qquad (D.4) $$

and

$$ -e_r \times \nabla_h Y_{\ell,m} = \left( e_{\theta} \frac{1}{\sin \theta} \frac{\partial}{\partial \phi} - e_{\phi} \frac{\partial}{\partial \theta} \right) Y_{\ell,m}, \qquad (D.5) $$

respectively. Here, $Y_{\ell,m}(\theta, \phi)$ is a spherical harmonic of degree $\ell$ and azimuthal order $m$, while $e_r$, $e_{\theta}$ and $e_{\phi}$ indicate the unit vectors in the $r$, $\theta$ and $\phi$ directions, respectively. Given that $\xi_h$ is expanded into the series, we may regard that $\lambda$ is fixed by

an average of the contribution of each term in the both components, specified by Eqs. (D.4) and (D.5). In the absence of rotation, the eigenfunction, $\xi_h$, is contributed by only one term of the spheroidal component, which satisfies

$$ \nabla_h \cdot \nabla_h Y_{\ell,m} = -\ell(\ell+1)Y_{\ell,m}. \qquad (D.6) $$

We can thus understand $\lambda = \ell(\ell+1)$ in this case. Although the toroidal component also contributes to $\xi_h$ when $s \neq 0$, $\nabla_h \cdot \xi_h$ has no sensitivity to the component (divergence free). From these considerations, we may interpret that the gradual decrease in $\lambda$ as $s$ increases is caused by the decaying contribution of the spheroidal components to $\xi_h$, which implies the growing contribution of the toroidal components.

## Appendix E: Differential rotation

In the description of stellar oscillations under the traditional approximation of rotation, it is usually assumed that the rotation rate, $v_{\text{rot}}$, inside the stars is uniform. However, if we assume that the wavelength of the oscillations is much shorter than the scale height of the equilibrium structure, which can be justified for the low-frequency oscillations of $\gamma$ Dor and SPB stars, we can show that the effect of weak differential rotation on the pulsation frequencies $v$ is negligibly small (cf. Dziembowski et al. 2007; Mathis 2009; Van Reeth et al. 2018).

Under the assumptions that $v_{\text{rot}}$ is a function of only the radius, $r$, and that the scale height of the variation in $v_{\text{rot}}$ is much longer than the wavelength of the oscillations, we can show that the form of the eigenmode condition, given by Eq. (1), remains the same. The only differences are found in that the oscillation frequency in the rotating frame, $v_{\text{co}}$, and the spin parameter, $s$, in Eq. (2) are not constant, but dependent on $r$ through the $r$-dependence of $v_{\text{rot}}$. We may thus formulate the (leading-order) change in the eigenfrequency, $\delta v$, due to the variation of $v_{\text{rot}}$ in the radial direction. Neglecting the frequency dependence of $r_1$ and $r_2$ in Eq. (1), we can derive

$$ \delta v = \frac{2}{s} \left[ 1 + \frac{ms}{2} - \left( 1 + \frac{1}{2} \frac{\partial \ln \lambda_{\ell,m}}{\partial \ln s} \right)^{-1} \right] \langle \delta v_{\text{rot}} \rangle. \qquad (E.1) $$

Here, $\delta v_{\text{rot}}$ means the deviation of $v_{\text{rot}}$ from some constant value, $v_{\text{rot},0}$, while $\langle f \rangle$ generally indicates the average of $f$ defined by

$$ \langle f \rangle = \frac{\int_{r_1}^{r_2} \frac{N}{r} f \, dr}{\int_{r_1}^{r_2} \frac{N}{r} dr} = \frac{P_0}{2\pi^2} \int_{r_1}^{r_2} \frac{N}{r} f \, dr. \qquad (E.2) $$

We note that $s$ and the partial derivative of $\lambda_{\ell,m}$ in Eq. (E.1) are evaluated for $v_{\text{rot}} = v_{\text{rot},0}$. Because the choice of $v_{\text{rot},0}$ is arbitrary (so far as $\delta v_{\text{rot}}$ remains small), we can particularly set $v_{\text{rot},0} = \langle v_{\text{rot}} \rangle$, for which we find $\langle \delta v_{\text{rot}} \rangle = 0$ and hence $\delta v = 0$ for any modes from Eq. (E.1). Thus, the oscillation frequency spectrum is not changed for any different rotation profiles if their averages $\langle v_{\text{rot}} \rangle$ are the same. In other words, the mode frequencies are sensitive to only the average $\langle v_{\text{rot}} \rangle$, and insensitive to any variation around it. We can thus interpret the estimate of the (uniform) rotation rate inferred from the oscillation frequencies as $\langle v_{\text{rot}} \rangle$. According to the evolutionary models of $\gamma$ Dor and SPB stars, the weight function $N/r$ generally gets larger in the deeper layers, takes a maximum just outside the convective core with a sharp peak (due to steep gradients in the distribution of the chemical compositions), and drops to essentially zero in the convective core. Therefore, the average $\langle v_{\text{rot}} \rangle$ is mainly weighted in the layers near but outside the convective core. The rotation frequency of slow rotators that is inferred from high-order $g$ modes is essentially the same average, $\langle v_{\text{rot}} \rangle$ (e.g. Goupil et al. 2013).