Neardegeneracy effects on the frequencies of rotationallysplit mixed modes in red giants
^{1} Université de Toulouse; UPSOMP; IRAP, 31013 Toulouse, France
email: sebastien.deheuvels@irap.omp.eu
^{2} CNRS; IRAP; 14 avenue Edouard Belin, 31400 Toulouse, France
^{3} Institut d’Astrophysique et de Géophysique de l’Université de Liège, Allée du 6 Août 17, 4000 Liège, Belgium
^{4} Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
^{5} Observatoire de Paris, LESIA, CNRS UMR 8109, 92195 Meudon, France
^{6} Department of Astronomy, Yale University, PO Box 208101, New Haven, CT 065208101, USA
Received: 14 March 2017
Accepted: 24 May 2017
Context. The Kepler space mission has made it possible to measure the rotational splittings of mixed modes in red giants, thereby providing an unprecedented opportunity to probe the internal rotation of these stars.
Aims. Asymmetries have been detected in the rotational multiplets of several red giants. This is unexpected since all the red giants whose rotation profiles have been measured thus far are found to rotate slowly, and low rotation, in principle, produces symmetrical multiplets. Our aim here is to explain these asymmetries and find a way of exploiting them to probe the internal rotation of red giants.
Methods. We show that in the cases where asymmetrical multiplets were detected, neardegeneracy effects are expected to occur, because of the combined effects of rotation and mode mixing. Such effects have not been taken into account so far. By using both perturbative and nonperturbative approaches, we show that neardegeneracy effects produce multiplet asymmetries that are very similar to the observations. We then propose and validate a method based on the perturbative approach to probe the internal rotation of red giants using multiplet asymmetries.
Results. We successfully apply our method to the asymmetrical l = 2 multiplets of the Kepler young red giant KIC 7341231 and obtain precise estimates of its mean rotation in the core and the envelope. The observed asymmetries are reproduced with a good statistical agreement, which confirms that neardegeneracy effects are very likely the cause of the detected multiplet asymmetries.
Conclusions. We expect neardegeneracy effects to be important for l = 2 mixed modes all along the red giant branch (RGB). For l = 1 modes, these effects can be neglected only at the base of the RGB. They must therefore be taken into account when interpreting rotational splittings and as shown here, they can bring valuable information about the internal rotation of red giants.
Key words: stars: oscillations / stars: rotation / stars: evolution / stars: individual: KIC 7341231
© ESO, 2017
1. Introduction
Rotation is known to have a large influence on the structure and the evolution of stars. In particular, it is expected to induce additional mixing of the chemical elements inside stars. Taking the effects of rotation into account is believed to be a key step to making progress in stellar modeling. However, the internal rotation profiles of stars and the way they are modified during their evolution remains uncertain. Rotationinduced processes of angular moment transport such as meridional circulation and shear instabilities (as they are currently modeled) are not efficient enough to account for the nearly solidbody rotation in the solar radiative interior that was found by helioseismologic analyses (Schou et al. 1998; Chaplin et al. 1999; EffDarwich & Korzennik 2013). Other types of processes are likely to operate. The main candidates invoked so far are internal gravity waves (Charbonnel & Talon 2005; Talon & Charbonnel 2008; Fuller et al. 2014) and magnetic fields that either have a fossil origin (e.g., Gough & McIntyre 1998; Spada et al. 2010; Maeder & Meynet 2014) or result from magnetohydrodynamic (MHD) instabilities (e.g. Spruit 1999; Rüdiger et al. 2015). However, the relative importance of these processes and the timescales over which they operate is still a matter of active debate.
Recent progresses in the seismology of subgiants and red giants have proved that it could significantly contribute to understanding angular momentum transport inside stars. Owing to the large core density of red giants, the frequencies of gravity modes in these stars are comparable to the frequencies of the acoustic modes stochastically excited in the convective envelope. This gives rise to mixed modes, which behave as gravity modes in the core and as pressure modes in the envelope. Such modes were initially found in stellar models by Dziembowski (1971), Scuflaire (1974), and Osaki (1975). They were first detected from the ground (Kjeldsen et al. 1995). More recently, with the advent of space missions CoRoT (Baglin et al. 2006) and Kepler (Borucki et al. 2010), mixed modes were detected in thousands of subgiants (Deheuvels et al. 2010; Campante et al. 2011) and red giants (Bedding et al. 2011; Mosser et al. 2011). The measurement of rotational splittings of mixed modes (Beck et al. 2012; Deheuvels et al. 2012; Mosser et al. 2012) has given direct insight into the internal rotation of red giants and their variations with evolution. In the subgiant phase, Deheuvels et al. (2014) showed that the core spins up and the envelope spins down. This was qualitatively expected, considering that the deepest layers below the Hburning shell contract, while the layers above expand. However, the core rotation rates measured with seismology are several orders of magnitude below the values predicted by theoretical models that include rotationinduced transport of angular momentum (Eggenberger et al. 2012; Marques et al. 2013). On the red giant branch (RGB), seismic measurements of the core rotation for several hundreds of Kepler targets led to the striking observation that the core of red giants in fact spins down (Mosser et al. 2012), which is also at odds with current theoretical models. These results are clear evidence that an efficient redistribution of angular momentum between the core and the envelope occurs in red giants, the origin of which remains unknown. More recently it has been shown that a similar conclusion can be drawn for intermediatemass stars in the coreHe burning phase, where an even milder radial differential rotation was observed, with coreenvelope ratios around two (Deheuvels et al. 2015). These novel observational constraints prompted new theoretical studies about the efficiency of angular momentum transport caused by internal gravity waves (Fuller et al. 2014), MHD instabilities (Rüdiger et al. 2015; Jouve et al. 2015), or mixed modes themselves (Belkacem et al. 2015b,a), but the question remains open.
So far, all the red giants whose rotation profile could be seismically measured have been found to be slow rotators; they rotate in a regime where the effects of the centrifugal force on stellar pulsations can be safely neglected and the rotation frequency remains small compared to the pulsation frequencies. In such cases, the effects of rotation can be treated as a firstorder perturbation to the equations of nonradial oscillations. For a sphericallysymmetric rotation profile Ω(r), the frequency shift for a mode of radial order n, degree l, and azimuthal order m can be expressed as (1)where the functions K_{n,l}(r) are the rotational kernels. According to Eq. (1), the components of a rotational multiplet are expected to be uniformly spaced, and thus symmetrical with respect to the central m = 0 component. While most rotational multiplets in red giants indeed show this symmetry, asymmetries have been reported in several Kepler red giants. Deheuvels et al. (2012) found significant asymmetries in the rotational multiplets of two neighboring l = 2 modes of the young red giant KIC 7341231. Asymmetries were also reported in the multiplets of several l = 1 modes in the spectra of the red giants KIC 5006817 (Beck et al. 2014) and KIC 4448777 (Di Mauro et al. 2016). These asymmetries need to be understood because they are at odds with what is expected from Eq. (1), which has been assumed by all the studies that have derived information about the internal rotation of red giants.
In this paper, we argue that for KIC 7341231 (and potentially also for the two other aforementioned targets), these asymmetries are caused by neardegeneracy effects, which occur because of the combined effects of rotation and mode mixing. So far, it has always been assumed that the effects of rotation on the frequencies of mixed modes could be estimated in the same manner as for regular modes. However, mode mixing occurs in red giants when p modes and g modes have frequencies close enough for them to couple through the evanescent region separating the two cavities. It is well known that when two modes have close frequencies, neardegeneracy effects may arise which modify the frequency corrections to be applied (e.g., Dziembowski & Goode 1992; Suárez et al. 2006). If the frequency spacing between mixed modes is of the order or smaller than the rotation rate, we expect neardegeneracy effects to become important. This issue has not been addressed so far for red giants. In this paper, we study this issue using KIC 7341231 as a test case. We show that neardegeneracy effects need to be taken into account in this star and we demonstrate that they produce multiplet asymmetries that can account for the observations.
The paper is organized as follows. In Sect. 2, we confirm and update the measurement of multiplet asymmetries in the oscillation spectrum of KIC 7341231 using the full Kepler dataset and we demonstrate that neardegeneracy effects cannot be neglected. In Sect. 3, we include the effects of neardegeneracy in the firstorder perturbative treatment of rotation in the oscillation equations. We show that this correction indeed produces asymmetrical multiplets that are in qualitative agreement with the observations, and we argue that the intensity of the asymmetries depends on the level of radial differential rotation and on the trapping of the modes. In Sect. 4, we use the nonperturbative oscillation code acor (Ouazzani et al. 2012) to validate the firstorder perturbative approach. In Sect. 5, we propose a new method to measure the average core and envelope rotation rates using the asymmetries in l = 2 rotational multiplets. This method is thoroughly tested and validated using simulated data, and then successfully applied to KIC 7341231 in Sect. 6. Section 7 is dedicated to a discussion of the results and we conclude in Sect. 8.
2. Asymmetry of rotational multiplets in Kepler data: the testcase of KIC 7341231
The star KIC 7341231 is a young red giant, which was observed with the Kepler spacecraft in short cadence (58.84876 s) from quarter Q5 through quarter Q17 of the mission. This represents 3.14 yr of nearly continuous data (duty cycle of 91%). Based on the data from quarters Q5 through Q8 (about one year of data), Deheuvels et al. (2012) could measure the rotational splittings of 17 l = 1 modes and 2 l = 2 modes. By using the Optimally Localized Averages (OLA) inversion method, they obtained a precise estimate of the mean core rotation (⟨ Ω ⟩ _{c}/ 2π = 710 ± 51 nHz). The authors showed that the measured splittings, however, could not be used to build an average kernel for the envelope that would efficiently suppress the core contribution, and they were thus only able to derive an upper limit for the envelope mean rotation rate (⟨ Ω ⟩ _{e}/ 2π< 150 ± 19 nHz). This limitation is inherent to the eigenfunctions of the detected modes, and the increased precision to which the rotational splittings could be measured with the complete Kepler dataset would not solve this problem. It might, however, be that the rotational splittings of additional modes have become measurable, which might help in canceling the core contribution in the envelope average kernel. This is out of the scope of the present study.
Fig. 1 Section of the oscillation spectrum of KIC 7341231 in the neighborhood of the 384μHz radial mode (light gray area). The dark gray area corresponds to a rotationallysplit l = 1 mode and the vertical arrows indicate the observed mode frequencies of two l = 2 mixed modes (see Table 1). The vertical dotted lines show the results of a fit of Eq. (24) to the observed mode frequencies. 

