Issue 
A&A
Volume 618, October 2018



Article Number  A109  
Number of page(s)  25  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201832777  
Published online  22 October 2018 
Period spacings in red giants^{⋆}
IV. Toward a complete description of the mixedmode pattern
LESIA, Observatoire de Paris, PSL Research University, CNRS, Université Pierre et Marie Curie, Université Paris Diderot, 92195 Meudon, France
email: benoit.mosser@obspm.fr
Received:
6
February
2018
Accepted:
21
July
2018
Context. Oscillation modes with a mixed character, as observed in evolved lowmass stars, are highly sensitive to the physical properties of the innermost regions. Measuring their properties is therefore extremely important to probe the core, but requires some care, due to the complexity of the mixedmode pattern.
Aims. The aim of this work is to provide a consistent description of the mixedmode pattern of lowmass stars, based on the asymptotic expansion. We also study the variation of the gravity offset ε_{g} with stellar evolution.
Methods. We revisit previous works about mixed modes in red giants and empirically test how period spacings, rotational splittings, mixedmode widths, and heights can be estimated in a consistent view, based on the properties of the mode inertia ratios.
Results. From the asymptotic fit of the mixedmode pattern of a large set of red giants at various evolutionary stages, we derive unbiased and precise asymptotic parameters. As the asymptotic expansion of gravity modes is verified with a precision close to the frequency resolution for stars on the red giant branch (10^{−4} in relative values), we can derive accurate values of the asymptotic parameters. We decipher the complex pattern in a rapidly rotating star, and explain how asymmetrical splittings can be inferred. We also revisit the stellar inclinations in two open clusters, NGC 6819 and NGC 6791: our results show that the stellar inclinations in these clusters do not have privileged orientation in the sky. The variation of the asymptotic gravity offset with stellar evolution is investigated in detail. We also derive generic properties that explain under which conditions mixed modes can be observed.
Key words: stars: evolution / stars: interiors / stars: oscillations
Full Table 1 is only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/618/A109
© ESO 2018
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Probing the cores of stars is difficult since, generally, stellar information arises from their photosphere. Fortunately, asteroseismology of evolved stars reveals stellar interiors in a unique and powerful way: gravity waves that propagate throughout the core couple with pressure waves and construct mixed modes that can be observed (Beck et al. 2011; Bedding et al. 2011; Benomar et al. 2014). The measurement of the global seismic properties of these mixed modes then carries unique information on the core structure (e.g., Montalbán et al. 2013; Lagarde et al. 2016; Bossini et al. 2015, 2017). Observations with the space missions CoRoT and Kepler have provided the measurement of the asymptotic period spacings (Mosser et al. 2012b; Vrard et al. 2016), of the differentialrotation profile in red giants (Beck et al. 2012; Deheuvels et al. 2014, 2015), and of the core rotation for about 300 stars analyzed by Mosser et al. (2012c).
Most of the previous studies are based on the measurement and analysis of global seismic parameters, such as the asymptotic large separation Δν and the asymptotic period spacings ΔΠ_{1} (e.g., Miglio et al. 2017). It is now time to access the properties of individual frequencies in red giants. Up to now, most of the studies (e.g., Baudin et al. 2012; Di Mauro et al. 2016) were limited to stars on the red giant branch (RGB). Two main reasons explain this restriction: first, the oscillation spectra benefit from a better relative frequency resolution for this evolutionary stage; second, the oscillation spectra remain simple, with rotational splittings smaller than period spacings. When stars evolve, these features become intricate, so that confusion is possible. For the most evolved stars, mixed modes are no longer observable (e.g., Baudin et al. 2012; Mosser et al. 2013; Stello et al. 2014).
The understanding of any complicated mixedmode oscillation pattern must be based on an unambiguous identification of the modes. Up to now, the most efficient method has relied on the use of the asymptotic expansion, completed by a clear description of the influence of rotation (Mosser et al. 2015). New insights on rotation were provided by the analysis depicted in Gehan et al. (2016), who have developed a methodology to measure rotational splittings in an automated way; Gehan et al. (2017 and 2018) showed how rapid rotation can be addressed efficiently. This efficiency derives from the use of stretched oscillation spectra.
In this work, we first examine in Sect. 2 how the different frequency spacings in the asymptotic mixedmode expansion can be expressed as a function of the mode inertia. New expressions are proposed for the mixedmode spacings and rotational splittings. Case studies are examined in Sect. 3 to test and validate these expressions. In Sect. 4, we take advantage of the precision of the fits to derive accurate asymptotic period spacings and gravity offsets; for the first time, we can exhibit the global evolution of these gravity offsets as a function of stellar evolution. New insights on the rotational splittings are proposed in Sect. 5; in particular, we show how the asymptotic expansion can be used to provide priors based upon physical assumptions for any fitting code used later in the analysis. Finally, we assess the conditions for observing mixed modes, based on global asymptotic parameters only (Sect. 6). Section 7 is devoted to our conclusions.
2. Mixedmode parameters
Following the work of Shibahashi (1979) and Unno et al. (1989), we derived asymptotic expansions of mixed modes for different seismic parameters: eigenfrequencies (Mosser et al. 2012b), period spacings (ChristensenDalsgaard 2012), rotational splittings (Goupil et al. 2013; Deheuvels et al. 2015), and mode widths and mode heights (Grosjean et al. 2014; Belkacem et al. 2015b,a; Mosser et al. 2017a). Here, we intend to revisit all these parameters that depict the mixedmode spectrum in order to provide a more precise and unified view.
2.1. Asymptotic expansion
The asymptotic expansion of mixed modes is an implicit relation between the phases θ_{p} and θ_{g} of the pressure and gravitywave contributions to the mixed modes, respectively. It reads(1)
where q is the coupling factor (Mosser et al. 2017b). The phases are related to the large separation Δν and the period spacing ΔΠ_{1}. The most convenient expressions of the phase refer respectively to the pure^{1} p and g mode spectra(2) (3)
where ν_{p} and ν_{g} are the asymptotic frequencies of pure pressure and gravity modes, respectively, and Δν_{p} is the frequency difference between the consecutive pure pressure radial modes with radial orders n_{p} and n_{p} + 1. In this work, we consider that the radial modes and pure dipole pressure modes obey the universal red giant oscillation pattern (Mosser et al. 2011b) and that the dipole gravity modes follow the asymptotic comblike pattern(4)
where ΔΠ_{1} is the period spacing and ε_{g} is the gravity offset.
Mosser et al. (2015) derived that the variation of the oscillation period P with the mixed radial order n writes(5)
A convenient way to write the parameter ζ is (Hekker & ChristensenDalsgaard 2017)(6)
where 𝒩 = Δν/(ν^{2}ΔΠ_{1}) is the density of gravity modes compared to pressure modes, in other words the number of mixed modes in a Δνwide interval. Compared to the original form presented in Mosser et al. (2015), the rapidly varying phase θ_{g} has been replaced by a function of θ_{p} that varies in a smooth way.
As demonstrated by Goupil et al. (2013) and used by subsequent work (Benomar et al. 2014; Deheuvels et al. 2015), the function ζ is connected to the inertia of mixed modes. Introducing the contributions of the envelope and of the core,(7)
and assuming that the envelope contribution of a mixed mode is similar to the inertia of the closest radial mode (I_{env} ≡ I_{np,0}), we find that the inertia of the dipole mode with mixed radial order n varies as(8)
For the sake of simplicity, we use hereafter the abridged notation I_{n} for the inertia of the dipole mixed modes and I_{0} for the closest radial modes, and follow the same convention for the mode heights and widths.
2.2. Seismic parameters
With ζ, we now intend to express the different seismic parameters.
2.2.1. Period spacing
Following ChristensenDalsgaard (2012) and Mosser et al. (2015), period spacings can be expressed as(9)
This expression is however ambiguous, since ζ may vary significantly between the periods P_{n+1} and P_{n} (> P_{n+1}). Therefore, we prefer to consider the expression resulting from the integration of Eq. (5)(10)
where we consider that the mixedmode radial order n is a continuous variable defined by dn = dτ/ΔΠ_{1}, where τ is the stretched period introduced by Mosser et al. (2015); i.e.,(11)
In fact, n takes consecutive integer values for each mixed mode. In this work, we use an estimate of n = n_{p} + n_{g} derived from the pressure and gravity radial orders; n_{p} is derived from the universal red giant oscillation pattern (Mosser et al. 2011b), whereas n_{g} is given by(12) (13)
where the correcting terms ±1/4 that depend on the evolutionary stage are justified in Sect. 4.4. They differ by 1/2, as depicted by the asymptotic relation (e.g., Tassoul 1980; Benomar et al. 2013). In red giants, the high density 𝒩 of mixed modes implies that n_{g} ≫ n_{p}, so that the mixedmode orders are negative.
From the definition of the stretched period, Eq. (10) reduces to(14)
This evident relation justifies the relevance of Eq. (10) instead of Eq. (9): using ⟨ζ⟩_{n} is necessarily more accurate than using ζ for computing period spacings.
2.2.2. Rotational splitting
As introduced by Goupil et al. (2013), the function ζ is used to express the mixedmode rotational splitting as a function of the mean rotational splittings related to pure gravity or pure pressure modes:(15)
As shown by subsequent works (e.g., Deheuvels et al. 2014; Di Mauro et al. 2016; Triana et al. 2017), it is difficult to derive from the observed rotational splittings more than these two mean quantities.
Again, we have to solve the ambiguity of the meaning of ζ in Eq. (15), since we can either consider the value^{2} ζ(ν_{n,0}), in the framework of the perturbation of the nonrotating frequency ν_{n,0}, or ζ(ν_{n,m}), considering that the inertia to be considered corresponds to the actual frequency ν_{n,m}. By analogy with the equation dealing with the period spacing, we propose to rewrite the rotational splitting δν_{rot} = ν_{n,m} − ν_{n,0}, in the limit case where the mean envelope rotation is negligible compared to the mean core rotation, as(16)
where δν_{rot,core} ≡ δν_{rot,g}. As for the radial order n in Eq. (10), we consider the azimuthal order m as a continuous variable varying from 0 to ±1. So, we have introduced two mean values of ζ,(17) (18)
to account for the period spacings and rotational splittings. The relevance of ⟨ζ⟩_{n} is already proven by Eqs. (10) and (11), whereas the relevance of ⟨ζ⟩_{m} has yet to be demonstrated. If we succeed, we will also have understood the relevance of the use of stretched periods for analyzing the mixed modes (Eq. (11)).
2.3. Mixedmode width, height, and amplitude
The work performed by the gas during one oscillation cycle is the same for all modes, associated with surface damping, when the radiative damping in the BruntVäisälä cavity is considered as negligible. Hence, Benomar et al. (2014) have estimated that the mode width of the mixed modes writes(19)
From this relation, we verify that mixed modes have smaller mode widths than radial modes. However, we recall that a family of stars behave differently, when mixed modes are depressed because of an extra damping in the radiative inner region (Mosser et al. 2012a, 2017a; García et al. 2014).
From Belkacem et al. (2015a) we also derive that the amplitude of a resolved dipole mixed mode is(20)
when the geometrical factor that conducts to a visibility of about 1.54 for red giant dipole modes (Mosser et al. 2012a, 2017a) is omitted. Such amplitudes correspond to similar heights for radial and dipole modes since A^{2} = πΓH/2. When, for nonresolved mixed modes, the width Γ_{n} is less than the frequency resolution δf_{res}, a dilution factor must be considered (Dupret et al. 2009). It expresses(21)
when radial modes are resolved, which is the common case.
2.4. Synthetic mixedmode pattern
The previous ingredients can be used to depict an oscillation pattern. Figure 1 shows the synthetic spectrum of a typical star on the low RGB, based on Eq. (16) for the rotational splittings, on Eq. (19) for the mode widths, and on Eq. (21) for the mode heights of unresolved modes. This spectrum resembles the description derived by Grosjean et al. (2014) from nonadiabatic computations, with a timedependent treatment of convection which provides the lifetimes of radial and nonradial mixed modes.
Fig. 1.
Description of the radial order closest to ν_{max} of the oscillation power spectrum of a typical RGB star. Frequencies, widths, and heights are estimated according to the function ζ. Quadupole and radial modes are plotted in gray, dipole mixed modes in dark blue (m = −1), light blue (m = 0), or purple (m = 1), respectively. 
3. Case studies
In this section, we use RGB stars showing clear oscillation spectra as case studies, in order to test the description of the mixedmode spacings, widths, heights, and rotational splittings, which were previously introduced. The first steps consist in identifying their oscillation spectra and in fitting as many dipole mixed modes as possible. One of the two stars considered here, KIC 6144777 was already investigated in many previous articles (e.g., Corsaro et al. 2015; García Saravia Ortiz de Montellano et al. 2018). The other one, KIC 3955033, was less studied since it shows a complicated mixedmode spectrum; it belongs to the list of red giants with period spacings automatically computed by Vrard et al. (2016). We used data downloaded from the KASOC site^{3}, processed using the Kepler pipeline developed by Jenkins et al. (2010), and corrected from outliers, occasional jumps, and drifts (see García et al. 2011, for details).
3.1. Identification of the mixed modes
The location of the mixed modes primarily relies on the firm identification of the pure pressuremode spectrum. The determination of the large separation Δν, first derived from the envelope autocorrelation function (Mosser & Appourchaux 2009), is based on the universal red giant oscillation pattern. This method provides the efficient identification of the radial modes and helps to locate the frequency ranges where mixed modes cannot be mistaken for radial or quadrupole modes. For ℓ = 1 modes, the secondorder asymptotic expansion writes(22)
where ε_{p} is the acoustic offset, n_{max} = ν_{max}/Δν − ε_{p}, and α = 0.076/n_{max}. The parameter d_{01} is function of the large separation, under the form A + B logΔν (where Δν is expressed in μHz), with A = 0.0553 and B = −0.036, as determined from the largescale analysis along the RGB conducted by Mosser et al. (2014). The accurate determination of d_{01} is crucial for the determination of the pure dipole pressure modes, hence for the determination of the minima of the function ζ. In that respect, the small modulation of the radialmode pattern induced by the soundspeed glitches (Miglio et al. 2010; Vrard et al. 2015) must be considered also. Therefore, we fit the actual position of the radial modes first, then use them to refine the pure pressure dipolemode frequencies according to(23)
The background parameters, derived as in Mosser et al. (2012a), are used to correct the granulation contribution in the frequency range around ν_{max}. Hence, mixed modes can be automatically identified in frequency ranges that have no radial and quadrupole modes when their heights are significantly above the background. The automatic selection of the modes relies on a statistical test: the heighttobackground ratio of the modes must be higher than a threshold level R_{p} in order to reject the null hypothesis to a low probability p. According to Appourchaux et al. (2006), the relation between R_{p} and p depends for longlived modes on the observation duration T_{obs} and on the width Δ_{ν} of the frequency range where a mode is expected. This relation expresses(24)
when expressed in noise unit. This situation applies here, since the precise identification of the mixedmode pattern is based on gravitydominated mixed modes. With 4year observations and the search of a couple of modes in a frequency range Δ_{ν} = Δν/𝒩, the threshold is typically 10 for a secure probability rejection at the 10^{−2} level. In practice, mixed modes with a heighttobackground ratio higher than 10 are used to initiate the fit. A lower threshold is enough for the final agreement, when the synthetic mixedmode pattern based on secure modes can be used to search for longlived mixed modes in narrow frequency ranges. We benefit from the fact that the asymptotic fit is precise and enables to search for thin modes in a frequency range Δ_{ν} narrower than 0.1 μHz. Therefore, a threshold of 7 is enough for rejecting the null hypothesis at the 1%level for these modes whose detection benefits from the information gained by larger peaks. The thin mode widths (Eq. (19)) are of great use to map the observed spectrum: a thin gravitydominated mixed mode must be found in the close vicinity, less than 4 times the mode width, of its expected position. For unresolved peaks, this condition is relaxed to 4 times the frequency resolution. The global seismic parameters of the gravity component are then derived from the methods described in Vrard et al. (2016) and Mosser et al. (2017b), with a leastsquare fit between the observed and asymptotic patterns.
At this stage, global seismic parameters are measured and mixed modes are identified, so that it is possible to measure their individual properties.
3.2. Individual fitting procedure
When fitting individually mixed modes, we aim at testing the validity of the asymptotic expression, but not at reaching the ultimate precision, which is the role of a dedicated fit of individual modes (e.g., Gaulme et al. 2009). Therefore, in order to simplify the fit, we supposed (and checked a posteriori) that all multiplets can be fitted independently. This is not the case in all red giant spectra, but it is verified for most stars on the early RGB or in the red clump.
From the asymptotic fit, we identify in the backgroundcorrected spectrum the power excess associated to each mode. Then, we determine the central frequency of the peak as the barycenter of the power excess. The height H and full width at half maximum Γ are simultaneously derived from the Lorentzian fit of the mode. We use Eqs. (19) and (21) as priors. Modes are fitted individually when the mode density is low, or simultaneously when the Lorentzians used as priors overlap.
The fitted spectrum and the seismic parameters of KIC 6144777, used as a first study case, are given in Fig. 2 and in Table A.1. We note the large agreement between the observed and asymptotic peaks. As in other stars showing a seismic signal with a high signaltonoise ratio (S/N), outliers with a heighttobackground value R higher than 7 are present. Their detection does not invalidate the method presented above: they correspond either to ℓ = 2 or 3 mixed modes, possibly also to ℓ = 4 modes, or to aliases (since the duty cycle is about 93%), or even to noise since the detection of 1 noisy peak with R ≥ 8 is expected in a 30μHz frequency range after 4 years of observation, assuming that the noise statistic is a χ^{2} with two degrees of freedom.
The quality of the fit is shown by the small residuals between the observed frequencies and the asymptotic fits (Fig. 3); we note that these residuals are comparable to the uncertainties, derived from Libbrecht (1992) or slightly larger when the quality of the fit may be affected by the high mode density. These residuals are of about the frequency resolution. When pressuredominated mixed modes are excluded, the standard deviation of the asymptotic fit is 11 nHz. This value represents 1.3 times the frequency resolution δf_{res}, or Δν/1000, or a relative precision at ν_{max} of about 10^{−4}. The quality of the fits is based on a small number of parameters: the radial mode frequencies, the mean location d_{01} of the expected pure pressure dipole modes, and four asymptotic parameters: the period spacing ΔΠ_{1}, the coupling factor q, gravitational offset ε_{g}, and the mean core rotational splitting δν_{rot}. Residuals reach maximum values near the pressuredominated mixed modes: there, deviations of about Δν/200 are observed, to be compared to the amplitudes of pressure glitches of about Δν/40 (Vrard et al. 2015). We suspect that these residuals are mostly due to the variation of the parameter d_{01} with frequency.
Fig. 2.
Fit of the oscillation pattern of the RGB star KIC 6144777, showing the pressure radial orders n_{p} from 9 to 12. The power spectrum density has been divided by the fit of the background. Radial and quadrupole modes are highlighted in red and green. The expected locations of dipole mixed modes are labelled with their mixed radial orders. When detected, mixed modes are highlighted in dark blue (m = −1), light blue (m = 0), or purple (m = 1). ℓ = 3 modes, which are also mixed, are located near the abscissa 0.22; extra peaks in the range [−0.2, −0.05] are mixed quadrupole modes. The gray dashed lines indicate the two thresholds used in this work, corresponding to heighttobackground ratios of 7 and 10. 
Fig. 3.
Relative residuals, multiplied by 1000, between the observed and asymptotic mixedmode frequencies in KIC 6144777. The color of the symbols indicates the azimuthal order: dark blue squares for m = −1, light blue diamonds for m = 0, or purple triangles for m = 1; 1σ uncertainties are also shown. The dashed line corresponds to a perfect fit. The dotted lines show the frequency resolution plus an extramodulation Δν(1 − ζ)/100, which is empirically used to define the quality of the fit. 
3.3. Relationships with ζ
With the identification of the mixedmode pattern, we aim to verify the relevance of the use of ⟨ζ⟩_{n} for the period spacings, to test the relevance of ⟨ζ⟩_{m} for the rotational splittings, and further test the predictions for the mode widths and heights.
3.3.1. Period spacings
Period spacings were fitted with different functions of ζ, according either to the integrated value ⟨ζ⟩_{n} (Eq. (10)) or to the arithmetical mean ζ = (ζ_{n} + ζ_{n+1})/2. The resulting plots are shown in Fig. 4. When ⟨ζ⟩_{n} is not used, one remarks that the ΔP(ζ) relation shows a modulation that results from the concavity of ζ. When ζ is close to unity for gravitydominated mixed modes, no modulation is seen; in the range [0.7, 0.9], where the function is convex, the period spacings are larger than predicted; below 0.7, where the function is concave, the period spacings are smaller than expected. The relation between ΔP and ⟨ζ⟩_{n} does not show such a modulation. Furthermore, the fit with ⟨ζ⟩_{n} is nearly linear, with residuals two times smaller. From this comparison, we confirm that the use of ⟨ζ⟩_{n} is preferable for fitting the period spacings.
Fig. 4.
Period spacings of the RGB star KIC 6144777. Top: plot as a function of the arithmetical mean value (ζ(ν_{n})+ζ(ν_{n + 1}))/2. Bottom: plot as a function of the mean value ⟨ζ⟩_{n}. The colors code the azimuthal orders, as in Fig. 2; the dashed line indicates the 1:1 relation; 1σ uncertainties on both the spacings and the mean values of ζ are indicated by vertical and horizontal rror bars. 
3.3.2. Rotational splittings
We performed similar test for the rotational splittings. We a priori excluded a dependence on ζ(ν_{n,0}), since we clearly observe asymmetrical splittings (see below, Sect. 5.1) that cannot be reproduced with ζ(ν_{n,0}). In fact, the rotation rate of KIC 6144777 is not important enough to observe any difference between the variations with either ζ or ⟨ζ⟩_{m}. We therefore performed the fit of the star KIC 3955033 (Fig. A.3), which shows a much more rapid rotation (Fig. 5). From the comparison of δν_{rot}(ζ) and δν_{rot}(⟨ζ⟩_{m}), we derive that this latter expression is more convenient since it provides a χ^{2} ten times smaller than when using ζ, associated with a much more precise estimate of the core rotation: δν_{rot,core} = 765 ± 10 nHz with ⟨ζ⟩_{m}, versus δν_{rot,core} = 730 ± 50 nHz with ζ. From this test, we conclude positively about the relevance of the use of ⟨ζ⟩_{m} for the rotational splittings.
Fig. 5.
Mean rotation splittings of the RGB star KIC 3955033. Top: plot as a function of ζ. Bottom: plot as a function of the mean value ⟨ζ⟩_{m}. The colors code the azimuthal orders; the dashed line indicates the 1:1 relation; 1σ uncertainties on both the splittings and the mean values of ζ are indicated by vertical and horizontal error bars. 
3.3.3. Widths, amplitudes, and heights
As expected from Eq. (19), the mixedmode width shows large variations: pressuredominated mixed modes have widths comparable to those of the radial modes, contrary to gravitydominated modes that are much thinner (Fig. 6, top panel). Figure 6 also shows the validity of Eq. (19), with the mixedmode width proportional to (1 − ζ), except for low values where the observations resolution hampers the measurement of very thin widths. The precision of the fit is limited by the stochastic excitation, especially for longlived peaks: the presence or absence of signal in a single frequency bin can modify the width in large proportion. This limit added to the limitation in frequency resolution does not allow us to test if small additional radiative damping affects the gravitydominated mixed modes (Dupret et al. 2009; Grosjean et al. 2014).
As shown by Mosser et al. (2015), Eq. (20) has a strong theoretical justification, since it expresses the conservation of energy: the sum of all the energy distributed in the mixed modes corresponds to the energy expected in the single pure pressure mode that should exist in absence of any coupling. So, our result is in line with the findings of Mosser et al. (2012a), who measured that, except for depressed modes, the observed total visibility of dipole modes matches the theoretical expectations.
Due to the stochastic nature of the excitation, the mode heights show a large spread (Fig. 7). Dips in the distributions occur when modes are not resolved. It is however clear in Fig. 7 bottom that the dipole mode heights follow the radial distribution according to the trend of Eq. (21). We note that all mixed modes associated with a given pressure radial order show a systematic behavior. For instance, all mixed modes of KIC 6144777 in the frequency range [137, 143 μHz] associated with the pressure mode n_{p} = 11 show lower amplitudes than expected from the Gaussian fit of the power excess. Such a behavior recalls us that the excitation of a mixed mode is due to its acoustic component.
Fig. 6.
Mode widths as a function of the cyclic frequency (top) or as a function of (1 − ζ)Γ_{0} (bottom), for the RGB star KIC 6144777. Radial modes are plotted with square symbols and dipole mixed modes with +; 1σ uncertainties on Γ are indicated by vertical error bars. The value 2δf_{res}/π (Eq. (21)) plotted as a dotted line is proportional to the 4year long frequency resolution. In the bottom plot, radial modes have been considered too, assuming they have ζ = 0 as pure pressure modes. The dashed line indicates the 1:1 relation. 
Fig. 7.
Mode heights as a function of the cyclic frequency (top) or as a function of the radial mode height, modified when the modes are not resolved, (bottom), for the RGB star KIC 6144777. Radial modes are plotted with square symbols and dipole mixed modes with +. The Fourier spectrum is plotted in red (green) for emphasizing the radial (quadrupole) modes. The dotdashed line provides the expected heights of dipole modes, under the assumption that the power excess mimics a Gaussian relation (dotted line). The dashed line indicates the 1:1 relation. 
3.4. Validation
From these two case studies and from other examples shown in the appendix, we can conclude that the asymptotic fits are relevant at all evolutionary stages, when the S/N is high enough. So, the equations developed in Sect. 2 allow us to depict the mixedmode spectrum with a very high accuracy, when the integrated values ⟨ζ⟩_{n} and ⟨ζ⟩_{m} are considered for the period spacings and the rotational splittings, respectively. Up to now, only red clump stars showing buoyancy glitches cannot be fitted with a single set of parameters.
4. Asymptotic period spacings and gravity offsets
In this section, we show how previous findings can be used to derive accurate period spacings. We also explore the variation of the gravity offsets ε_{g} with stellar evolution. These studies rely on the determination of the puregravity mode pattern.
4.1. Observations
Our analysis was conducted over 372 red giants at various evolutionary stages, mainly from Mosser et al. (2014) and Vrard et al. (2016), with stars also considered in Beck et al. (2012), Kallinger et al. (2012), Deheuvels et al. (2014), and Corsaro et al. (2015). Data were obtained as for the two stars considered in Sect. 3. When available, effective temperatures are from APOGEE spectra (Albareti et al. 2017). Selection criteria are mainly based upon the noise level, with Kepler magnitudes brighter than 12 on the low RGB or 14 for more evolved stars. Following the method exposed in Sect. 3.1, we need data with a S/N high enough to allow the identification of gravitydominated mixed modes. When such modes are too few, measurements are impossible. This condition induces a selection bias, specifically addressed in Sect. 6.
The 372 stars that were analyzed are shown in a seismic diagram (Fig. 8). We considered stars from the low RGB (Fig. A.1) to more evolved RGB stars (Fig. A.2). The spectrum of the evolved RGB star KIC 2443903 (Fig. A.4) corresponds to a case near the limit of visibility of gravitydominated mixed modes, with a mode density 𝒩 ≃ 22.4 close to the limit value above which the detection is impossible (Sect. 6). The fitting process for red clump stars can be achieved only when the amplitude of the buoyancy glitches remains limited (Fig. A.5); the same limitation appears in the secondary red clump (Fig. A.6). In fact, except for redclump stars with large buoyancy glitches (Cunha et al. 2015; Mosser et al. 2015), the asymptotic expansion provides a relevant fit. We could then obtain precise measurements of the asymptotic gravity parameters in Eq. (4) and of their uncertainties for a large number of stars. We must report one exception: KIC 3216736 is the only RGB star of our sample where we found buoyancy glitches and could not provide a relevant fit of the spectrum, but only an échelle diagram based on stretched periods (Fig. 9). Since we have tested more than 160 stars on the RGB, with a systematic approach, we can conclude that the most common case on the RGB is the absence of buoyancy glitches, as expected theoretically (Cunha et al. 2015).
Characterizing the sample we studied in terms of bias is difficult. Apart from the RGB stars that were already studied in detail in previous works, we have mostly treated the stars with increasing KIC numbers. This systematic method implies that we did not introduce any further bias compared to the Kepler sample of red giants. Considering a high enough S/N, which is almost equivalent to select bright stars in the red giant domain, is not supposed to introduce biases either. Contrary to many previous studies, we are not limited to stars showing rotational splittings smaller than the confusion limit (). However, the presence of a strong cutoff (Sect. 6) limits the sample, when gravitydominated mixed modes disappear. Redclump stars with strong buoyancy glitches are absent in our data set since the fitting process requires then to account for the extramodulation, which can be quite large (about ΔΠ_{1}/10). When mixed modes are depressed, the low heighttobackground ratio of the mixed modes allows the measurement of ΔΠ_{1} (Mosser et al. 2017a) but is not enough for fitting the pattern. Both cases deserve specific care beyond the scope of this work.
Fig. 8.
Seismic diagram of the 372 red giants studied in this work, with used as a proxy for the luminosity. The color codes the stellar mass. Stars on the RGB are plotted with + symbols, red clump stars with ⋄. The dotted and dashed lines indicate the limit of the visibility of mixed modes for RGB and clump stars, respectively, as defined by Eq. (31). Evolution tracks on the RGB, computed with MESA for solar metallicity (Gehan et al. 2018), are shown for the stellar masses 1.0, 1.3, 1.6, and 1.9 M_{⊙}. The error box indicates the typical 1σ uncertainties on T_{eff} and ν_{max}. 
Fig. 9.
Stretchedperiod échelle diagram of KIC 3216736, the only RGB star in our sample showing buoyancy glitches. The spectrum is simple, since only m = 0 dipole mixed modes are present, but shows a largeperiod modulation instead of the expected vertical alignment. Modes plotted in light blue are pressure dominated; extra peaks that do not follow the global trend are either ℓ = 3 modes or ℓ = 2 mixed modes. Red figures indicate the radial orders of the radial modes. For clarity, the power spectrum density is also plotted twice, top to tail. 
4.2. Pure gravity modes
The identification of the mixed modes depicted in Sect. 3.1 allows us to retrieve the periods of the pure gravity modes and to infer global asymptotic parameters of the gravity components. We compute these periods from the mixedmode frequencies ν, using Eqs. (1) and (2),(25)
Close to each radial mode, when θ_{p} varies from values less than but close to π/2 to values higher than but close to −π/2, the atan correcting term introduces an offset of −ΔΠ_{1}, which in fact allows to relate the (𝒩 + 1) mixed modes in the Δνwide interval to 𝒩 only gravity modes. In order to use all mixed modes, including the m = 1 components, we corrected first the rotational splittings, using Eq. (16) in order to obtain ν values that are corrected from the rotational splitting.
From the periods of the gravity modes 1/ν_{g}, we could then derive the asymptotic parameters ΔΠ_{1} and ε_{g}, assuming the firstorder asymptotic expression for pure gravity modes (Eq. (4)). In practice, a first estimate of ΔΠ_{1} derived from the formalism of Mosser et al. (2015) and Vrard et al. (2016) is used in Eq. (25), then iterated with a leastsquare fit of the linear variation of the gravity modes (Eq. (4)).
4.3. Asymptotic period spacings
Up to now, measurements of ΔΠ_{1} considering that ε_{g} is a free parameter have been obtained for a few stars only (Buysschaert et al. 2016; Hekker et al. 2018, for 3 and 22 observed stars, respectively). The offset ε_{g} being arbitrarily fixed, Mosser et al. (2012b and 2014) reported a very high precision for the period spacings, of typically 0.1 s for stars on the RGB and 0.3 s in the red clump. However, owing to the choice of ε_{g} = 0, their period spacings were slightly affected by a bias of about a fraction of . The values reported by Vrard et al. (2016), free of any hypothesis on ε_{g}, are not biased but show uncertainties typically five to fifteen times higher than the new values. Their comparison with our data confirms the absence of systematic offsets (Fig. 10). So, the new method ensures accuracy, in the sense that the measurements of ΔΠ_{1} are now free of any hypothesis on ε_{g} and prove that the asymptotic expansion for gravity modes (Eq. (4)) is relevant. The relative accuracy we obtained for the period spacings, assuming Eq. (4), is shown in Fig. 11. The median relative accuracies on the RGB and in the red clump are similar, of about 7 × 10^{−4}. They translate, respectively, into 0.06 s on the RGB and 0.22 s in the red clump; Hekker et al. (2018) reach a similar precision.
Fig. 10.
Comparison of the asymptotic period spacings with previous values. Light blue triangles show the bias in period spacings computed under the assumption ε_{g} = 0 (Mosser et al. 2014), whereas dark blue diamonds are free of this hypothesis (Vrard et al. 2016). 
4.4. Gravity offsets ε_{g}
We could measure ε_{g} for a large set of stars. We however have to face the indetermination of ε_{g} modulo 1: we simply assume that ε_{g} is in the range [−0.5, 0.5]. The ε_{g} values computed for the set of stars presented in the paper is given in Table 1 and plotted in Fig. 12, where an histogram is also given. The complete table is available at the CDS. Uncertainties on ε_{g} are small, related to the uncertainties in ΔΠ_{1} by(26)
This relation comes from the derivative of Eq. (4). As a result, the median uncertainties are of about 0.08 on the RGB and 0.06 in the red clump.
We noticed that the median value of ε_{g} on the RGB is in fact close to 1/4, which is the expected asymptotic value in absence of stratification below the convection zone (Provost & Berthomieu 1986), derived from the contribution ℓ/2 − ε_{as} with ℓ = 1 and ε_{as} = 1/4. Hence, we inferred that the degeneracy on the determination of ε_{g} is removed. We then noted a slight decrease in ε_{g} when stars evolve on the RGB, with an accumulation of values close to 0 for redclump stars. Hekker et al. (2018) reported values of ε_{g} in the range [−0.2, 0.5] for 21 stars on the RGB, but did not identify the accumulation of values in the range [0.20, 0.35] predicted by Takata (2016a) for stars on the low RGB. Our measurements fully confirm this prediction. From a check of their data set, we interpret the differences in ε_{g} as resulting from less precise gravity spacings when large rotational splittings apparently modify the period spacings. As made clear by the recent theoretical developments of the asymptotic expansion (Takata 2016b, a), the accurate measurement of the leadingorder term ΔΠ_{1} is necessary to provide reliable estimates of ε_{g}.
We can study the variation of ε_{g} along stellar evolution. On the RGB, the asymptotic expansion predicts ε_{g} = 1/4 − ϑ (Provost & Berthomieu 1986), where ϑ is a measure of the stratification just below the convection zone. From this dependence, we can infer that the term ϑ is certainly very small for most stars on the low RGB. Higher values are suspected for evolved RGB stars, but with too few stars to firmly conclude, whereas lower values are seen for evolved RGB. We checked that the change of regime of ε_{g} is not associated with the luminosity bump since it occurs for more evolved stars than our sample (Khan et al. 2018). In the red clump, the ϑ correction seems important, on the order of 0.3, with a larger spread than observed on the RGB.
An extended study of ε_{g} can now be performed to use this parameter as a probe of the stratification occurring in the radiative region. This study is however beyond the scope of this work.
Period spacings and gravity offsets.
Fig. 12.
Left: variation of ε_{g} with Δν, with the same style as Fig. 11. The horizontal dark gray domain corresponds to the expected range predicted for RGB stars by Takata (2016a), whereas the dotdashed line shows the value ε_{g,as} = 1/4 derived from the asymptotic expansion (Provost & Berthomieu 1986). Uncertainties on ε_{g} are indicated by vertical lines; uncertainties on Δν are smaller than the symbol size. Right: histograms of the distributions of ε_{g} on the RGB (blue curve) and in the red clump (red curve). The dotdashed line and the gray domain have the same meaning as indicated above. 
5. Rotation
The fits based on the function ζ also allow us to analyze rotational splittings in detail.
5.1. Splitting asymmetry
Recently, asymmetries in the rotational splittings were reported by Deheuvels et al. (2017), as the signature of the combined effects of rotation and mode mixing. Using both perturbative and nonperturbative approaches, they computed neardegeneracy effects and could fit the data. In fact, the asymptotic development of mixed modes also describes the combined effects of rotation and mode mixing, so that the rotational splittings based on ⟨ζ⟩_{m} (Eqs. (16) and (18)) are not symmetric. Inversely, the symmetrical rotational splitting based on ζ (Eq. (15)) does not reproduce the observed asymmetry. Hence, observing asymmetrical triplets is a way to prove the relevance of the use of ⟨ζ⟩_{m} instead of ζ.
Observing the asymmetry is challenging but possible for stars with a rapid rotation rate. As explained by Gehan et al. (2017, 2018), rapid rotation means δν_{rot} ≥ Δν/𝒩 for seismology. This rotation is however very slow in terms of interior structure, so that the formalism developed by Goupil et al. (2013) and Deheuvels et al. (2014), summarized by Eq. (16), remains relevant. It simplifies the study, as shown by Ouazzani et al. (2013) who treated the case where rotational splittings can be as large as Δν. We fitted the mixedmode spectrum of KIC 3955033 with both the symmetrical and asymmetrical splitting. At high radial order n_{p}, it is hard to distinguish them. At low orders, when the rotational splittings exceed the mixedmode spacings, the symmetrical splittings fail whereas the asymmetrical one provides a consistent solution along the whole spectrum. The radial order n_{p} = 8 is shown in Fig. 13, the whole spectrum is shown in Fig. A.3.
Fig. 13.
Splitting asymmetry at low frequency in KIC 3955033. Each dipole mixed mode of the spectrum is labelled with its radial and azimuthal orders. The rotational splittings of the radial orders from −142 to −140, plotted with diamonds, do no match the function ζ. Only the multiplet with n = −141 is complete: the m = +1 splitting is much larger than the m = −1 splitting; the colored regions indicate the ranges over which the function ζ is integrated for the components of the multiplet n = −141. The dashed lines indicate heighttobackground values of 7 and 10. 
5.2. Surface rotation
For stars on the low RGB, surface rotation can be inferred from the rotational splittings (Eq. (15)). The measurement is however difficult, since it results from an extrapolation at ζ = 0, when values are mostly obtained above ζ = 0.6 only (Fig. 5). The highest level of precision, hence the use of ⟨ζ⟩_{m} instead of ζ, is required for deriving a correct estimate of the surface rotation. The case of KIC 3955033 is illustrative, with a negative surface rotation when using ζ; the use of ⟨ζ⟩_{m} provides a null value (5 ± 20 nHz). This case also confirms the general situation shown by previous works (Goupil et al. 2013; Di Mauro et al. 2016; Triana et al. 2017): deriving surface rotation can be achieved for the low RGB only.
5.3. Stellar inclination
From its ability to fit the gravitydominated modes that carry useful information, the asymptotic fit can be used to derive the stellar inclination i. The amplitude of the m = 0 component of the dipole multiplet is proportional to sin^{2} i whereas the sum of the amplitudes of the m = ±1 mode is proportional to cos^{2} i. From Eq. (20), a correction factor of 1/(1 − ζ) should be applied on the amplitudes: its differential effect is however much below the precision one can get on i.
We tested our results on a set of stars for which the inclinations measured with other methods have been obtained. We checked that our results are relevant, with a precision limited by the uncertainties on the amplitude measurements. In order to avoid bias, we consider only peaks with a heighttobackground ratio larger than 8. Nevertheless, we noted that the stochastic excitation of the modes induces a small bias for large inclinations. Equatoron inclinations, near 90^{°}, cannot be retrieved precisely, with measurements reduced toward the range 70–80^{°}. As a consequence, they are rare in our analysis. However, many stars show inclinations that, according to the uncertainties, are compatible with equatoron measurement, so that the bias does not affect the following analysis.
We measured inclinations of red giants in the open clusters NGC 6819 observed by Kepler (e.g., Basu et al. 2011; Stello et al. 2011; Miglio et al. 2012). We selected the stars that exhibit mixed modes and could fit 20 mixedmode patterns with the asymptotic expansion. In one case, the asymptotic fit is impossible, due to a low S/N. In two other cases, different possible solutions exist, based either on different period spacings, or on different rotational splittings, but without any ambiguity for the inclination measurement: when two peaks dominate per period spacing, the inclination is necessarily high, whereas it is low when one single peak only is present. We completed this list with other NGC 6819 members listed in Handberg et al. (2017) and could fit two additional stars, which incidentally show a large inclination. Results for the inclinations and rotational splittings are given in Table 2. As shown in Fig. 14, the distribution of the stellar inclinations mimics the sin i relation expected for random inclinations, except near 90^{°}, due to bias mentioned above. A similar test performed on the open cluster NGC 6791 reaches the same conclusion.
Low stellar inclinations in NGC 6819 and 6791 were measured by Corsaro et al. (2017), using a Bayesian analysis, from which aligned spins were inferred. Our measurements however contradict their claim, as shown in Fig. 15 for NGC 6819. In fact, our measurements compared to theirs agree for high inclinations, whereas they mostly disagree for low inclinations. Comparison with the asymptotic fits shows that their Bayesian rotational splittings are most often overestimated and that the related inclinations are most often underestimated.
In fact, the asymptotic fit can be used as a prior for the Bayesian fit. It indicates that the rotational splitting is derived from the thin gravitydominated mixed modes, with narrow widths (Eq. (19)) and average rotational splittings slightly inferior to the mean core rotation δν_{rot} (Eq. (16)). Mixed modes at low pressure radial orders, with frequencies much below ν_{max}, are especially informative, since previous work has shown that their radial mode widths, hence their mixedmode widths according to Eq. (19), are the thinnest possible (Fig. 5 of Vrard et al. 2018). Figure 1 of the supplementary material of Corsaro et al. (2017) provides an explanation of the discrepant Bayesian values. Their fit of the star KIC 5112373 in NGC 6819 provides nearly uniform large mode widths, relevant for the pressuredominated mixed modes but much too high for gravitydominated modes, in contradiction with the physical variation indicated by Eq. (19). As a consequence, their fit assumes that all the power is concentrated in the m = 0 mode; the resulting stellar inclination is 20 ± 8^{°}. We show the asymptotic solution of KIC 5112373 in Fig. 16, with thin gravity dominated mixed modes and the clear identification of triplets. Since m = ±1 modes are observed all along the spectrum, our solution for the inclination is larger, about 47 ± 18^{°}.
We provide another example with the star KIC 2437976, a NGC 6791 member. As shown in Fig. 17, rotational splittings are explained in a consistent way with thin unresolved gravitydominated mixed modes and a rotation rate rapid enough to ensure that close modes do not belong to the same multiplets. All peaks can be explained by the m = ±1 modes. In practice, m = 0 modes are absent, so that this star has necessarily an inclination close to 90^{°}, whereas Corsaro et al. (2017) measured i ≃ 0^{°}. We conclude that some of the low inclinations reported in Corsaro et al. (2017) are incompatible with the analysis presented here. It seems that the difference is due to a too low range of the linewidth priors in the Bayesian analysis, which favors a solution with a low inclination angle and a high splitting. As a result, stellar spins in old open clusters are neither aligned nor quasi parallel to the line of sight. Our study emphasizes a major role for the asymptotic analysis: providing relevant estimates of all features of the mixedmode pattern, including mode widths.
Asymptotic and rotational parameters in NGC 6819.
Fig. 14.
Histogram of the inclinations measured in NGC 6819. The dashed line indicates the sin i distribution. 
Fig. 15.
Comparison of the inclinations of the spin axis of the stars in NGC 6819. Inclinations measured by Corsaro et al. (2017) are plotted on the xaxis, while inclinations derived from the asymptotic analysis are on the yaxis; 1σ uncertainties are indicated by vertical and horizontal error bars. The red symbol shows the case where no asymptotic fit could be found, and the green ones to cases without Bayesian fit. The dashed line corresponds to the 1:1 relation. 
Fig. 16.
Fit of the mixed modes corresponding to n_{p} = 7 in KIC 5112373 (NGC 6819 member). The color codes the azimuthal order: m = +1 in purple, m = −1 in blue. The gray dashed lines indicate the two thresholds used in this work, corresponding to heighttobackground ratios of 7 and 10. Contrary to the analysis conducted by Corsaro et al. (2017), modes with m = ±1 are clearly identified. 
Fig. 17.
Fit of the mixed modes corresponding to n_{p} = 9 in KIC 2437976 (NGC 6791 member). The color codes the azimuthal order: m = +1 in purple, m = −1 in blue. The location of m = 0 modes is indicated in light blue, but none shows a large height for this star seen equatoron. The gray dashed lines indicate the two thresholds used in this work, corresponding to heighttobackground ratios of 7 and 10. Many peaks above the threshold value 5.5 that rejects the null hypothesis at the 5%level follow the mixedmode pattern. 
6. Observability of the mixed modes
All the information derived from mixed modes relies on their observability. The properties of the function ζ can be used to assess under which conditions mixed modes can be actually observed. To achieve this, we investigate first the domain where pressuredominated mixed modes are observed, then the condition for observing gravitydominated mixed modes.
6.1. Pressuredominated mixed modes
We can define the frequency range where mixed modes are pressuredominated (pm) from the full width at half minimum of the ζ function. So, these modes cover a range, expressed in terms of the pressure phase θ_{p} (Eq. (6)), verifying(27)
under the assumption that q is small, which is verified for all stars except at the transition between subgiants and red giants (Mosser et al. 2017b). When expressed in frequency and compared to the large separation, this condition corresponds to a frequency range surrounding each pure pressure modes with a width δν_{pm} defined by(28)
The variations in q and 𝒩 explain the narrowing of the region with pressuredominated mixed modes when stars evolve on the RGB. An example is shown in Fig. A.4. The expression of δν_{pm} also shows that redclump stars, with larger q show pressuredominated mixed modes in a broader region than RGB stars.
6.2. Visible gravitydominated mixed modes
The nondilution of the mode height expressed by Eq. (21) can be used to define a criterion of visibility of the gravitydominated (gm) mixed modes. So, they are clearly visible when they show heights similar to those of the pressure modes (H_{n} = H_{0}), hence when Γ_{0}(1 − ζ) ≥ 2δf_{res}/π (Eq. (19)). This condition translates into(29)
Except at the transition from subgiants to red giants, where mixed modes are unambiguously visible (Benomar et al. 2013; Deheuvels et al. 2014), the terms q^{2} are negligible, so that modes are clearly visible if(30)
This condition for observing gravitydominated mixed modes has many consequences:

