Issue 
A&A
Volume 589, May 2016



Article Number  A93  
Number of page(s)  26  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201527967  
Published online  19 April 2016 
Measuring the extent of convective cores in lowmass stars using Kepler data: toward a calibration of core overshooting
^{1} Université de Toulouse, UPSOMP, IRAP, 31400 Toulouse, France
email: sebastien.deheuvels@irap.omp.eu
^{2} CNRS, IRAP, 14 avenue Édouard Belin, 31400 Toulouse, France
^{3} Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, 4150762 Porto, Portugal
^{4} Centro de Astrofísica e Faculdade de Ciências, Universidade do Porto, Rua das Estrelas, 4150762 Porto, Portugal
^{5} Stellar Astrophysics Centre, Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
^{6} LESIA, CNRS UMR 8109, Université Pierre et Marie Curie, Université Denis Diderot, Observatoire de Paris, 92195 Meudon, France
^{7} Observatoire de Paris, GEPI, CNRS UMR 8111, 92195 Meudon, France
^{8} Institut de Physique de Rennes, Université de Rennes 1, CNRS UMR 6251, 35042 Rennes, France
^{9} Univ. ParisSud, Institut d’Astrophysique Spatiale, UMR 8617, CNRS, Bâtiment 121, 91405 Orsay Cedex, France
Received: 15 December 2015
Accepted: 1 March 2016
Context. Our poor understanding of the boundaries of convective cores generates large uncertainties on the extent of these cores and thus on stellar ages. The detection and precise characterization of solarlike oscillations in hundreds of mainsequence stars by CoRoT and Kepler has given the opportunity to revisit this problem.
Aims. Our aim is to use asteroseismology to consistently measure the extent of convective cores in a sample of mainsequence stars whose masses lie around the mass limit for having a convective core.
Methods. We first tested and validated a seismic diagnostic that was proposed to probe the extent of convective cores in a modeldependent way using the socalled r_{010} ratios, which are built with l = 0 and l = 1 modes. We applied this procedure to 24 lowmass stars chosen among Kepler targets to optimize the efficiency of this diagnostic. For this purpose, we computed grids of stellar models with both the Cesam2k and mesa evolution codes, where the extensions of convective cores were modeled either by an instantaneous mixing or as a diffusion process.
Results. We found that 10 stars in our sample are in fact subgiants. Among the other targets, were able to unambiguously detect convective cores in eight stars, and we obtained seismic measurements of the extent of the mixed core in these targets with a good agreement between the Cesam2k and mesa codes. By performing optimizations using the LevenbergMarquardt algorithm, we then obtained estimates of the amount of extra mixing beyond the core that is required in Cesam2k to reproduce seismic observations for these eight stars, and we showed that this can be used to propose a calibration of this quantity. This calibration depends on the prescription chosen for the extra mixing, but we found that it should also be valid for the code mesa, provided the same prescription is used.
Conclusions. This study constitutes a first step toward calibrating the extension of convective cores in lowmass stars, which will help reduce the uncertainties on the ages of these stars.
Key words: asteroseismology / convection / stars: evolution / stars: interiors
© ESO, 2016
1. Introduction
The extent of chemically mixed regions associated to stellar convective cores is notoriously uncertain. Several physical processes that remain challenging to describe theoretically are known to extend convective cores beyond the theoretical Schwarzschild limit. The most often cited among them is core overshooting. According to Schwarzschild’s criterion, the boundary of a convective core corresponds to the layer above which upwardmoving convective blobs are braked. However, this criterion neglects the inertia of the ascending blobs, which are expected to penetrate the radiative zone over a certain distance (overshoot). The theoretical complexity of this phenomenon is illustrated well by the large number of developments that were proposed to describe it (e.g., Saslaw & Schwarzschild 1965; Shaviv & Salpeter 1971; Roxburgh 1978; Zahn 1991 to quote only a few) and by the diversity of the predicted distances d_{ov} over which convective eddies are expected to overshoot in the stable region (predicted values for d_{ov} range from 0 to 2 H_{P}, where H_{P} is the local pressure scale height). Current numerical simulations of overshooting are encouraging, but they are still far from reproducing the very high turbulence of stellar convection and cannot be used yet to obtain reliable prescriptions for core overshooting (see Dintrans 2009 for a review). Another complication arises from the fact that convective cores can also be extended due to rotationallyinduced mixing (see Maeder 2009 and references therein). As a result, it remains to be determined (1) over which distance convective cores are extended; (2) what the temperature stratification is like in these core extensions; and (3) how chemical elements are mixed in these regions.
Since convective cores constitute reservoirs for nuclear reactions, the uncertainty on their sizes generates significant uncertainties on stellar ages, especially near the end of the main sequence (MS). For instance, Lebreton et al. (2014) estimated that an extension of convective cores over a typical distance of 0.2 H_{P} can generate errors on stellar ages as large as 30% at the turnoff. It also affects the isochrones that have turnoff masses above ~ 1.1 M_{⊙}, hence the ages of rather young clusters.
To account for the combined effects of core overshooting and rotational mixing, 1D stellar models often consider an ad hoc extra mixing at the edge of the convective core, which is either modeled as an instantaneous mixing (simple extension of the mixed core) or as a diffusion process (Ventura et al. 1998), i.e., as a noninstantaneous mixing (see Noels et al. 2010 for a review). In both cases, the extent of the extra mixing (usually known as the overshooting distance d_{ov} even though overshooting may not be the only mechanism at work) depends on one free parameter. These models are clearly overly simplistic, but current observations have not yet permitted more complex models to be constrained. The overshooting distance has been observationally constrained by fitting isochrones to the colormagnitude diagrams of open clusters (e.g., Maeder & Mermilliod 1981; VandenBerg et al. 2006) and by performing calibrations using eclipsing binaries (e.g., Claret 2007; Stancliffe et al. 2015). These studies typically pointed toward an instantaneous mixing over a distance d_{ov} ~ 0.2H_{P} (where H_{P} is the local pressure scale height) with rather large startostar variations. The case of lowmass stars (typically M ≲ 1.5 M_{⊙}) is known to be problematic within this formalism. Indeed, for stars with small convective cores, the overshooting region becomes unrealistically large because H_{P}(r) → ∞ when r goes to zero. This has prompted several authors to consider an overshoot parameter α_{ov} that increases with stellar mass in the approximate mass range 1.1 M_{⊙} ≲ M ≲ 1.5 M_{⊙} (e.g., Pietrinferni et al. 2004; Bressan et al. 2012). In these cases, an ad hoc linear increase of α_{ov} as a function of M was chosen, with some success in reproducing the turnoff of clusters with turnoff masses around 1.3 M_{⊙} (Pietrinferni et al. 2004). The problem remains, however, poorly constrained in this range of mass, and the use of eclipsing binary systems for this purpose is unfortunately of little help (Valle et al. 2016).
Recently, constraints on the extent of the extra mixing beyond convective cores have been obtained from asteroseismology. Sharp variations in the mean molecular weight profile at the boundary of the mixed core create a glitch to which oscillation modes are sensitive, which can be used to measure the extent of the mixed region associated to convective cores. This approach has been successfully applied to solarlike pulsators in the main sequence (Deheuvels et al. 2010b; Goupil et al. 2011; Silva Aguirre et al. 2013; Guenther et al. 2014; Appourchaux et al. 2015), in the subgiant phase (Deheuvels & Michel 2010, 2011), and to several mainsequence B stars (e.g., Degroote et al. 2010; Neiner et al. 2012; Moravveji et al. 2015). All these studies report the need for extended convective cores and confirmed the great potential of asteroseismology to measure this extension. However, we are still lacking consistent seismic studies of larger samples of stars, which are needed to better understand how the overshooting distance varies with stellar parameters.
In this paper, we took advantage of the detection of solarlike oscillations in hundreds of solarlike pulsators with an unprecedented level of precision by the space mission Kepler (Borucki et al. 2010) to consistently measure the extent of the convective core in a larger sample of stars. We have focused on stars whose masses lie around the mass limit for having a convective core (M ≳ 1.1 M_{⊙} at solar metallicity). For these stars, a large part of the core luminosity comes from the burning of ^{3}He outside of equilibrium. Core overshooting can considerably increase the abundance of ^{3}He in the core and therefore also the core luminosity, size, and lifetime (Roxburgh 1985; Deheuvels et al. 2010b). For instance an instantaneous overshooting over a distance of 0.1 H_{P} in a 1.3M_{⊙} star generates an increase of as much as 50% in the convective core radius during the main sequence^{1}. As a consequence, these stars are particularly good tracers of the existence and amount of core overshooting.
It has been shown in previous studies that the small separations built with l = 0 and l = 1 modes are particularly sensitive to the structure of the core (Provost et al. 2005; Deheuvels et al. 2010b; Silva Aguirre et al. 2011) and that their ratios r_{010} to the large separations are nearly insensitive to the socalled nearsurface effects (Roxburgh & Vorontsov 2003). In Sect. 2, we show that this diagnostic can be used to obtain a modeldependent estimate of the extent of the mixed core by building a grid of models with the evolution code Cesam2k (Morel & Lebreton 2008). We then select a subsample of 24 solarlike pulsators among Kepler targets that are the most likely to provide constraints on the amount of core overshooting based on the results of our grid of models, and we extract their mode frequencies from their oscillation spectra in Sect. 3. In Sect. 4, we compare the observed ratios r_{010} to those of two grids of models computed with Cesam2k and mesa (Paxton et al. 2011). We consistently detect convective cores in eight of the selected targets and we obtain measurements of the extent of the mixed core in these stars. In Sect. 5 we show how our results can be used to calibrate the prescription used in the code Cesam2k, and we address the question whether such a calibration can be adapted in mesa.
2. Estimating the core size with seismology
2.1. Asteroseismic diagnostics
A sharp gradient of the mean molecular weight μ builds up at the boundary of the homogeneous convective core, which induces rapid variations in the sound speed profile, and even makes it discontinuous in the case of a growing core without microscopic diffusion. It is well known that such a glitch in c(r) adds an oscillatory modulation to the expression of the mode frequencies as a function of the radial order. The period of this modulation is directly related to the depth of the glitch (Gough 1990). This is not specific to the boundary of convective cores, and such acoustic glitches can also be produced by the base of convective envelopes or the helium ionization regions.
When the period of this oscillation is smaller than the frequency range of the observed frequencies, the acoustic depth of the glitch can be estimated in a modelindependent way. It has been recently shown that the depth of the second helium ionization zone and the base of the convective envelope can be estimated with such a diagnostic (Lebreton & Goupil 2012; Mazumdar et al. 2014). Unfortunately, the glitch caused by convective cores induces a longerperiod oscillation and only a fraction of the period can be observed. This makes it more difficult to obtain modelindependent information about the boundary of convective cores. Cunha & Brandão (2011) and Brandão et al. (2014) showed that the amplitude of the sound speed discontinuity at the core edge may be recovered in some favorable cases. It is not clear whether a modelindependent estimate of the extent of the mixed core can be obtained.
However, it has been shown by several studies that a modeldependent measurement of the core size can be obtained through seismology. Combinations of mode frequencies built with l = 0 and l = 1 modes are well suited for this type of study because they are particularly sensitive to the core structure (Provost et al. 2005; Deheuvels et al. 2010b). Roxburgh & Vorontsov (2003) advised to use the fivepoint separations d_{01} and d_{10} defined as They showed that the ratios between these small separations and the large separations constructed as where Δν_{l}(n) = ν_{l,n}−ν_{l,n−1} are largely insensitive to the structure of the outer layers, which makes them almost immune to the socalled nearsurface effects. These ratios, referred to as r_{010} when combined together, have been used e.g., to estimate the depth of the convective envelope and the second helium ionization zone in the Sun (Roxburgh 2009) or to establish the existence of a convective core in a Kepler target (Silva Aguirre et al. 2013). We note that Cunha & Metcalfe (2007) proposed to use a combination of frequencies using modes of degrees up to 3 (dr_{0213}), which can interestingly be related to the intensity of the sound speed jump at the edge of growing cores. However, l = 3 modes have low amplitudes in stars other than the Sun, and although several detections of such modes have been obtained (e.g., Deheuvels et al. 2010a; Metcalfe et al. 2010), it remains exceptional to reliably estimate their frequencies over several consecutive radial orders. In this study, we have tested and used the diagnostic based on the r_{010} ratios.
Figure 1 shows the behavior of the r_{010} ratios for a model of 1.2 M_{⊙} evolved from the zeroage main sequence (ZAMS) to the beginning of the subgiant phase. The ratios are represented only in the frequency range where modes are expected to be observed, i.e., over about 12 radial orders around the frequency of maximum power of the oscillations ν_{max}. As mentioned above, only a fraction of the period of the oscillation induced by the edge of the core can be observed, and the r_{010} ratios can in fact be well approximated by secondorder polynomials throughout the MS, as can be seen in Fig. 1.
Several studies have shown that the slope and mean value of r_{010} ratios are a good indicator of the size of the mixed core (Popielski & Dziembowski 2005; Deheuvels et al. 2010b; Silva Aguirre et al. 2011). However, these previous studies either focused on a particular star or worked with models that share the same physical properties other than the mixing at the edge of the core. We know that several other parameters, such as the abundance of heavy elements, have a significant impact on the size of the convective core. We here aimed at testing the efficiency of this diagnostic tool.
Fig. 1 Variations in the ratio r_{010} around ν_{max} as a function of frequency for models of 1.2 M_{⊙} from the ZAMS (dark blue) to the beginning of the post main sequence (dark red). The dashed lines correspond to fits of 2nd order polynomials. 
2.2. Testing the diagnostic of r_{010} ratios
2.2.1. Description of the grid
To determine in which circumstances the extent of the core can be estimated with the r_{010} ratios, we computed a grid of models using the stellar evolution code Cesam2k (Morel & Lebreton 2008).
We used the OPAL 2005 equation of state and opacity tables as described in Lebreton et al. (2008). The nuclear reaction rates were computed using the NACRE compilation (Angulo et al. 1999) except for the ^{14}N(p,γ)^{15}O reaction where we adopted the revised LUNA rate (Formicola et al. 2004). The atmosphere was described by Eddington’s gray law. We assumed the classical solar mixture of heavy elements of Asplund et al. (2009) (hereafter AGSS09). Convection was treated using the CanutoGoldmanMazzitelli (CGM) formalism (Canuto et al. 1996). This description involves a free parameter, the mixing length, which is taken as a fraction α_{CGM} of the pressure scale height H_{P}. We here assumed a value of α_{CGM} calibrated on the Sun (α_{⊙} = 0.64, Samadi et al. 2006).
To account for the physical processes that could increase the size of convective cores, we considered an instantaneous mixing beyond convective cores over a distance d_{ov} taken as a fraction α_{ov} of the pressure scale height H_{P}. The free parameter α_{ov} is as often referred to as the overshoot parameter. In order to avoid the overshooting region from unrealistically extending over a distance as large as the core itself, Cesam2k models define the overshooting distance as (5)where r_{s} is the Schwarzschild limit of the core. We note that this is the case during most of the main sequence for stars with masses ≲1.5 M_{⊙}, as shown by Fig. 2. We have imposed the adiabatic temperature gradient in the overshoot region.
Fig. 2 Variations in the pressure scale height H_{P} (dashed line) and the radius of the extended convective core R_{c} (solid blue line) with age for a 1.3 M_{⊙} Cesam2k model with solar metallicity, a solarcalibrated value for the mixing length, Y_{0} = 0.26, and α_{ov} = 0.1. The gray solid line indicates the Schwarzschild limit. 
Microscopic diffusion is known to increase the abundance in heavy elements in the core as the star evolves, and thus to increase the size of convective cores. In this section, microscopic diffusion is not included in the models so that the core extension imposed by the overshoot parameter α_{ov} can be partly attributed to its effects. The contribution from microscopic diffusion is addressed in Sect. 4.
The grid was computed with masses ranging from 0.9 to 1.5 M_{⊙} (step 0.05 M_{⊙}), metallicities from −0.4 to 0.4 dex (step 0.1 dex), and two values of the initial helium abundance (0.26 or 0.30). Models were computed for values of α_{ov} ranging from 0 to 0.3 (step 0.05). For each evolutionary sequence, the mode frequencies were computed with the oscillation code losc (Scuflaire et al. 2008) for about 60 models between the ZAMS and the beginning of the subgiant phase. We stopped the evolution as soon as mixed modes appear around ν_{max}, because these modes cause brutal variations in the r_{010} ratios and prevent them from being directly used as a diagnostic for the core size.
Fig. 3 Left: evolutionary tracks of stellar models of 1.2 M_{⊙} in the (a_{1},a_{0}) plane for different amounts of core overshooting: α_{ov} = 0 (gray), 0.1 (blue), 0.15 (cyan), 0.2 (green), 0.25 (red), and 0.3 (purple). Full (resp. dashed) lines indicate that the model has a convective (resp. radiative) core. Right: variations in the size of the convective core as a function of age for the same models. 
For each of the models along the evolutionary tracks, we fitted 2nd order polynomials of the type (6)to the r_{010} ratios. The parameters β, γ_{1}, and γ_{2} were chosen to ensure that P(ν) is a sum of orthogonal polynomials for each model. The fits were performed in the approximate frequency range where modes are expected to be observed, i.e., about 12 orders around the frequency of maximum power of oscillations ν_{max}. This latter frequency was estimated for stellar models by assuming that it scales as the acoustic cutoff frequency. This assumption, which is the basis of the socalled seismic scaling relations, was observationally verified to work at the level of a few percent at least (Stello et al. 2008; Huber et al. 2011; Silva Aguirre et al. 2012), and is gaining theoretical support (Belkacem et al. 2011). We note that during most of the MS, the r_{010} ratios vary roughly linearly with frequency in the range of observed frequencies, so that the coefficient a_{2} of the fit is negligible.
2.2.2. Evolutionary tracks in the (a_{1},a_{0}) plane
Before commenting on the results of the grid, we show as an example the evolutionary tracks in the (a_{1},a_{0}) plane (slope versus mean value) of 1.2M_{⊙} models for different amounts of core overshooting (Fig. 3a). For comparison, the variations in the size of the convective core for the same models are shown as a function of age in Fig. 3b. As mentioned by Silva Aguirre et al. (2011), the trajectory of models in the (a_{1},a_{0}) plane depends in a complex way on the evolutionary stage, the size of the convective core, and the amplitude of the glitch in the sound speed. However, we can still broadly understand it. At the beginning of the MS, the stars with different α_{ov} start roughly at the same point in the (a_{1},a_{0}) plane (bottom right corner in Fig. 3a). Indeed, the μgradient at the edge of the core has not had time to build up yet, so the r_{010} ratios are still nearly independent from the size of the convective core. As the star evolves, the glitch in the sound speed profile builds up, which causes the amplitude of the oscillations of the r_{010} ratios to increase. Therefore both the mean value a_{0} and the absolute value of the slope  a_{1}  of the ratios increase. But also, as the star evolves, its ν_{max} frequency decreases. As a result, the range of observable frequencies shifts to a different part of the oscillation produced by the glitch. As can be seen in Fig. 1, when stars reach the end of the MS, the r_{010} ratios lie around a maximum of this oscillation, which results in a decrease of the absolute value of the slope  a_{1} . For postmain sequence stars, the mean slope a_{1} even becomes positive. This explains why the evolutionary tracks of models in the (a_{1},a_{0}) plane are vaguely circular, as can be seen in Fig. 3.
For models with larger amounts of overshooting, the convective core is larger. Therefore, the period of the oscillation caused by the glitch is shorter and the absolute mean slope  a_{1}  of the r_{010} ratios is larger. As a result, stars with larger α_{ov} are shifted to the left in the (a_{1},a_{0}) plane, and they draw larger circles. This confirms previous statements that the position in the (a_{1},a_{0}) plane is discriminant for the size of the core if all parameters other than α_{ov} are fixed.
2.2.3. Results of the grid
Fig. 4 Location of models in the (a_{1},a_{0}) plane at fixed Δν. Colors indicate the amount of core overshooting: α_{ov} = 0 (gray), 0.1 (blue), 0.15 (cyan), 0.2 (green), 0.25 (red), 0.3 (magenta). Open squares indicate models with a convective core, and crosses, models with radiative cores. The black open circles indicate models that are in the postmainsequence (X_{c}< 10^{2}). 
Global parameters of the selected targets.
When solarlike oscillations are detected in a star, it is usually straightforward to estimate the mean large separation of its acoustic modes Δν. We thus chose to show the results of the grid at fixed values of Δν. This time, each evolutionary sequence of our grid is represented as a dot in the (a_{1},a_{0}) plane, provided its large separation matches the chosen value of Δν at some point along the evolution.
Figure 4 shows the location of the models in the (a_{1},a_{0}) plane for four values of Δν: 110, 95, 70, and 65 μHz. For Δν = 110 μHz (top left plot), there is a relative degeneracy of the models in the (a_{1},a_{0}) plane. This can be understood because only lowmass unevolved stars reach such a high value of Δν. Highermass stars begin the MS with a lower Δν, and this quantity further decreases as the star evolves^{2}. For this reason, few of the stars with Δν = 110 μHz have a convective core. And those that have one are still close to the ZAMS, so the μgradient has not had time to build up yet and the r_{010} ratio still does not feel it. The diagnostic is thus less efficient for Δν ≳ 110 μHz.
For lower values of Δν, different populations are represented: (1) evolved lowmass stars (in the PoMS for the lowest masses) and (2) MS highermass stars. In these cases, Fig. 4 clearly shows that the location of a model in the (a_{1},a_{0}) plane can be used to estimate:

the evolutionary state: as mentioned before, when stars leave theMS, the mean slope a_{1} of the ratios increases and becomes positive.As a result, PoMS models occupy a place in the (a_{1},a_{0}) plane that isincreasingly distinct from that of MS models, as the largeseparation decreases. This opens the possibility to determine theevolutionary status of a star from its location in the (a_{1},a_{0}) plane;

the existence and the size of the convective core: for stars with large separations below ~95 μHz, models with α_{ov} = 0, 0.1, 0.2, and 0.3 occupy distinct regions in the (a_{1},a_{0}) plane, which suggests that it should be possible to measure the size of the mixed core by using the location of the star in this plane.
We stress that the effects of metallicity on the size of the core are here taken into account in a very conservative way, since the models of Fig. 4 include a wide range of metallicities (−0.4 to 0.4 dex). In practice, the metallicity of an observed star is usually known with a much better accuracy if spectroscopic measurements are available. We thus conclude that the r_{010} ratios are in principle an efficient tool to measure the size of convective cores, provided the observed star is evolved enough to have developed a glitch in the sound speed at the edge of the core.
3. Extracting the r_{010} ratios from Kepler targets
3.1. Selection of targets
Based on the tests performed on stellar models in Sect. 2.2, we established a set of criteria to select Kepler targets for which the r_{010} ratios should provide a good diagnostic for the core structure. We selected stars for which