Open with DEXTER 
Extracted parameters for the components of two l = 2 multiplets in avoided crossing in the spectrum of KIC 7341231.
Deheuvels et al. (2012) reported clear asymmetries for two closelyspaced l = 2 mixed modes in the spectrum of KIC 7341231 using the oneyearlong observation available at that time. For this reason, these two modes were at that time excluded from the dataset used for rotation inversions. This is unsatisfactory since these modes clearly carry information on the internal rotation of the star. Additionally, since the coupling between the pmode and gmode cavities is weaker for l = 2 modes, these modes could provide a more precise estimate of the envelope rotation. This study is aimed at understanding these asymmetries in rotational multiplets and we used KIC 7341231 as a testcase.
We updated the measurements of multiplet asymmetries in this star using the fulllength Kepler data (quarters Q5 through Q17). Figure 1 shows a section of the power spectrum of the star in the neighborhood of the two l = 2 modes. The m = (−2,0,2) components are the most prominent because the inclination angle of the star is close to 90°; only these components were detected by Deheuvels et al. (2012) with the observation available at that time. When using the full Kepler data, we find that the m = + 1 components of both l = 2 multiplets also stand out significantly above the noise. We extracted the mode parameters (frequency, height, width) of all the components of the multiplets using the fitting technique described in Deheuvels et al. (2012). The results are summarized in Table 1. The detected components of the two multiplets, indicated by arrows in Fig. 1, show clear asymmetries.
The asymmetry of the multiplets can be conveniently quantified by computing the quantity (2)With this definition, δ_{asym} = 0 for a symmetric multiplet, while  δ_{asym}  = 1 if either ω_{− m} or ω_{m} overlaps with the ω_{0} mode. A positive value for δ_{asym} means that the splittings of the m> 0 components are larger than those of the m< 0 components (and conversely for δ_{asym}< 0). For the two l = 2 mixed modes of KIC 7341231, we found δ_{asym} ~ −0.24 ± 0.03 (red mode in Fig. 1) and δ_{asym} ~ 0.14 ± 0.01 (blue mode in Fig. 1) using the m = ± 2 components.
The two multiplets that show asymmetries in the spectrum of KIC 7341231 have close frequencies, which means that neardegeneracy effects may arise, as mentioned in Sect. 1. Using Table 1, we find that the frequency spacing between the m = 0 components of the two multiplets is 2.24 ± 0.02 μHz. The core rotation rate of KIC 7341231 as obtained by Deheuvels et al. (2012) thus amounts to ~32% of the frequency spacing between the two modes, which is not negligible. Consequently, we expect neardegeneracy effects to arise and they need to be taken into account. In the following, we show that neardegeneracy effects can indeed produce multiplet asymmetries very similar to the observed ones. For the sake of completeness, in Sect. 7 we mention other possible causes for multiplet asymmetries, which we find less likely for KIC 734123.
3. Firstorder perturbative approach including neardegeneracy
We tested the influence of neardegeneracy effects on rotational multiplets by including them in a firstorder perturbative treatment of rotation. This has already been studied in detail by previous authors (e.g., Dziembowski & Goode 1992; Suárez et al. 2006; Ouazzani & Goupil 2012; see also Goupil 2009, for a review) so we briefly recall the procedure and refer to these authors for more details.
3.1. Nondegenerate frequencies
The oscillation equations including the firstorder effects of rotation can be written as (3)where ω corresponds to the mode eigenfrequency and ξ is the associated eigenfunction. The firstorder correction ℒ_{1} to the nonrotating operator ℒ_{0} includes the effect of the advection relative to the inertial frame and the contribution from the Coriolis force. The two operators are given by For simplicity, in the following the mode eigenfunctions are normalized by the mode inertias, so that ⟨ ξ_{a}  ξ_{a} ⟩ = 1, where the inner product is defined as (6)and L^{2} = l(l + 1).
In the general case without neardegeneracy, the firstorder correction ω_{1} to the mode eigenfrequency is given by the variational principle, and therefore (7)If we further assume a spherically symmetric rotation profile Ω(r), Eq. (7) reduces to (8)where the rotational kernel K(r) is given by (9)Equation (8) shows that without neardegeneracy, the splittings vary linearly with m, so that the components of a rotational multiplet are expected to be uniformly spaced, and thus symmetrical with respect to the central m = 0 component.
3.2. Neardegenerate frequencies
We now consider the case of two modes with same degree l, same azimuthal order m, and neardegenerate frequencies, that is,  ω_{0,a}−ω_{0,b}  ~ Ω, where the subscripts a and b refer to the two modes. In this case, the modes are coupled and corrections need to be included in the perturbative calculation of mode frequencies. To firstorder, the eigenfunctions of neardegenerate modes can be written as (10)As described in Appendix A, the perturbed eigenfrequencies are then given by the following expression (11)where ω_{a} and ω_{b} correspond to the firstorder perturbed frequencies of modes a and b when neardegeneracy effects are ignored; for example, for mode a(12)The quantity ω_{1,ab} represents the coupling between the two modes. For a spherically symmetric rotation profile, we have (13)where (14)Equation (11) shows that if  ω_{a}−ω_{b}  ≫ ω_{1,ab}, then we obtain ω_{±} ∈ {ω_{a},ω_{b}}, that is, we recover the eigenfrequencies of the nondegenerate case.
As shown in Appendix A, we can also obtain the values of A and B in Eq. (10) for each component of both multiplets. We can thus estimate the perturbations to the mode eigenfunctions caused by neardegeneracy effects (for nondegenerate modes, changes in the eigenfunctions perturb the mode frequencies only at second order).
3.3. Test on stellar models
We tested the effects of neardegeneracy on the rotational multiplets using stellar models. For this purpose, we used the bestfit stellar model obtained for KIC 7341231 by Deheuvels et al. (2012; model B). This model was computed with the evolutionary code cesam2k (Morel 1997) and optimized to match the observed frequencies of the l = 1 mixed modes for this star. We refer the reader to Deheuvels et al. (2012) for information about the input physics of this model and the optimization procedure that was adopted. One interesting feature of this model is that it also roughly reproduces the frequencies of the two l = 2 mixed modes that are undergoing an avoided crossing between 380 and 384 μHz (even though these modes were not included in the fitting procedure followed by Deheuvels et al. 2012). This model can therefore be used to investigate the origin of the observed asymmetries in these multiplets.
We also had to assume a rotation profile for the star. For flexibility reasons, we chose an analytic rotation profile of the type (15)where Ω_{c} and Ω_{e} correspond to the core and surface rotation, respectively. We here took Ω_{c}/ 2π = 710 nHz and Ω_{e}/ 2π = 150 nHz (i.e., a coreenvelope contrast of ~ 5), as per the results obtained by Deheuvels et al. (2012) for KIC 7341231. The influence of the level of radial differential rotation is studied in Sect. 3.5. We assumed a smooth transition between the core and surface rotation rates (d = 0.06) located inside the evanescent region (r_{0} = 0.2).
Fig. 2 Variations in the frequencies of two l = 2 rotational multiplets during an avoided crossing. The sequence of models B0 through B8 was used, assuming the rotation profile given by Eq. (15) with a coreenvelope contrast of five. The frequencies of rotational multiplets were computed using the firstorder pertubative approach including neardegeneracy effects (open blue circles), and using the nonperturbative oscillation code acor (solid and dashed black curves). The filled red circles correspond to the frequencies obtained with the approximate Eq. (24) for model B4. 

Open with DEXTER 
Fig. 3 Variations in the asymmetry of two l = 2 rotational multiplets during an avoided crossing, computed using Eq. (2) for m = ± 2 components. The symbols are identical to those of Fig. 2. 

Open with DEXTER 
To study the effects of neardegeneracy on the mode frequencies during the avoided crossing between the two l = 2 multiplets, we recomputed the evolution of our bestfit model of KIC 734123, but stopped the evolution at a time t_{0} just before the avoided crossing (the corresponding model is hereafter referred to as model B0). We then resumed the evolution with small time steps (1 Myr) in order to span the entire avoided crossing. This sequence of models is labeled B0 through B8 and referred to as such in subsequent discussions. For each model in the sequence, we computed the unperturbed mode frequencies (ω_{0,a} and ω_{0,b}) and eigenfunctions (ξ_{0,a} and ξ_{0,b}) of the two l = 2 modes that are in avoided crossing. We then calculated ω_{1,a} and ω_{1,b} using Eq. (8), and ω_{1,ab} using Eq. (13). The perturbed mode frequencies were then obtained using Eq. (11). The variations in the mode frequencies for our sequence of models are plotted as a function of time in Fig. 2 (blue open circles). We note that the frequencies of m = 0 modes in fact correspond to the frequencies of unperturbed modes. Before the avoided crossing, the lowestfrequency multiplet has a gmode behavior and the eigenfrequencies of its components thus increase owing to the core contraction. The highestfrequency multiplet behaves as a p mode and the frequencies of its components are decreasing owing to the increase of the stellar radius. When these two multiplets have comparable frequencies, the modes become mixed and the multiplets progressively exchange natures.
We found that the frequencies of the two l = 2 multiplets do indeed show significant levels of asymmetry during the avoided crossing. We quantified these asymmetries by using Eq. (2) for m = ± 2 components in the same way as was done with the observations. The results are shown in Fig. 3. As expected, the asymmetries of the multiplets are maximal when the modes that are bumping each other are at their smallest frequency difference, that is, when the modes are the most mixed. Away from the avoided crossing, the asymmetries decrease and eventually become negligible.
The asymmetries obtained with our firstorder perturbative approach including neardegeneracy effects show encouraging similarities with the ones observed in the spectrum of KIC 7341231 using Kepler data. Indeed, the higherfrequency mode is found to have a positive δ_{asym} during the avoided crossing, while the lowerfrequency multiplet has a negative δ_{asym}. This is in agreement with the observed asymmetries (see Sect. 2). The orders of magnitude of the asymmetries are also comparable, although the ones found here are smaller. This difference is discussed in detail in Sect. 5.
3.4. Mode eigenfunctions
As mentioned in Sect. 3.2, we were also able to estimate the eigenfunctions of the components of both multiplets for our sequence of models. These eigenfunctions can provide valuable insight about the trapping of each component during the avoided crossing. In Fig. 4, we show the integrand of the mode inertia for all the mcomponents of the higherfrequency l = 2 multiplet at the time during the avoided crossing when the asymmetry is the strongest (t = t_{0} + 3 Myr). The radial displacements were normalized to their surface value so that the contribution from the pmode cavity is nearly identical for all the components. As can be seen in Fig. 4, the contribution of the core to the mode inertia varies by a factor of approximately ten between the m = −2 and the m = + 2 components. This shows that during the avoided crossing, the mcomponents of an l = 2 multiplet are trapped differently inside the star. This sheds new light on the asymmetries of multiplets. Indeed, within a multiplet, the components that have a more glike behavior mainly probe the core rotation, while the plike components are more sensitive to the envelope rotation. This suggests that there is a link between the internal rotation of the star and the asymmetries in the rotational multiplets of mixed modes, which is addressed in the following section.
Fig. 4 Integrand of the mode inertia for the mcomponents of the higherfrequency l = 2 multiplet at the time during the avoided crossing where it shows the strongest asymmetry (t = t_{0} + 3 Myr). The radial and horizontal displacements (ξ_{r} and ξ_{h}) were computed using a perturbative approach as explained in Sect. 3 (gray dashed lines) and with the code acor (solid black lines). 

