Highprecision acoustic helium signatures in 18 lowmass lowluminosity red giants
Analysis from more than four years of Kepler observations^{⋆}
^{1} Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain
^{2} Universidad de La Laguna, Departamento de Astrofísica, 38206 La Laguna, Tenerife, Spain
^{3} Laboratoire AIM, CEA/DSM – CNRS – Univ. Paris Diderot – IRFU/SAp, Centre de Saclay, 91191 GifsurYvette Cedex, France
email: enrico.corsaro@cea.fr
^{4} Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
Received: 18 February 2015
Accepted: 13 April 2015
Context. Highprecision frequencies of acoustic modes in red giant stars are now available thanks to the long observing length and high quality of the light curves provided by the NASA Kepler mission, thus allowing the interior of evolved cool lowmass stars to be probed with an unprecedented level of detail.
Aims. We characterize the acoustic signature of the helium second ionization zone in a sample of 18 lowmass lowluminosity red giants by exploiting new modefrequency measurements derived from more than four years of Kepler observations.
Methods. We analyzed the second frequency differences of radial acoustic modes in all the stars of the sample by using the Bayesian code Diamonds.
Results. We find clear acoustic glitches due to the signature of helium second ionization in all the stars of the sample. We could measure the acoustic depth and the characteristic width of the acoustic glitches with a precision level on average around ~2% and ~8%, respectively. We find good agreement with theoretical predictions and existing measurements from the literature. Finally, we derive the amplitude of the glitch signal at ν_{max} for the second differences and for the frequencies with an average precision of ~6%, obtaining values in the range 0.14−0.24 μHz and 0.08−0.33 μHz, respectively, which can be used to investigate the helium abundance in the stars.
Key words: asteroseismology / stars: oscillations / stars: latetype / stars: interiors / methods: statistical / methods: data analysis
Appendix A is available in electronic form at http://www.aanda.org
© ESO, 2015
1. Introduction
The socalled acoustic glitches are regions of sharpstructure variation located in the interior of stars and caused by the presence of a change in the energy transport from radiative to convective, by a rapid variation in the chemical composition, or by ionization zones of chemical elements, such as hydrogen and helium. As originally predicted for the Sun (e.g., Vorontsov 1988; Gough 1990), these regions produce tiny and regular variations in the frequency of the acoustic (p) modes that can be detected by direct measurement of the characteristic large frequency separation, namely the frequency separation between modes having the same angular degree.
By studying the glitch signature in the Sun, it was possible to measure the acoustic position of the base of the convective zone and of the helium second ionization (He ii) zone. The signature can also be used to provide estimates of the helium abundance in the envelope and the extent of the overshooting (e.g., Basu & Antia 1995; Basu 1997; Monteiro & Thompson 2005; ChristensenDalsgaard et al. 2011). The same effect was expected to be observed in distant stars (e.g., Monteiro et al. 2000; Mazumdar & Antia 2001; Ballot et al. 2004). Thanks to the advent of the CoRoT (Baglin et al. 2006) and Kepler space missions (Borucki et al. 2010; Koch et al. 2010), an outstanding number of highquality photometric observations for thousands of stars has been released. This yielded the frequency shifts caused by the acoustic glitches to be discovered and analyzed in many lowmass, mainsequence, subgiant and redgiant stars (RGs) (Mosser et al. 2010; Miglio et al. 2010; Mazumdar et al. 2012, 2014; Verma et al. 2014), allowing the positions of the base of the convective zone, of the He ii zone in mainsequence and subgiant stars, and of the He ii zone in the case of the red giants to be constrained.
The asteroseismology of red giant stars, in particular, has led to several important breakthroughs in the stellar physics of lowmass stars in recent years (e.g. Beck et al. 2011; Mosser et al. 2011a; Bedding et al. 2011; Beck et al. 2012; Deheuvels et al. 2012). The characterization of the glitch signatures is able to provide tighter constraints on the chemical composition and the internal stratification of the star, and potentially allows retrieving helium abundances in distant stars, essential for population studies (e.g., see Broomhall et al. 2014, hereafter B14, and references therein). More recent studies focusing on these evolved cool stars have analyzed the glitches due to the He ii zone for more than a hundred targets observed by Kepler (Vrard et al. 2014), and thoroughly investigated the properties of the signature from a theoretical point of view (B14, see also ChristensenDalsgaard et al. 2014, for more discussion).
The recent availability of Kepler datasets spanning more than four years of nearly continuous observations, coupled with the development of new computational advances in asteroseismic data analysis (e.g., see Corsaro & De Ridder 2014; Corsaro et al. 2015, hereafter C15), enables the study of the acoustic glitch signatures in red giants with an unprecedented level of detail.
In this paper, we report on the evidence of clear acoustic glitch signatures due to the He ii zone in the sample of 19 red giants recently investigated by C15, so we fully characterize the oscillatory signal by means of a Bayesian approach.
2. Data analysis
As noticed by C15, the lowmass lowluminosity red giants (LRGs) are candidates that are well suited to testing stellar structure models and stellar evolution theory. The lessevolved stage in the red giant branch (RGB) of the stellar evolution for the LRGs implies the highest frequency of maximum power ν_{max} for a red giant (between 100 and 200 μHz), hence a broader power excess caused by the oscillations and, consequently, a larger number of radial orders observed (in general between six and nine). By having a larger number of high signaltonoise ratio pmode frequencies available, one is thus able to constrain the signature of the glitches more efficiently.
In this work we analyze the sample of LRGs studied by C15, who fitted and characterized their full oscillation spectrum using Kepler observations from Q0 till Q17.1, a total of ~1470 days, with a frequency resolution of δν_{bin} ≃ 0.008 μHz. The stars have ν_{max} values ranging from 110 to ~190 μHz and masses in the interval 1−2 M_{⊙}. We adopt the highprecision individual frequency measurements from C15, and follow the theoretical approach of B14.
In the present analysis we refer to the first (frequency) difference as the large frequency separation of a given angular degree, which is a function of the frequency in the power spectral density of the star, . For a radial order n, is thus defined as (1)where ν_{n,ℓ} is the central frequency of the mode with angular degree ℓ and radial order n. In addition, we compute the second (frequency) difference (see, e.g., Gough 1990), , defined for a single radial order as (2)We fit the acoustic glitch signatures with the model introduced by Houdek & Gough (2007) and used by B14 for RGs, defined as (3)with ω_{n,ℓ} ≡ 2πν_{n,ℓ} and Δ_{2}ω_{n,ℓ} ≡ 2πΔ_{2}ν_{n,ℓ}, A a dimensionless amplitude of the signature signal, τ_{He} the acoustic depth of the He ii zone, b its characteristic width, and φ and c are constant phase shift and offset, respectively, of the oscillatory signal. Following the arguments discussed by B14, we apply the fit to the second differences only because they are less prone to additional varying components, such as hydrogen ionization and nonadiabatic processes, and to the general frequency dependence of the large separation caused by the development of the secondorder term of the asymptotic relation (Mosser et al. 2011b). The second differences are at the same time available in a reasonably high number of measurements (two less than the total number of modes obtained for a given angular degree), still allowing the model parameters to be constrained without leading to degeneracies in the solutions.
Despite the possible presence of the oscillatory component in modes of angular degree ℓ> 0, we point out that only radial mode frequencies are used for the final fit. The reason behind this choice is mainly the need to exploit pure pmode character oscillations (see also B14), which in the case of RGs, are only available in the form of ℓ = 0 modes. This is because the coupling occurring between p modes of angular degree ℓ > 0 and g modes arising from the radiative interior can hamper the asymptotic behavior of the corresponding modes by producing socalled mixed modes, whose frequencies deviate from the expected position of a pure p mode oscillation (e.g., Beck et al. 2011).
We perform all the fits following a Bayesian approach using Diamonds (Corsaro & De Ridder 2014), hence exploiting a nestedsampling Monte Carlo method to perform the inference and estimate the free parameters of the model given by Eq. (3) from their individual marginal probability distributions (see Corsaro & De Ridder 2014, for more details). The configuring parameters of Diamonds (following the definitions by Corsaro & De Ridder 2014) used for all the computations are initial enlargement fraction 1.0 ≤ f_{0} ≤ 1.7, shrinking rate α = 0.020.03, number of live points N_{live} = 1000, number of clusters 1 ≤ N_{clust} ≤ 4, number of total drawing attempts M_{attempts} = 10^{4}, number of nested iterations before the first clustering M_{init} = 1000, and number of nested iterations with the same clustering M_{same} = 50.
For this analysis we adopted a normal likelihood function, such as the one used by Corsaro et al. (2013), which takes the uncertainties into account, with correlations included, on the measurements of the second differences. This assumes that the residuals arising from the difference between predicted and measured second differences are Gaussiandistributed. Since we did not have any initial guesses available from the literature for the given stars for the set up of priors, we used uniform (i.e., flat) prior probability distributions for all the free parameters of Eq. (3), with lower and upper boundaries for each parameter range obtained by comparison with existing measurements of the acoustic depths derived by Miglio et al. (2010), Mazumdar et al. (2014), and the theoretical results by B14 in the observed range of ν_{max}. The choice of uniform priors also yielded a faster computation with Diamonds, as already discussed by Corsaro & De Ridder (2014) and C15.
Following the discussion by Ballot et al. (2004) and B14, we computed the acoustic radius of the He ii zones, t_{He II}, since it represents a quasi unbiased estimator of the acoustic position of the glitch. This is done by using the mean large frequency separation, ⟨ Δν ⟩, obtained from the radial mode frequencies provided by C15, giving the total acoustic radius of the star, T = (2 ⟨ Δν ⟩ )^{1}, hence the acoustic radius of the He ii zone, t_{He II} ≡ T−τ_{He II}.
Finally, to provide measurements that can be used to model the helium abundance in the envelope of the stars, following B14, we extracted the amplitude of the signal at ν_{max} from Eq. (3), obtaining (4)with ω_{max} = 2πν_{max}, and ν_{max} derived from the background fit done by C15. Following Verma et al. (2014), we also derived the amplitude of the signal in the frequencies, A_{He}, given as (5)where ⟨ Δν ⟩ is the same mean large frequency separation used to calculate t_{He II}, and τ_{He II} is the same acoustic depth used in Eq. (3). The adoption of A_{He} to retrieve the helium content in the envelope was preferred since this value is not influenced by the location of the glitch (see, e.g., Mazumdar & Antia 2001; Mazumdar et al. 2014; Verma et al. 2014). For clarity, we stress that A_{max} is derived from a Bayesian approach by using the same sampling of the posterior probability distribution obtained by Diamonds for the free parameters of the glitch model (see also Corsaro & De Ridder 2014, Fig. 7, for an analogous case presenting the sampling from Diamonds for a combination of different inferred parameters). The parameters t_{He II} and A_{He} simply follow from their definitions as presented above, by using both the value ⟨ Δν ⟩ , computed from the radial mode frequencies presented by C15 for each star of our sample, and the estimated model parameters of Eq. (3) (see Sect. 3 for more details).
3. Results
The results for and for the star KIC 12008916 are shown in Fig. 1 (top and bottom panels, respectively) and can be found in Appendix A for all the other LRGs, together with the tables with the individual measurements of the radial angular frequencies ω_{n,0}, and of the corresponding second angular frequency differences used in this work. We discarded KIC 10123207 from the fit because of the low number of available measurements (foursecond differences only, one less than the minimum required to fit the model given by Eq. (3)). To help the reader visualize the presence of the oscillatory signal in the first differences of the angular degrees ℓ = 0,2,3 and in the second differences of the angular degrees ℓ = 2,3, we included lowdegree (3−4) polynomial fits. The 1σ uncertainties on the first and second differences derived from a standard error propagation of the uncertainties of the individual mode frequencies, following Eqs. (1) and (2), respectively, are overlaid in each plot, though they are not visible in most of the cases because they are smaller than the size of the symbols used for the measurements. The uncertainties on all the measurements are listed in the corresponding tables in Appendix A for completeness.
Fig. 1 Acoustic glitches for KIC 12008916. Top panel: the first difference, , Eq. (1). Blue squares represent values computed from ℓ = 0 modes, while green triangles from ℓ = 2 modes, with polynomial fits (dashed lines with same color as symbols) overlaid to visualize the oscillatory trend. The dashed vertical line marks the position of ν_{max} for a reference. Bottom panel: the second difference , Eq. (2), with the same symbol description as for the top panel. The solid red line indicates the fit to the ℓ = 0 measurements given by Eq. (3) with the estimated parameters listed in Table 1, as derived by Diamonds. The horizontal dotted line denotes the offset level c/ (2π), which is useful for visualizing the amplitude of the signature. The inset shows a zoom in of one of the measurements to visualize the precision level of the fit. 