the mean large separation is below 110 μHz, so that the diagnostic tool is efficient;

no mixed modes are contaminating the r_{010} ratios;

a long enough data set is available, so that a good precision can be attained in the estimates of the parameters a_{i}. Even with 9 months of Kepler data, the r_{010} ratios of a target studied by Silva Aguirre et al. (2013) were contaminated by a spurious increase in the lowsignaltonoise part of the spectrum. To avoid these features that might bias our estimates of the a_{i} parameters, we selected only stars that were observed for at least 9 months;

the observed modes are narrow enough: we excluded F stars, whose modes are too wide to unambiguously distinguish the l = 1 ridge from the l = 0 and l = 2 ridges in an échelle diagram. We note that Bayesian methods have been proposed to identify the degree of the modes and extract the mode frequencies even in these cases (e.g., Benomar et al. 2009). However, this type of analysis requires dedicated works, which can be undertaken as an interesting followup of this work to explore the sizes of convective cores in highermass stars.
We applied these criteria to the solarlike pulsators whose global parameters were determined by Chaplin et al. (2014) and obtained a list of 24 targets, which are given in Table 1. Most of these stars were also observed spectroscopically from the ground, which yielded estimates of the effective temperature and of the surface metallicity. When available, these measurements are specified in Table 1.
Fig. 5 Ratios r_{010} computed for KIC 6106415 (left) and KIC 12258514 (right) using the mode frequencies extracted from the Kepler oscillation spectra (see text). The colored dashed lines correspond to 2ndorder polynomial fits to the observed ratios using either the raw covariance matrix (gray lines) or the covariance matrix modified through truncated SVD (see Sect. 3.3). 
3.2. Extraction of the mode frequencies
The mode frequencies of 13 out of the 24 selected targets were already extracted from Kepler observations by Appourchaux et al. (2012). However, this study was performed with nine months of Kepler data, whereas at current time almost three years of data are available in the most favorable cases. We thus decided to reanalyze all the targets of the selected sample using the full Kepler data sets available (until Q16) to date. For this purpose, we used a maximum likelihood estimation (MLE) method in the same way as previously applied to CoRoT and Kepler targets (e.g., Appourchaux et al. 2008; Deheuvels et al. 2010b). For each star, we adjusted Lorentzian profiles to all the modes simultaneously (global fits). We here neglected the rotational splitting of the modes and fitted only one component for each multiplet of degree l and radial order n. Since the stars of the sample are expected to be slow rotators, the rotational multiplets should be approximately symmetrical with respect to their m = 0 component. As a result, we expect negligible bias due to rotation in our estimates of the mode frequencies. We stress that in this work, we were only interested in estimating the a_{i} parameters of a polynomial fit to the r_{010} ratios of the observed stars. As a result, we did not seek to estimate the frequencies of lower signaltonoise modes around the edges of the frequency range of observed modes. We obtained estimates of the mode parameters over 9 to 15 overtones for the 24 targets. The results are given in Tables A.1 to A.6 in Appendix A. Our results are in good agreement with those obtained by Appourchaux et al. (2012) for the targets that are among our sample. We indeed found that 31% (resp. 8%, 3%) of the fitted mode frequencies agree within 1 (resp. 2, 3) σ with the results of Appourchaux et al. (2012), which is close to what is statistically expected.
We used the estimated mode frequencies to evaluate the global seismic parameters of the selected targets. A linear regression of the frequencies of l = 0 modes as a function of the radial order n provided an estimate of the mean large separation Δν. The obtained values are given in Table 1. We then performed a gaussian fit to the mode amplitudes as a function of frequency. The central frequency of the fitted Gaussian provides an estimate of the frequency of maximum power of the oscillations ν_{max} (see Table 1). Seismic scaling relations were then used to relate the global seismic parameters Δν and ν_{max}, and the effective temperature T_{eff} to the stellar mass and radius. The underlying assumption behind seismic scaling relations was already mentioned in Sect. 2.2.1. Whenever it was available, we used the spectroscopic T_{eff} obtained by Bruntt et al. (2012). For the three stars of the sample that were not observed by Bruntt et al. (2012), we used a photometric estimate obtained from the recipe proposed by Pinsonneault et al. (2012), which was applied to the griz photometry available from the Kepler input catalogue (KIC). We thus obtained stellar masses ranging from 0.94 to 1.39 M_{⊙} (see Table 1). We note that for all the stars for which both spectroscopic and photometric estimates of T_{eff} were available, the agreement on the stellar masses obtained with both sets of T_{eff} is excellent (below 1σ for all stars except one at 1.7σ).
3.3. Polynomial fit to r_{010} ratios
We used the fitted mode frequencies listed in Tables A.1 to A.6 in Appendix A to compute the r_{010} ratios of all the stars of the sample. Two representative examples are shown in Fig. 5. KIC6106415 (left plot) is still in a phase where the r_{010} ratios are roughly linear in the range of observed frequencies, while the ratios of KIC12258514 have a more parabolic shape. As predicted by stellar models, we found that the observed ratios are well reproduced by 2nd degree polynomials. For several targets, the r_{010} ratios deviate from a mere parabola because of a shortperiod oscillation around the parabolic general trend. This is expected and corresponds to the signature of the base of the convective envelope. In this work, the polynomial fit that we applied to the r_{010} ratios filters out this contribution. This is to our advantage here since we are merely interested in probing the core properties in this study. However, we stress that these signatures of the bottom of the convective envelope can potentially yield precious modelindependent constraints on the stellar structure (Mazumdar et al. 2014) and deserve further investigation. The dip in the profile of the adiabatic index Γ_{1} corresponding to the region of second ionization of helium can also create a shortperiod oscillation in seismic indexes, however the r_{010} ratios are almost insensitive to these shallow regions and the amplitude of the corresponding oscillation is expected to be negligible.
To fit polynomials to the observed r_{010} ratios, one needs to take into account the high level of correlation between the data points. Indeed, each mode frequency is used by several data points. The covariance matrix between linear combinations of the mode frequencies (e.g., between the d_{01} and d_{10} separations as defined by Eqs. (1) and (2)) can easily be computed analytically, but it is much harder for the r_{010} ratios because of the division by the large separations. We therefore resorted to Monte Carlo simulations using the observed mode frequencies and their associated error bars to estimate the covariance matrix C for each star. This approach supposes that the errors in the mode frequency estimates are normally distributed, which has been shown to be a valid approximation (Benomar et al. 2009), except for low signaltonoiseratio modes, which we have excluded here^{3}. The optimal parameters a_{0}, a_{1}, and a_{2} of the polynomial described in Eq. (6) were then obtained by a leastsquare minimization of the residuals weighted by the coefficients of the inverse W of the covariance matrix, as described in Appendix B. This type of fitting is now applied routinely to fit stellar models constrained by combinations of mode frequencies (e.g., Silva Aguirre et al. 2013; Lebreton & Goupil 2014). However, when applied directly to our simple case of a polynomial fit of the r_{010}, we obtained poor fits to the observed ratios (see gray dashed lines in Fig. 5).
After careful inspection of the results, we found that the covariance matrix C is in fact illconditioned, with a conditioning of the order of 10^{5} or 10^{6}. As a result, the covariance matrixes are nearly noninvertible, which explains the poor agreement obtained by direct fitting. This property is not specific to our particular case, and we expect any covariance matrix built with combinations of frequencies to show similar behavior as the number of points increases. The conditioning of matrix C increases as the number of modes involved in the combinations of frequencies increases, which explains why the problem is so obvious for the r_{010} ratios, but with a large enough number of points, it also arises for threepoint separations. To remedy this problem, we applied truncated SVD to the covariance matrix as explained in Appendix B. We found that suppressing the 5 smallest eigenvalues of matrix C is generally enough to obtain satisfactory fits to the observed r_{010} ratios (red dashed lines in Fig. 5).
4. Measuring the size of mixed cores in Kepler targets
Since the r_{010} ratios have been shown to efficiently cancel out the contribution from the outer layers (Roxburgh & Vorontsov 2003), the observed ratios could be directly compared to those of models. We thus compared the observed r_{010} ratios to those of two grids of models: the one computed with Cesam2k, which was described in Sect. 2.2, and a second equivalent grid that was built with the evolutionary code mesa (Paxton et al. 2011, 2013), which is described below. Obviously, these grids are too coarse to provide in themselves statistically reliable estimates of the stellar parameters, and in particular of the amount of core overshooting. However, based on the tests performed in Sect. 2.2, these grids can be used to identify stars with a convective core and obtain a rough estimate of the extension of the mixed core in these stars. As a second step presented in Sect. 5, these estimates were refined using a more sophisticated optimization procedure.
Fig. 6 Location in the (a_{1},a_{0}) plane (star symbols and black error bars) of the stars of the sample that were found to be on the MS with a convective core in this study. Models that reproduce the observed large separation, the spectroscopic estimate of metallicity, and the stellar mass derived from scaling laws within 3σ errors are overplotted. The symbols have the same meaning as in Fig. 4. 
Fig. 7 Location in the (a_{1},a_{0}) plane of the MS stars for which the presence of a convective core is uncertain (note that the evolutionary status of KIC9410862 is in fact uncertain). The symbols have the same meaning as in Fig. 6. 
Fig. 8 Location in the (a_{1},a_{0}) plane of the first six PoMS stars of the sample. The symbols have the same meaning as in Fig. 6. 
Fig. 9 Location in the (a_{1},a_{0}) plane of the last four PoMS stars of the sample. The symbols have the same meaning as in Fig. 6. 
4.1. CESAM2K models
For each star of the sample, we selected among the grid described in Sect. 2.2.1 the models that have a surface metallicity within 3σ of the spectroscopic [Fe/H] (all metallicities were included in the cases where no spectroscopic measurement was available), and a stellar mass within 3σ of the estimate obtained from scaling laws (see Table 1). Among the selected evolutionary sequences, we retained only the models whose mean large separations bracket the observed Δν. We note that for both models and observations, the mean value of Δν was estimated using only the modes below ν_{max} so that the corresponding large separations are only slightly affected by nearsurface effects.
For the selected models, we fitted polynomials to the r_{010} ratio as defined by Eq. (6). For this purpose, we used the same modes and the same values of β, γ_{1}, and γ_{2} (see Eq. (6)) as those found from the observations, so that the parameters a_{i} of the models can be directly compared to the observed ones. Since the models that we retained do not exactly match the observed large separation, we performed an interpolation to obtain the parameters a_{i} that correspond exactly to the observed Δν. This process was repeated for all the stars of the sample. Figure 6 through 9 show the location of the selected models and the observations in the (a_{1},a_{0}) plane.
The first comforting observation is that all the observed stars occupy a place in the (a_{1},a_{0}) plane that is populated by models. This shows that in all cases, there exist models that simultaneously reproduce the observed trend of the r_{010} ratio and the other global observational constraints.
Secondly, as anticipated in the previous section, the evolutionary status of the observed stars can be unambiguously established in most cases using the diagnostic from the r_{010} ratios. For 13 stars of the sample, the profile of the r_{010} ratio is only compatible with MS models, the PoMS models lying at least several σ away in the (a_{1},a_{0}) plane (see Figs. 6 and 7). Conversely, 10 stars are clearly in the PoMS phase judging by their location in the (a_{1},a_{0}) plane (see Figs. 8 and 9). We stress that it was not obvious at first sight that these 10 stars are in the subgiant phase. Indeed, the PoMS status of solarlike pulsators is generally established by the presence of mixed modes in their oscillation spectrum. However, at the beginning of the subgiant phase, the lowest order g modes have not yet reached the frequency range of observed modes and such a diagnostic cannot be applied. It is the case for these 10 stars of the sample, and we here showed that the general trend of the r_{010} ratios is a powerful diagnostic for the evolutionary status in this case. The evolutionary status remains ambiguous only for one star of the sample, KIC9410862, which is either at the end of the main sequence or at the beginning of the subgiant phase (Fig. 7).
Among the 13 MS targets, eight have values of the parameters a_{0} and a_{1} that can be reproduced only by models that have a convective core. The locations of these stars in the (a_{1},a_{0}) plane are shown in Fig. 6. As predicted in Sect. 2.2, we were able to use the position in the (a_{1},a_{0}) plane of the stars that have a convective core to obtain an estimate of the amount of core overshooting. Interestingly, the eight stars draw a quite consistent picture of the extension of convective cores in lowmass stars.

