Angular momentum redistribution by mixed modes
in evolved lowmass stars
Key Words.:
Waves  Stars: oscillations  Stars: interiors  Stars: rotationThe detection of mixed modes in subgiants and red giants by the CoRoT and Kepler spaceborne missions allows us to investigate the internal structure of evolved lowmass stars, from the end of the main sequence to the central heliumburning phase. In particular, the measurement of the mean core rotation rate as a function of the evolution places stringent constraints on the physical mechanisms responsible for the angular momentum redistribution in stars. It showed that the current stellar evolution codes including the modelling of rotation fail to reproduce the observations. An additional physical process that efficiently extracts angular momentum from the core is thus necessary.
Our aim is to assess the ability of mixed modes to do this. To this end, we developed a formalism that provides a modelling of the wave fluxes in both the mean angular momentum and the mean energy equations in a companion paper. In this article, mode amplitudes are modelled based on recent asteroseismic observations, and a quantitative estimate of the angular momentum transfer is obtained. This is performed for a benchmark model of 1.3 at three evolutionary stages, representative of the evolved pulsating stars observed by CoRoT and Kepler.
We show that mixed modes extract angular momentum from the innermost regions of subgiants and red giants. However, this transport of angular momentum from the core is unlikely to counterbalance the effect of the core contraction in subgiants and early red giants. In contrast, for more evolved red giants, mixed modes are found efficient enough to balance and exceed the effect of the core contraction, in particular in the hydrogenburning shell. Our results thus indicate that mixed modes are a promising candidate to explain the observed spindown of the core of evolved red giants, but that an other mechanism is to be invoked for subgiants and early red giants.
1 Introduction
The CoRoT (Baglin et al. 2006a, b; Michel et al. 2008) and Kepler (Borucki et al. 2010) spaceborne missions provided a wealth of observed stars exhibiting solarlike oscillations, from the mainsequence to the redgiant phases (see Chaplin & Miglio 2013, for a review). A large number of those stars are lowmass evolved stars showing a rich spectrum of mixed modes, which behave as acoustic modes in the stellar envelope and as gravity modes in the core. They are therefore detectable at the stellar surface while yielding information on the innermost regions (e.g., Dziembowski 1971; Scuflaire 1974; Aizenman et al. 1977; Dziembowski et al. 2001; Dupret et al. 2009).
While mixed modes allowed unveiling the structure of the red giants core (Bedding et al. 2011; Mosser et al. 2011, 2014), the detection of associated rotational splittings enabled measuring the mean core rotation of subgiant and red giant stars (Beck et al. 2012; Deheuvels et al. 2012; Mosser et al. 2012b; Deheuvels et al. 2014). It turns out that the core of lowmass red giant stars slows down during evolution. Since local conservation of angular momentum would imply a spinup of the core, these observations indicate a strong transfer of angular momentum from the inner to the outer layers of these stars. However, current stellar evolutionary rotating models that take the transport of angular momentum into account fail to reproduce the observations. Meridional circulation or shear instabilities are not efficient enough to slow down the redgiant cores (e.g., Eggenberger et al. 2012; Marques et al. 2013; Ceillier et al. 2013), and recent attempts to include the effect of propagative internal gravity waves (Fuller et al. 2014) or magnetic fields using the TaylorSpruit dynamo formalism (Cantiello et al. 2014) failed to solve the problem. Nonetheless, for subgiants and early red giants, Rüdiger et al. (2015) have recently shown that magnetorotational instabilities of a toroidal magnetic field could explain the angular momentum redistribution in subgiants and early red giants (see also Maeder & Meynet 2014).
Hence, an additional physical mechanism that efficiently extracts angular momentum from the core of red giants seems to be needed. Nonradial modes in subgiants and red giants are potential candidates. Indeed, prograde and retrograde mixed modes are differentially damped in the presence of rotation, allowing a net transport of angular momentum (e.g., Ando 1986; Lee & Saio 1993). Our aim is thus to investigate the ability of mixed modes to transfer angular momentum from the inner radiative interior to the outer convective layers.
In the first article of this series (hereafter Paper I), we used the transformed Eulerian mean (TEM) formalism to account for the effect of wave driving through both the wave momentum flux in the mean angular momentum equation and the wave heat flux in the mean energy equation. The wave field was then modelled using the asymptotic, quasiadiabatic, and slowrotation approximations, valid in the dense radiative layers of evolved stars. This allowed us to obtain an explicit expression of the wave fluxes appearing in the mean angular momentum and energy equations. This article, the second of this series, is dedicated to a quantitative determination of the angular momentum transported by mixed modes, which requires the knowledge of mode amplitudes. Recent CoRoT and Kepler observations (Baudin et al. 2011; Mosser et al. 2012a; Samadi et al. 2012) allow us to determine mixedmode amplitudes and subsequently to determine the efficiency of angular momentum transport by mixed modes.
The paper is organized as follows: Sect. 2 recalls the formalism developed in Paper I, while Sect. 3 presents the modelling of the equilibrium models and related oscillations. In Sect. 4, based on the observations, we describe the modelling of mode amplitudes. Finally, the estimate of the angular momentum transport efficiency is presented in Sect. 5, and Sect. 6 is dedicated to discussions and conclusions.
2 Theoretical formalism
In this section, we recall the formalism developed in Paper I as well as the related assumptions for both the mean flow and the wave fluxes appearing in the mean equations.
2.1 Mean flow equations
Following Paper I, we adopt the transformed Eulerian mean formalism (TEM). In addition, we assume shellular rotation (see Maeder 2009, for an extensive discussion) so that at dominant order the mean angular momentum equation is
(1) 
where is the rotation angular frequency that depends only on the radius, is the density, the overbar and the symbol denotes the horizontal average (i.e. azimuthal and meridional). The fluxes are defined by
(2)  
(3)  
(4) 
where , the overbar is the azimuthal average, is the residual meridional velocity, is the radial eddy viscosity, is the specific entropy, and the prime denotes perturbations associated with the nonradial oscillations, so that , , are the azimuthal, radial, and meridional component of the wave velocity field, and the wave Eulerian perturbation of entropy. We have also introduced the Lagrangian derivative .
Under the same assumptions as for Eq. (1), the mean entropy equation is
(5) 
with
(6) 
and
(7) 
where is the temperature, is the nuclear energy generation rate and is the thermal conductivity (see, e.g., Sect. 6 of Mathis & Zahn 2004, and also Maeder & Zahn 1998). This expresses energy conservation on level surfaces.
To quantify the effect of the waves on the rotation profile, we further need to provide a modelling of the wave field^{1}^{1}1In this article, by wave we denote global mixed modes..
2.2 Wave flux modelling
Our aim is to determine the angular momentum redistribution induced by mixed modes in the radiative region of lowmass evolved stars. Therefore, the wave fluxes in Eqs. (1) and (5) must be modelled. This has been done in Paper I, so that for the momentum equation we have
(8)  
where the coefficients are given by Eqs. (A.25) to (A.28) of Paper I, and (with the modal frequency). The amplitude corresponds to the amplitude of a mode of a given angular degree, , and azimuthal order, . We note that the amplitude is considered statistically constant in time since we deal with solarlike oscillations that result from a balance between mode driving and damping (see Samadi 2011, for details).
For the entropy equation, we obtained
(9) 
where
(10) 
with the luminosity, the buoyancy frequency, and the actual and adiabatic temperature gradients, respectively, and the radial component of the eigenfunction for a mode of angular degree and azimuthal order . The overbars and angular brackets have been dropped for ease of notation.
We also recall that in the course of deriving Eqs. (8) and (9) several simplifying approximations were performed, namely:

The quasiadiabatic approach: It consists of neglecting the difference between adiabatic and nonadiabatic eigenfunctions in the full wave equations. This approximation is valid provided the local thermal timescale is much longer than the modal period. This is justified in the radiative region of evolved lowmass stars.
Consequently, our formalism is valid in the dense radiative core of evolved slowly rotating stars.
3 Computing selected stellar models
We have computed an evolutionary sequence of a benchmark model of a 1.3 M evolved lowmass star representative of CoRoT and Kepler observations. We selected three models at three postmain sequence stages of this evolutionary sequence, as discussed below. The location of the models in the HertzsprungRussell (HR) diagram is shown in Fig. 1, and their characteristics are summarised in Table 1.
Model  Age (Myr)  (Hz)  (Hz)  

M0  4401  2.47  4.95  39.5  666 
M1  4590  3.55  6.78  22.9  339 
M2  4791  9.51  38.3  5.24  48.8 
3.1 Equilibrium models and related eigenfunctions and eigenfrequencies
We used the stellar evolution code CESTAM (Marques et al. 2013) to compute the equilibrium models. The atmosphere was computed assuming a grey Eddington approximation. Convection was included according to Canuto et al. (1996), with a mixinglength parameter . The initial chemical composition follows Asplund et al. (2005), with a helium mass fraction of and a metallicity of . We used the OPAL equation of state (Rogers et al. 1996) and opacities (Iglesias & Rogers 1996), complemented, at K, by the Alexander & Ferguson (1994) opacities. We adopted the NACRE nuclear reaction rates from Angulo et al. (1999) except for the reaction, where we used the reaction rates given in Imbriani et al. (2004).
3.2 Rotation profiles
We consider that rotation has negligible effects on the equilibrium structure of the star and that it can be considered as a perturbation for the oscillations. This approximation is consistent with the formalism developed in Paper I.
The rotation profiles are therefore determined by computing synthetic profiles that match the observations. The synthetic profiles should have a core that rotates faster than the envelope, with a transition region located at the Hburning shell in the case of red giant stars (models M1 and M2). The width of the transition region is related to the width of the Hburning shell. We consequently used the following profile
(11) 
where and are the surface and central angular velocities, respectively, and and are the radius and width of the burning shell. was chosen to reproduce the mean core rotation rate of evolved lowmass stars as derived from Kepler observations (see Mosser et al. 2012b; Deheuvels et al. 2014), and the surface rotation rate () was obtained so as to conserve the total angular momentum from the subgiant stage (model M0).
For the subgiant model (model M0), we also used Eq. (11), but with both and given by the radius of the central convective zone at the end of the main sequence (just before it disappears). In this case, the central and surface rotation ( and , respectively) were chosen to reproduce the results of Deheuvels et al. (2012). The values of the parameters for models M0, M1, and M2 are summarised in Table 2, and the resulting rotation profiles are displayed in Fig. 2.
4 Modelling mode amplitudes
There are mainly two approaches to compute mode amplitudes. The first is based on a full nonadiabatic computation including a timedependent treatment of convection. This procedure is timeconsuming, however, and still suffers from theoretical uncertainties (see Dupret et al. 2009; Grosjean et al. 2014, for details). The second is based on recent CoRoT and Kepler observations that allowed us to measure the amplitudes of solarlike modes in a wide variety of stars in different evolutionary stages. This wealth of observations made it possible to establish scaling relations that provide mode amplitudes versus global stellar parameters (e.g., Mosser et al. 2012a; Samadi et al. 2012). In this paper, we adopt the latter approach.
4.1 Radial mode amplitudes
We first express the amplitudes of radial modes. The surface velocity of the radial modes, at the frequency of maximum height () in the oscillation power spectrum, follows the scaling relation (see Belkacem & Samadi 2013, for details)
(12) 
where is the effective temperature, is the large separation, and the symbols stand for the solar values. Here we adopt the values K, Hz, and Hz.
To account for the frequency dependence of radial mode amplitudes, we follow the work of Mosser et al. (2010) , who used a Gaussian envelope centred on with a fullwidth at halfmaximum , that is,
(13) 
with . Based on Kepler observations, we take Hz for M0, Hz for M1, and Hz for M2 for the following.
4.2 Nonradial mode amplitudes
The next step is to consider the amplitudes of nonradial modes. The ratio between radial and nonradial mode amplitudes of comparable frequencies can be written as (see Benomar et al. 2014, for details)
(14) 
where is the mode mass, is the eigenfunction, is the mass profile, is the total radius, and is the mode linewidth that is related to the mode damping rate by .
For high angular degrees, the radiative work of the oscillations in the innermost regions dominantes the work done in the uppermost layers (see Grosjean et al. 2014, for details). It follows that the mode amplitudes of high angular degrees become very small due to strong radiative damping. We thus define a cutoff angular degree, , below which we assume that the work of nonradial modes equals the work of radial modes
(15) 
where the product between the mode mass and the mode damping rate equals the work integral (e.g., Dupret et al. 2009; Benomar et al. 2014; Grosjean et al. 2014). Now, using Eq. (15) together with Eq. (14) leads to
(16) 
For , we assume that the radiative damping in the innermost region becomes dominant over the damping of the upper layers so that mode amplitudes become negligible.
The cutoff angular degree is obtained by assuming that the work done over an oscillation period through the effect of radiative losses in the innermost layers becomes of the same order of magnitude as the work done in the uppermost layers. We obtain for M0, for M1, and for M2 (see Appendix A for details). Note that this approach for the modelling of mode amplitudes provides a lower limit to the amount of angular momentum transported by mixed modes since we neglect the contribution of angular degrees higher than . However, this is enough to determine whether mixed modes are able to slow down the rotating core of red giants.
5 Quantitative estimate of angular momentum extracted by mixed modes
Our aim is to estimate the rate of angular momentum transported by mixed modes for the benchmark models M0, M1, and M2. To this end, as in Paper I, we introduce the notation
(17) 
where the righthand side is given by Eq. (8). Using the approach presented in Sect. 4 to derive the amplitudes, we compute the wave fluxes in the righthand side of Eqs. (1) and (5) using Eqs. (8) and (9).
The contribution of the mixed modes to the mean angular momentum equation is shown in Fig. 4 for models M0, M1, and M2. The main effect of mixed modes is to slow down the core of the stars in the region near the maximum of the buoyancy frequency that also corresponds to the edge of the fast rotating core. For the redgiant models, this region is the hydrogenburning shell. The decrease shown towards the upper radiative layers is the result of both a lower rotation rate and a lower value of the buoyancy frequency. This behaviour is similar in the three models. In addition, the amount of angular momentum extracted by mixed modes hardly depends on the rotation gradient between the core and the envelope (the sensitivity of the results to the parameter in Eq. (11) remains weak).
The amount of angular momentum extracted by mixed modes increases with the evolutionary stage of the star. This effect is the result of several factors. First, as shown in Fig. 3, mode amplitudes increase from models M0 to M2 and thus more energy is available to transport angular momentum. Moreover, the buoyancy frequency also increases and so does the radial wave number in the gravitymode cavity. As the energy exchanges between modes and the background is proportional to the radial wave number, the amount of angular momentum extracted also increases. Finally, the number of mixed modes between two pdominated modes significantly increases between M0 and M2.
Furthermore, the angular momentum is mainly carried away by gdominated modes of the highest angular degrees. Indeed, these modes are very efficiently trapped in the core and undergo a larger radiative damping (due to the larger radial wave number). The predominance of high angular degree modes implies that is sensitive to our determination of . As mentioned before, the estimate of the amount of angular momentum extracted in the radiative regions is thus a lower limit. In other words, adding the contribution of angular degrees higher than would strengthen these trends.
To proceed and assess the efficiency of the transport of angular momentum for our benchmark models, we compare two timescales, namely:

the timescale associated with the efficiency of the transport of angular momentum by mixed modes, defined as
(18) 
the timescale associated with the contraction of the star, defined as
(19) where .
The two timescales are shown in Fig. 5 (top panel). In all the radiative regions of models M0 and M1, the extraction of angular momentum remains negligible compared to the star contraction. In contrast, for model M2, the timescale of angular momentum extraction is of the same order of magnitude as the contraction timescale and even lower in the hydrogenshellburning region, for which the contraction is maximal (see Fig. 5, bottom panel). However, in the very centre or in the upper radiative region, the mixed mode extraction of angular momentum is again negligible.
Therefore, we conclude that the extraction of angular momentum by mixed modes is negligible in subgiants and early red giants, but becomes strong in the hydrogenburning shell in stars higher on the redgiant branch. In these cases, mixed modes are able to counterbalance the spinup due to the star contraction and can thus enforce a spindown in those layers. We note, however, that the exact location in the HR diagram of the transition between inefficient and efficient extraction by mixed modes is likely to depend on the mass and internal physics of the models.
6 Concluding remarks
Based on the theoretical formalism developed in paper I, we computed the amount of angular momentum extracted by mixed modes in evolved lowmass stars. To this end, we considered benchmark models of a subgiant (model M0), an early red giant (model M1), and a more evolved red giant (model M2) for which we computed realistic mode amplitudes based on recent CoRoT and Kepler observations. We found that mixed modes extract angular momentum in the innermost region of subgiants and red giants. This extraction of angular momentum is found negligible for the subgiant and lower RGB models, whereas it is efficient enough to counterbalance the effect of the stellar contraction for the most evolved model.
For this evolved star, mixed modes efficiently transport the angular momentum in the hydrogenburning shell (that corresponds to the maximum of the buoyancy frequency) but becomes less efficient in the very central layers. We note, however, that mixed modes can be considered as an external torque for meridional circulation (e.g., Mathis et al. 2013). Therefore, such a localized spindown in the hydrogenburning shell is likely to feed a strong redistribution of angular momentum in the central layers by meridional circulation. This would nevertheless deserve a quantitative estimate.
Indeed, evolutionary calculations including the effect of mixed modes coupled with meridional circulation and shear instabilities are the next step for a definitive conclusion on the ability of mixed modes to slow down the cores of evolved red giants. Modelling of mode amplitudes would require further investigation since we assumed that the contribution of highangular degree modes would be negligible due to their high radiative damping. Therefore, we believe that our calculations provide a lower limit. Full nonadiabatic computations would provide more accurate results, however.
We conclude that, for subgiants and early red giants, our results indicate that extraction of angular momentum by mixed modes cannot explain the current observations (Deheuvels et al. 2014). Other physical mechanisms such as internal gravity waves (Talon & Charbonnel 2008; Fuller et al. 2014) or magnetic fields (Rüdiger et al. 2015) should also be present in these stars in order to extract enough angular momentum from the core. In contrast, as demonstrated in this work, the transport of angular momentum by mixed modes is a serious explanation for the observed spin down of the core of evolved stars on the redgiant branch. At least, given the order of magnitude of the induced angular momentum extraction, it seems difficult to neglect this mechanism in the future.
Finally, we note that an extension of this work to central helium burning stars would be desirable in the future. Indeed, even if the mode amplitudes are similar for stars on the vertical branch and on the clump for the same (e.g., Mosser et al. 2012a), the maximum of buoyancy frequency is lower and shifted toward higher radius for clump stars (e.g., Montalbán et al. 2013). Therefore, we can expect a lower efficiency of the transport by mixed modes but this needs to be definitively quantified.
Acknowledgements.
We acknowledge the ANR (Agence Nationale de la Recherche, France) program IDEE (nANR12BS050008) “Interaction Des Étoiles et des Exoplanètes” as well as financial support from ”Programme National de Physique Stellaire” (PNPS) of CNRS/INSU, France. RMO aknowledges funding for the Stellar Astrophysics Centre, provided by The Danish National Research Foundation, and funding for the ASTERISK project (ASTERoseismic Investigations with SONG and Kepler) provided by the European Research Council (Grant agreement no.: 267864).References
 Aizenman et al. (1977) Aizenman, M., Smeyers, P., & Weigert, A. 1977, A&A, 58, 41
 Alexander & Ferguson (1994) Alexander, D. R. & Ferguson, J. W. 1994, ApJ, 437, 879
 Ando (1986) Ando, H. 1986, A&A, 163, 97
 Angulo et al. (1999) Angulo, C., Arnould, M., Rayet, M., et al. 1999, Nuclear Physics A, 656, 3
 Appourchaux et al. (2012) Appourchaux, T., Benomar, O., Gruberbauer, M., et al. 2012, A&A, 537, A134
 Asplund et al. (2005) Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Astronomical Society of the Pacific Conference Series, Vol. 336, Cosmic Abundances as Records of Stellar Evolution and Nucleosynthesis, ed. T. G. Barnes III & F. N. Bash, 25–+
 Baglin et al. (2006a) Baglin, A., Auvergne, M., Barge, P., et al. 2006a, in ESA Special Publication, Vol. 1306, ESA Special Publication, ed. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, 33
 Baglin et al. (2006b) Baglin, A., Auvergne, M., Boisnard, L., et al. 2006b, in COSPAR Meeting, Vol. 36, 36th COSPAR Scientific Assembly, 3749
 Baudin et al. (2011) Baudin, F., Barban, C., Belkacem, K., et al. 2011, A&A, 535, C1
 Beck et al. (2012) Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55
 Bedding et al. (2011) Bedding, T. R., Mosser, B., Huber, D., et al. 2011, Nature, 471, 608
 Belkacem et al. (2012) Belkacem, K., Dupret, M. A., Baudin, F., et al. 2012, A&A, 540, L7
 Belkacem & Samadi (2013) Belkacem, K. & Samadi, R. 2013, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 865, Lecture Notes in Physics, Berlin Springer Verlag, ed. M. Goupil, K. Belkacem, C. Neiner, F. Lignières, & J. J. Green, 179
 Benomar et al. (2014) Benomar, O., Belkacem, K., Bedding, T. R., et al. 2014, ApJ, 781, L29
 Borucki et al. (2010) Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977
 Cantiello et al. (2014) Cantiello, M., Mankovich, C., Bildsten, L., ChristensenDalsgaard, J., & Paxton, B. 2014, ApJ, 788, 93
 Canuto et al. (1996) Canuto, V. M., Goldman, I., & Mazzitelli, I. 1996, ApJ, 473, 550
 Ceillier et al. (2013) Ceillier, T., Eggenberger, P., García, R. A., & Mathis, S. 2013, A&A, 555, A54
 Chaplin & Miglio (2013) Chaplin, W. J. & Miglio, A. 2013, ARA&A, 51, 353
 ChristensenDalsgaard (2008) ChristensenDalsgaard, J. 2008, Ap&SS, 316, 113
 ChristensenDalsgaard (2011) ChristensenDalsgaard, J. 2011, ADIPLS: Aarhus Adiabatic Oscillation Package (ADIPACK), astrophysics Source Code Library
 Deheuvels et al. (2014) Deheuvels, S., Doğan, G., Goupil, M. J., et al. 2014, A&A, 564, A27
 Deheuvels et al. (2012) Deheuvels, S., García, R. A., Chaplin, W. J., et al. 2012, ApJ, 756, 19
 Dupret et al. (2009) Dupret, M.A., Belkacem, K., Samadi, R., et al. 2009, A&A, 506, 57
 Dziembowski (1971) Dziembowski, W. A. 1971, Acta Astron., 21, 289
 Dziembowski et al. (2001) Dziembowski, W. A., Gough, D. O., Houdek, G., & Sienkiewicz, R. 2001, MNRAS, 328, 601
 Eggenberger et al. (2012) Eggenberger, P., Montalbán, J., & Miglio, A. 2012, A&A, 544, L4
 Fuller et al. (2014) Fuller, J., Lecoanet, D., Cantiello, M., & Brown, B. 2014, ApJ, 796, 17
 Godart et al. (2009) Godart, M., Noels, A., Dupret, M.A., & Lebreton, Y. 2009, MNRAS, 396, 1833
 Goupil et al. (2013) Goupil, M. J., Mosser, B., Marques, J. P., et al. 2013, A&A, 549, A75
 Grosjean et al. (2014) Grosjean, M., Dupret, M.A., Belkacem, K., et al. 2014, A&A, 572, A11
 Iglesias & Rogers (1996) Iglesias, C. A. & Rogers, F. J. 1996, ApJ, 464, 943
 Imbriani et al. (2004) Imbriani, G., Costantini, H., Formicola, A., et al. 2004, A&A, 420, 625
 Lee & Saio (1993) Lee, U. & Saio, H. 1993, MNRAS, 261, 415
 Maeder (2009) Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (Astronomy and Astrophysics Library, Volume . ISBN 9783540769484. Springer Berlin Heidelberg, 2009)
 Maeder & Meynet (2014) Maeder, A. & Meynet, G. 2014, ApJ, 793, 123
 Maeder & Zahn (1998) Maeder, A. & Zahn, J.P. 1998, A&A, 334, 1000
 Marques et al. (2013) Marques, J. P., Goupil, M. J., Lebreton, Y., et al. 2013, A&A, 549, A74
 Mathis et al. (2013) Mathis, S., Decressin, T., Eggenberger, P., & Charbonnel, C. 2013, A&A, 558, A11
 Mathis & Zahn (2004) Mathis, S. & Zahn, J.P. 2004, A&A, 425, 229
 Michel et al. (2008) Michel, E., Baglin, A., Auvergne, M., et al. 2008, Science, 322, 558
 Montalbán et al. (2013) Montalbán, J., Miglio, A., Noels, A., et al. 2013, ApJ, 766, 118
 Mosser et al. (2011) Mosser, B., Barban, C., Montalbán, J., et al. 2011, A&A, 532, A86
 Mosser et al. (2010) Mosser, B., Belkacem, K., Goupil, M., et al. 2010, A&A, 517, A22
 Mosser et al. (2014) Mosser, B., Benomar, O., Belkacem, K., et al. 2014, A&A, 572, L5
 Mosser et al. (2012a) Mosser, B., Elsworth, Y., Hekker, S., et al. 2012a, A&A, 537, A30
 Mosser et al. (2012b) Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012b, A&A, 548, A10
 Rogers et al. (1996) Rogers, F. J., Swenson, F. J., & Iglesias, C. A. 1996, ApJ, 456, 902
 Rüdiger et al. (2015) Rüdiger, G., Gellert, M., Spada, F., & Tereshin, I. 2015, A&A, 573, A80
 Samadi (2011) Samadi, R. 2011, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 832, Lecture Notes in Physics, Berlin Springer Verlag, ed. J.P. Rozelot & C. Neiner, 305
 Samadi et al. (2012) Samadi, R., Belkacem, K., Dupret, M.A., et al. 2012, A&A, 543, A120
 Scuflaire (1974) Scuflaire, R. 1974, A&A, 36, 107
 Talon & Charbonnel (2008) Talon, S. & Charbonnel, C. 2008, A&A, 482, 597
Appendix A Cutoff angular degrees
The cutoff angular degree, as introduced in Sect. 4.2, can be defined as the angular degree for which the work performed over an oscillation period through the effect of radiative losses in the innermost layers becomes of the same order of magnitude as for the work performed in the uppermost layers. More precisely, is obtained at through the relation
(20) 
The lefthand side of Eq. (20) can be identified with a radial mode of comparable frequency for which the mode mass is computed as described in Sect. 3 and the associated damping rate is obtained using the scaling relation as provided by Belkacem et al. (2012) (see also Baudin et al. 2011; Appourchaux et al. 2012). It gives Hz for M0, Hz for M1, and Hz for M2.
The righthand side of Eq. (20) can be identified with a gravity dominated mode for which the mode mass is computed following Sect. 3 and the damping rate is obtained using the asymptotic relation (e.g. Dziembowski et al. 2001; Godart et al. 2009)
(21) 
where is the pressure, is the luminosity, is the base of the convective region, , is the temperature gradient and its adiabatic counterpart.