It can be fulfilled only if the definition of the right term is ensured, which requires a frequency resolution low enough compared to the radial mode width. With Γ_{0} in the range [100, 150 nHz], the observation must last 50–75 days at least. In fact, mixed modes were observable with CoRoT runs lasting about 150 days (Mosser et al. 2011a), but are hardly observable with K2 80day time series (Stello et al. 2017).

When stars evolve on the RGB, the decrease in q and increase in 𝒩 contribute to the narrowing of observable modes. Mixed modes are more easily visible in the red clump, owing to larger q values. This criterion is implicitly used by Elsworth et al. (2017) for their determination of the evolutionary state of redgiant stars.

All mixed modes are clearly visible when the condition expressed by Eq. (30) is always met, that is when 𝒩 ≤ q(πΓ_{0}/2δf_{res} − 1). This condition is met for subgiants, on the lower RGB, and for secondaryclump stars (Mosser et al. 2014).

No mixed mode can be observed when the condition is so drastic that only pressuredominated mixed modes can be observed. The combination of the conditions expressed by Eqs. (27) and (30) yields the limit of visibility of gravitydominated mixed modes, expressed by a condition on the mixedmode density (31)
In the conditions of observation of Kepler, with typical parameters defined as in Mosser et al. (2017a), this limit corresponds to a mode density 𝒩 of about 25, for RGB and clump stars, over which no gravitydominated mixed modes can be identified. This theoretical estimate is observed in practice, with a few exceptions with larger 𝒩 (Fig. 8). On the RGB, observation of mixed modes with Kepler is limited to Δν ≥ 6 μHz, whereas the limit is around 3 μHz for clump stars. As a consequence, visible mixed modes in an oscillation spectrum with Δν in the range [3, 6 μHz] most often indicate a redclump star. Incidentally, the location of the RGB bump was recently identified by Khan et al. (2018) in the range [5, 6 μHz], depending on the stellar mass and metallicity. This means that sounding the bump with mixed modes will be very difficult, if not impossible.
7. Conclusion
The asymptotic analysis allows us to depict the whole properties of the mixedmode spectrum in a consistent way. Period spacings, rotational splittings, mode widths, and mode heights, all depend on the mode inertia so that all are related to the parameter ζ. We could derive interesting properties:

The asymptotic fit of the mixed modes proves to be precise and unbiased. Its precision for the RGB stars is so high that the asymptotic expansion of gravity modes can be validated when buoyancy glitches are absent. This ensures the delivery of accurate asymptotic parameters ΔΠ_{1}, q, and ε_{g}. We found only one RGB star with such buoyancy glitches; on the contrary, buoyancy glitches are often present in redclump stars.

The period spacings and rotational splittings are better estimated with integrated values of the function ζ. The use of these mean values ⟨ζ⟩_{n} and ⟨ζ⟩_{m} is useful for evolved RGB stars and is mandatory for stars with intricate splittings and spacings. Using the stretched period (Mosser et al. 2015) is in fact equivalent.

The gravity asymptotic parameters ΔΠ_{1} and ε_{g} can now be accurately determined, with typical accuracy of respectively 0.06 s and 0.1 on the RGB, and 0.22 s and 0.08 in the red clump. This opens the way to a fruitful dialogue with theoretical developments (Takata 2006, 2016b,a) and modeling (e.g., Bossini et al. 2015; Cunha et al. 2015).

We have made clear that observing mixed modes in evolved red giants requires an observation duration longer than ≃100 days. However, gravitydominated mixed modes are no longer observable when the stars are more evolved than Δν ≃ 6 μHz on the RGB, or Δν ≃ 3 μHz in the red clump. These thresholds are indicative values: the natural spread of the seismic parameters with respect to their mean values explain slight differences.