All the targets require an extended core compared to the classical Schwarzschild criterion. Indeed, all the stars that have a convective core lie several σ away from models computed without overshooting.

None of the targets were found to be consistent with a core overshooting above α_{ov} = 0.2.

The only target which is consistent with a core overshooting around α_{ov} = 0.2 (KIC7206837) corresponds to the highestmass star of the sample (1.54 ± 0.09 M_{⊙} according to seismic scaling relations). This raises the question of a potential mass dependence of the amount of core overshooting as implemented in the evolution code Cesam2k, which is addressed in more details in Sect. 5.
We stress that seismology provides information about the size of the mixed core at the current age of the star. The amounts of overshooting that are quoted above are those required so that the evolution code Cesam2k produces cores with an appropriate size. One should be careful that the values that were obtained for α_{ov} hold only for the prescription of core overshooting that is implemented in Cesam2k and they should not be directly applied to other codes. We discuss this point in details in Sect. 5.
A more relevant result to quote is the extent of the mixed core obtained from seismic constraints. To determine this for each of the stars for which a convective core was detected, we selected a subset of five optimal models from the grid of models, defined as those that minimize the quantity (7)where the correspond to the N observables used to constrain the models, namely the effective temperature T_{eff}, the surface metallicity (Z/X) (if available), the asteroseismic log g, and the parameters a_{0} and a_{1} of the 2nd order polynomial fit of the observed r_{010} ratio. The σ_{i} are the measurement errors, and the are the values corresponding to the observables computed from the models. We note that the observables can be regarded as independent (since we fitted a sum of orthogonal polynomials to the observed ratios) so that Eq. (7) holds. For each star, the fractional mass of the convective core M_{c}/M_{⋆} for the five best models is shown in Fig. 10 (blue squares for Cesam2k models). We note that the spreads in M_{c}/M_{⋆} observed in Fig. 10 cannot be interpreted as uncertainties on this quantity. Indeed, to estimate proper uncertainties one should have chosen the set of optimal models based on the variations of the χ^{2} function compared to the lowest value of χ^{2} in the grid (Δχ^{2} = 1, 4, and 9 provide 1, 2, and 3σ errors, respectively) but the grid computed here is too coarse to make such an approach possible^{4}.
4.2. MESA models
As mentioned above, we have also computed a second grid of models with the evolution code mesa (Paxton et al. 2011, 2013).
The mesa models were computed using the OPAL 2005 equation of state from the tables of Rogers & Nayfonov (2002), which are completed at lower temperature by the tables of Saumon et al. (1995). mesa opacity tables are constructed by combining radiative opacities with the electron conduction opacities from Cassisi et al. (2007). Radiative opacities are taken from Ferguson et al. (2005) for 2.7 < log T< 4.5 and OPAL opacities Iglesias & Rogers (1993, 1996) for 3.75 < log T< 8.7. The low temperature opacities of Ferguson et al. (2005) include the effects of molecules and grains on the radiative opacity. The nuclear reaction rates module from mesa contains the rates computed by Caughlan & Fowler (1988) and Angulo et al. (1999) (NACRE), with preference given to the NACRE rates when available. The atmosphere was described as Hopf’s gray law. We used the solar mixture from Grevesse & Noels (1993). Convection was treated using the classical mixinglength theory (MLT, BöhmVitense 1958) with a fixed mixing length parameter α_{MLT} = 1.9, which corresponds to a solar calibration (Paxton et al. 2011).
Core overshooting is included and described as a diffusive process, following Herwig (2000). For this purpose, an extra diffusion is added at the edge of the core, with a coefficient (8)where D_{0} is the MLTderived diffusion coefficient near the Schwarzschild boundary, H_{P} is the pressure scale height at this location, and f is the adjustable overshooting parameter. To avoid unrealistically large extensions of convective cores, the current version of mesa uses a modified value for the pressure scale height, defined as (9)in the case where the mixing length ℓ_{MLT} = α_{MLT}H_{P} becomes larger than the Schwarzschild limit r_{s} of the core. This prescription is different from the one adopted in the Cesam2k code. When using the same prescription for core overshooting (instantaneous or diffusive) and the same overshooting parameter at the boundary of small convective cores, the approach followed by mesa is expected to yield core extensions that are smaller by a factor α_{MLT} compared to the extensions produced with the Cesam2k approach in the saturated regime.
Gravitational settling and chemical diffusion are taken into account by solving the equations of Burgers (1969) using the method and diffusion coefficients of Thoul et al. (1994).
Fig. 10 Fractional mass of the convective core for the eight stars that were found to have a convective core in this study. For each star, the open symbols correspond to the core size of the five models of the two grids (blue squares for Cesam2k models, red circles for mesa models) that yield the lowest values of χ^{2} as defined by Eq. (7). The filled squares give the core sizes obtained from a LevenbergMarquardt optimization and the evolution code Cesam2k (see Sect. 5.1). 
Fig. 11 Location in the (a_{1},a_{0}) plane of four stars of the sample compared to the location of models computed with the evolution code mesa. Symbols are the same as in Fig. 6, except for the colors, which indicate diffusive overshooting parameters of: f = 0.004 (gray), 0.010 (blue), 0.016 (cyan), 0.022 (green), 0.028 (red), 0.035 (magenta). 
For each star of the sample, we performed the same model selection as was done with Cesam2k models, and for each selected model we fitted 2ndorder polynomials to the r_{010} ratios in the same way as described in Sect. 4.1. This allowed us to compare the location of the observed stars in the (a_{1},a_{0}) plane to that of mesa models. Figure 11 shows the results obtained for four stars of the sample, which are representative of the different cases identified in Sect. 4.1: KIC8228742 and KIC7206837 are in the MS and have a convective core, KIC6106415 is in the MS but has no convective core, and KIC6933899 is in the PoMS.
The mesa grid agrees with the Cesam2k grid on the evolutionary status of all the stars of the sample. The star KIC94110862, whose evolutionary status was uncertain based on Cesam2k models, was found to be more consistent with models shortly after the end of the MS using the mesa grid. Additionally, the eight stars identified as having a convective core with the Cesam2k grid were also found to have one with the mesa grid. The locations of two of these stars in the (a_{1},a_{0}) plane are shown in the upper panels of Fig. 11. It is clear that the extension can be estimated from the a_{0} and a_{1} parameters, as was claimed in Sect. 4.1. Interestingly, all the conclusions reached with Cesam2k models about the amount of overshooting that is required are confirmed. The eight stars with convective cores all require an extended core with overshooting parameters ranging from 0.010 to 0.035, and the star that requires the largest amount of overshooting corresponds to the highestmass stars of the sample (KIC7206837) as was found in Sect. 4.1.
Obviously, the overshooting parameters obtained from the mesa models are not directly comparable to those found from the Cesam2k grid because a diffusive overshooting was chosen in mesa models. A more detailed comparison is provided in Sect. 5.2, but we can already compare directly the absolute sizes of the extended cores found with both evolution codes. For all the stars that have a convective core, we selected the five models of the mesa grid that minimize the χ^{2} function as defined by Eq. (7). The fractional mass of the mixed core M_{c}/M_{⋆} in these models is shown in Fig. 10. Interestingly, there is a quite good agreement on the size of the extended cores obtained with both evolution codes, in spite of the different prescriptions for core overshooting. This is further indication that the seismic diagnostic based on r_{010} ratios can provide a measurement of the size of the mixed core mostly independently of the input physics, as was already suggested by Silva Aguirre et al. (2011).
4.3. Instantaneous vs. diffusive mixing beyond convective cores
In this study, we have chosen to adopt two different prescriptions for core overshooting, an instantaneous overshooting (Cesam2k models) and a diffusive overshooting (mesa models), with the aim of comparing the two most frequently used prescriptions to Kepler data. It is interesting to address the question whether we can distinguish between these two types of mixing beyond convective cores using r_{010} ratios. The mesa code offers the possibility to test this since both treatments have been implemented. We computed a 1.3M_{⊙}mesa model including diffusive overshooting with a parameter f = 0.020, which we evolved until X_{c} has dropped to 0.2 (chosen arbitrarily). We also computed a mesa model including a step overshooting with α_{ov} = 0.22 and a slightly higher mass (1.31 M_{⊙}) evolved until it has the same large separation as the diffusiveovershooting model. We found that both models are undistinguishable from an observational point of view (within typical observational errors), and they also share a very similar behavior of the r_{010} ratios, as is shown in Fig. 12. This shows that the seismic diagnostic based on the r_{010} ratios is unfortunately not capable of distinguishing between the two scenarios regarding the nature of the extra mixing beyond the core.
Fig. 12 Profile of r_{010} for two mesa models: one with a diffusive overshooting (f = 0.02) and a 1.3M_{⊙} mass (blue squares) and the other with a step overshooting (α_{ov} = 0.22) and a 1.31M_{⊙} mass (red circles). Both models have the same mean large separation. 
5. Toward a calibration of core overshooting for lowmass stars
In Sect. 4, we were able to measure the sizes of mixed cores in eight lowmass stars using seismology. The question is then how these results can be used to estimate the efficiency of the extra mixing beyond convective cores. Answering this question is not straightforward. One could consider simply comparing the convective core masses obtained in Sect. 4 to the convective core masses that would be obtained with identical stellar parameters but no mixing beyond the core. This is, however, inapplicable in practice because increasing the size of the convective core at the beginning of the main sequence has large subsequent effects on its composition and evolution. For stars in the mass range that we considered here, the main effect is that the abundance of ^{3}He in the core increases, which increases its luminosity, and thus also its size because the Schwarzschild radius increases. For instance, extending the convective core of a 1.3M_{⊙} star over 10% of the Schwarzschild radius in fact results in an increase of the core radius of as much as 50% during the main sequence. Another consequence is that the lifetime of small convective cores can be dramatically extended (see Roxburgh 1985; Deheuvels et al. 2010b). For instance, a stellar model of KIC62245718 computed with the same stellar parameters as those found in Sect. 4 but without including any extra mixing beyond the core has lost its convective core at current age.
It therefore seems that the problem of the efficiency of convective core extensions cannot be studied independently from the evolution of the star, even though asteroseismology only tells us about the size of the mixed core at current age. We thus chose to estimate the efficiency of the extra mixing beyond the core by adjusting the overshooting parameter (α_{ov} for instantaneous mixing or f for diffusive mixing) considered constant throughout the evolution, so that stellar models have the right convective core size at current age. As mentioned in the introduction, so far we have had to model convective core extensions using such simplistic parametric models because we lack observational constraints that would justify using more complex models. Our aim in this section is to search for correlations between the efficiency of overshooting and properties of stellar interiors, which might eventually give us better insight on the physical processes that are responsible for core extensions, and lead us to prefer more realistic modelings of this phenomenon. On the shorter term, this type of study can enable us to propose a calibration of the overshooting parameter, which can later be used in 1D stellar models.
5.1. Calibration of core overshooting in CESAM2K
To calibrate core overshooting in Cesam2k models, we needed to obtain more quantitative estimates of the amounts of core overshooting that are required for the stars of the sample.
5.1.1. Stars with a convective core
We performed optimizations for the eight stars that were found to have a convective core in Sect. 4. For this purpose, we used the LevenbergMarquardt algorithm, which is an appealing alternative to gridsearch minimization when the number of free parameters is large. This algorithm combines the low sensitivity to initial guesses of the gradient search method and the rapidity of convergence of the NewtonRaphson method. Its use has first been suggested for the purpose of stellar modeling by Miglio & Montalbán (2005). The main drawback of such an optimization technique is the risk to converge toward a secondary minimum of the cost function if the initial guesses are to far from the optimum set of parameters. In our particular case, this risk is minimized since we used the best models of the grid computed in Sect. 4.1 as initial guesses.
To find optimal models, we minimized the quantity χ^{2} as defined in Eq. (7). We used the same observables as those listed in Sect. 4.2, to which we added the frequency of the lowestorder observed radial mode. This observable is preferred to the observed mean large separation because of its lower dependence on the structure of the outer layers. We note that the a_{2} parameter of the 2nd order polynomial fit of the observed r_{010} ratio was here included as a constraint. This parameter becomes constraining for evolved stars, for which the observed r_{010} ratios depart from a simple linear relation (see Fig. 1). To reproduce these observables five parameters were left free: the stellar mass, age, initial helium abundance Y_{i}, initial metallicity (Z/X)_{i}, and the parameter of core overshooting α_{ov}. We imposed a lower limit of 0.24 for Y_{i} in order to exclude models with initial helium abundances significantly below the standard big bang nucleosynthesis (SBBN) values of Y_{0} = 0.248 ± 0.007 (Steigman 2010). To limit the number of free parameters, we kept the mixing length fixed to α_{CGM} = 0.64, which was obtained from a solar calibration. As a consequence, the fit that we performed has two degrees of freedom and a reduced value was thus obtained by dividing the regular χ^{2} by two. For each star, two types of optimizations were performed, one where the effects of microscopic diffusion are neglected, and another that includes these effects following the formalism of Burgers (1969). This procedure enabled us to test the influence of microscopic diffusion on the amount of core overshooting that is required. As mentioned above, diffusion increases the abundance of heavy elements in the core and thus the opacity, which results in an increase in the size of the convective core. We therefore expected to require less core overshooting when microscopic diffusion is included. Since Cesam2k does not include the computation of radiative accelerations of chemical elements, their effect was neglected in this study. Radiative levitation acts against gravitational settling in the interior of stars with masses above about 1.2 M_{⊙}, so our models including microscopic diffusion likely overestimate the sinking of heavy elements in this mass range. We thus expect our models computed with microscopic diffusion and stellar masses above 1.2 M_{⊙} to provide us with an upper limit to the effects of diffusion, in particular on the sizes of convective cores.
Fitted parameters and characteristics of the bestfit models obtained for the stars of the sample that have a convective core.
Fig. 13 Amount of core overshooting found for the stars of the sample that have a convective core as a function of the fitted stellar mass (left) and as a function of the fitted initial metallicity (right). Blue squares indicate models computed without microscopic diffusion and gray circles, models where microscopic diffusion is included following Burgers (1969). The vertical arrows indicate upper limits of α_{ov} (see Sect. 5.1.2). 
The parameters of the bestfit models are given in Table 2. The quoted error bars were obtained as the diagonal coefficients of the inverse of the Hessian matrix. The results confirm that the amount of core overshooting can be well constrained by using the parameters a_{i}. The values obtained for α_{ov} range from 0.07 to 0.18 in the case without diffusion, which is in good agreement with the results of the grids of models (Sect. 4.1). As foreseen, the models that include microscopic diffusion require lower amounts of core overshooting to reproduce the seismic observations, with values ranging from 0.05 to 0.15. However, our results show that the effects of diffusion cannot account in themselves for the entire extension of convective cores since core overshooting was required for all eight stars of the sample. We note that for several stars of the sample, the fitted value of the initial helium abundance Y_{i} coincides with the lower limit of 0.24 that we have imposed to avoid subSBBN helium abundances. Similar results have been found in several studies where seismic modelings were performed (e.g., Metcalfe et al. 2014; Silva Aguirre et al. 2015). This is potentially the consequence of the wellknown correlation between stellar mass and helium abundance (Lebreton & Goupil 2012). For these stars, we have performed additional fits imposing a higher limit to Y_{i} (0.26) and found results that agree within 1σ errors with the values quoted in Table 2 (in particular, we found very little difference in the sizes of convective cores, which is our main interest here). The optimizations also provided estimates of the stellar mass, which are given in Table 2. The agreement with estimates from scaling laws is quite good (below 1.3σ for all the stars). We note that KIC12009504 was already modeled by Silva Aguirre et al. (2013) who already found that this stars possesses a convective core that extends beyond the Schwarzschild limit. Our results for this star are in good agreement with those of Silva Aguirre et al. (2013).
The values of for some of our fits are significantly larger than 1, which in principle indicates either disagreements between models and observations, or underestimated error bars for the observables. Table 2 gives the level of agreement with observations for each fitted parameter normalized by observational 1σ errors. It shows that a very good level of agreement is reached for the a_{0} and a_{1} parameters, as was expected based on the results of Sect. 4. On the contrary, disagreements above the 3σ level arise for the a_{2} parameter. This occurs mainly for stars where the r_{010} ratios vary nearly linearly with frequency, so that the a_{2} coefficient is small. In this case, the observational estimate of a_{2} can be altered by the shortperiod oscillation that arises because of the glitch at the base of the convection zone (see Sect. 3.3). Disagreements above the 2σ level also arise for the effective temperature and the surface metallicity. We note that we have used the error bars of Bruntt et al. (2012) for these quantities, which have been deemed somewhat underestimated in previous studies (Silva Aguirre et al. 2013). This might at least partly explain this disagreement. Also, we note that the agreement with the observed surface metallicities improves when including microscopic diffusion in models.
Using our optimizations, we could also obtain estimates of the total size of the mixed core in the eight stars. Since the size of the core is not a fitted parameter, the optimization algorithm does not directly provide error bars on these obtained values. However, they can be deduced from the relation (10)where the b_{j} terms correspond to the P free parameters and the derivatives can be evaluated with the models used to compute the Hessian matrix. The fractional masses of the convective cores for the eight stars are plotted along with their error bars in Fig. 10.
The refined estimates of the amount of core overshooting in the eight stars that have a convective core enabled us to test correlations between the overshooting parameter α_{ov} and other stellar parameters. Figure 13a shows the obtained values of α_{ov} as a function of the stellar mass. We observe that there seems to be a tendency of core overshooting to increase with stellar mass in this mass range. This tendency is less clear for the models where microscopic diffusion was included (gray circles), but we still found in this case that the three less massive stars of the sample require less core overshooting that the five more massive ones. Clearly more data points are required to be more conclusive, but if such a tendency is confirmed, then an empirical law could be derived and implemented in the Cesam2k code in order to better model the extent of mixed cores for stars in this mass range. We also note that we have found no apparent dependency of the amount of core overshooting required with stellar metallicity (see Fig. 13b).
5.1.2. Stars without a convective core
Information about core overshooting can also be drawn from stars that have no convective core but lie just below the mass limit for having one. Indeed, above a certain amount of core overshooting, the models all develop a convective core and the profile of the r_{010} becomes at odds with the observations. So these targets can be used to obtain an upper limit to the amount of core overshooting. For these targets, we performed optimizations using the LevenbergMarquardt algorithm as before, except that we fixed the parameter of core overshooting to predefined values ranging from 0 to 0.3. The result of this procedure is shown as an example for the case of KIC 10516096. For α_{ov} = 0, the fit converges toward a PoMS model with a mass of 1.12 M_{⊙}, an age of about 6.4 Gyr, a metallicity of (Z/X) = 0.0229 and no convective core. For , the fits converge toward roughly the same model. The only difference between the bestfit models is that the initial convective core survives longer for higher values of α_{ov} (about 1 Gyr for α_{ov} = 0.15 compared to 30 Myr for α_{ov} = 0). However, even with α_{ov} = 0.15 the convective core vanishes long before the end of the MS and its effect on the core structure has been washed out by the age of 6.4 Gyr. On the contrary, for α_{ov} = 0.2 this model keeps a convective core until the end of the MS. As a result, the duration of the MS is extended and by the time the model reaches the observed large separation, it is still in the MS with a convective core and the r_{010} ratio of this model is in poor agreement with the observations. Therefore, to decrease the χ^{2}, the fit converges toward a model with higher metallicity (Z/X = 0.0281) for which the convective core vanishes before the end of the MS. However, this latter model is in less good agreement with the observations as can be seen in Fig. 14. We thus obtained an upper limit of the overshooting parameter of about 0.19 for this star (value of α_{ov} above which the obtained χ^{2} is larger than min(χ^{2}) + 9). Similar results were found for one other PoMS star (KIC 6933899) and three MS stars (KIC 6106415, KIC 6116048, and KIC 8394589). For all these stars, the agreement deteriorates for an upper limit 0.16 <α_{lim}< 0.20. These constraints were added as vertical arrows in Fig. 13. Unfortunately, they are too loose to confirm the tendency of α_{ov} to increase with mass that was found in Sect. 5.1.1.
Fig. 14 Value of the χ^{2} of the bestfit model as a function of the (fixed) amount of core overshooting for the star KIC 10516096. 
Characteristics of the bestfit models obtained for KIC5184732 when modifying the chosen input physics.
5.1.3. Sensitivity to input physics
We here briefly address the question of the sensitivity of our results to some of the choices of the model input physics. We focused on one star (KIC5184732) chosen arbitrarily among the stars which was found to have a convective core and was modeled in Sect. 5.1.1. We performed additional optimizations of this target modifying each time one assumption on the model input physics. We note that the influence of microscopic diffusion, in particular on the size of the mixed core, was already addressed in Sect. 5.1.1. We did not expect the measurement of the mixed core size to be modified because we have confirmed in this study that its inference is mostly independent of the model physics. However, the amount of overshooting required to produce the appropriate core size at current age does depend on the input physics.
Equation of state
Our reference Cesam2k models were computed using the OPAL05 equation of state (Rogers & Nayfonov 2002). To estimate uncertainties linked to the choice of EoS, we performed a new optimization for the target KIC5184732 using the OPAL01 EoS instead. As can be seen in Table 3, the fitted parameters all lie within 1σ errors of the results obtained with the OPAL05 EoS.
Nuclear reaction rates
We also calculated models of KIC5184732 conserving the NACRE nuclear reaction rate for the ^{14}N(p,γ)^{15}O reaction instead of the revised rate obtained from the LUNA facility (Formicola et al. 2004), which was used in Sect. 5.1. Table 3 gives the obtained fitted parameters. The amount of overshooting that is required to produce a mixed core with the appropriate size is significantly reduced. This is understandable since the NACRE cross section for the ^{14}N(p,γ)^{15}O reaction was about 30% higher than the revised LUNA rate. As a consequence, models computed with this previous cross section have a larger luminosity in the core, and thus a larger mixed core. The other fitted parameters are little modified compared to the reference fit. In particular, the size of the mixed core is unchanged, within statistical errors.
Solar mixture
We adopted the solar mixture of AGSS09 for which (Z/X_{⊙}) = 0.0181 in our reference models in Sect. 5.1. We here explored the impact of considering instead the solar mixture of Grevesse & Noels (1993; GN93), for which (Z/X_{⊙}) = 0.0244. As can be seen from Table 3, this new optimization converged toward a solution with roughly the same abundance of heavy elements as in the reference fit using AGSS09. As a result, the fitted parameters are very similar to the reference case.
5.2. Applicability to MESA models
We now address the question whether the prescription obtained for the Cesam2k code in Sect. 5.1 can be applied to mesa models. We found in Sect. 4.3 that an instantaneous overshooting with α_{ov} is roughly equivalent to a diffusive overshooting with f ~ α_{ov}/ 10, as was already pointed out in several studies before (e.g., Noels et al. 2010). At first sight this correspondence leads to believe that the mesa models require core extensions larger than the Cesam2k models. For instance, for the target KIC7206837, an instantaneous overshooting with α_{ov} = 0.18 was found necessary with Cesam2k (see Table 2), while mesa models required a diffusive overshoot parameter of f = 0.035, which would translate into α_{ov} ≈ 0.35 according to the established correspondence. However, we have mentioned in Sect. 4.2 that when using the same prescription for core overshooting (instantaneous or diffusive) and the same overshooting parameter, mesa yields core extensions that are smaller by a factor α_{MLT} compared to the extensions produced with Cesam2k. Since mesa models were computed with α_{MLT} = 1.9, the overshooting parameters obtained with mesa should be divided by a factor 1.9 to be compared to the Cesam2k overshoot parameters. By doing this, we find that the diffusive overshooting parameter of f = 0.035 obtained with mesa for KIC7206837 is equivalent to an instantaneous overshooting with α_{ov} = 0.35/1.9 ≈ 0.18 using the Cesam2k formalism, which is in agreement with the value of α_{ov} obtained with Cesam2k for this star.
To push further the comparison between Cesam2k and mesa in terms of convective core size, we checked that if the exact same formalism is used for core overshooting (and therefore the same prescription for small convective cores), the two codes provide similar sizes for the extended convective cores. For this purpose, we evolved a 1.3M_{⊙} model with both Cesam2k and mesa, either without or with overshooting. In the latter case we used an instantaneous overshooting with α_{ov} = 0.1 in both codes, and redefined in mesa the overshooting distance for small convective cores using Eq. (5) instead of Eq. (9). As shown by Fig. 15, the variations in the core size with age are very similar for Cesam2k (solid lines) and for mesa (dashed lines), both in the case without overshooting (black lines) and in the case with overshooting (red lines). The only slight differences occur right after the exhaustion of the initial ^{12}C in the core, whose burning outside of equilibrium creates the sharp peak in the core size between 15 and 25 Myr, and at the end of the main sequence, whose duration is slightly different in the two codes because of small differences in the input physics.
Fig. 15 Variations in the size of the convective core with age for a 1.3 M_{⊙} model without overshooting (black lines) and with α_{ov} = 0.1 (red lines). The gray lines indicate the Schwarzschild limit for the case with overshooting. Cesam2k models are shown as solid lines, while mesa models are represented by dashed lines. 
We thus conclude that the prescription for the overshooting parameter as a function of stellar mass obtained with Cesam2k models should also be applicable to mesa models, provided the exact same formalism is considered for core overshooting. Consistency tests such as the one presented above should be performed before applying this prescription to other stellar evolution codes.
6. Conclusion
The main result of this paper is the detection of a convective core in eight mainsequence solarlike pulsators observed with the Kepler space mission, and the asteroseismic measurement of the extent of the core in these stars.
For this purpose, we tested the seismic diagnostic for the size of the core based on the r_{010} ratios, which had been successfully applied to isolated targets before (e.g., Silva Aguirre et al. 2013) but whose general validity had not been addressed. By computing a grid of stellar models with varying mass, age, helium abundance, metallicity, and core overshooting, we established that the slope and mean value of the r_{010} ratios can be used to estimate (1) whether the star has left the main sequence or not; (2) whether it has a convective core or not; and (3) the extent of the convective core if the star possesses one. The efficiency of this diagnostic stems from the presence of a sharp μgradient at the boundary of the mixed core, which adds an oscillatory component to the r_{010} ratios. Since unevolved stars have not yet built up such a μgradient, the diagnostic is ineffective on these targets.
Based on this, we selected a subset of 24 G and lateF solarlike pulsators among Kepler targets, avoiding too unevolved stars. We extracted the oscillation mode frequencies of these stars using the complete Kepler data set (nearly four years) and fitted secondorder polynomials to the observed r_{010} ratios. At this occasion, we realized that the covariance matrix of the observables is very illconditioned, which in some cases leads the fit astray. We therefore resorted to truncated SVD to solve the problem. This issue should be kept in mind, as it can be suspected to occur in any seismic modeling where combinations of mode frequencies are used as observables, as is now frequently done (e.g., Lebreton & Goupil 2014).
By comparing the slope and mean value of the r_{010} ratios of the 24 selected targets to those of a grid of models computed with the Cesam2k code, we were able to establish that