Open with DEXTER 
We find that all the stars we analyzed have clear acoustic glitches due to the signature of the He ii zone in ℓ = 0 and 2 modes up to the second differences. We can also see the presence of acoustic glitches in ℓ = 3 modes for most of the stars thanks to at least four different frequency measurements that are available. As mentioned in Sect. 2, we find that the measurements for modes having angular degrees ℓ = 2,3 often deviate from those of the radial oscillations (see, e.g. Figs., A.6, A.8, A.10, and A.13). As also indicated by C15 for the case of the mode linewidths of the ℓ = 2 modes, this different behavior relies on the presence of both mixed quadrupole modes and rotational split components. When using an individual Lorentzian profile to fit the frequency region containing the oscillation peak, either an ℓ = 2 or 3 mode, as done by C15, the effects mentioned before can significantly change the measured frequency of the peak. A reliable treatment of the mixed modes and of the rotational split components for ℓ = 2,3 modes is, however, difficult due to the high proximity of the individual peaks.
The model fit to the acoustic glitch signatures of KIC 12008916 is shown in the bottom panel of Fig. 1 for the case of , and similarly for the other stars in Appendix A. All the estimated parameters of Eq. (3) are provided in Table 1. The inset in the bottom panel of Fig. 1 provides a closer view of one of the measurements for visualizing the precisionlevel achieved in the fit. In particular, we find that the model given by Eq. (3) yields a remarkable fit quality for most of the stars, with average uncertainties of ~2% for τ_{He II} and ~8% for b. Following the analysis presented by Corsaro et al. (2013), we obtained the weighted Gaussian rms of the residuals, σ_{rms} (listed in Table 1 as well).
To computate the rms, we adopted the weights , σ_{i}, which are the uncertainties on the second frequency difference coming from those reported in Tables A.1. The quantity σ_{rms} provides additional information to the reader because it allows the quality of the fits to be compared between different stars and the precision achieved on the individual fits to the given uncertainties of the data points to be related. We note that for all the fits presented in this work, the values for σ_{rms} are remarkably low, ranging from 10^{2} down to 10^{3}μHz in the best cases, thus in many cases reaching the same order of the precision level obtained on the individual frequencies of the radial modes. For a reference to the reader, in Table 1 we also provide the values for the total acoustic radius T, with its 1σ standard deviation and the values of ν_{max} obtained by C15.
The stars KIC 8475025, KIC 9145955, KIC 10200377, and KIC 11913545 each show a component at high frequency that is not properly predicted by the adopted model. This mainly relies on some residual frequency dependence of the second differences that becomes more pronounced towards the wings of the region containing the oscillations. However, we note that the measurements at higher frequencies all have larger error bars (up to 10 times) with respect to the others, because of the larger linewidths of the peaks occurring at high frequency (see C15 for more details). The fits derived, except for KIC 8475025, are therefore not significantly affected by the measurements of the second differences falling at highfrequency, whereas they are almost entirely constrained by those close to ν_{max}. This is inspected by refitting the glitch model without the highest frequency measurement of the second difference (showing the deviating behavior), hence noticing that the new estimated free parameters of the model lie well within the uncertainties of those reported in Table 1. For KIC 8475025, however, we find that the fit is unstable owing to the large deviation (more than 0.2 μHz) of the second difference measurement falling at the highest frequency (see Fig. A.7). This is because the measured oscillation frequency of the highest frequency radial mode is likely to be affected by additional sources, such as mixed modes and rotational split components arising from the neighbor ℓ = 2 mode, which are enhanced by the large mode linewidths (see C15 for more details). To stabilize the fit for this star and provide estimates of the model parameters that are comparable to the other stars in the sample, we have therefore chosen to discard the last measurement of the second difference for this particular target. In the case of KIC 8366239, KIC 9267654, and KIC 10200377, the highestfrequency values are marked as unreliable, according to the Bayesian peak significance test done by C15, although they were included in the fit since they do not produce any significant change in the results for the same reasons discussed above.
The measurements of the acoustic radius t_{He II}, the amplitude of the signal A_{max} from Eq. (4), and the corresponding characteristic width b, are shown in Fig. 2 (top, middle, and bottom panels, respectively) for all the stars in the sample. We note that while the model parameters (Eq. (3)) and their corresponding 68% Bayesian credible intervals are estimated by means of Diamonds (see Table 1 and Corsaro & De Ridder 2014 for more details on the derivation of the Bayesian uncertainties), the uncertainties for the additional parameters t_{He II}, A_{max}, and A_{He} were obtained in a subsequent step. In particular, for A_{max}, we used the same sampling of the posterior probability distribution as obtained by Diamonds, so we have derived the median and the corresponding 68% Bayesian credible intervals directly from the marginal probability distributions of A_{max}. For t_{He II} and A_{He}, the uncertainties follow from those of the acoustic depth through the definition of the acoustic radius, and by a rescaling of the uncertainties on A_{max} through Eq. (5), respectively. All the resulting values are listed in Table 1 as well. We note that the precision obtained on our measurements of the acoustic radii of the He ii zones is about ten times higher than obtained by Miglio et al. (2010) using CoRoT data. In addition, all the values match those predicted by B14 along the entire range of ν_{max} investigated, showing a clear increasing trend towards lower ν_{max}, as expected for more evolved stages of the evolution in the RGB. The derived amplitudes in frequency, A_{He}, are within the range 0.08−0.33 μHz, and are varying from star to star with uncertainties on average around ~6%, thus opening the possibility of studying the He abundance by direct comparison with stellar models.
Fig. 2 Acoustic radius of the He ii zone (top panel) and corresponding amplitude of the oscillatory signal at ν_{max} (middle panel), and characteristic width b (bottom panel), as a function of ν_{max} for all the stars of the sample. The 68% Bayesian uncertainties listed in Table 1 are overlaid for all the measurements. 