Open with DEXTER 
3.5. Link between multiplet asymmetry and differential rotation
The link between differential rotation and multiplet asymmetry can be simply understood based on Fig. 2. This figure shows that the avoided crossings between the mcomponents of the modes occur at slightly different epochs. The m = + 2 modes are the first to bump with each other, followed by the m = + 1 modes, and so on. The time delay between the avoided crossings of the different mcomponents is caused by the difference between the splittings of gmodes and those of pmodes. For the rotation profile that we assumed here, the splittings of g modes are larger than those of p modes. As a result, the m = + 2 modes are the first ones to become close enough to each other for the coupling between p and gmodes to start modifying their frequencies, as can be seen in Fig. 2. The link with multiplet asymmetries is clear. For instance, at t−t_{0} = 3 Myr, in the lowerfrequency multiplet, the m< 0 components are still behaving as g modes, while the m> 0 components start to resemble p modes because they are already coupling with the components of the higherfrequency mode. The splittings of the m< 0 components are therefore larger than those of the m> 0 components, thus producing a negative asymmetry (see Fig. 3). Conversely, in the higherfrequency multiplet at that same time, the m< 0 components are still plike, while the m> 0 components have already taken a glike behavior. This produces a positive asymmetry, as confirmed by Fig. 3.
If this interpretation is correct, we expect the asymmetries to increase as the degree of differential rotation increases. Conversely, we expect the asymmetries to vanish when the p modes and the g modes have identical rotational splittings. For l = 2 modes, this occurs when ⟨ Ω ⟩ _{c} = 6/5 ⟨ Ω ⟩ _{e}, that is, when the rotation is nearly solidbody like^{1}. When the envelope is spinning faster than the core, we expect the asymmetries to increase again as the envelopetocore ratio increases. To check this, we modified the surface rotation Ω_{e} in the input rotation profile given by Eq. (15) to produce coreenvelope contrasts of 20, 10, 5, 1, and 0.5. For each of these rotation profiles, the mode frequencies were computed as before using Eq. (11) over the course of the avoided crossing. The asymmetries were then calculated using Eq. (2). We plotted the degree of asymmetry of both multiplets as a function of time in Fig. 5. To guide the eye in this figure, we performed a cubic spline interpolation of δ_{asym}(t−t_{0}) for each rotation profile. When the core spins faster than the envelope, we observe that the asymmetry of the multiplets increases as the coretoenvelope rotation ratio increases, as predicted above. The asymmetries are at their smallest in the case of a solidbody rotation. When the envelope spins faster than the core, the asymmetries increase again as the envelopetocore rotation ratio increases, as we had foreseen. This clearly confirms that differential rotation is indeed the cause of multiplet asymmetries during avoided crossings. This opens the interesting possibility of using l = 2 multiplets to probe the internal rotation of red giants.
Fig. 5 Degree of asymmetry δ_{asym} of the two l = 2 multiplets during the avoided crossing. Different rotation profiles were considered, with coreenvelope contrasts (Ω_{c}/ Ω_{e}) of 20 (purple circles), 10 (red triangle), 5 (blue squares), 1 (green downwardtriangle), and 0.5 (cyan stars). The higherfrequency multiplet is represented by filled symbols and the lowerfrequency multiplet by empty symbols. To guide the eye, we overplotted cubic spline interpolations δ_{asym}(t−t_{0}) for each rotation profile — solid (resp. dashed) lines for higherfrequency (resp. lowerfrequency) multiplet. 

Open with DEXTER 
3.6. Link between multiplet asymmetry and mode trapping
Figure 5 also shows that the development of the asymmetry is not simultaneous for both multiplets during the avoided crossing. Two different cases can be distignusihed.

When ⟨ Ω ⟩ _{c}> 6/5 ⟨ Ω ⟩ _{e}, the higherfrequency multiplet reaches its maximum asymmetry before the lowerfrequency multiplet does (in terms of the absolute value). This can be understood as follows. Before the avoided crossing, the higherfrequency multiplet has a plike behavior. At the beginning of the avoided crossing, the m = + 2 component becomes more glike and its splitting increases because of the fast core rotation. As a result, the numerator of δ_{asym} in Eq. (2) increases. Since the other components still have a plike behavior, their rotational splittings are sensitive mainly to the envelope rotation and thus the denominator of δ_{asym} remains small, leading to a large degree of asymmetry. On the other hand, the lowerfrequency multiplet initially has a glike behavior. When the avoided crossing begins, its m = + 2 component takes on a plike character and its splitting decreases because of the slow envelope rotation. The numerator of δ_{asym} becomes increasingly negative but since the other components are still gdominated, their splittings are sensitive to the core rotation and thus the denominator of δ_{asym} remains large. As a result,  δ_{asym}  is initially smaller than for the higherfrequency multiplet. Near the end of the avoided crossing, the situation is reversed and  δ_{asym}  is larger for the lowerfrequency multiplet. It is important to point this behavior out because it shows that the asymmetries of rotational multiplets carry information about how far the star is in the avoided crossing process, and thus on the trapping of the two l = 2 modes. We make use of this property in Sect. 5 when developing a procedure to infer core and envelope rotation rates from the frequencies of l = 2 mixed modes.

When ⟨ Ω ⟩ _{c}< 6/5 ⟨ Ω ⟩ _{e}, the splittings of l = 2p modes are larger than those of g modes, and the evolution of multiplet asymmetries during the avoided crossing is the opposite of the previous case. We can indeed see in Fig. 5 that for ⟨ Ω ⟩ _{c}/ ⟨ Ω ⟩ _{e} = 0.5, the lowerfrequency multiplet reaches its maximum of asymmetry (in absolute value) earlier than the higherfrequency multiplet.
4. Nonperturbative approach
Before tackling the question of probing the internal rotation of giants using l = 2 rotational multiplets (see Sect. 5), we tested the validity of the firstorder perturbative approach that includes neardegeneracy effects developed in Sect. 3. For this purpose, we used the oscillation code acor (Adiabatic Computations of Oscillations including Rotation, Ouazzani et al. 2012), which is one of only two codes that take into account the effects of rotation in a nonperturbative manner (the other being the Twodimensional Oscillation Programme, a.k.a. top, developed by Reese et al. 2006). These codes were in fact designed to compute the pulsations of rapidlyrotating and therefore significantly distorted stars.
acor solves the equations of hydrodynamics perturbed by Eulerian fluctuations. Rotation is included by taking into account both the effects of the centrifugal distortion and those of the Coriolis acceleration. The numerical method is based on a spectral multidomain method that expands the angular dependence of pulsation modes on sphericalharmonic series. The radial differentiation is done by means of a sophisticated finite differences method which is accurate up to the fifth order in terms of the radial resolution (Scuflaire et al. 2008). The code has been validated by comparison with the results of TOP (Reese et al. 2006) for polytropic models and the agreement between the two codes was found to be excellent (Ouazzani et al. 2012).
We used acor to compute the mode frequencies for the sequence of models B0 through B8, which follows the avoided crossing between the two l = 2 multiplets. We used the same rotation profile as in Sect. 3 (Eq. (15) with a coreenvelope contrast of five). The results are overplotted in Fig. 2 (black solid and dashed lines). The agreement with the frequencies obtained with the firstorder perturbative approach is excellent.
We also calculated multiplet asymmetries δ_{asym} by inserting acor mode frequencies into Eq. (2). As can be seen in Fig. 3, the results match the asymmetries obtained with firstorder perturbed frequencies in the vicinity of the avoided crossing (2 Myr <t−t_{0}< 5 Myr) very well. Outside this range, differences appear between the two approaches. For Myr, the asymmetry of the lowestfrequency multiplet increases again, as can be seen in Fig. 3. At this point, this mode still behaves mainly as a p mode, but it is already in the process of bumping with a gravitydominated mode at lower frequency. Similarly, at t = t_{0}, the asymmetry of the largestfrequency multiplet is negative, which is caused by the bumping with another gravitydominated mode at larger frequency. In comparison, the asymmetries predicted by the firstorder perturbative approach vanish outside the avoided crossing because the two modes under consideration are no longer neardegenerate and adjacent l = 2 modes were not included. This discrepancy could be resolved by including the two l = 2 adjacent modes in the firstorder perturbative approach; however, this is not necessary since l = 2 mixed modes can be detected only when they have the most mixed character, that is, only in the close vicinity of the maximum of an avoided crossing.
Finally, we also computed the eigenfunctions of all components of the two multiplets along the sequence of models using acor. The radial and horizontal displacements were used to compute the integrand of the mode inertia for the higherfrequency l = 2 multiplet at t = t_{0} + 3 Myr, and the results were overplotted in Fig. 4. The agreement with the eigenfunctions obtained with the firstorder perturbative approach is once again excellent.
These comparisons clearly validate the use of a firstorder perturbative approach taking neardegeneracy into account to estimate the effects of rotation on the frequencies of mixed modes during avoided crossings. In particular, they confirm that secondorder effects of rotation on the mode frequencies can indeed be safely neglected here.
5. Probing internal rotation using l = 2 multiplets
We have shown in Sect. 3 that the asymmetry of l = 2 multiplets depends directly on the amount of differential rotation. In this section, we investigate whether the frequencies of l = 2 asymmetric multiplets can be used to measure mean core rotation ⟨ Ω ⟩ _{c} and envelope rotation ⟨ Ω ⟩ _{e} in red giants. For this purpose, the perturbative expression of the frequency shifts caused by rotation given by Eq. (11) is convenient because it opens up the possibility for applying inversion methods to recover the internal rotation profile.
5.1. A simplified expression for firstorder perturbation of mode frequencies
We separated the firstorder perturbation of mode eigenfrequencies into contributions from the core and the envelope. For instance, we write for mode a (which we here considered to correspond to the lowerfrequency multiplet) (16)where and The kernel K_{a}(r) is obtained from Eq. (9). The radius r_{c} was defined as the outer turning point of the gravity waves (N(r_{c}) = ω_{0,a}, where N(r) is the BruntVäisälä frequency), so that ⟨ Ω ⟩ _{c} is an average of the rotation rate over the gmode cavity. Correspondingly, ⟨ Ω ⟩ _{e} is an average of the angular velocity in the layers above r_{c}. It is essentially a measure of the mean rotation rate in the pmode cavity because the evanescent zone (region between r_{c} and the inner turning point of the pmode cavity) contributes very little to this average. We also introduced β_{c,b} and β_{e,b} for mode b (corresponding to the higherfrequency multiplet) with expressions analogous to Eq. (17) and (18). We note that r_{c} varies little with mode frequency because of the sharp decrease of the BruntVäisälä frequency at the edges of the gmode cavity. We thus took the same value for both modes. Likewise, the values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} are very similar whether we use the kernel of mode a or that of mode b in Eqs. (19) and (20). In the following discussion, we will thus refer to these quantities without mentioning which kernel was used to compute them. The values of β_{c} and β_{e} were computed for modes a and b for our sequence of models spanning the avoided crossing between the two l = 2 modes. The results are shown in Fig. 6. One can clearly see that mode a (lowerfrequency multiplet) is initially sensitive to core rotation (gdominated) and its sensitivity gradually shifts to the envelope (pdominated), while mode b (higherfrequency multiplet) has the opposite behavior.
Fig. 6 Variations in the parameters β_{c}, β_{e}, γ_{c}, and γ_{e} (defined in Eqs. (17), (18), (22), and (23)) along the course of the avoided crossing between the two l = 2 multiplets for our reference sequence of models B0 through B8. Mode a (resp. b) corresponds to the lower (resp. higher) frequency multiplet. Core quantities are represented by open circles and envelope quantities by open squares. Interpolations using cubic splines are overplotted (blue dashed lines for core parameters and red dotdash lines for envelope parameters). 