10 of these targets are in the post main sequence and therefore do not possess convective cores;

13 targets are in the main sequence (the evolutionary status of the remaining target is uncertain) and among them eight stars have a convective core;

the convective cores of these eight targets extend beyond the classical Schwarzschild boundary.
Interestingly, identical conclusions were reached using a similar grid of models computed with the mesa code. We were able to obtain measurements of the extent of the convective cores of the eight targets that possess one, with a good agreement between the values obtained with Cesam2k and mesa. We also produced precise estimates of the stellar parameters of these eight stars that we obtained through seismic modelings. Consequently, these stars are ideal targets to test and potentially calibrate theoretical models of physical processes that could be responsible for the extension of convective cores, such as core overshooting itself or rotational mixing. Before realistic models of these processes are available, the results obtained in this paper can be used to calibrate the simple parametric models of convective core extensions that are included in most 1D stellar evolution codes.
We addressed this question using the code Cesam2k, in which cores are extended over a fraction α_{ov} of either the pressure scale height H_{P}, or the radius of the core in the sense of the Schwarzschild limit if it is smaller than H_{P}. We were able to efficiently constrain α_{ov} for the eight stars, obtaining values ranging from 0.07 to 0.18. We showed that microscopic diffusion is responsible for only a small fraction of the core extension. Interestingly, we observed a tendency of α_{ov} to increase with stellar mass, which opens the possibility to derive an empirical law for α_{ov}(M) in the mass range of observed targets (), and thus to a calibration of what is usually referred to as core overshooting, but in fact encompasses the effects of all nonstandard processes that extend convective cores. One must be careful that such a calibration necessarily depends on the prescription chosen to model the extension of convective cores in 1D stellar models. We can also suspect that it depends on the evolution code itself. We have however verified in this study that the sizes of convective cores produced by the code mesa are very similar to those produced by the code Cesam2k, provided the same prescription for core overshooting is adopted.
This study thus constitutes a first step towards the calibration of the extension of convective cores in lowmass stars. Constraints on the extent of the convective cores of more stars will be required to confirm and enrich our results. In that respect, the plato mission (Rauer et al. 2014), which was recently selected by ESA, will be particularly helpful. Reciprocally, obtaining a calibration of the distance over which convective cores extend will reduce the uncertainties on stellar ages, which will be useful to stellar physics in general, and in particular to the plato mission, for which the precise determination of stellar ages is crucial.
When fitting the modes following a Bayesian approach coupled with a Markov chain Monte Carlo algorithm, the covariance matrix can be estimated without having to assume normally distributed errors (e.g., Davies et al. 2016).
Proper uncertainties on the core sizes are obtained from optimizations in Sect. 5.1.
Acknowledgments
The authors wish to thank the anonymous referee for suggestions that helped clarify the paper. This work was performed using HPC resources from CALMIP (Grant 2015P1435). We acknowledge support from the Centre National d’Études Spatiales (CNES, France). I.M.B. and M.S.C. are supported by Fundação para a Ciência e a Tecnologia (FCT) through the Investigador FCT contract of reference IF/00894/2012 and POPH/FSE (EC) by FEDER funding through the program COMPETE. Funds for this work were provided also by the FCT research grant UID/FIS/04434/2013 and by EC, under FP7, through the project FP7SPACE2012312844. 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). V.S.A. acknowledges support from VILLUM FONDEN (research grant 10118).
References
 Angulo, C., Arnould, M., Rayet, M., et al. 1999, Nucl. Phys. A, 656, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Appourchaux, T., Michel, E., Auvergne, M., et al. 2008, A&A, 488, 705 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Appourchaux, T., Chaplin, W. J., García, R. A., et al. 2012, A&A, 543, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Appourchaux, T., Antia, H. M., Ball, W., et al. 2015, A&A, 582, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Belkacem, K., Goupil, M. J., Dupret, M. A., et al. 2011, A&A, 530, A142 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Benomar, O., Appourchaux, T., & Baudin, F. 2009, A&A, 506, 15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 BöhmVitense, E. 1958, Z. Astrophys., 46, 108 [Google Scholar]
 Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Brandão, I. M., Cunha, M. S., & ChristensenDalsgaard, J. 2014, MNRAS, 438, 1751 [NASA ADS] [CrossRef] [Google Scholar]
 Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Bruntt, H., Basu, S., Smalley, B., et al. 2012, MNRAS, 423, 122 [NASA ADS] [CrossRef] [Google Scholar]
 Burgers, J. M. 1969, Flow Equations for Composite Gases (Academic Press) [Google Scholar]
 Canuto, V. M., Goldman, I., & Mazzitelli, I. 1996, ApJ, 473, 550 [NASA ADS] [CrossRef] [Google Scholar]
 Cassisi, S., Potekhin, A. Y., Pietrinferni, A., Catelan, M., & Salaris, M. 2007, ApJ, 661, 1094 [NASA ADS] [CrossRef] [Google Scholar]
 Caughlan, G. R., & Fowler, W. A. 1988, Atomic Data and Nuclear Data Tables, 40, 283 [Google Scholar]
 Chaplin, W. J., Basu, S., Huber, D., et al. 2014, ApJS, 210, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Claret, A. 2007, A&A, 475, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cunha, M. S., & Brandão, I. M. 2011, A&A, 529, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cunha, M. S., & Metcalfe, T. S. 2007, ApJ, 666, 413 [NASA ADS] [CrossRef] [Google Scholar]
 Davies, G. R., Silva Aguirre, V., Bedding, T. R., et al. 2016, MNRAS, 456, 2183 [NASA ADS] [CrossRef] [Google Scholar]
 Degroote, P., Aerts, C., Baglin, A., et al. 2010, Nature, 464, 259 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Deheuvels, S., & Michel, E. 2010, Astron. Nachr., 331, 929 [NASA ADS] [CrossRef] [Google Scholar]
 Deheuvels, S., & Michel, E. 2011, A&A, 535, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Bruntt, H., Michel, E., et al. 2010a, A&A, 515, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Michel, E., Goupil, M. J., et al. 2010b, A&A, 514, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dintrans, B. 2009, Comm. Asteroseismol., 158, 45 [Google Scholar]
 Ferguson, J. W., Alexander, D. R., Allard, F., et al. 2005, ApJ, 623, 585 [NASA ADS] [CrossRef] [Google Scholar]
 Formicola, A., Imbriani, G., Costantini, H., et al. 2004, Phys. Lett. B, 591, 61 [NASA ADS] [CrossRef] [Google Scholar]
 Gough, D. O. 1990, in Progress of Seismology of the Sun and Stars (Berlin: Springer Verlag), eds. Y. Osaki, & H. Shibahashi, Lect. Not. Phys., 367, 283 [Google Scholar]
 Goupil, M. J., Lebreton, Y., Marques, J. P., et al. 2011, J. Phys. Conf. Ser., 271, 012032 [Google Scholar]
 Grevesse, N., & Noels, A. 1993, in Origin and Evolution of the Elements, eds. N. Prantzos, E. VangioniFlam, & M. Casse (Cambridge University Press), 15 [Google Scholar]
 Guenther, D. B., Demarque, P., & Gruberbauer, M. 2014, ApJ, 787, 164 [NASA ADS] [CrossRef] [Google Scholar]
 Herwig, F. 2000, A&A, 360, 952 [NASA ADS] [Google Scholar]
 Huber, D., Bedding, T. R., Stello, D., et al. 2011, ApJ, 743, 143 [NASA ADS] [CrossRef] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1993, ApJ, 412, 752 [NASA ADS] [CrossRef] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Lebreton, Y., & Goupil, M. J. 2012, A&A, 544, L13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lebreton, Y., & Goupil, M. J. 2014, A&A, 569, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lebreton, Y., Monteiro, M. J. P. F. G., Montalbán, J., et al. 2008, Ap&SS, 316, 1 [Google Scholar]
 Lebreton, Y., Goupil, M. J., & Montalbán, J. 2014, in EAS Publ. Ser., 65, 99 [Google Scholar]
 Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars, Astron. and Astrophys. Lib. (Berlin, Heidelberg: Springer Verlag) [Google Scholar]
 Maeder, A., & Mermilliod, J. C. 1981, A&A, 93, 136 [NASA ADS] [Google Scholar]
 Mazumdar, A., Monteiro, M. J. P. F. G., Ballot, J., et al. 2014, ApJ, 782, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Metcalfe, T. S., Monteiro, M. J. P. F. G., Thompson, M. J., et al. 2010, ApJ, 723, 1583 [NASA ADS] [CrossRef] [Google Scholar]
 Metcalfe, T. S., Creevey, O. L., Doğan, G., et al. 2014, ApJS, 214, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Miglio, A., & Montalbán, J. 2005, A&A, 441, 615 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Moravveji, E., Aerts, C., Papics, P. I., Andres Triana, S., & Vandoren, B. 2015, A&A, 580, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Morel, P., & Lebreton, Y. 2008, Ap&SS, 316, 61 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Neiner, C., Mathis, S., Saio, H., et al. 2012, A&A, 539, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Noels, A., Montalban, J., Miglio, A., Godart, M., & Ventura, P. 2010, Ap&SS, 328, 227 [Google Scholar]
 Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [Google Scholar]
 Paxton, B., Cantiello, M., Arras, P., et al. 2013, ApJS, 208, 4 [NASA ADS] [CrossRef] [Google Scholar]
 Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 168 [NASA ADS] [CrossRef] [Google Scholar]
 Pinsonneault, M. H., An, D., MolendaŻakowicz, J., et al. 2012, ApJS, 199, 30 [NASA ADS] [CrossRef] [Google Scholar]
 Popielski, B. L., & Dziembowski, W. A. 2005, Acta Astron., 55, 177 [NASA ADS] [Google Scholar]
 Provost, J., Berthomieu, G., Bigot, L., & Morel, P. 2005, A&A, 432, 225 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rauer, H., Catala, C., Aerts, C., et al. 2014, Exp. Astron., 38, 249 [NASA ADS] [CrossRef] [Google Scholar]
 Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064 [Google Scholar]
 Roxburgh, I. W. 1978, A&A, 65, 281 [NASA ADS] [Google Scholar]
 Roxburgh, I. W. 1985, Sol. Phys., 100, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Roxburgh, I. W. 2009, A&A, 493, 185 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Roxburgh, I. W., & Vorontsov, S. V. 2003, A&A, 411, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Samadi, R., Kupka, F., Goupil, M. J., Lebreton, Y., & van’t VeerMenneret, C. 2006, A&A, 445, 233 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Saslaw, W. C., & Schwarzschild, M. 1965, ApJ, 142, 1468 [NASA ADS] [CrossRef] [Google Scholar]
 Saumon, D., Chabrier, G., & van Horn, H. M. 1995, ApJS, 99, 713 [NASA ADS] [CrossRef] [Google Scholar]
 Scuflaire, R., Montalbán, J., Théado, S., et al. 2008, Ap&SS, 316, 149 [NASA ADS] [CrossRef] [Google Scholar]
 Shaviv, G., & Salpeter, E. E. 1971, ApJ, 165, 171 [NASA ADS] [CrossRef] [Google Scholar]
 Silva Aguirre, V., Ballot, J., Serenelli, A. M., & Weiss, A. 2011, A&A, 529, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Silva Aguirre, V., Casagrande, L., Basu, S., et al. 2012, ApJ, 757, 99 [NASA ADS] [CrossRef] [Google Scholar]
 Silva Aguirre, V., Basu, S., Brandão, I. M., et al. 2013, ApJ, 769, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Silva Aguirre, V., Davies, G. R., Basu, S., et al. 2015, MNRAS, 452, 2127 [NASA ADS] [CrossRef] [Google Scholar]
 Stancliffe, R. J., Fossati, L., Passy, J.C., & Schneider, F. R. N. 2015, A&A, 575, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Steigman, G. 2010, J. Cosmol. Astropart. Phys., 4, 29 [Google Scholar]
 Stello, D., Bruntt, H., Preston, H., & Buzasi, D. 2008, ApJ, 674, L53 [NASA ADS] [CrossRef] [Google Scholar]
 Thoul, A. A., Bahcall, J. N., & Loeb, A. 1994, ApJ, 421, 828 [NASA ADS] [CrossRef] [Google Scholar]
 Valle, G., Dell’Omodarme, M., Prada Moroni, P. G., & Degl’Innocenti, S. 2016, A&A, 587, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 VandenBerg, D. A., Bergbusch, P. A., & Dowler, P. D. 2006, ApJS, 162, 375 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Ventura, P., Zeppieri, A., Mazzitelli, I., & D’Antona, F. 1998, A&A, 334, 953 [NASA ADS] [Google Scholar]
 Zahn, J. 1991, A&A, 252, 179 [NASA ADS] [Google Scholar]