We have demonstrated the nonalignment of the rotation axis of the stars belonging to the old open clusters NGC 6791 and NGC 6819. These results contradict previous findings by Corsaro et al. (2017) and illustrate how useful the asymptotic fit will be in the future when used to define priors to any Bayesian or other type of fit of mixed modes.
Acknowledgments
We thank the entire Kepler team, whose efforts made these results possible. BM warmly thanks Yvonne Elsworth, James Kuszlewicz and Masao Takata for their comments on the draft submitted for internal review on the website of the Kepler Asteroseismic Science Operations Center. We acknowledge financial support from the Programme National de Physique Stellaire (CNRS/INSU). BM acknowledges the support of the International Space Institute (ISSI) for the program AsteroSTEP (Asteroseismology of STEllar Populations).
References
 Albareti, F. D., Allende Prieto, C., Almeida, A., et al. 2017, ApJS, 233, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Appourchaux, T., Berthomieu, G., Michel, E., et al. 2006, in The CoRoT Mission PreLaunch Status – Stellar Seismology and Planet Finding, eds. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, ESA SP, 1306, 377 [NASA ADS] [Google Scholar]
 Basu, S., Grundahl, F., Stello, D., et al. 2011, ApJ, 729, L10 [NASA ADS] [CrossRef] [Google Scholar]
 Baudin, F., Barban, C., Goupil, M. J., et al. 2012, A&A, 538, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Beck, P. G., Bedding, T. R., Mosser, B., et al. 2011, Science, 332, 205 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Bedding, T. R., Mosser, B., Huber, D., et al. 2011, Nature, 471, 608 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Belkacem, K., Marques, J. P., Goupil, M. J., et al. 2015a, A&A, 579, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Belkacem, K., Marques, J. P., Goupil, M. J., et al. 2015b, A&A, 579, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Benomar, O., Bedding, T. R., Mosser, B., et al. 2013, ApJ, 767, 158 [NASA ADS] [CrossRef] [Google Scholar]
 Benomar, O., Belkacem, K., Bedding, T. R., et al. 2014, ApJ, 781, L29 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Bossini, D., Miglio, A., Salaris, M., et al. 2015, MNRAS, 453, 2290 [NASA ADS] [CrossRef] [Google Scholar]
 Bossini, D., Miglio, A., Salaris, M., et al. 2017, MNRAS, 469, 4718 [NASA ADS] [CrossRef] [Google Scholar]
 Buysschaert, B., Beck, P. G., Corsaro, E., et al. 2016, A&A, 588, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 ChristensenDalsgaard, J. 2012, in Progress in Solar/Stellar Physics with Helio and Asteroseismology, eds. H. Shibahashi, M. Takata, & A. E. LynasGray, ASP Conf. Ser., 462, 503 [NASA ADS] [Google Scholar]
 Corsaro, E., De Ridder, J., & García, R. A. 2015, A&A, 579, A83 [NASA ADS] [CrossRef] [EDP Sciences] [PubMed] [Google Scholar]
 Corsaro, E., Lee, Y.N., García, R. A., et al. 2017, Nat. Astron., 1, 0064 [NASA ADS] [CrossRef] [Google Scholar]
 Corsaro, E., De Ridder, J., & García, R. A. 2018, A&A, 612, C2 [CrossRef] [EDP Sciences] [Google Scholar]
 Cunha, M. S., Stello, D., Avelino, P. P., ChristensenDalsgaard, J., & Townsend, R. H. D. 2015, ApJ, 805, 127 [Google Scholar]
 Deheuvels, S., Doğan, G., Goupil, M. J., et al. 2014, A&A, 564, A27 [Google Scholar]
 Deheuvels, S., Ballot, J., Beck, P. G., et al. 2015, A&A, 580, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Ouazzani, R. M., & Basu, S. 2017, A&A, 605, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Di Mauro, M. P., Ventura, R., Cardini, D., et al. 2016, ApJ, 817, 65 [NASA ADS] [CrossRef] [Google Scholar]
 Dupret, M., Belkacem, K., Samadi, R., et al. 2009, A&A, 506, 57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Elsworth, Y., Hekker, S., Basu, S., & Davies, G. R. 2017, MNRAS, 466, 3344 [NASA ADS] [CrossRef] [Google Scholar]
 García, R. A., Hekker, S., Stello, D., et al. 2011, MNRAS, 414, L6 [NASA ADS] [CrossRef] [Google Scholar]
 García, R. A., Pérez Hernández, F., Benomar, O., et al. 2014, A&A, 563, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 García Saravia Ortiz de Montellano, A., Hekker, S., & Themeßl, N. 2018, MNRAS, 476, 1470 [NASA ADS] [CrossRef] [Google Scholar]
 Gaulme, P., Appourchaux, T., & Boumier, P. 2009, A&A, 506, 7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gehan, C., Mosser, B., & Michel, E. 2016, ArXiv eprints [arXiv:1612.05414] [Google Scholar]
 Gehan, C., Mosser, B., & Michel, E. 2017, Eur. Phys. J. Web Conf., 160, 04005 [CrossRef] [Google Scholar]
 Gehan, C., Mosser, B., Michel, E., Samadi, R., & Kallinger, T. 2018, A&A, 616, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Goupil, M. J., Mosser, B., Marques, J. P., et al. 2013, A&A, 549, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Grosjean, M., Dupret, M.A., Belkacem, K., et al. 2014, A&A, 572, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Handberg, R., Brogaard, K., Miglio, A., et al. 2017, MNRAS, 472, 979 [NASA ADS] [CrossRef] [Google Scholar]
 Hekker, S., & ChristensenDalsgaard, J. 2017, A&ARv, 25, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hekker, S., Elsworth, Y., & Angelou, G. C. 2018, A&A, 610, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jenkins, J. M., Caldwell, D. A., Chandrasekaran, H., et al. 2010, ApJ, 713, L87 [NASA ADS] [CrossRef] [Google Scholar]
 Kallinger, T., Hekker, S., Mosser, B., et al. 2012, A&A, 541, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Khan, S., Hall, O. J., Miglio, A., et al. 2018, ApJ, 859, 156 [NASA ADS] [CrossRef] [Google Scholar]
 Lagarde, N., Bossini, D., Miglio, A., Vrard, M., & Mosser, B. 2016, MNRAS, 457, L59 [NASA ADS] [CrossRef] [Google Scholar]
 Libbrecht, K. G. 1992, ApJ, 387, 712 [Google Scholar]
 Miglio, A., Montalbán, J., Carrier, F., et al. 2010, A&A, 520, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Miglio, A., Brogaard, K., Stello, D., et al. 2012, MNRAS, 419, 2077 [NASA ADS] [CrossRef] [Google Scholar]
 Miglio, A., Chiappini, C., Mosser, B., et al. 2017, Astron. Nachr., 338, 644 [NASA ADS] [CrossRef] [Google Scholar]
 Montalbán, J., Miglio, A., Noels, A., et al. 2013, ApJ, 766, 118 [NASA ADS] [CrossRef] [Google Scholar]
 Mosser, B., & Appourchaux, T. 2009, A&A, 508, 877 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Barban, C., Montalbán, J., et al. 2011a, A&A, 532, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Belkacem, K., Goupil, M., et al. 2011b, A&A, 525, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Elsworth, Y., Hekker, S., et al. 2012a, A&A, 537, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012b, A&A, 540, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012c, A&A, 548, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Dziembowski, W. A., Belkacem, K., et al. 2013, A&A, 559, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Benomar, O., Belkacem, K., et al. 2014, A&A, 572, L5 [NASA ADS] [CrossRef] [EDP Sciences] [PubMed] [Google Scholar]
 Mosser, B., Vrard, M., Belkacem, K., Deheuvels, S., & Goupil, M. J. 2015, A&A, 584, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Belkacem, K., Pinçon, C., et al. 2017a, A&A, 598, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Pinçon, C., Belkacem, K., Takata, M., & Vrard, M. 2017b, A&A, 600, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ouazzani, R.M., Goupil, M. J., Dupret, M.A., & Marques, J. P. 2013, A&A, 554, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Provost, J., & Berthomieu, G. 1986, A&A, 165, 218 [NASA ADS] [Google Scholar]
 Shibahashi, H. 1979, PASJ, 31, 87 [NASA ADS] [Google Scholar]
 Stello, D., Meibom, S., Gilliland, R. L., et al. 2011, ApJ, 739, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Stello, D., Compton, D. L., Bedding, T. R., et al. 2014, ApJ, 788, L10 [NASA ADS] [CrossRef] [Google Scholar]
 Stello, D., Zinn, J., Elsworth, Y., et al. 2017, ApJ, 835, 83 [NASA ADS] [CrossRef] [Google Scholar]
 Takata, M. 2006, PASJ, 58, 893 [Google Scholar]
 Takata, M. 2016a, PASJ, 68, 109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Takata, M. 2016b, PASJ, 68, 91 [NASA ADS] [CrossRef] [Google Scholar]
 Tassoul, M. 1980, ApJS, 43, 469 [NASA ADS] [CrossRef] [Google Scholar]
 Triana, S. A., Corsaro, E., De Ridder, J., et al. 2017, A&A, 602, A62 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial Oscillations of Stars (Tokyo: University of Tokyo Press) [Google Scholar]
 Vrard, M., Mosser, B., Barban, C., et al. 2015, A&A, 579, A84 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vrard, M., Mosser, B., & Samadi, R. 2016, A&A, 588, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Vrard, M., Kallinger, T., Mosser, B., et al. 2018, A&A, 616, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A:
Seismic parameters
We used KIC 6144777 as a case study (Fig. 2). Table A.1 provides the fit of its radial dipole mixed modes. Our results are in agreement with those published by Corsaro et al. (2015) and derive a similar number of modes (about 100), but also show differences:

The determination of the frequencies in Corsaro et al. (2015) can be as precise as 0.3 nHz. This precision of about δf_{res}/30 was corrected into about δf_{res}/10 in their corrigendum (Corsaro et al. 2018), which remains surprisingly good; the frequencies we obtain are given with a precision that is at best about half the frequency resolution (≃4 nHz).

Their mode widths are quite different and, most often, larger than ours;

Heights also differ, which can come from a different treatment of the time series.
The potential of the comparison between methods based on different principles is very high: coupling the physics of the asymptotic expansion and the power of a pure numerical approach is the next step for delivering duly identified mixed modes.
The échelle diagrams of the stars mentioned in the main text are also presented:

KIC 10272858 lies on the low part of the RGB (Fig. A.1);

KIC 11353313 is on the RGB (Fig. A.2);

KIC 3955033 is a RGB star with a rapid core rotation (Fig. A.3); its frequencies are given in Table A.2;

KIC 2443903 is more evolved on the RGB, at the limit of detection of mixed modes (Fig. A.4);