Open with DEXTER 
4. Conclusions
By exploiting the set of individual mode frequencies extracted by C15 for a sample of 19 LRGs with a precision level up to 10^{3} μHz, we computed the first differences, Eq. (1), and the second differences, Eq. (2), for presenting the results on the acoustic glitches of these stars. In this work, we have shown that the acoustic glitches are remarkably clear for all the red giants of the sample, and for both ℓ = 0 and ℓ = 2 modes, up to the second difference (where five to seven different measurements are available for each star, except KIC 10123207, which instead has only four and was not considered in the analysis), with many cases involving ℓ = 3 modes as well.
By adopting the model for the second differences introduced by Houdek & Gough (2007), Eq. (3), we have extracted the acoustic depth, the characteristic width, and amplitude of the signal of all the He ii zones of the stars in our sample (see Table 1) with an unprecedented level of detail for red giant stars (on average ~2% for the acoustic depths, ~8% for the characteristic widths, and ~6% for the amplitudes of the glitch signal in both the second difference and frequency), improved by about one order of magnitude with respect to existing measurements of acoustic depths from the literature. These values, reflecting the highprecision obtained on the individual frequency measurements of the radial modes, confirm the theoretical predictions done by B14 in the same range of ν_{max}. We also stress that the glitch model given by Houdek & Gough (2007) is able to predict the observations very exhaustively for most of the stars (well within the given uncertainties of the measurements for most of the data points available, as shown in Fig. 1 and in the other
figures in the Appendix). This is also supported by our computation of σ_{rms}, listed in Table 1, which is at the same precision level as given by the measurements of the second difference for most of the stars analyzed.
Finally, the set of values for A_{max} and A_{He} derived in this work, where the latter are not influenced by the position of the glitch, coupled with the high precision achieved, will be useful for investigating the helium content in the envelope of the stars and possibly contributing to study scenarios of helium enrichment in lowmass stars (e.g., see Gratton et al. 2012, and references therein).
Acknowledgments
E.C. is funded by the European Community’s Seventh Framework Programme (FP7/20072013) under grant agreement N° 312844 (SPACEINN). The research leading to these results has received funding from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007–2013) ERC grant agreement N° 227224 (PROSPERITY), from the Fund for Scientific Research of Flanders (G.0728.11). E.C. thanks A.M. Broomhall and A. Miglio for useful discussions.
References
 Baglin, A., Michel, E., Auvergne, M., & COROT Team. 2006, in Proc. of SOHO 18/GONG 2006/HELAS I, Beyond the spherical Sun, ESA SP, 624, 34 [Google Scholar]
 Ballot, J., TurckChièze, S., & García, R. A. 2004, A&A, 423, 1051 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Basu, S. 1997, MNRAS, 288, 572 [NASA ADS] [Google Scholar]
 Basu, S., & Antia, H. M. 1995, MNRAS, 276, 1402 [NASA ADS] [Google Scholar]
 Beck, P. G., Bedding, T. R., Mosser, B., et al. 2011, Science, 332, 205 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Bedding, T. R., Mosser, B., Huber, D., et al. 2011, Nature, 471, 608 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Broomhall, A.M., Miglio, A., Montalbán, J., et al. 2014, MNRAS, 440, 1828 (B14) [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J., Monteiro, M. J. P. F. G., Rempel, M., & Thompson, M. J. 2011, MNRAS, 414, 1158 [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J., Silva Aguirre, V., Elsworth, Y., & Hekker, S. 2014, MNRAS, 445, 3685 [NASA ADS] [CrossRef] [Google Scholar]
 Corsaro, E., & De Ridder, J. 2014, A&A, 571, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Corsaro, E., Fröhlich, H.E., Bonanno, A., et al. 2013, MNRAS, 430, 2313 [NASA ADS] [CrossRef] [Google Scholar]
 Corsaro, E., De Ridder, J., & García, R. A. 2015, A&A, in press DOI: 10.1051/00046361/201525895 [Google Scholar]
 Deheuvels, S., García, R. A., Chaplin, W. J., et al. 2012, ApJ, 756, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Gough, D. O. 1990, in Lect. Notes Phys., 367, Progress of Seismology of the Sun and Stars, eds. Y. Osaki, & H. Shibahashi (Berlin: Springer Verlag), 283 [Google Scholar]
 Gratton, R. G., Carretta, E., & Bragaglia, A. 2012, A&ARv, 20, 50 [NASA ADS] [CrossRef] [Google Scholar]
 Houdek, G., & Gough, D. O. 2007, MNRAS, 375, 861 [NASA ADS] [CrossRef] [Google Scholar]
 Koch, D. G., Borucki, W. J., Basri, G., et al. 2010, ApJ, 713, L79 [NASA ADS] [CrossRef] [Google Scholar]
 Mazumdar, A., & Antia, H. M. 2001, A&A, 368, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mazumdar, A., Michel, E., Antia, H. M., & Deheuvels, S. 2012, A&A, 540, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mazumdar, A., Monteiro, M. J. P. F. G., Ballot, J., et al. 2014, ApJ, 782, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Miglio, A., Montalbán, J., Carrier, F., et al. 2010, A&A, 520, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Monteiro, M. J. P. F. G., & Thompson, M. J. 2005, MNRAS, 361, 1187 [NASA ADS] [CrossRef] [Google Scholar]
 Monteiro, M. J. P. F. G., ChristensenDalsgaard, J., & Thompson, M. J. 2000, MNRAS, 316, 165 [NASA ADS] [CrossRef] [Google Scholar]
 Mosser, B., Belkacem, K., Goupil, M.J., et al. 2010, A&A, 517, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Barban, C., Montalbán, J., et al. 2011a, A&A, 532, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Belkacem, K., Goupil, M. J., et al. 2011b, A&A, 525, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verma, K., Faria, J. P., Antia, H. M., et al. 2014, ApJ, 790, 138 [NASA ADS] [CrossRef] [Google Scholar]
 Vorontsov, S. V. 1988, in Advances in Helio and Asteroseismology, eds. J. ChristensenDalsgaard & S. Frandsen, IAU Symp., 123, 151 [Google Scholar]
 Vrard, M., Mosser, B., & Barban, C. 2014, in SF2A2014: Proc. Annual meeting of the French Society of Astronomy and Astrophysics, eds. J. Ballet, F. Martins, F. Bournaud, R. Monier, & C. Reylé, 531 [Google Scholar]
Online material
Appendix A: Results for the fitting the acoustic signatures
All the individual angular frequency measurements and corresponding second differences of the radial modes are listed in the Tables A.1. The results for the first differences and second differences of all the LRGs are shown in the Figs. A.1 to A.17, similar to those provided for KIC 12008916 in Fig. 1.
Fig. A.1 Same description as in Fig. 1 but for KIC 3744043, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. 

Open with DEXTER 
Fig. A.2 Same description as in Fig. 1 but for KIC 6117517, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER 
Fig. A.3 Same description as in Fig. 1 but for KIC 6144777, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. 

Open with DEXTER 
Fig. A.4 Same description as in Fig. 1 but for KIC 7060732. 

Open with DEXTER 
Fig. A.5 Same description as in Fig. 1 but for KIC 7619745. 

Open with DEXTER 
Fig. A.6 Same description as in Fig. 1 but for KIC 8366239. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER 
Fig. A.7 Same description as in Fig. 1 but for KIC 8475025. 

Open with DEXTER 
Fig. A.8 Same description as in Fig. 1 but for KIC 8718745, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. 

Open with DEXTER 
Fig. A.9 Same description as in Fig. 1 but for KIC 9145955, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER 
Fig. A.10 Same description as in Fig. 1 but for KIC 9267654. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER 
Fig. A.11 Same description as in Fig. 1 but for KIC 9475697, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER 
Fig. A.12 Same description as in Fig. 1 but for KIC 9882316. 

Open with DEXTER 
Fig. A.13 Same description as in Fig. 1 but for KIC 10200377, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER 
Fig. A.14 Same description as in Fig. 1 but for KIC 10257278, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER 
Fig. A.15 Same description as in Fig. 1 but for KIC 11353313, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER 
Fig. A.16 Same description as in Fig. 1 but for KIC 11913545. 

Open with DEXTER 
Fig. A.17 Same description as in Fig. 1 but for KIC 11968334, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER 
Radial angular frequencies ω_{n,0} and second angular frequency differences Δ_{2}ω_{n,0} for all the stars of the sample.
All Tables
Radial angular frequencies ω_{n,0} and second angular frequency differences Δ_{2}ω_{n,0} for all the stars of the sample.
All Figures
Fig. 1 Acoustic glitches for KIC 12008916. Top panel: the first difference, , Eq. (1). Blue squares represent values computed from ℓ = 0 modes, while green triangles from ℓ = 2 modes, with polynomial fits (dashed lines with same color as symbols) overlaid to visualize the oscillatory trend. The dashed vertical line marks the position of ν_{max} for a reference. Bottom panel: the second difference , Eq. (2), with the same symbol description as for the top panel. The solid red line indicates the fit to the ℓ = 0 measurements given by Eq. (3) with the estimated parameters listed in Table 1, as derived by Diamonds. The horizontal dotted line denotes the offset level c/ (2π), which is useful for visualizing the amplitude of the signature. The inset shows a zoom in of one of the measurements to visualize the precision level of the fit. 

Open with DEXTER  
In the text 
Fig. 2 Acoustic radius of the He ii zone (top panel) and corresponding amplitude of the oscillatory signal at ν_{max} (middle panel), and characteristic width b (bottom panel), as a function of ν_{max} for all the stars of the sample. The 68% Bayesian uncertainties listed in Table 1 are overlaid for all the measurements. 

Open with DEXTER  
In the text 
Fig. A.1 Same description as in Fig. 1 but for KIC 3744043, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. 

Open with DEXTER  
In the text 
Fig. A.2 Same description as in Fig. 1 but for KIC 6117517, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER  
In the text 
Fig. A.3 Same description as in Fig. 1 but for KIC 6144777, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. 

Open with DEXTER  
In the text 
Fig. A.4 Same description as in Fig. 1 but for KIC 7060732. 

Open with DEXTER  
In the text 
Fig. A.5 Same description as in Fig. 1 but for KIC 7619745. 

Open with DEXTER  
In the text 
Fig. A.6 Same description as in Fig. 1 but for KIC 8366239. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER  
In the text 
Fig. A.7 Same description as in Fig. 1 but for KIC 8475025. 

Open with DEXTER  
In the text 
Fig. A.8 Same description as in Fig. 1 but for KIC 8718745, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. 

Open with DEXTER  
In the text 
Fig. A.9 Same description as in Fig. 1 but for KIC 9145955, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER  
In the text 
Fig. A.10 Same description as in Fig. 1 but for KIC 9267654. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER  
In the text 
Fig. A.11 Same description as in Fig. 1 but for KIC 9475697, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER  
In the text 
Fig. A.12 Same description as in Fig. 1 but for KIC 9882316. 

Open with DEXTER  
In the text 
Fig. A.13 Same description as in Fig. 1 but for KIC 10200377, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER  
In the text 
Fig. A.14 Same description as in Fig. 1 but for KIC 10257278, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER  
In the text 
Fig. A.15 Same description as in Fig. 1 but for KIC 11353313, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER  
In the text 
Fig. A.16 Same description as in Fig. 1 but for KIC 11913545. 

Open with DEXTER  
In the text 
Fig. A.17 Same description as in Fig. 1 but for KIC 11968334, with yellow star from ℓ = 3 modes and corresponding polynomial fit with same color. Open symbols represent measurements that used modes with detection probability under the threshold suggested by C15. 

Open with DEXTER  
In the text 