Appendix A: Mode parameters of the selected targets
Tables A.1 through A.6 give the fitted frequencies of the oscillation modes in the 24 selected stars among Kepler targets (see Sect. 3).
Mode parameters of KIC8394589, KIC9098294, KIC9410862, and KIC6225718.
Mode parameters of KIC10454113, KIC6106415, KIC10963065, and KIC6116048.
Mode parameters of KIC5184732, KIC3656476, KIC7296438, and KIC4914923.
Mode parameters of KIC12009504, KIC8938364, KIC7680114, and KIC10516096.
Mode parameters of KIC7206837, KIC8176564, KIC8694723, and KIC12258514.
Mode parameters of KIC6933899, KIC11244118, KIC7510397, and KIC8228742.
Appendix B: Polynomial fit to the observed r_{010} ratios
We denote as y = (y_{1},...,y_{n}) the set of n observed r_{010} ratios, and as x = (x_{1},...,x_{n}) the corresponding frequencies. We would like to fit a 2ndorder polynomial to the functional y(x). Since we would like to interpret individual coefficients of the polynomial regression, we need them to be independent from one another. We are thus required to use orthogonal polynomials. For this purpose, we fit polynomials of the type (B.1)to the observed ratios y, where β, γ_{1}, and γ_{2} need to be determined so as to ensure the independence of the a_{k} coefficients. The values taken by the polynomial P(ν) at each point of vector x can be written in a vectorial way as (B.2)where a = (a_{0},a_{1},a_{2}), and J is an n × 3 matrix such that In our case, the observables y are combinations of mode frequencies and are thus highly correlated. We denote as C the covariance matrix of the observables. At first, let us assume that C is invertible and denote its inverse as W. In this case, the χ^{2} function to be minimized can be written as (B.6)where the exponent T indicates matrix transposition. The gradient of this function is given by (B.7)and the optimal set of coefficients is then (B.8)where corresponds to the error matrix for the coefficients a_{k}. To ensurethat these coefficients are uncorrelated, we must thus require that the nondiagonal coefficients of the error matrix vanish. This yields the set of equations which can be solved to yield β, γ_{1}, and γ_{2}.
When fitting polynomials to the r_{010} ratios of models, the errors on the mode frequencies are assumed to be zero, so the covariance matrix and its inverse W are the identity, which simplifies the set of Eqs. (B.9) to (B.11). In particular, β reduces the mean of the vector of frequencies x.
For polynomial fits to the observations, one complication arises. As mentioned in Sect. 3.3, the determinant of the covariance matrix C of the observed r_{010} ratios is vanishingly small. As a result, C is almost noninvertible and numerical problems arise when trying to invert it. To remedy this, we resorted to a truncated SVD (singular value decomposition) approach. We computed the eigenvalues of matrix C, further denoted , and its eigenvectors denoted as matrix P, where P_{k} is the eigenvector corresponding to the kth eigenvalue. The covariance matrix C can thus be rewritten as P^{T}DP, where D is the diagonal matrix whose diagonal coefficients are the λ_{i}.
The conditioning of matrix C, given by the ratio λ_{1}/λ_{n}, is very large (of the order of 10^{6}). We thus truncate the n−m smallest eigenvalues in order to improve the conditioning. For this purpose, the covariance matrix is approximated by the matrix , where is a diagonal matrix whose diagonal coefficients are the m largest λ_{i} and is an (m,n) matrix composed of the m eigenvectors corresponding to the retained eigenvalues. For each star, we suppressed eigenvalues until the impact of suppressing an eigenvalue on the results of the polynomial fit is negligible. We found that suppressing n−m = 5 eigenvalues is generally enough to ensure the latter condition.
The χ^{2} function to be minimized becomes (B.12)where . In this case, the values of β, γ_{1}, γ_{2} required to ensure the independence of the a_{k} are given by Eqs. (B.9) to (B.11) where the coefficients of matrix W are to be replaced by those of matrix . The bestfit coefficients a_{k} are obtained by doing the same thing with Eq. (B.8).
All Tables
Fitted parameters and characteristics of the bestfit models obtained for the stars of the sample that have a convective core.
Characteristics of the bestfit models obtained for KIC5184732 when modifying the chosen input physics.
All Figures
Fig. 1 Variations in the ratio r_{010} around ν_{max} as a function of frequency for models of 1.2 M_{⊙} from the ZAMS (dark blue) to the beginning of the post main sequence (dark red). The dashed lines correspond to fits of 2nd order polynomials. 