KIC 1723700 is in the red clump star (Fig. A.5);

and KIC 1725190 is a secondary red clump star (Fig. A.6).
Oscillation pattern of the RGB star KIC 6144777.
Oscillation pattern of the RGB star KIC 3955033.
Fig. A.1.
Fit of the oscillation pattern of the low RGB star KIC 10272858, at the limit of validity of the asymptotic pattern. Owing to the small radial orders, small shifts are seen between observed and asymptotic spectra. Same style as Fig. 2, but ℓ = 3 modes appear near the abscissa 0.28. 
Fig. A.3.
Fit of the oscillation pattern of the RGB star KIC 3955033. The overlap of mixed modes with different mixedmode orders is the signature of the rapid core rotation. The second rotation crossing, where all components of the multiplets overlap (Gehan et al. 2017), occurs at the mixedorder n = −122. Same style as Fig. 2. 
Fig. A.4.
Fit of the oscillation pattern of the evolved RGB star KIC 2443903, near the limit of capability of identification, with a large crowding due to the high mode density. The second rotation crossing, where all m components apparently coincide, occurs at n = −189 (with an abscissa ≃0.1 and n_{p} = 9); the third crossing, where m = 1 components apparently coincide with m = 0 inbetween, occurs at n = −233 (with an abscissa ≃0.25 and n_{p} = 7). Same style as Fig. 2. Note that the modes with large heights at an abscissa ≃ 0.2 are ℓ = 3 modes. 
Fig. A.5.
Fit of the oscillation pattern of the redclump star KIC 1723700. Buoyancy glitches explain the small shifts between observed and asymptotic spectra but do no hamper the mode identification. Same style as Fig. 2. 
Fig. A.6.
Fit of the oscillation pattern of the secondaryclump star KIC 1725190. Buoyancy glitches explain the small shifts between observed and asymptotic spectra but do no hamper the mode identification. Same style as Fig. 2. 
Appendix B:
Stars in open clusters
All stars studied by Corsaro et al. (2017) were investigated. The fitting process is challenging, due to the dim magnitudes of such dim stars in open clusters. However, the combination of all pressure radial orders near ν_{max} provides in most cases an unambiguous fit, and at least a few mixedmode radial orders provide clear splittings.