Open with DEXTER 
The coupling term ω_{1,ab} was also rewritten as (21)where and the function K_{ab}(r) is defined by Eq. (13). The values of γ_{c} and γ_{e} during the avoided crossing are shown in Fig. 6. As described in Appendix B, a simplified expression can be obtained for ω_{1,ab}, which shows that the coupling term vanishes in the particular case of a rotation profile for which the splittings of pure g modes are identical to those of pure p modes. In this case, we recover the eigenfrequencies of the nondegenerate case (see Eq. (11)) and the multiplets remain symmetric, in agreement with what was found in Sect. 3.5.
As a result, the firstorder perturbed mode frequencies can be expressed as (24)where we have defined and similar expressions hold for and δβ_{e} as well.
We thus obtained approximate expressions for the perturbed mode frequencies involving only four parameters (⟨ Ω ⟩ _{c}, ⟨ Ω ⟩ _{e}, ω_{0,a}, and ω_{0,b}) and the rotational kernels of the chosen reference model. This gives us the possibility of estimating the mean core and envelope rotation using l = 2 modes by fitting Eq. (24) to the observed mode frequencies. We note that β_{c} and β_{e} can in fact be approximately estimated from seismic observables without the need of a reference stellar model, as was shown by Goupil et al. (2013; see also Mosser et al. 2015). If approximate expressions relating γ_{c} and γ_{e} to observable quantities can be derived, then Eq. (24) could be used to estimate ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} directly from the observed mode frequencies. This is, however, out of the scope of the present work, and instead we compute β_{c}, β_{e}, γ_{c}, and γ_{e} from reference stellar models in this study.
5.2. Preliminary test
To test this approach, we used the eigenfrequencies obtained with acor for a model in the close neighborhood of the avoided crossing between the l = 2 multiplets (model B4) as artificial data, assuming the same rotation profile as in Sect. 3 (Eq. (15) with a coreenvelope contrast of five). In this preliminary test, we assumed that all the components of the multiplets are detected (ten frequencies) and that the frequencies are free of noise. We used the same model (model B4) as the reference model to perform inversions of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e}. We used the nonrotating version of acor to compute the unperturbed eigenfunctions of model B4, which were then inserted into Eqs. (17), (18), (22), and (23) to obtain the values of β_{c}, β_{e} for both modes, as well as γ_{c}, γ_{e}. We then performed an iterative fit of Eq. (24) to the eigenfrequencies obtained with acor, using ⟨ Ω ⟩ _{c}, ⟨ Ω ⟩ _{e}, ω_{0,a}, and ω_{0,b} as free parameters of the fit. The optimal values obtained for the core and envelope rotation rates (⟨ Ω ⟩ _{c,fit}/ 2π = 716 nHz and ⟨ Ω ⟩ _{e,fit}/ 2π = 154 nHz) are very close to the true values of these quantities computed using Eqs. (19) and (20) for our reference rotation profile (⟨ Ω ⟩ _{c,mod}/ 2π = 710 nHz, and ⟨ Ω ⟩ _{e,mod}/ 2π = 154 nHz). The frequencies obtained using the bestfit parameters and Eq. (24) are overplotted in Fig. 2 (red filled circles at t = t_{0} + 4 Myr). The agreement with the frequencies computed with acor is better than 10 nHz, that is, below the precision to which rotational splittings can be measured with the four years of Kepler observations.
These results are encouraging, but they do not guarantee that we can recover ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} using the observed frequencies of l = 2 multiplets. Indeed, in this preliminary test, we considered the ideal case where the reference model is “perfect”. If we perform a similar fit, but use the values of β_{c}, β_{e}, γ_{c}, and γ_{e} computed with a different model (say for instance model B5), no satisfactory fit to the frequencies of acor can be found and the recovered values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} are incorrect. The reason for this is that the trapping of the modes differs between the chosen reference model and the model that was used to produce rotational multiplets.
The need for a reference model that produces a correct trapping of the modes was already emphasized in works that performed rotation inversions using dipolar modes. However, in the latter case, a large number of l = 1 modes can be detected, and these have a wide variety of natures ranging from gdominated to pdominated. Deheuvels et al. (2012) showed that all the nonrotating models that provide satisfactory matches to the frequencies of the m = 0 components of the observed l = 1 modes predict a similar trapping of the dipolar modes, quantified by the parameter ζ defined as the ratio between the mode inertia in the gmode cavity and the total mode inertia, that is, (27)Goupil et al. (2013) showed that the trapping of l = 1 modes can even be estimated reliably from the mode frequencies themselves, independently of stellar models, using the WentzelKramersBrillouin (WKB) approximation. For l = 2 modes the situation seems much more complex. Deheuvels et al. (2012) mentioned that the models that satisfactorily reproduce the frequencies of l = 1 modes predict different mode trappings for the l = 2 mixed modes, which means that these models cannot be directly used to perform inversions using l = 2 modes. Besides, the coupling between the pmode and gmode cavities is weaker for l = 2 modes than for l = 1 modes because the evanescent region that separates them is wider. We can thus detect l = 2 modes only when their frequencies are very close to those of pure p modes. Outside these frequency ranges, the l = 2 modes have inertias that are too large for the modes to have detectable heights in the oscillation spectra, even with the longest Kepler datasets (see, e.g., Grosjean et al. 2014). This makes it much harder to precisely estimate the asymptotic properties of l = 2 pure g modes, which are required if we want to apply the method of Goupil et al. (2013) to l = 2 modes.
5.3. New fitting procedure
In Sect. 3.5, we have shown that the asymmetries of the rotational multiplets carry information not only about the internal rotation, but also about how far the star is in the avoided crossing process, and thus on the trapping of the l = 2 modes. This could yield a solution to the problem mentioned above. We explored the possibility of also recovering the trapping of the l = 2 modes (expressed by the ζ parameter) when fitting Eq. (24) to the data. The idea is to allow the age of the reference model to vary slightly in our fitting procedure by considering the time elapsed since the beginning of the avoided crossing, Δt ≡ t−t_{0}, as an additional free parameter. We note that because of the much weaker coupling, the timescale over which the trapping of l = 2 modes is modified is small compared to the equivalent quantity for l = 1 modes. As a result, the small age adjustment that is made in this fit has a negligible impact on the trapping of l = 1 modes in the model.
With this new fitting procedure, the value of Δt changes at each iteration, and in principle the rotational kernels need to be recomputed every time, which is numerically very costly. In fact the only quantities that are needed are β_{c}, β_{e}, γ_{c}, and γ_{e}, whose variations during the avoided crossing are quite smooth (see Fig. 6). We thus chose to perform a cubicspline interpolation of these quantities using our sequence of models B0 through B8. The results are overplotted in Fig. 6. For several values of Δt, we compared the interpolated values of β_{c}, β_{e}, γ_{c}, and γ_{e} to the ones computed using the rotational kernels and found agreements at a level of 1% or better. Inserting the interpolated parameters into Eq. (24), this 1% discrepancy translates into frequency differences at most of the order of 1 nHz, which is small compared to the precision to which the mode frequencies can be measured from the observations. We thus used interpolated values of β_{c}, β_{e}, γ_{c}, and γ_{e} in the fits presented in the subsequent sections.
For clarity, we briefly summarize the steps of the fitting procedure:

1.
Compute a reference model that satisfactorily reproduces theobservational constraints of the star under study. This modelshould in particular reproduce the frequencies of the observedl = 1 modes and roughly match the frequencies of thel = 2 modes that are undergoing an avoided crossing.

2.
Recompute the evolution of the reference model until a time t_{0} just before the avoided crossing.

3.
Restart the evolution with small timesteps to span the entire avoided crossing.

4.
For each model of the sequence, compute the values of the β and γ parameters using Eqs. (17), (18), (22), and (23).

5.
Perform cubic spline interpolations of the variations in the β and γ parameters with time so that they can be estimated at any time during the avoided crossing.

6.
Perform an iterative fit of Eq. (24) to the observed frequencies of the component of the two l = 2 multiplets, using ⟨ Ω ⟩ _{c}, ⟨ Ω ⟩ _{e}, ω_{0,a}, and ω_{0,b} as free parameters. At each iteration of the fit, the β and γ parameters corresponding to the value of Δt are calculated using the interpolations.
Values of the parameters ⟨ Ω ⟩ _{c}, ⟨ Ω ⟩ _{e}, Δt, ζ_{0,a}, and ζ_{0,b} obtained when fitting Eq. (24) to the mode frequencies of acor for models A, C, D, and E, assuming a rotation profile following Eq. (15) with a core rotation ⟨ Ω ⟩ _{c}/ 2π = 710 nHz and coreenvelope contrasts of 100, 20, 10, 5, and 2.
5.4. Test of the method with noisefree artificial data
To test the fitting procedure described above, we again generated artificial data (frequencies of l = 2 rotational multiplets in avoided crossing) using other models than model B. We conveniently used models A, C, D, and E obtained by Deheuvels et al. (2012). All these models were found to reproduce the observed frequencies of the m = 0 components of the observed l = 1 mixed modes of KIC 7341231 quite satisfactorily. They also match the observed atmospheric parameters of the star. The main difference between these models is the metallicity, which ranges from − 0.75 dex (model A) to − 1.75 dex (model E), reflecting the poorly constrained metallicity of the star. All these models were used by Deheuvels et al. (2012) to perform inversions of the rotation profile of KIC 7341231, and were found to yield almost identical results. They all include an l = 2 avoided crossing between two l = 2 modes around 380 μHz, but they predict different mode trappings for these two l = 2 mixed modes. For models A, C, D, and E, we computed the frequencies of the two l = 2 rotational multiplets with acor, assuming the same rotation profile as before (Eq. (15) with a core rotation of 710 nHz) and we considered coreenvelope contrasts of 100, 20, 10, 5, and 2.
For each of these scenarios, we fitted Eq. (24) to the frequencies of acor, considering ⟨ Ω ⟩ _{c}, ⟨ Ω ⟩ _{e}, ω_{0,a}, ω_{0,b}, and Δt as free parameters. As described above, for each iteration of the fit we interpolated through the values of β_{c}, β_{e}, γ_{c}, and γ_{e} of the sequence of models B0 though B8. The inversions of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} are thus performed using a different model than the ones used to generate the artificial data. The results are given in Table 2. The first observation is that the recovered mean core and envelope rotation rates agree very well with the actual values computed with the input rotation profile and Eqs. (19) and (20), regardless of the chosen coreenvelope ratio and the model that is used to generate artificial data.
The optimal value of Δt can also be used to obtain an estimate of ζ_{0,a} and ζ_{0,b}, which quantify the trapping of the m = 0 components of the two l = 2 multiplets. Indeed, we can compute ζ_{0,a} and ζ_{0,b} for the models B0 through B8 and perform cubic spline interpolations of these quantities in the same way as was done for the β and γ parameters. The values of ζ_{0,a} and ζ_{0,b} can then be obtained at the fitted time Δt. The corresponding values are given in Table 2. They can be compared to the actual values of ζ_{0,a} and ζ_{0,b} computed by injecting the eigenfunctions of the model that was used to generate data into Eq. (27). As can be seen in Table 2, the agreement between the recovered values of ζ_{0,a} and ζ_{0,b} and the actual ones computed from the model is excellent. This confirms that we can indeed determine the trapping of the modes using the frequencies of two l = 2 rotational multiplets in avoided crossing.
5.5. Test of the method under realistic conditions
Next, we tested our method with more realistic artificial data by including noise and by taking into account that the observed l = 2 multiplets are generally incomplete. We performed Monte Carlo simulations using models A, C, D, and E and considered different rotation profiles. In the following, we show only the results obtained with model D (the other models give undistinguishable results), assuming a rotation profile with ⟨ Ω ⟩ _{c}/ 2π = 710 nHz as before and a coreenvelope contrast of 5. The influence of a change in this parameter is mentioned when relevant.
5.5.1. Effects of adding noise on mode frequencies
We first assumed uncertainties of 15 nHz for the frequencies of the components of rotation multiplets. This value corresponds to the typical error bar resulting from the extraction of mode frequencies in the oscillation spectra of red giants observed during four years with Kepler (see for instance Table 1). We added a realization of Gaussian noise with a standard deviation of 15 nHz to the mode frequencies computed with acor and fitted Eq. (24) to the resulting set of frequencies. This operation was then repeated 1000 times. The values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} recovered from our Monte Carlo simulation for this case are shown in Fig. 7. Using the distribution of the fitted parameters, we derived ⟨ Ω ⟩ _{c}/ 2π = 721 ± 18 nHz and ⟨ Ω ⟩ _{e}/ 2π = 150 ± 16 nHz, which is in close agreement with the actual values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} computed using Eqs. (19) and (20) (see Table 2). It is clear from Fig. 7 that even with noisy data one can efficiently recover the average core and envelope rotation rates.
We repeated this procedure using different sets of initial parameters for the iterative fits. For instance, we tried starting from a solidbody rotation or with profiles such that ⟨ Ω ⟩ _{e}> ⟨ Ω ⟩ _{c}. We also varied the initial value of Δt. In all cases, the obtained results are undistinguishable from the ones mentioned above, which shows that the procedure has little dependence on the chosen initial parameters. We note, however, that an ambiguity may arise in the determination of Δt. As was mentioned in Sect. 3.5, a larger  δ_{asym}  for the higherfrequency multiplet means that the star is either at the beginning of the avoided crossing with a core rotating faster than the envelope, or near the end of the avoided crossing with an envelope rotating faster the core (see Fig. 5). This potential degeneracy was not found when using models A, C, D, or E to generate artificial data, but it arose when using model B3. When choosing ⟨ Ω ⟩ _{c}> ⟨ Ω ⟩ _{e} in the initial parameters of the fit, we recovered the correct mean core and envelope rotation rates with Δt = 3 Myr, as expected. However, when starting with ⟨ Ω ⟩ _{e}> ⟨ Ω ⟩ _{c}, the fit converged to a spurious solution with a fast rotating envelope near the end of the avoided crossing (⟨ Ω ⟩ _{c}/ 2π = 194 nHz and ⟨ Ω ⟩ _{e}/ 2π = 621 nHz, with Δt = 5.1 Myr). However, if this potential ambiguity arises, it should be easy to resolve, since the two scenarios provide very different splittings for dipolar modes. In all the red giants studied so far using Kepler data, the rotational splittings of l = 1 modes safely ruled out the possibility of an envelope rotating faster than the core.
Fig. 7 Values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} obtained by fitting Eq. (24) to the mode frequencies of model D computed with acor, assuming a rotation profile with ⟨ Ω ⟩ _{c}/ 2π = 710 nHz and ⟨ Ω ⟩ _{c}/ ⟨ Ω ⟩ _{e} = 5. A random noise with an rms of 15 nHz was added to the mode frequencies and 1000 iterations were performed. The quantity β_{c} ⟨ Ω ⟩ _{c} + β_{e} ⟨ Ω ⟩ _{e} is conserved along the blue lines, where β_{c} and β_{e} correspond to mode a (dashed line) or mode b (dotdashed line). The red cross indicates the actual values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} computed with Eqs. (19) and (20). 

Open with DEXTER 
Fig. 8 Distribution of the values of ⟨ Ω ⟩ _{c} (left) and ⟨ Ω ⟩ _{e} (right) obtained by fitting Eq. (24) to the mode frequencies of acor for model D, assuming a rotation profile with ⟨ Ω ⟩ _{c}/ 2π = 710 nHz and ⟨ Ω ⟩ _{c}/ ⟨ Ω ⟩ _{e} = 5. A random noise with an rms of 15 nHz (top), 30 nHz (middle), and 50 nHz (bottom) was added to the mode frequencies. The colors correspond to the different scenarios that were tested: ideal case where all the components are detected (black), inclination angle of 90° (purple) or 30° (red), and case of an overlap with the closest radial mode (blue). The vertical dashed lines indicate the actual values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} computed with Eqs. (19) and (20). 