In the text 
Fig. 2 Variations in the pressure scale height H_{P} (dashed line) and the radius of the extended convective core R_{c} (solid blue line) with age for a 1.3 M_{⊙} Cesam2k model with solar metallicity, a solarcalibrated value for the mixing length, Y_{0} = 0.26, and α_{ov} = 0.1. The gray solid line indicates the Schwarzschild limit. 

In the text 
Fig. 3 Left: evolutionary tracks of stellar models of 1.2 M_{⊙} in the (a_{1},a_{0}) plane for different amounts of core overshooting: α_{ov} = 0 (gray), 0.1 (blue), 0.15 (cyan), 0.2 (green), 0.25 (red), and 0.3 (purple). Full (resp. dashed) lines indicate that the model has a convective (resp. radiative) core. Right: variations in the size of the convective core as a function of age for the same models. 

In the text 
Fig. 4 Location of models in the (a_{1},a_{0}) plane at fixed Δν. Colors indicate the amount of core overshooting: α_{ov} = 0 (gray), 0.1 (blue), 0.15 (cyan), 0.2 (green), 0.25 (red), 0.3 (magenta). Open squares indicate models with a convective core, and crosses, models with radiative cores. The black open circles indicate models that are in the postmainsequence (X_{c}< 10^{2}). 