Figure B.1 provides the asymptotic fit of KIC 5024476, member of the open cluster NGC 6819 observed by Kepler. We note that m = ±1 modes are clearly identified and derive a stellar inclination i = 79 ± 11^{°} for this star. This result is in disagreement with Corsaro et al. (2017) who found an inclination i = 20 ± 7^{°}.

Similar conclusions are reached for KIC 2437976 (Fig. B.2), member of the open cluster NGC 6791. Corsaro et al. (2017) found an inclination i = 0 ± 10^{°}, despite the fact m = 1 modes are clearly identified and indicate i = 76 ± 14^{°}.
Fig. B.1.
Fit of the oscillation pattern of the RGB star KIC 5024476, member of the open cluster NGC 6819. The dim magnitude of the cluster stars explains the low S/N. However, unambiguous doublets are identified all along the mixedmode spectrum; m = 0 modes are mostly absent and m = 1 modes dominate the mixedmode spectrum, so that a nearly poleon inclination is not possible. Same style as Fig. 2. 
Fig. B.2.
Fit of the oscillation pattern of the RGB star KIC 2437976, member of the open cluster NGC 6791. The dim magnitude of the cluster stars explains the low S/N. The identification at radial order 9, supported by the radial orders 8 and 10, is unambiguously conclusive: m = 0 modes are mostly absent and m = 1 modes dominate the mixedmode spectrum, so that a poleon inclination is not possible. Same style as Fig. 2. 
Oscillation pattern of the redclump star KIC 5024476 in NGC 6819.
Oscillation pattern of the RGB star KIC 2437976 in NGC 6791.
All Tables
All Figures
Fig. 1.
Description of the radial order closest to ν_{max} of the oscillation power spectrum of a typical RGB star. Frequencies, widths, and heights are estimated according to the function ζ. Quadupole and radial modes are plotted in gray, dipole mixed modes in dark blue (m = −1), light blue (m = 0), or purple (m = 1), respectively. 

In the text 
Fig. 2.
Fit of the oscillation pattern of the RGB star KIC 6144777, showing the pressure radial orders n_{p} from 9 to 12. The power spectrum density has been divided by the fit of the background. Radial and quadrupole modes are highlighted in red and green. The expected locations of dipole mixed modes are labelled with their mixed radial orders. When detected, mixed modes are highlighted in dark blue (m = −1), light blue (m = 0), or purple (m = 1). ℓ = 3 modes, which are also mixed, are located near the abscissa 0.22; extra peaks in the range [−0.2, −0.05] are mixed quadrupole modes. The gray dashed lines indicate the two thresholds used in this work, corresponding to heighttobackground ratios of 7 and 10. 