Open with DEXTER 
We observe from Fig. 7 that the recovered values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} are highly correlated. This can be understood because the average rotation ⟨ Ω ⟩ over the whole interior, which can be written as (28)can be recovered from the observed modes frequencies much more precisely than ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} individually. As a result, the quantity β_{c} ⟨ Ω ⟩ _{c} + β_{e} ⟨ Ω ⟩ _{e} is expected to be roughly conserved for all the iterations of the Monte Carlo simulation. The conservation of this quantity in the (⟨ Ω ⟩ _{e}, ⟨ Ω ⟩ _{c}) plane defines a straight line, which is overplotted in Fig. 7 and accounts well for the correlation between the fitted values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e}.
To study the influence of the uncertainties of mode frequencies on the performance of our method, we performed Monte Carlo simulations considering error bars of 30 and 50 nHz. The distributions of the recovered values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} are shown in Fig. 8 (black histograms). The error bars of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} increase to ~35 nHz for σ_{ν} = 30 nHz and to ~65 nHz for σ_{ν} = 50 nHz. In the latter case, the distributions of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} become slightly bimodal, with a small secondary peak at ⟨ Ω ⟩ _{c}/ 2π ~ 480 nHz and ⟨ Ω ⟩ _{e}/ 2π ~ 360 nHz. These spurious solutions are obtained along with extremum values of Δt (either 0 Myr or 8 Myr), that is, away from the avoided crossing. These solutions therefore have negligible multiplet asymmetries, and they correspond to the solutions that would be obtained if neardegeneracy effects were ignored. The fact that the fits sometimes converge to this spurious solution in our noisiest scenario can be understood because the frequency uncertainty (50 nHz) becomes comparable to the asymmetry of the multiplets (the numerator of δ_{asym} is around 135 nHz for the rotation profile that we have chosen for these simulations). If we use rotation profiles with larger coreenvelope contrasts, the asymmetry of rotational multiplets during the avoided crossing increases and the bimodality of the recovered rotation rates vanishes, even for σ_{ν} = 50 nHz.
In our Monte Carlo simulations, the parameter Δt was also recovered with very little scatter (3.63 ± 0.02 for σ_{ν} = 15 nHz). As was done before, we used the optimal value of Δt to estimate ζ_{0,a} and ζ_{0,b} for each simulation. The distribution obtained for ζ_{0,a} is shown in Fig. 9 (the distribution of ζ_{0,b} is qualitatively similar). The recovered values are in close agreement with the actual values of the ζ parameters computed by inserting the eigenfunctions of model D into Eq. (27) (we obtained ζ_{0,a} = 0.481 ± 0.008 and ζ_{0,a} = 0.545 ± 0.008 for σ_{ν} = 15 nHz). This confirms that the trapping of the modes can be precisely recovered, even when the effects of noise are included.
5.5.2. Effects of missing components
In addition to the effects of the noise, we must take into account the fact that in general not all the components of the l = 2 multiplets can be detected. This is caused mainly by the inclination of the star, but also by the possible overlap of certain components with the adjacent radial or mixed dipolar modes. To determine the effects of missing components on the fitting procedure, we tested the following scenarios:
 1.
Ideal case (all mcomponents are detected),
 2.
Inclination angle of 90° (m = ± 1 components missing),
 3.
Inclination angle of ~30° (m = ± 2 components missing),