In the text 
Fig. 5 Ratios r_{010} computed for KIC 6106415 (left) and KIC 12258514 (right) using the mode frequencies extracted from the Kepler oscillation spectra (see text). The colored dashed lines correspond to 2ndorder polynomial fits to the observed ratios using either the raw covariance matrix (gray lines) or the covariance matrix modified through truncated SVD (see Sect. 3.3). 

In the text 
Fig. 6 Location in the (a_{1},a_{0}) plane (star symbols and black error bars) of the stars of the sample that were found to be on the MS with a convective core in this study. Models that reproduce the observed large separation, the spectroscopic estimate of metallicity, and the stellar mass derived from scaling laws within 3σ errors are overplotted. The symbols have the same meaning as in Fig. 4. 

In the text 
Fig. 7 Location in the (a_{1},a_{0}) plane of the MS stars for which the presence of a convective core is uncertain (note that the evolutionary status of KIC9410862 is in fact uncertain). The symbols have the same meaning as in Fig. 6. 

In the text 
Fig. 8 Location in the (a_{1},a_{0}) plane of the first six PoMS stars of the sample. The symbols have the same meaning as in Fig. 6. 

In the text 
Fig. 9 Location in the (a_{1},a_{0}) plane of the last four PoMS stars of the sample. The symbols have the same meaning as in Fig. 6. 

In the text 
Fig. 10 Fractional mass of the convective core for the eight stars that were found to have a convective core in this study. For each star, the open symbols correspond to the core size of the five models of the two grids (blue squares for Cesam2k models, red circles for mesa models) that yield the lowest values of χ^{2} as defined by Eq. (7). The filled squares give the core sizes obtained from a LevenbergMarquardt optimization and the evolution code Cesam2k (see Sect. 5.1). 

In the text 
Fig. 11 Location in the (a_{1},a_{0}) plane of four stars of the sample compared to the location of models computed with the evolution code mesa. Symbols are the same as in Fig. 6, except for the colors, which indicate diffusive overshooting parameters of: f = 0.004 (gray), 0.010 (blue), 0.016 (cyan), 0.022 (green), 0.028 (red), 0.035 (magenta). 

In the text 
Fig. 12 Profile of r_{010} for two mesa models: one with a diffusive overshooting (f = 0.02) and a 1.3M_{⊙} mass (blue squares) and the other with a step overshooting (α_{ov} = 0.22) and a 1.31M_{⊙} mass (red circles). Both models have the same mean large separation. 

In the text 
Fig. 13 Amount of core overshooting found for the stars of the sample that have a convective core as a function of the fitted stellar mass (left) and as a function of the fitted initial metallicity (right). Blue squares indicate models computed without microscopic diffusion and gray circles, models where microscopic diffusion is included following Burgers (1969). The vertical arrows indicate upper limits of α_{ov} (see Sect. 5.1.2). 

In the text 
Fig. 14 Value of the χ^{2} of the bestfit model as a function of the (fixed) amount of core overshooting for the star KIC 10516096. 

In the text 
Fig. 15 Variations in the size of the convective core with age for a 1.3 M_{⊙} model without overshooting (black lines) and with α_{ov} = 0.1 (red lines). The gray lines indicate the Schwarzschild limit for the case with overshooting. Cesam2k models are shown as solid lines, while mesa models are represented by dashed lines. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.