In the text 
Fig. 3.
Relative residuals, multiplied by 1000, between the observed and asymptotic mixedmode frequencies in KIC 6144777. The color of the symbols indicates the azimuthal order: dark blue squares for m = −1, light blue diamonds for m = 0, or purple triangles for m = 1; 1σ uncertainties are also shown. The dashed line corresponds to a perfect fit. The dotted lines show the frequency resolution plus an extramodulation Δν(1 − ζ)/100, which is empirically used to define the quality of the fit. 

In the text 
Fig. 4.
Period spacings of the RGB star KIC 6144777. Top: plot as a function of the arithmetical mean value (ζ(ν_{n})+ζ(ν_{n + 1}))/2. Bottom: plot as a function of the mean value ⟨ζ⟩_{n}. The colors code the azimuthal orders, as in Fig. 2; the dashed line indicates the 1:1 relation; 1σ uncertainties on both the spacings and the mean values of ζ are indicated by vertical and horizontal rror bars. 

In the text 
Fig. 5.
Mean rotation splittings of the RGB star KIC 3955033. Top: plot as a function of ζ. Bottom: plot as a function of the mean value ⟨ζ⟩_{m}. The colors code the azimuthal orders; the dashed line indicates the 1:1 relation; 1σ uncertainties on both the splittings and the mean values of ζ are indicated by vertical and horizontal error bars. 

In the text 
Fig. 6.
Mode widths as a function of the cyclic frequency (top) or as a function of (1 − ζ)Γ_{0} (bottom), for the RGB star KIC 6144777. Radial modes are plotted with square symbols and dipole mixed modes with +; 1σ uncertainties on Γ are indicated by vertical error bars. The value 2δf_{res}/π (Eq. (21)) plotted as a dotted line is proportional to the 4year long frequency resolution. In the bottom plot, radial modes have been considered too, assuming they have ζ = 0 as pure pressure modes. The dashed line indicates the 1:1 relation. 

In the text 
Fig. 7.
Mode heights as a function of the cyclic frequency (top) or as a function of the radial mode height, modified when the modes are not resolved, (bottom), for the RGB star KIC 6144777. Radial modes are plotted with square symbols and dipole mixed modes with +. The Fourier spectrum is plotted in red (green) for emphasizing the radial (quadrupole) modes. The dotdashed line provides the expected heights of dipole modes, under the assumption that the power excess mimics a Gaussian relation (dotted line). The dashed line indicates the 1:1 relation. 

In the text 
Fig. 8.
Seismic diagram of the 372 red giants studied in this work, with used as a proxy for the luminosity. The color codes the stellar mass. Stars on the RGB are plotted with + symbols, red clump stars with ⋄. The dotted and dashed lines indicate the limit of the visibility of mixed modes for RGB and clump stars, respectively, as defined by Eq. (31). Evolution tracks on the RGB, computed with MESA for solar metallicity (Gehan et al. 2018), are shown for the stellar masses 1.0, 1.3, 1.6, and 1.9 M_{⊙}. The error box indicates the typical 1σ uncertainties on T_{eff} and ν_{max}. 

In the text 
Fig. 9.
Stretchedperiod échelle diagram of KIC 3216736, the only RGB star in our sample showing buoyancy glitches. The spectrum is simple, since only m = 0 dipole mixed modes are present, but shows a largeperiod modulation instead of the expected vertical alignment. Modes plotted in light blue are pressure dominated; extra peaks that do not follow the global trend are either ℓ = 3 modes or ℓ = 2 mixed modes. Red figures indicate the radial orders of the radial modes. For clarity, the power spectrum density is also plotted twice, top to tail. 

In the text 
Fig. 10.
Comparison of the asymptotic period spacings with previous values. Light blue triangles show the bias in period spacings computed under the assumption ε_{g} = 0 (Mosser et al. 2014), whereas dark blue diamonds are free of this hypothesis (Vrard et al. 2016). 

In the text 
Fig. 11.
Relative precision of the asymptotic period spacings. Same style as in Fig. 8. 

In the text 
Fig. 12.
Left: variation of ε_{g} with Δν, with the same style as Fig. 11. The horizontal dark gray domain corresponds to the expected range predicted for RGB stars by Takata (2016a), whereas the dotdashed line shows the value ε_{g,as} = 1/4 derived from the asymptotic expansion (Provost & Berthomieu 1986). Uncertainties on ε_{g} are indicated by vertical lines; uncertainties on Δν are smaller than the symbol size. Right: histograms of the distributions of ε_{g} on the RGB (blue curve) and in the red clump (red curve). The dotdashed line and the gray domain have the same meaning as indicated above. 

In the text 
Fig. 13.
Splitting asymmetry at low frequency in KIC 3955033. Each dipole mixed mode of the spectrum is labelled with its radial and azimuthal orders. The rotational splittings of the radial orders from −142 to −140, plotted with diamonds, do no match the function ζ. Only the multiplet with n = −141 is complete: the m = +1 splitting is much larger than the m = −1 splitting; the colored regions indicate the ranges over which the function ζ is integrated for the components of the multiplet n = −141. The dashed lines indicate heighttobackground values of 7 and 10. 

In the text 
Fig. 14.
Histogram of the inclinations measured in NGC 6819. The dashed line indicates the sin i distribution. 

In the text 
Fig. 15.
Comparison of the inclinations of the spin axis of the stars in NGC 6819. Inclinations measured by Corsaro et al. (2017) are plotted on the xaxis, while inclinations derived from the asymptotic analysis are on the yaxis; 1σ uncertainties are indicated by vertical and horizontal error bars. The red symbol shows the case where no asymptotic fit could be found, and the green ones to cases without Bayesian fit. The dashed line corresponds to the 1:1 relation. 

In the text 
Fig. 16.
Fit of the mixed modes corresponding to n_{p} = 7 in KIC 5112373 (NGC 6819 member). The color codes the azimuthal order: m = +1 in purple, m = −1 in blue. The gray dashed lines indicate the two thresholds used in this work, corresponding to heighttobackground ratios of 7 and 10. Contrary to the analysis conducted by Corsaro et al. (2017), modes with m = ±1 are clearly identified. 

In the text 
Fig. 17.
Fit of the mixed modes corresponding to n_{p} = 9 in KIC 2437976 (NGC 6791 member). The color codes the azimuthal order: m = +1 in purple, m = −1 in blue. The location of m = 0 modes is indicated in light blue, but none shows a large height for this star seen equatoron. The gray dashed lines indicate the two thresholds used in this work, corresponding to heighttobackground ratios of 7 and 10. Many peaks above the threshold value 5.5 that rejects the null hypothesis at the 5%level follow the mixedmode pattern. 

In the text 
Fig. A.1.
Fit of the oscillation pattern of the low RGB star KIC 10272858, at the limit of validity of the asymptotic pattern. Owing to the small radial orders, small shifts are seen between observed and asymptotic spectra. Same style as Fig. 2, but ℓ = 3 modes appear near the abscissa 0.28. 

In the text 
Fig. A.2.
Fit of the oscillation pattern of the RGB star KIC 11353313. Same style as Fig. 2. 

In the text 
Fig. A.3.
Fit of the oscillation pattern of the RGB star KIC 3955033. The overlap of mixed modes with different mixedmode orders is the signature of the rapid core rotation. The second rotation crossing, where all components of the multiplets overlap (Gehan et al. 2017), occurs at the mixedorder n = −122. Same style as Fig. 2. 

In the text 
Fig. A.4.
Fit of the oscillation pattern of the evolved RGB star KIC 2443903, near the limit of capability of identification, with a large crowding due to the high mode density. The second rotation crossing, where all m components apparently coincide, occurs at n = −189 (with an abscissa ≃0.1 and n_{p} = 9); the third crossing, where m = 1 components apparently coincide with m = 0 inbetween, occurs at n = −233 (with an abscissa ≃0.25 and n_{p} = 7). Same style as Fig. 2. Note that the modes with large heights at an abscissa ≃ 0.2 are ℓ = 3 modes. 

In the text 
Fig. A.5.
Fit of the oscillation pattern of the redclump star KIC 1723700. Buoyancy glitches explain the small shifts between observed and asymptotic spectra but do no hamper the mode identification. Same style as Fig. 2. 

In the text 
Fig. A.6.
Fit of the oscillation pattern of the secondaryclump star KIC 1725190. Buoyancy glitches explain the small shifts between observed and asymptotic spectra but do no hamper the mode identification. Same style as Fig. 2. 

In the text 
Fig. B.1.
Fit of the oscillation pattern of the RGB star KIC 5024476, member of the open cluster NGC 6819. The dim magnitude of the cluster stars explains the low S/N. However, unambiguous doublets are identified all along the mixedmode spectrum; m = 0 modes are mostly absent and m = 1 modes dominate the mixedmode spectrum, so that a nearly poleon inclination is not possible. Same style as Fig. 2. 

In the text 
Fig. B.2.
Fit of the oscillation pattern of the RGB star KIC 2437976, member of the open cluster NGC 6791. The dim magnitude of the cluster stars explains the low S/N. The identification at radial order 9, supported by the radial orders 8 and 10, is unambiguously conclusive: m = 0 modes are mostly absent and m = 1 modes dominate the mixedmode spectrum, so that a poleon inclination is not possible. Same style as Fig. 2. 

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.