Overlap of a part of the higherfrequency multiplet with the closest radial mode (m = + 1, + 2 components of this multiplet missing).
Fig. 9 Distribution of the values of ζ_{0,a} recovered when fitting Eq. (24) to the mode frequencies of acor for model D, assuming a rotation profile with ⟨ Ω ⟩ _{c}/ 2π = 710 nHz and ⟨ Ω ⟩ _{c}/ ⟨ Ω ⟩ _{e} = 5. The lines and colors have the same meaning as in Fig. 8. The vertical dashed lines indicate the actual value of ζ_{0,a} computed with the eigenfunctions of the model and Eq. (27).
Open with DEXTER
Scenario 1 was already tested in Sect. 5.5.1. We repeated the same Monte Carlo simulations for scenarios 2, 3, and 4. The distributions of the obtained values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} are shown in Fig. 8. For scenarios 2 and 4, the missing components result in a slight increase in the error bars of the recovered core and envelope rotation rates, but overall the performance of the method remains quite comparable to the ideal case where all the components are detected. In the case of a low inclination angle (scenario 3), the performance of the method is much poorer, with the bimodality of the distributions of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} appearing even with σ_{ν} = 15 nHz. This is expected since in the absence of the m = ± 2 components, the asymmetry can only be detected through the m = ± 1 components and it is much smaller. With the rotation profile that we have chosen here, ω_{m = −1} + ω_{m = 1}−2ω_{m = 0} is around 35 nHz, which is comparable to the uncertainties of the mode frequencies. This explains the large increase of the secondary peak in the distributions for scenario 3 as σ_{ν} increases. As mentioned in Sect. 5.5.1, this secondary peak corresponds to the cases where the fit converged toward Δt = 0 Myr or Δt = 8 Myr, that is, toward a configuration without multiplet asymmetries. These cases are clearly identified in Fig. 9 because they form sharp peaks at ζ_{0,a} ~ 0.2 and ζ_{0,a} ~ 0.96, which correspond to Δt = 8 Myr or Δt = 0 Myr, respectively.
We therefore conclude that provided the inclination angle of the star is large enough to detect m = ± 2 components (i ≳ 50°), the full set of Kepler data (four years of observations) should be sufficient to reliably estimate the mean core and envelope rotation rates using only two l = 2 modes undergoing an avoided crossing.
6. Application to KIC 7341231
The young red giant KIC 7341231 is a favorable case for study considering the results of our Monte Carlo simulation. The inclination angle of the star is high (i = 85° ± 5°, Deheuvels et al. 2012) so that the m = ± 2 components of the l = 2 multiplets are clearly visible. Compared to scenario 2 which was tested in Sect. 5.5.2, we also have measurements for the frequencies of the m = + 1 components, even though the associated error bars are larger than for the other components. The error bars of the m = −2,0, and 2 components are all around or below 15 nHz. Consequently, we expected to be able to reliably extract ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} for this star using our method.
We fitted Eq. (24) to the observed frequencies of the two l = 2 multiplets which are listed in Table 1. The fit converged to a core rotation of ⟨ Ω ⟩ _{c}/ 2π = 771 ± 13 nHz and an envelope rotation of ⟨ Ω ⟩ _{e}/ 2π = 45 ± 12 nHz. The values obtained for the other fitted parameters are given in Table 3. The quoted error bars were obtained by performing Monte Carlo simulations using the observed mode frequencies and their uncertainties. We repeated the fitting procedure using a wide range of initial parameters. We found that if ⟨ Ω ⟩ _{c}> ⟨ Ω ⟩ _{e} in the set of initial parameters, we always recover the solution that is given in Table 3, regardless of the initial guess for Δt. As mentioned in Sect. 5.5, we found that if the fit is started with initial rotation rates such that ⟨ Ω ⟩ _{e}> ⟨ Ω ⟩ _{c}, it converges to a solution with the envelope rotating much faster than the core (⟨ Ω ⟩ _{c}/ 2π = 67 ± 14 nHz and ⟨ Ω ⟩ _{e}/ 2π = 675 ± 12 nHz, with Δt = 3.98 ± 0.02). This alternate solution can be clearly excluded for KIC 7341231 because if it were the case, the glike dipolar modes would have smaller splittings than the plike ones, contrary to what is observed.
The mode frequencies produced by the set of bestfit parameters are overplotted in Fig. 1 (vertical dotted lines). They are in very close agreement with the observed mode frequencies. The agreement can be qualitatively assessed by computing the χ^{2} using the observed frequencies and error bars listed in Table 1. For the bestfit parameters, we obtained χ^{2} = 4.8. This translates into a reduced χ^{2} of 1.6 (the fit involves five free parameters for eight observables, thus yielding a number of degrees of freedom of three), which indicates a statistically good agreement. With this procedure we were also able to predict the location of the m = −1 components of the multiplets, which could not be detected from the observed oscillation spectrum. For the lowerfrequency multiplet, the theoretical location of the m = −1 component coincides with a slight excess of power, while nothing can be seen in the power spectrum around the frequency predicted for the m = −1 component of the higherfrequency multiplet.
7. Discussion
7.1. Other possible sources of multiplets asymmetry
As we have shown in this paper, neardegeneracy effects are very likely the origin of the multiplet asymmetries in KIC 7341231. Conditions that cause neardegeneracy of modes are expected to occur in this star, and our treatment of these effects shows that they produce asymmetries similar to those that are observed. Additionally, no asymmetries were found in the l = 1 multiplets of KIC 7341231, which implies that whatever phenomenon is producing multiplet asymmetries in this star does not affect dipolar modes. This can be accounted for by neardegeneracy effects, because the evanescent region for l = 1 modes is thinner, which produces a much stronger coupling between the gmode and pmode cavities. As a result, the frequency separation between dipolar mixed modes is larger than for quadrupolar modes and neardegeneracy effects are much weaker. Based on Deheuvels et al. (2012), the smallest frequency spacing between l = 1 modes in the spectrum of KIC 7341231 (7.64 ± 0.05 μHz) represents more than ten times the core rotation rate of the star, which confirms that neardegeneracyeffects should be small for dipolar modes in this star. For the sake of completeness, we briefly mention other potential sources of asymmetries.
Asymmetries in rotational multiplets can be the signature of a rotation rate that is fast enough to produce nonnegligible effects of the centrifugal force and secondorder effects of the Coriolis force. In this case, higherorder perturbations need to be included in the oscillation equations, which break the symmetry of rotational multiplets (e.g., Dziembowski & Goode 1992; Suárez et al. 2006). However, the rotation rate estimated for KIC 7341231 by Deheuvels et al. (2012) is much too low for secondorder effects of rotation to significantly contribute to rotation splittings. Besides, no asymmetries were detected in the l = 1 multiplets of KIC 7341231, contrary to what would be expected in the presence of fast rotation.
It is also known that latitudinal differential rotation can produce asymmetric rotational multiplets. To first order, the frequency shifts caused by rotation can be written as (29)Since the kernels K_{n,l,m}(r,θ) depend on the azimuthal order m, the splittings of the different m components, usually computed as (ω_{− m}−ω_{m})/2m, differ. However, it can be readily shown from the expression of the kernels that K_{n,l,−m} = K_{n,l,m}, so that δω_{n,l,−m} = −δω_{n,l,m} (e.g., Aerts et al. 2010). As a result, the asymmetry as defined by δ_{asym} in Eq. (2) is expected to cancel at the first order in the presence of latitudinal differential rotation. Secondorder effects of rotation may on the other hand produce nonzero values of δ_{asym}. In order to rule out latitudinal differential rotation as the cause of the observed asymmetries, we performed nonperturbative calculations with acor. The rotation law was chosen to be differential in latitude only, of the form Ω(θ) = Ω_{eq}(1−ΔΩ/Ω_{eq}cosθ^{2}), where Ω_{eq} and Ω_{pol} are the equatorial and polar rotation, respectively, and ΔΩ ≡ Ω_{eq}−Ω_{pol}. We chose Ω_{eq}/ 2π = 710 nHz and a latitudinal differential rotation of ΔΩ/Ω_{eq} = 0.36, which corresponds to the solar value. With this rotation profile, we obtained asymmetries of the order of δ_{asym} ~ 10^{3}, that is, much lower than the observed values. We thus conclude that latitudinal differential rotation cannot account on its own for the observed asymmetries.
Finally, multiplet asymmetries can also be the signature of an internal magnetic field. The possibility of a buried magnetic field inclined with respect to the rotation axis can be excluded for KIC 7341231, since such a field would produce a magnetic splitting of each of the 2l + 1 components in each rotational multiplets (Dicke 1982; Dziembowski & Goode 1985), and this is not observed. However, an internal field aligned with the rotation axis would not change the number of components, but would alter the fine structure pattern of the multiplets, potentially producing asymmetries. Establishing the impact of such a field on the rotational multiplets would require a dedicated study and is out of the scope of the present paper. We can, however, remark that a buried magnetic field would produce more distortion in the gdominated multiplets than in the pdominated ones (Hasan & ChristensenDalsgaard 1992), which is the opposite of what is observed in KIC 7341231 and KIC 5006817. Also, it is hard to imagine how an internal magnetic field would produce asymmetries in quadrupolar multiplets but not in dipolar ones in KIC 7341231.
7.2. Neardegeneracy effects in red giants
In this paper, we have focused on the case of the young red giant KIC 7341231. The internal rotation of this star is not peculiar compared to other similar targets, thus neardegeneracy effects are expected to be important for many other giants. In fact, since the intensity of the coupling between the pmode and gmode cavities decreases during the evolution (e.g., Mosser et al. 2017), we anticipate that neardegeneracy effects become even larger for giants more evolved than KIC 7341231, which lies at the base of the RGB. Therefore, these effects should always be considered when using the rotational multiplets of quadrupolar modes to probe the internal rotation of red giants.
For dipolar modes, neardegeneracy effects are smaller than those for quadrupolar modes, as explained in Sect. 7.1. However, they need to be considered for giants slightly more evolved than KIC 7341231. This is the case for KIC 5006817, for which asymmetrical l = 1 multiplets have been reported (Beck et al. 2014). It is straightforward to show that neardegeneracy effects are expected in this star. The closest frequency separation between l = 1 mixed modes in this star is ~0.8 μHz, which is even smaller than the core rotation rate of 0.93 μHz inferred for this star by Beck et al. (2014). In fact, the asymmetries detected in this star show a lot of similarities with the ones produced by neardegeneracy effects: they arise only for the modes that show the highest level of mixing (see Fig. 6 of Beck et al. 2014) and it is clear from Fig. 5 of Beck et al. (2014) that for all pairs of asymmetric multiplets, the lowestfrequency multiplet has a negative asymmetry, while the highestfrequency multiplet has a positive asymmetry, as we have found in Sect. 3. In the case of KIC 4448777, for which Di Mauro et al. (2016) reported asymmetries in the multiplets of several l = 1 modes, the authors found a mean core rotation of 0.75 μHz, which represents about 40% of the closest frequency spacing between l = 1 mixed modes. We thus expect neardegeneracy effects to be nonnegligible in this star also, although a dedicated study would be necessary to determine whether they are the cause of the reported asymmetries.
8. Conclusion
In this paper we investigated the origin of the asymmetries that have been detected in the rotational multiplets of mixed modes in red giants. We showed that in both red giants in which asymmetrical multiplets have been reported (KIC 7341231, Deheuvels et al. 2012; KIC 5006817, Beck et al. 2014), neardegeneracy effects are expected to arise because of the small frequency spacing between mixed modes. Such effects had not been included in previous studies. Using the young red giant KIC 7341231 as a test case, we accounted for the effects of neardegeneracy in the framework of a firstorder perturbative approach and showed that these effects indeed produce multiplet asymmetries with very similar features compared to the observations. We also validated this perturbative approach using the code acor (Ouazzani et al. 2012), which solves the oscillation equations including the effects of rotation in a nonperturbative manner.
We showed that the asymmetries of rotational multiplets arising for neardegenerate mixed modes are linked to differential rotation between the core (gmode cavity) and the envelope (pmode cavity). Indeed, for neardegenerate mixed modes, the components of the rotational multiplets are trapped differently inside the star. Each component is thus sensitive to an average of the internal rotation with different weights given to the core and the envelope. By exploiting this property, we developed a method to measure the mean core and envelope rotation using the rotational multiplets of neardegenerate mixed modes. This method was then tested using Monte Carlo simulations on artificial data, where we included the effects of the uncertainties of mode frequencies and considered different scenarios for the mode visibilities within the multiplets. We showed that, provided the dataset is long enough for the asymmetries to be detected at a statistically significant level, our method can be used to successfully and precisely recover the mean rotation rates of the core and envelope using the frequencies of only two l = 2 neardegenerate rotational multiplets. Higherinclination angles are more favorable, because then the m = ± 2 components can be detected, and they are the ones that show the largest asymmetries.
Next, we successfully applied our method to the young red giant KIC 7341231. We obtained a core rotation rate of ⟨ Ω ⟩ _{c}/ 2π = 771 ± 13 nHz for this star. This measurement is in good agreement with the results of Deheuvels et al. (2012), who had obtained ⟨ Ω ⟩ _{c}/ 2π = 710 ± 51 nHz by applying the OLA inversion method to the dipolar modes. The improved precision on the measurement of ⟨ Ω ⟩ _{c} is likely related to the longer time series that was used here compared to Deheuvels et al. (2012), who had only one year of data at their disposal. This result shows that we are now able to measure the core rotation rate using quadrupolar mixed modes, which had not been possible earlier. We intend to investigate in a future work whether this measurement, combined with the one obtained with dipolar modes, can be used to place constraints on the latitudinal differential rotation of the core – we know that modes of different degree are sensitive to different latitudes of the stellar rotation profile.
Interestingly, our method also enabled us to provide an estimate of the envelope rotation of KIC 7341231, which had not been possible for this star using the dipolar modes. We obtained ⟨ Ω ⟩ _{e}/ 2π = 45 ± 12 nHz, which is in agreement with the upper limit given by Deheuvels et al. (2012). We thus obtained an estimate of the coreenvelope contrast of 17 ± 5 for KIC 7341231. This is an improvement compared to the results of Deheuvels et al. (2012) who could only provide a lower limit of five for the coreenvelope rotation ratio. This measurement could be used to improve our knowledge of the timescale over which angular momentum is transported in red giants (Eggenberger et al. 2017).
Based on our results, we expect neardegeneracy effects to be important for quadrupolar mixed modes all along the redgiant branch. For dipolar mixed modes, these effects are negligible at the base of the redgiant branch, but must be taken into account for more evolved giants. We have shown in this paper that rotational multiplets affected by neardegeneracy effects require special treatment and that they give us information about the internal rotation of red giants.
Acknowledgments
We would like to thank T. Appourchaux who initially suggested that we look into the role of rotation on avoided crossings. We are also thankful to M. A. Dupret, J. Ballot, M. J. Goupil, and B. Mosser for fruitful discussions that were very helpful. We also thank the anonymous referee for suggestions that improved the clarity of the paper. S.D. also wishes to thank the KITP at UCSB for hosting the research program “Asteroseismology in the Space Age”, in which this work started. S.D. acknowledges support from the PNPS under the grant “Rotation interne et magnétisme des sousgéantes et géantes Kepler” and from the Centre National d’Études Spatiales (CNES). Funding for the Stellar Astrophysics Centre is provided by The Danish National Research Foundation (Grant agreement No.: DNRF106). The research is supported by the ASTERISK project (ASTERoseismic Investigations with SONG and Kepler) funded by the European Research Council (Grant agreement No.: 267864). R.M.O. acknowledges financial support by a CNES postdoctoral fellowship. S.B. acknowledges NASA grant NNX16AI09G and NSF grant AST1514676.
References
 Aerts, C., ChristensenDalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology (Springer Science + Business Media BV) [Google Scholar]
 Baglin, A., Auvergne, M., Barge, P., et al. 2006, in ESA SP, 1306, eds. M. Fridlund, A. Baglin, J. Lochard, & L. Conroy, 33 [Google Scholar]
 Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Beck, P. G., Hambleton, K., Vos, J., et al. 2014, A&A, 564, A36 [NASA ADS] [CrossRef] [EDP Sciences] [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]
 Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Campante, T. L., Handberg, R., Mathur, S., et al. 2011, A&A, 534, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chaplin, W. J., ChristensenDalsgaard, J., Elsworth, Y., et al. 1999, MNRAS, 308, 405 [NASA ADS] [CrossRef] [Google Scholar]
 Charbonnel, C., & Talon, S. 2005, Science, 309, 2189 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Deheuvels, S., Bruntt, H., Michel, E., et al. 2010, A&A, 515, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., García, R. A., Chaplin, W. J., et al. 2012, ApJ, 756, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Deheuvels, S., Doğan, G., Goupil, M. J., et al. 2014, A&A, 564, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Ballot, J., Beck, P. G., et al. 2015, A&A, 580, A96 [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]
 Dicke, R. H. 1982, Nature, 300, 693 [NASA ADS] [CrossRef] [Google Scholar]
 Dziembowski, W. A. 1971, Acta Astron., 21, 289 [NASA ADS] [Google Scholar]
 Dziembowski, W., & Goode, P. R. 1985, ApJ, 296, L27 [NASA ADS] [CrossRef] [Google Scholar]
 Dziembowski, W. A., & Goode, P. R. 1992, ApJ, 394, 670 [NASA ADS] [CrossRef] [Google Scholar]
 EffDarwich, A., & Korzennik, S. G. 2013, Sol. Phys., 287, 43 [NASA ADS] [CrossRef] [Google Scholar]
 Eggenberger, P., Montalbán, J., & Miglio, A. 2012, A&A, 544, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eggenberger, P., Lagarde, N., Miglio, A., et al. 2017, A&A, 599, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fuller, J., Lecoanet, D., Cantiello, M., & Brown, B. 2014, ApJ, 796, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Gough, D. O., & McIntyre, M. E. 1998, Nature, 394, 755 [NASA ADS] [CrossRef] [Google Scholar]
 Goupil, M. J. 2009, in The Rotation of Sun and Stars, eds. J.P. Rozelot, & C. Neiner (Berlin: Springer Verlag), Lect. Notes Phys., 765, 45 [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]
 Hasan, S. S., & ChristensenDalsgaard, J. 1992, ApJ, 396, 311 [NASA ADS] [CrossRef] [Google Scholar]
 Jouve, L., Gastine, T., & Lignières, F. 2015, A&A, 575, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kjeldsen, H., Bedding, T. R., Viskum, M., & Frandsen, S. 1995, AJ, 109, 1313 [NASA ADS] [CrossRef] [Google Scholar]
 Maeder, A., & Meynet, G. 2014, ApJ, 793, 123 [NASA ADS] [CrossRef] [Google Scholar]
 Marques, J. P., Goupil, M. J., Lebreton, Y., et al. 2013, A&A, 549, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Morel, P. 1997, A&AS, 124, 597 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Mosser, B., Barban, C., Montalbán, J., et al. 2011, A&A, 532, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012, A&A, 548, A10 [NASA ADS] [CrossRef] [EDP Sciences] [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., Pinçon, C., Belkacem, K., Takata, M., & Vrard, M. 2017, A&A, 600, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Osaki, J. 1975, PASJ, 27, 237 [NASA ADS] [Google Scholar]
 Ouazzani, R.M., & Goupil, M.J. 2012, A&A, 542, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ouazzani, R.M., Dupret, M.A., & Reese, D. R. 2012, A&A, 547, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reese, D., Lignières, F., & Rieutord, M. 2006, A&A, 455, 621 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rüdiger, G., Gellert, M., Spada, F., & Tereshin, I. 2015, A&A, 573, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Schou, J., Antia, H. M., Basu, S., et al. 1998, ApJ, 505, 390 [NASA ADS] [CrossRef] [Google Scholar]
 Scuflaire, R. 1974, A&A, 36, 107 [NASA ADS] [Google Scholar]
 Scuflaire, R., Montalbán, J., Théado, S., et al. 2008, Ap&SS, 316, 149 [NASA ADS] [CrossRef] [Google Scholar]
 Spada, F., Lanzafame, A. C., & Lanza, A. F. 2010, MNRAS, 404, 641 [NASA ADS] [CrossRef] [Google Scholar]
 Spruit, H. C. 1999, A&A, 349, 189 [NASA ADS] [Google Scholar]
 Suárez, J. C., Goupil, M. J., & Morel, P. 2006, A&A, 449, 673 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Talon, S., & Charbonnel, C. 2008, A&A, 482, 597 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Firstorder perturbative approach including neardegeneracy effects
Let us consider the case of two modes (indicated by subscripts a and b) with neardegenerate eigenfrequencies, that is,  ω_{0,a}−ω_{0,b}  ~ Ω. In this case, the eigenfunctions of the neardegenerate modes can be written as ξ = Aξ_{0,a} + Bξ_{0,b}. By injecting this expression into Eq. (3) and taking the inner product with ξ_{0,a} and ξ_{0,b} successively, the following system is obtained (A.1)where we have defined (A.2)To first order, one can write , where ω_{a} corresponds to the firstorder perturbed frequency of mode a in the absence of neardegeneracy (see Eq. (12)) and a similar expression holds for mode b. Equation (A.1) can thus be rewritten as (A.3)Solutions to Eq. (A.3) can exist only if (A.4)This yields the following expression for the frequencies of the neardegenerate modes (A.5)Equation (A.3) can also be used to obtain a relation between A and B, and therefore to estimate the eigenfunctions of the perturbed modes.
Appendix B: An approximate expression for ω_{1,ab}
We here obtain an approximate expression for the coupling term ω_{1,ab}, starting from Eq. (21), in order to interpret its behavior as a function of the rotation profile. For this purpose, we introduce the quantity (B.1)where K_{ab} is defined by Eq. (14). Taking into account the fact that the unperturbed eigenfunctions of modes a and b are orthogonal, Eq. (B.1) can be rewritten (B.2)In the gmode cavity, the horizontal displacement is much larger than the vertical displacement for all modes, so that the crossterms ξ_{r}ξ_{h} can be safely neglected in this region. In the pmode cavity, these crossterms dominate, but since the horizontal displacement is much smaller that the radial displacement, the overall contribution of the pmode cavity to γ is smaller than the contribution of the gmode cavity. We can thus approximately write(B.3)In the gmode cavity, considering that ξ_{h} ≫ ξ_{r}, the quantity γ_{c} (defined by Eq. (22)) can be approximated by Hence (B.6)Consequently, Eq. (21) can be rewritten as (B.7)In the particular case of a rotation profile for which the rotational splitting of pure g modes is identical to that of pure p modes, one has ⟨ Ω ⟩ _{c} = L^{2}/ (L^{2}−1) ⟨ Ω ⟩ _{e}, and it can be seen from Eq. (B.7) that ω_{1,ab} vanishes.
All Tables
Extracted parameters for the components of two l = 2 multiplets in avoided crossing in the spectrum of KIC 7341231.
Values of the parameters ⟨ Ω ⟩ _{c}, ⟨ Ω ⟩ _{e}, Δt, ζ_{0,a}, and ζ_{0,b} obtained when fitting Eq. (24) to the mode frequencies of acor for models A, C, D, and E, assuming a rotation profile following Eq. (15) with a core rotation ⟨ Ω ⟩ _{c}/ 2π = 710 nHz and coreenvelope contrasts of 100, 20, 10, 5, and 2.
All Figures
Fig. 1 Section of the oscillation spectrum of KIC 7341231 in the neighborhood of the 384μHz radial mode (light gray area). The dark gray area corresponds to a rotationallysplit l = 1 mode and the vertical arrows indicate the observed mode frequencies of two l = 2 mixed modes (see Table 1). The vertical dotted lines show the results of a fit of Eq. (24) to the observed mode frequencies. 

Open with DEXTER  
In the text 
Fig. 2 Variations in the frequencies of two l = 2 rotational multiplets during an avoided crossing. The sequence of models B0 through B8 was used, assuming the rotation profile given by Eq. (15) with a coreenvelope contrast of five. The frequencies of rotational multiplets were computed using the firstorder pertubative approach including neardegeneracy effects (open blue circles), and using the nonperturbative oscillation code acor (solid and dashed black curves). The filled red circles correspond to the frequencies obtained with the approximate Eq. (24) for model B4. 

Open with DEXTER  
In the text 
Fig. 3 Variations in the asymmetry of two l = 2 rotational multiplets during an avoided crossing, computed using Eq. (2) for m = ± 2 components. The symbols are identical to those of Fig. 2. 

Open with DEXTER  
In the text 
Fig. 4 Integrand of the mode inertia for the mcomponents of the higherfrequency l = 2 multiplet at the time during the avoided crossing where it shows the strongest asymmetry (t = t_{0} + 3 Myr). The radial and horizontal displacements (ξ_{r} and ξ_{h}) were computed using a perturbative approach as explained in Sect. 3 (gray dashed lines) and with the code acor (solid black lines). 

Open with DEXTER  
In the text 
Fig. 5 Degree of asymmetry δ_{asym} of the two l = 2 multiplets during the avoided crossing. Different rotation profiles were considered, with coreenvelope contrasts (Ω_{c}/ Ω_{e}) of 20 (purple circles), 10 (red triangle), 5 (blue squares), 1 (green downwardtriangle), and 0.5 (cyan stars). The higherfrequency multiplet is represented by filled symbols and the lowerfrequency multiplet by empty symbols. To guide the eye, we overplotted cubic spline interpolations δ_{asym}(t−t_{0}) for each rotation profile — solid (resp. dashed) lines for higherfrequency (resp. lowerfrequency) multiplet. 

Open with DEXTER  
In the text 
Fig. 6 Variations in the parameters β_{c}, β_{e}, γ_{c}, and γ_{e} (defined in Eqs. (17), (18), (22), and (23)) along the course of the avoided crossing between the two l = 2 multiplets for our reference sequence of models B0 through B8. Mode a (resp. b) corresponds to the lower (resp. higher) frequency multiplet. Core quantities are represented by open circles and envelope quantities by open squares. Interpolations using cubic splines are overplotted (blue dashed lines for core parameters and red dotdash lines for envelope parameters). 

Open with DEXTER  
In the text 
Fig. 7 Values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} obtained by fitting Eq. (24) to the mode frequencies of model D computed with acor, assuming a rotation profile with ⟨ Ω ⟩ _{c}/ 2π = 710 nHz and ⟨ Ω ⟩ _{c}/ ⟨ Ω ⟩ _{e} = 5. A random noise with an rms of 15 nHz was added to the mode frequencies and 1000 iterations were performed. The quantity β_{c} ⟨ Ω ⟩ _{c} + β_{e} ⟨ Ω ⟩ _{e} is conserved along the blue lines, where β_{c} and β_{e} correspond to mode a (dashed line) or mode b (dotdashed line). The red cross indicates the actual values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} computed with Eqs. (19) and (20). 

Open with DEXTER  
In the text 
Fig. 8 Distribution of the values of ⟨ Ω ⟩ _{c} (left) and ⟨ Ω ⟩ _{e} (right) obtained by fitting Eq. (24) to the mode frequencies of acor for model D, assuming a rotation profile with ⟨ Ω ⟩ _{c}/ 2π = 710 nHz and ⟨ Ω ⟩ _{c}/ ⟨ Ω ⟩ _{e} = 5. A random noise with an rms of 15 nHz (top), 30 nHz (middle), and 50 nHz (bottom) was added to the mode frequencies. The colors correspond to the different scenarios that were tested: ideal case where all the components are detected (black), inclination angle of 90° (purple) or 30° (red), and case of an overlap with the closest radial mode (blue). The vertical dashed lines indicate the actual values of ⟨ Ω ⟩ _{c} and ⟨ Ω ⟩ _{e} computed with Eqs. (19) and (20). 

Open with DEXTER  
In the text 
Fig. 9 Distribution of the values of ζ_{0,a} recovered when fitting Eq. (24) to the mode frequencies of acor for model D, assuming a rotation profile with ⟨ Ω ⟩ _{c}/ 2π = 710 nHz and ⟨ Ω ⟩ _{c}/ ⟨ Ω ⟩ _{e} = 5. The lines and colors have the same meaning as in Fig. 8. The vertical dashed lines indicate the actual value of ζ_{0,a} computed with the eigenfunctions of the model and Eq. (27). 

Open with DEXTER  
In the text 