Issue 
A&A
Volume 591, July 2016



Article Number  A99  
Number of page(s)  10  
Section  Stellar structure and evolution  
DOI  https://doi.org/10.1051/00046361/201628311  
Published online  21 June 2016 
Asteroseismology of 19 lowluminosity red giant stars from Kepler
^{1} Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain
email: fph@iac.es
^{2} Universidad de La Laguna, Dpto. Astrofísica, 38206 La Laguna, Tenerife, Spain
^{3} Laboratoire AIM, CEA/DRFCNRS, Université Paris 7 Diderot, IRFU/SAp, Centre de Saclay, 91191 GifsurYvette, France
^{4} Royal Observatory of Belgium, Ringlaan 3, 1180 Brussels, Belgium
^{5} Institute of Astronomy, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium
Received: 15 February 2016
Accepted: 28 April 2016
Context. Frequencies of acoustic and mixed modes in red giant stars are now determined with high precision thanks to the long continuous observations provided by the NASA’s Keplermission. Here we consider the eigenfrequencies of nineteen lowluminosity red giant stars selected in a recent publication for a detailed peakbagging analysis.
Aims. Our objective is to obtain stellar parameters by using individual mode frequencies and spectroscopic information.
Methods. We use a forward modelling technique based on a minimization procedure combining the frequencies of the pmodes, the period spacing of the dipolar modes, and the spectroscopic data.
Results. Consistent results between the forward modelling technique and values derived from the seismic scaling relations are found but the errors derived using the former technique are lower. The average error for log g is 0.002 dex, compared to 0.011 dex from the frequency of maximum power, ν_{max}, and 0.10 dex from the spectroscopic analysis. Relative errors in the masses and radii are on average 2% and 0.5% respectively, compared to 3% and 2% derived from the scaling relations. No reliable determination of the initial helium abundances and the mixing length parameters could be made. Finally, for our grid of models with given input physics, we found that lowmass stars require higher values of the overshooting parameter.
Key words: asteroseismology / stars: general / stars: interiors
© ESO, 2016
1. Introduction
The Keplermission from NASA (Borucki et al. 2009) and its recent version K2 (Howell et al. 2014) are providing individual eigenfrequencies for a huge number of stars, including thousands of red giants (e.g. Stello et al. 2013, 2015; Chaplin et al. 2015). These data allow us to determine accurate stellar properties that help to constrain stellar evolution models (Metcalfe et al. 2010; Mathur et al. 2012) and improve the determination of the properties of the exoplanets they might host (Guillot & Havel 2011; Huber et al. 2013). A review of the progress with Keplerin the field of the red giant stars can be found in ChristensenDalsgaard (2012) and García & Stello (2015).
Obtaining stellar properties from pulsation spectra can be done with a variety of techniques such as forward modelling that seeks those models whose frequencies best match the observed ones, or inverse modelling that consists of some (usually linearised) relationship between the frequencies and the stellar structure. Intermediate approaches take advantage of some analytical or asymptotic approximation to form frequency combinations aimed at isolating some aspect of the stellar structure. In this paper we use a forward modelling approach and compare the results with those obtained from simple scaling relations and acoustic helium signature fits, the socalled acoustic glitches (e.g. Vorontsov 1988; Gough 1990). A similar comparison for mainsequence and subgiants stars was done by Metcalfe et al. (2014). They found that the uncertainties in the masses and radii were improved by a factor of three when individual frequencies were fitted compared to the use of empirical scaling relations.
Target stars and their observed parameters.
The interest in making such a comparison is worth emphasizing: on the one hand, the forward technique is modeldependent, introducing systematic errors into the stellar parameters, but the structural models and frequency computations can be done with uptodate physics. On the other hand, the scaling relations have wide observational support (see Bedding et al. 2003; Stello et al. 2008; and Mosser et al. 2013 work specifically on red giants), but they are not fully understood from a theoretical point of view. This mainly concerns the relation between the frequency at maximum power, ν_{max}, and the acoustical cutoff frequency, ν_{cut}, which is supposed to justify the scaling relation , where g is the surface gravity and T_{eff} the effective temperature. Although some theoretical work has been done in an effort to understand this scaling relation (Belkacem et al. 2013), neither its extent or accuracy is clear. In fact Jiménez et al. (2015) were able to measure the cutoff frequency for several stars and found a better agreement between ν_{max} and ν_{cut} than is suggested by Belkacem et al. (2013). Thus both techniques should be regarded as complementary, and confronting their results enables us to gain confidence in the values and uncertainties of the stellar parameters derived. The same is true when comparing results derived from forward modelling and a glitch fit, the latter being modelindependent but making use of an asymptotic approximation.
The splitting of mixed modes show that redgiant cores rotate faster than their convective envelope (Beck et al. 2012; Deheuvels et al. 2012) and opens up the possibility of probing their internal rotation rates (Deheuvels et al. 2014, 2015; Di Mauro et al. 2016). In fact, for our target stars Corsaro et al. (2015a) determined the frequencies of many mixed dipolar modes including a high number of rotational splittings from which one can extract information on the internal rotation of this kind of lowluminosity red giant stars. We will analyse such information in a separate paper.
2. Observations and models
2.1. Observational data
The target stars considered in the present work and analysed by Corsaro et al. (2015a) are listed in Table 1. They are lowmass, lowluminosity red giant stars (specifically only stars with ν_{max}> 110 μHz were considered) observed by Keplerover more than four years and were selected because of their good S/N and the availability of gravity period spacing measurements from Mosser et al. (2012). Mode frequencies were obtained using the Bayesian tool Diamonds (Corsaro & De Ridder 2014). Here one needs to implement a model for the expected frequency pattern, and while the ℓ = 1 mixed modes were fitted individually, for each ℓ = 2 and ℓ = 3 peak only a single Lorentzian profile was used (see Corsaro et al. 2012 for more details). This means that the ℓ = 2 frequencies can be affected by the presence of mixed quadrupole modes and related rotationally split components. Unfortunately it is not possible to disentangle quadrupole mixed modes and rotational split components because of their high density in the ℓ = 2 frequency region.
We include frequencies from the nineteen stars analysed by Corsaro et al. (2015a) but exclude modes identified with a probability lower than 99%, as suggested by their peak significance test. The large separations Δν given in Table 1 are originally from Mosser et al. (2012) whereas the frequencies of maximum power, ν_{max}, are a byproduct of the peak bagging carried out by Corsaro et al. (2015a). On the other hand, values of effective temperature (T_{eff}), surface gravity (log g), and surface metallicity (Z/X) were taken from the APOKASC Catalogue (Pinsonneault et al. 2014) except for stars D, Q and S for which spectroscopic data were not available, and photometric values of T_{eff} from Pinsonneault et al. (2012) were used as input parameters.
2.2. Stellar code and pulsation
Model fitting is based on a grid of stellar models evolved from the premain sequence with the MESA code (Paxton et al. 2011), version number 7184. Models were computed with the OPAL opacities (Iglesias & Rogers 1996) and GS98 metallicity mixture (Grevesse & Sauval 1998). Microscopic diffusion of elements was included; otherwise the standard MESA input physics was used. The choice of the metallicity mixture and its observed Z/X value for the Sun was decided since as opposed to recent spectroscopic estimations (Asplund et al. 2005, 2009), it gives a good agreement between solar standard models and helioseismic observations (see e.g. Basu & Antia 2008).
The starting grid is composed of evolution sequences with masses (M) from 0.85 M_{⊙} to 2.0 M_{⊙}, initial helium abundances (Y_{ini}) from 0.25 to 0.32, initial metalicities (Z_{ini}) from 0.0025 to 0.04, and mixing length parameters (α) from 1.5 to 2.2. The density of the original grid was increased to ensure that at least two values of every parameter were found in the solutions obtained after adding random noise to the data. This was checked for all the stars. Specifically, the final mass step in the range 0.9 ≤ M ≤ 1.5 is 0.025 M_{⊙} and the metallicity step for Z ≤ 0.025 is 0.0025. We do not find it necessary to increase the original steps of 0.01 and 0.1 for Y_{ini} and α, respectively.
No overshooting was considered in this global grid, but for each star, after the other stellar parameters were fixed, new models were computed using the exponential prescription of Herwig (2000). Here the particle spreading in the overshoot region can be described as a diffusion process with a diffusion coefficient D_{ov} given by (1)where z is the distance from the edge of the convection zone, H_{v} the velocity scale height of the overshooting convective elements at the edge of the convection zone, and H_{p} the pressure scale height at the same point. For the free parameter f_{ov}, we have considered values from 0 to 0.03 in steps of 0.0006. The same f_{ov} value was used for the core and the envelope and throughout the evolutionary sequence. Although, as specifically implemented in MESA, the formulation corresponds to that introduced by Herwig (2000) to investigate the overshooting on AGB stars, it is in fact a simplified version of the formulation given by Freytag et al. (1996), who analysed the envelope overshooting in solarlike stars, main sequence A type stars and white dwarfs. Recently, Moravveji et al. (2015) carried out a seismic analysis of core overshooting in a main sequence B star, obtaining satisfactory results compared to other overshooting formulations.
Frequencies were computed with ADIPLS code (ChristensenDalsgaard 2008). The code uses the adiabatic approximation and neglects the interaction between convection and oscillations. Mode degrees from ℓ = 0 to ℓ = 3 were considered.
For a typical evolutionary sequence in the initial grid, we save between 100 and 200 models, from the subgiant phase to the red giant phase with an upper radius of about 8 R_{⊙}. Owing to the very rapid change in the dynamical time scale of the models, such grids are too coarse in the time steps. Nevertheless, as detailed in Sect. 4, we have checked that interpolations between models provide estimations of the pmode frequencies, the period spacing of the ℓ = 1 modes, and the stellar parameters with errors much lower than the observational ones. We have not attempted to fit single mixed modes; hence this procedure is safe and consumes relatively less time.
3. Fitting procedure
We have considered a χ^{2} minimization method including simultaneously mode frequencies and spectroscopic data. Specifically we minimize the function (2)Regarding the spectroscopic parameters, when available we have included the effective temperature (T_{eff}) the surface gravity (log g) and the surface metallicity (Z/X); namely (3)where δT_{eff}, δ(Z/X) and δg correspond to differences between the observations and the models whereas σ_{Teff}, σ_{ZX} and σ_{g} are their respective observational errors.
The other three terms in Eq. (2)are determined from the mode frequencies. The term is aimed at minimizing the mean density through a term related to the large separation. In principle this term is not necessary in the minimization since the same information can be included in the term with the frequency differences. However, owing to the socalled surface effects not considered in the model and frequency computations, there will be some discrepancies between the large separation Δν of the models and that of the observations. For instance, for the Sun, using a frequency interval around ν_{max} representative of our set of stars, we obtain for the mean difference between adjacent radial modes a value of Δν = 134.8 μHz for the observations and Δν = 136.0 μHz for a solar model; that is a difference of 0.9%. Although ultimately a solar calibration could be performed, hopefully cancelling some of the uncertainties, it seems better to fix the relevant constant in such a way that the discrepancy in the solar case is removed as far as possible.
Introducing the dynamical time, t_{dyn} = (R^{3}/GM)^{1/2}, and the dimensionless frequencies σ_{nl} given by ω_{nl} = t_{dyn}σ_{nl}, the relative frequency differences between models and observations for radial modes can be expressed as (4)Hence, one might expect that fitting the frequency differences for radial oscillations to a function of frequency, namely (5)where A_{0} is a constant, P_{i} a Legendre polynomial of the order i, and x corresponds to 1 /ω linearly scaled to the interval [−1,1 ], the constant term A_{0} will be close to zero for models with the correct mean density. In what follows a value of k = 2, corresponding to a parabolic function, has been adopted. We have checked that when values of k = 3,4 are considered, the results for A_{0} are the same to within the errors.
In practice, even for the best model, unknown surface effects will introduce a term in F(σ), possibly including a constant that should be translated into the form of an uncertainty in the determination of A_{0}. In fact, when considering frequency differences between the Sun and a solar model, and limiting again the range of radial orders around ν_{max} to that typically observed in the stars considered here, we obtain A_{0} = 0.001. In this way, the error introduced in the solar mean density is almost an order of magnitude smaller than the discrepancy derived from the large separation computed as a simple average of frequency differences. It is worth mentioning that, had we used the full range of known radial modes for the Sun, the surface term could be isolated by taking into account that at very low frequencies such terms tend to zero.
For the minimization procedure we define the quantity (6)where A_{00} is an offset caused by the simplified physics used in the model and frequency computations. In principle one can take the solar value, A_{00} = 0.001, but we did not find this offset completely satisfactory for our target stars and we discuss it later in this paper. On the other hand σ_{A0} is the error associated to A_{0} and in principle could also be taken as the discrepancy σ_{A0} = 0.001 found for the Sun. We note that this value is at least one order of magnitude higher than the formal error found in a typical fit to Eq. (5): hence, considering this higher uncertainty is our main reason for dealing with the terms and separately. In other words, and as suggested by Eq. (4), will be intended for a minimization of the dimensionless frequencies. Its computation is detailed in the next paragraph.
The term corresponds to the frequency differences of the ℓ = 0,2,3 modes after removing a smooth function of frequency. The surface term is computed using only radial oscillations in a similar way to Eq. (5)except that, as suggested by our tests (see Sect. 4 below), the frequency differences were scaled with the dimensionless energy I_{nl}, defined as in Aerts et al. (2010), (7)Here we have also adopted a value of k = 2. Tests for some stars show that higher values do not significantly reduce the minimum value of χ^{2} while using k = 1 gives substantially higher values.
Then, in the minimization procedure we consider radial as well as ℓ = 2,3 modes. The corresponding function to be minimized is (8)where j runs for all the modes with degrees ℓ = 0,2,3, N is the number of modes considered in the fit, and σ_{ωj} the relative error in the frequency ω_{j}. We note that the polynomial function subtracted in Eq. (8) includes the constant coefficient, since a similar term (but with a very different uncertainty) was already included in .
As noted in Sect. 2.1, the observed modes include only one ℓ = 2 and one ℓ = 3 peak in every Δν interval. To mimic the observations, at least to a first approximation, for each observed peak with ℓ = 2,3 we have taken an average of all the eigenfrequencies with the same degree and within a frequency interval of 2 μHz around the mode with the lowest dimensionless energy and the correct asymptotic radial order, weighted by where is an interpolation of the dimensionless energy of the ℓ = 0 modes to the frequencies of the nonradial modes. A value of is assigned to the observed mode in Eq. (8). In practice, for modes with degree ℓ = 3 this is basically equivalent to searching for the corresponding “pure” pmode, but for modes with degree ℓ = 2, an average between two modes is often required.
The term corresponds to differences in the period spacing of the ℓ = 1 modes, ΔΠ_{1}. However for the sake of rapidity and robustness in the computations, we have used a simpler parameter, based on the work by Jiang & ChristensenDalsgaard (2014). The main simplification compared to that work is that we fix the values of the hypothetical ℓ = 1 pure pmodes by using the frequencies of ℓ = 0 and 2 modes, only fitting the coupling parameter and the period spacing that are assumed to be the same in the whole frequency range. We are not interested in using accurate equations for obtaining precise values of the period spacings but rather in using the same simple fit for the observations and the models.
The detailed computation is as follows. First we compute the small separation, δν_{02} between adjacent ℓ = 0 and ℓ = 2 modes as an average over all available pairs. We then estimate the frequencies of the hypothetical pure ℓ = 1 modes with the asymptotic equation as . Afterwards a first estimate of the period spacing is obtained with a linear fit of the dipolar period spacings, Π_{j} = 1 /ν_{j}−1 /ν_{j + 1} to a second order polynomial function of the variable (9)where the bar denotes average of two consecutive dipolar modes, j and j + 1, q_{0} is a guessed initial value for the coupling parameter q, and where is the closest pure pmode to ν_{j}. The period spacing ΔΠ_{1} corresponds to the zerothorder polynomial^{1}. This value is used as an initial guess for a nonlinear fit. Introducing the global parameters ω_{g} = 2π^{2}/ ΔΠ_{1} and ω_{p} = 2Δν, and for every dipolar mode with angular frequency ω_{j} = 2πν_{j} the functions f_{1j} = q^{2} + 2/3(πΔν_{j}/ Δν)^{2}, , and f_{3j} = 1 + f_{2j}/ω_{g}, the period spacings are fitted to (10)where ω_{g} and q are the coefficients to be determined.
The corresponding minimization function is given by (11)where σ_{ΔΠ1} is the uncertainty in the period spacing. The formal errors resulting from the fit are too small and, if they were to be adopted, the minimization procedure based on Eq. (2)would overweight the period spacing compared to the pmode frequencies or the spectroscopic parameters. This could cause the undesirable effect that many global and envelope structural properties become determined mainly by models with the correct ΔΠ_{1}. To prevent such a situation, we have used the following criteria. First, we consider a high value of the uncertainty, namely σ_{ΔΠ1} = 10σ_{A0}, to ensure that the other terms in Eq. (2)will give the proper global stellar parameters (more precisely, they will fix the model input parameters M, Y_{ini}, Z_{ini} and α). Then, we use a grid of models with different values of the overshooting parameter f_{ov} while introducing in Eq. (11)the formal uncertainties for the period spacing, but still imposing the condition σ_{ΔΠ1} ≥ σ_{A0}. We have checked that this worked as expected in the sense that some properties of the best fitted models, such as the core size or the age, are changed in the second iteration whereas others, such as the depth of the second helium ionization zone, remain almost unchanged.
4. Testing the methodology
As noted earlier, our minimization procedure is based on a grid that is too coarse to obtain the bestfitted models. Hence we interpolate their parameters and frequencies to a finer grid in ages. We then search for the minimum χ^{2} corresponding to every evolutionary sequence. We have found that the values of χ^{2} for our final models (which are explicitly computed from the premain sequence to the fitted age at the end of the minimization procedure) agree with the interpolated ones. As an example, Fig. 1 shows χ^{2} for an evolutionary sequence in the original grid with parameters matching those of KIC 003744043 (letter A in Tables 1 and 2). We show values for all the terms in Eq. (2): χ^{2} (black), (blue), (red), (green), and (magenta). The open circles correspond to models in the grid, the solid lines to values interpolated, and the full circles to the best fitted model. Since the example shown is for the best evolutionary sequence (within errors), all the terms have their minima at very close ages. To avoid misinterpretations, we note that this figure does not provide an indication of the uncertainty in the age since other parameter combinations can also give similar good fits but at different ages.
Although Fig. 1 illustrates that the time step in our grid is fine enough, it can also be seen that the dependence of with age is not completely smooth, even in this very short range of ages. Including the individual frequencies of the mixed modes in the analysis with the goal of improving the results, might result in more irregular χ^{2} functions than those shown in Fig. 1 and hence further tests would be required to validate the procedure.
Fig. 1 Values of χ^{2} as function of age for an evolutionary sequence with stellar parameters that match those of KIC 003744043 within errors (letter A in Tables 1 and 2). The black open circles are computed values for models and frequencies, whereas the continuous black line corresponds to parameters and frequencies interpolated to intermediate ages. The filled black circle is the χ^{2} value derived from the bestfit model explicitly computed at the end of the procedure to exactly match the interpolated age. Also shown are values of (blue), (red), (green), and (magenta). Here the solid lines and the full and open circles have the same meaning as before. 
Let us consider the offset A_{00} introduced in Eq. (6). As stated above, for the Sun we obtain A_{00} = 0.001 but our tests show that using that number for the red giants gives rise to positive values of the frequency differences ω(observations)−ω(best model) in the low frequency range, at least for some of the stars. Since the missed surface effects are expected to overestimate the theoretical frequencies, it seems natural to impose the condition that the frequency differences will be negative for the whole frequency range, at least to within the dispersion of the minimization procedure. For our sample of stars, a higher than solar value is required with a minimum offset of A_{00} = 0.003 to 0.004. We have finally taken A_{00} = 0.004 and σ_{A0} = 0.001 in Eq. (6). Alternatively, had we taken A_{00} = 0 but explicitly imposed the condition that the frequency differences cannot be positive, we would have found mean densities for our best models that differed from those reported here by a factor of 0.9986 ± 0.0015. Since the source of error for the mean density is mainly that of A_{0}, this dispersion is consistent with the input value of σ_{A0} = 0.001 adopted as the estimate of the uncertainties in A_{0}. We are aware that the constant shift depends on the physics of the models used (including the eigenfrequency computations), and also on the range of radial order spanned by the observed frequencies. Nevertheless, the former are the same and the latter very similar for all of our target stars. A comparison with a calibrated scaling relationship will be given in Sect. 5. In Fig. 2 we show the frequency differences between our best models and the observations for all the radial modes and the full set of stars. As can be seen, when plotted against the normalized frequency ν/ν_{max} they are similar for all the stars.
Fig. 2 Frequency differences between the observations and the best models, [ ω(observations)−ω(best model) ] /ω (red points), and associated dimensionless energies in arbitrary units (blue squares) for radial oscillations. For a better comparison between modes in different stars, the horizontal axis gives the frequencies normalized to the frequency of maximum power. All nineteen stars are included. The black line is a secondorder polynomial fit. 
Fig. 3 Frequency differences between the observations and the best models after subtracting the surface effects. Red points are for ℓ = 0, blue points for ℓ = 2 and green points for ℓ = 3. All nineteen stars are included. 
Figure 3 shows the residuals, that is, the relative frequency differences between the best models and the observations after subtracting the surface effects, for modes with degrees ℓ = 0,2,3 and all the stars. The mean residuals for modes with degrees ℓ = 0 and ℓ = 3 is 0.0002 whereas for modes with ℓ = 2 we obtain 0.0004 . These figures can be compared with the observational errors of 0.00008, 0.00010, and 0.00007 for modes with degrees ℓ = 0,2 and 3 respectively. Thus for the ℓ = 0 and 3 our residuals are about two times higher than the observational errors whereas they are about four times higher for the ℓ = 2. The higher value for the ℓ = 2 modes is probably caused by the single peak considered in the fit to the observed spectrum to what is regularly a pair of quadrupolar modes of mixed character (prior to any rotational splitting consideration). The way we have dealt with this issue (see the paragraph after Eq. (8)) is therefore not entirely satisfactory.
Figure 2 shows that the surface terms missing in the theoretical computations give rise to a smoother frequency dependence in the relative frequency differences (red points) than that of the dimensionless energies (blue points) at the lowest frequencies, indicating that the discrepancies between the theory and the observations should involve layers below the upper turning points of the lowest observed radial frequencies. It also suggests computing by using a simple polynomial fit to δω/ω rather than to Iδω/ω. In that case the surface term will be computed using Eq. (5)while Eq. (8)should be replaced by (12)where is an interpolation of the dimensionless energy of the ℓ = 0 modes to the frequencies of the nonradial modes. If a secondorder polynomial is considered, as in our reference case, we find residuals of 0.0003 for the ℓ = 0 and 3 modes and 0.0006 for the ℓ = 2, which is about 1.4 times higher than when Eq. (8)is used. We thus decided to continue using Eq. (7)for subtracting the surface effects.
Fig. 4 Period spacing obtained in this work compared to that in Corsaro et al. (2015a). Red points are values with errors lower than 0.1%. The continuous line is x = y and the dashed one y = x−2.5 s. 
As noted above, the procedure used here for determining the period spacing is a simplified version of that in Jiang & ChristensenDalsgaard (2014), and the resulting equation for the period spacing, Eq. (10), is highly inaccurate compared to the observational errors. Specifically we found that for the best models the individual periods deviate on average some 2.5 s from their fitted values. Nevertheless we expect a similar behaviour for the observed modes and, hence, most of the differences will be cancelled out, provided the observed and theoretical ΔΠ_{1} values are computed in the same way. In fact the average value for the standard deviations between the individual periods and their fitted values is 4 s for the observations, 1.5 times higher than the one derived from the best fitted models. Figure 4 shows the period spacing for the nineteen stars computed with Eq. (10)against the ΔΠ_{1} values given by Corsaro et al. (2015a based on the works by Mosser et al. 2012, 2014). Red points corresponds to those ΔΠ_{1} values determined here with an error lower than 0.1% The main difference is that our period spacings are 2.5 s smaller. The blue point on the right with a period of ΔΠ_{1} = 88 s according to Corsaro et al. (2015a) corresponds to KIC 008366239 for which we only used six modes in the fit.
5. Results
5.1. Individual stars
Parameters of the best models.
Table 2 summarizes the results for all the stars in our sample. To estimate the errors in the output parameters we add normally distributed errors to the observed frequencies, the coefficient A_{0}, and the spectroscopic parameters and search for the model with the minimum χ^{2} in every realization. In this way we estimate mean and 1σ values for χ^{2} and the stellar parameters. The last line in Table 2 shows mean σ values. Additionally we have redone all the analysis with a grid with a density along the mass axis of half the original. Compared to the full grid, a dispersion of 3% in the mass is obtained as the mean for the 19 stars. Since other parameters are correlated with the masses, they also change. In particular we obtain a mean dispersion of 1% in the radii, 9% in the luminosity, 0.96 Gyr in the age and 177 s in τ_{HeII}. These values are higher than those given in Table 2 but they should be regarded as upper limits to the uncertainties introduced by the limited number of models in the grid.
Table 3 gives the correlation matrix for which elements have been computed as averages of the linear Pearson correlation coefficients for all the stars in our dataset. Results are limited to the first step in the minimization procedure where only models without overshooting are considered and five parameters are changed: M, Y_{ini}, Z_{ini}, α, and age. The high correlation between M and Y has also been found in seismic analysis of main sequence stars, e.g. Metcalfe et al. (2014). There is also a very high correlation between the age and Z_{ini}. Had we used the mean density rather than age as the fifth parameter, the correlation coefficients on the last column would become very close to zero.
Fig. 5 Minimum values of χ^{2} as a function of M – changing Y_{ini}, Z_{ini} and α – for all the stars. Points with error bars correspond to χ^{2} values with 1σ uncertainties. The letters indicated in the legend correspond to those in Tables 1 and 2. 
Figure 5 shows χ^{2} as a function of M. For each mass we search for the minimum χ^{2} value by changing the remaining parameters (Y_{ini}, Z_{ini} and α). Models without overshooting were considered here. Figure 6 to 8 are similar but show Z_{ini}, Y_{ini}, and α respectively. In general the mass and the initial metallicity are well determined but the initial helium abundance and the mixing length parameter are not.
Correlation matrix for the grid of models that have no overshooting, corresponding to the first step in the minimization procedure.
Fig. 6 Minimum values of χ^{2} corresponding to every Z_{ini} for all the stars. Points with error bars correspond to χ^{2} values with 1σ uncertainties. The letters indicated in the legend correspond to those in Tables 1 and 2. 
Fig. 7 Minimum values of χ^{2} corresponding to every Y_{ini} for all the stars. Points with error bars correspond to χ^{2} values with 1σ uncertainties. The letters indicated in the legend correspond to those in Tables 1 and 2. 
Fig. 8 Minimum values of χ^{2} corresponding to every α for all the stars. Points with error bars correspond to χ^{2} values with 1σ uncertainties. The letters indicated in the legend correspond to those in Tables 1 and 2. 
Fig. 9 Values of χ^{2} against the overshooting parameter f_{ov}. The points with error bars correspond to χ^{2} values with 1σ uncertainties. The letters indicated in the legend correspond to those in Tables 1 and 2. The large circles correspond to models with M< 1 M_{⊙} and the large squares are for models with M> 1.39 M_{⊙}. 
As noted above, once the parameters M,Z_{ini},Y_{ini},α are determined by a χ^{2} minimization, we considered a grid of models with overshooting where ΔΠ_{1} is properly included in the minimization function. Figure 9 shows the corresponding χ^{2} values as a function of the overshooting parameter defined in Eq. (1). Although f_{ov} is not always well determined we note that in some cases an upper or lower limit can be inferred.
5.2. Overall results
Fig. 10 Values of log g obtained in this work plotted against observed values. The blue points correspond to the spectroscopic values and were used in the minimization procedure while the vertical coordinate of the red points was obtained by using the observed values of ν_{max} and T_{eff}, assuming a relationship . The observed values of ν_{max} were not used in the minimization procedure. The continuous black line corresponds to x = y. 
As might be expected, whereas the input and output values for T_{eff} are basically the same, values of log g are improved once the asteroseismic information is included. Specifically, the values of T_{eff} obtained in the minimization procedure are 75 ± 81 K higher than the input ones, which is consistent with the input and output errors of 108 K and 56 K respectively.
On the other hand, we obtain a mean error of 0.0018 dex for the output values of log g whereas the input spectroscopic errors are 0.1 dex on average. In Fig. 10 we compare both values of log g for all the target stars. The results are consistent and we have not found any bias, the dispersion between both data being 0.09 dex. Furthermore, in Fig. 10 we compare the values of log g obtained here with those derived from ν_{max} and T_{eff}, assuming a relationship of the form and calibrated with the Sun. Values of ν_{max} and T_{eff} and their errors are taken from Table 1. The resulting errors in log g are on average 0.01 dex, about five times higher than that from the minimization procedure. In any case, it seems clear from Fig. 10 that the agreement is much better in this instance, the mean difference being of 0.001 dex with a dispersion of 0.006 dex. We recall that we have not used the values of ν_{max} as input parameters, so the log g values obtained from the minimization procedure and those from ν_{max} are observationally linked only through T_{eff}. Hence, this comparison proves the consistency between both methods and indicates that the formal output error of 0.002 dex for log g obtained from the minimization procedure could be realistic.
The output values for the surface Z/X found here was lower than the input spectroscopic values by 0.002 ± 0.007. This can be compared with the average input and output errors of 0.0035 and 0.002 respectively. Hence, we do not find any bias but the dispersion is about twice as high as expected.
Fig. 11 Comparison between the masses obtained in the minimization procedure and those derived from simple scaling relations. The horizontal coordinate corresponds to the masses of the best models and the vertical coordinate gives the asteroseismic values derived from the scaling relationships. 
Fig. 12 Comparison between the radii obtained in the minimization procedure and those derived from simple scaling relationships. The horizontal coordinate corresponds to the radii of the best models while the vertical coordinate gives the asteroseismic values derived from the scaling relationships. 
The values of M and R determined here can be compared to those derived from the scaling relationships. Because for most of the stars we have used the spectroscopic values of T_{eff} from the APOKASC Catalogue (Pinsonneault et al. 2014) whereas Mosser et al. (2012) used photometric determinations that are on average lower by 145 K, the masses and radii derived here and those reported by Mosser et al. (2012) are also systematically shifted. For a consistent comparison we have computed values of M and R from the scaling relationships given by Mosser et al. (2013) but using values of T_{eff}, Δν and ν_{max} from Table 1. Results for individual stars are given in the last two columns of Table 2, and in Figs. 11 and 12 we compare the values of M and R found here with those derived from the scaling relationships. The masses derived from the minimization procedure are on average 0.006 times lower than those derived from the scaling relationships with a dispersion of 0.040. On the other hand the average relative errors are 0.023 and 0.030 for the minimization procedure and the scaling relationships respectively. Thus, both methods give consistent results, including their error determinations. For the radii we found relative differences between the two methods of 0.002 on average with a dispersion of 0.015. The formal relative errors are on average of 0.005 and 0.017 for the minimization procedure and the scaling relationships respectively. Hence, both methods again prove to be consistent.
It is also possible to compare the outputs from the forward modelling with those derived from the acoustic glitches, as in Fig. 13 where we show the acoustic depth of the second helium ionization zone, τ_{HeII}, versus log g. Blue points are for the forward modelling used in this work while the red points correspond to the values derived by Corsaro et al. (2015b) from a nonlinear fit of the second differences of the radial oscillations to a model introduced by Houdek & Gough (2007). As noted by Broomhall et al. (2014), τ_{HeII} depends mainly on the dynamical state of the star and to a much lesser extent on the He abundance; hence, the simple relationship between τ_{HeII} and log g shown in Fig. 13 for the results from the forward modelling. Values derived from the acoustic glitches are modelindependent but the fact that they do not reproduce such a relationship seems to indicate that errors in τ_{HeII} are higher than reported. Perhaps the problem arises because the number of measurements available is not much higher than the five free parameters of the theoretical equation considered. A test of different models and approaches to fit glitches could be considered in the future.
Fig. 13 Acoustic depth of the HeII ionization zone versus log g. Blue points are for the best models and red points the values derived by Corsaro et al. (2015b) from a fit to the acoustic glitches. 
Fig. 14 Overshooting parameter versus mass. Red and blue points follows the same criterion as in Fig. 4, the red points corresponding to stars where ΔΠ_{1} was obtained with a formal error lower than 0.1%. 
Finally in Fig. 14 we show the exponential overshooting parameter f_{ov} (see Eq. (1)) versus the mass, M. It seems that there is some correlation between both parameters, lower masses corresponding to higher f_{ov} values. To gain and idea of the relevance of such relation, in Fig. 9 we identify the three stars with masses M< 1 M_{⊙} (large circles) and the three stars with the highest masses (large squares). It seems clear that for all of our stars with M< 1 M_{⊙}, values f_{ov}> 0.025 give rise to a decrease in χ^{2} in such a way that lower f_{ov} should be rejected for these stars. In a similar way, Fig. 9 also indicates that for the more massive stars, the highest values of f_{ov} in the grid must be excluded. This seems more significant for KIC 009882316 (letter L in Fig. 9) which actually corresponds to the only massive star with a a good determination of ΔΠ_{1} (the rightmost red point in Fig. 14). Of course, this conclusion is modeldependent and the physical implications are hard to extract.
6. Conclusions
We have used a forward modelling technique for obtaining stellar parameters of nineteen lowmass, lowluminosity red giant stars for which highly accurate frequencies are available from the Kepler observations. In this first paper we have limited the work to the pmode frequencies and the period spacing of the ℓ = 1 modes. The relative frequency differences between our best models and the observations, once the surface effects are removed, are on average about 0.0002, twice as high as the observational errors, for modes with degrees ℓ = 0 and 3 and about 0.0004, four times higher than the observational errors, for the quadrupolar modes. The fact that these latter modes are worse fitted is probably caused by the regularly mixed nature of the ℓ = 2 which is hard to deal with properly, both observationally and theoretically.
The use of the pmode eigenfrequencies and the spectroscopic values of T_{eff} and surface metallicities have allowed us to determine the masses and radii of the stars with uncertainties of 2% and 0.5% respectively. These figures can be compared with the 3% and 2% uncertainties derived from the scaling relationships, that use only the general parameters Δν, ν_{max}, and T_{eff}. The consistency between both methods gives confidence in the individual values and the estimated errors of the stellar parameters reported in Table 2. However, it should be noted that the forward modelling is not free of systematic errors due to the input physics and the methodology used. We have not attempted to estimate such uncertainties, but as a guide, for mainsequence and subgiant stars, Chaplin et al. (2014) estimated errors of 3.7% and 1.3% for the mass and radius respectively due to these factors. Given the agreement we have found between the two methods, these figures seem rather high for our set of red giants. However this does not guarantee that other parameters such as the age or τ_{HeII} would be affected by the input physics and, hence, their errors can be higher than those given in Table 2.
On the other hand, other input parameters such as the initial helium abundance, Y_{ini}, and the mixing length parameter, α, could not be unambiguously determined. In principle one might think that owing to the highly accurate frequency measurements which allow us to clearly detect the presence of the glitch signatures caused by the second helium ionization zone (Corsaro et al. 2015b), the helium abundance could implicitly be constrained by the forward frequency comparison between models and observations. However, one should take into account that isolating Y_{ini} from other parameters is not a simple issue. In fact, in a related analysis for the solar case Pérez Hernández & ChristensenDalsgaard (1994) found that changes in Y_{ini} and the specific entropy of the adiabatic convection zone (parametrized by the mixing length parameter α) are highly correlated. It was only that for the Sun the known depth of the convection zone fixed the specific entropy that Y_{ini} could be determined with low uncertainties by using the acoustic signatures of the HeII zone.
We have found a correlation between overshooting and mass as shown in Fig. 14. However we are not claiming any general physical implication for it. First, our grid of models with overshooting was limited to a second step in the search for the best models, and was only introduced after the other parameters were fixed. We think this is reasonable for M and Z_{ini} but perhaps other parameters that were not well determined, such as Y_{ini} and α, should not be fixed in this second step. Second, the result is modeldependent and in particular only the prescription introduced by Herwig (2000) was considered. Also, a change in the opacity tables or the metallicity mixture can give rise to different results. Such improvements, of course, would increase by some orders of magnitudes the number of models to be considered.
To derive Eq. (9), we start from Eq. (31) in Jiang & ChristensenDalsgaard (2014) and note that δx can be expressed in our notation as πΔν_{j}/ Δν. Also ΔΠ here is Π_{j} and ω_{g} = 2π^{2}/ ΔΠ_{1}. Then if Δν is assumed to be known and a guess q_{0} for q is taken, this equation can be written as Π_{j} ≃ ΔΠ_{1}(1 + C ΔΠ_{1}y_{j})^{1}, where C is a known constant.
Acknowledgments
This paper made use of the IAC Supercomputing facility HTCondor (http://research.cs.wisc.edu/htcondor/), partly financed by the Ministry of Economy and Competitiveness with FEDER funds, code IACA133E2493. E.C. is funded by the European Community’s Seventh Framework Programme (FP7/20072013) under grant agreement n°312844 (SPACEINN). R.A.G. acknowledges received funding from the CNES GOLF and PLATO grants at CEA and the ANR (Agence Nationale de la Recherche, France) program IDEE (n ANR12BS050008) “Interaction Des Étoiles et des Exoplanètes”.
References
 Aerts, C., ChristensenDalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology (Springer) [Google Scholar]
 Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Cosmic Abundances as Records of Stellar Evolution and Nucleosynthesis, eds. T. G. Barnes, III, & F. N. Bash, ASP Conf. Ser., 336, 25 [Google Scholar]
 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Basu, S., & Antia, H. M. 2008, Phys. Rep., 457, 217 [NASA ADS] [CrossRef] [Google Scholar]
 Beck, P. G., Montalban, J., Kallinger, T., et al. 2012, Nature, 481, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Bedding, T., Kiss, L., & Kjeldsen, H. 2003, in IAU Joint Discussion, 12 [Google Scholar]
 Belkacem, K., Samadi, R., Mosser, B., Goupil, M.J., & Ludwig, H.G. 2013, in ASP Conf. Ser. 479, eds. H. Shibahashi, & A. E. LynasGray, 61 [Google Scholar]
 Borucki, W., Koch, D., Batalha, N., et al. 2009, in IAU Symp., 253, 289 [Google Scholar]
 Broomhall, A.M., Miglio, A., Montalbán, J., et al. 2014, MNRAS, 440, 1828 [NASA ADS] [CrossRef] [Google Scholar]
 Chaplin, W. J., Basu, S., Huber, D., et al. 2014, ApJS, 210, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Chaplin, W. J., Lund, M. N., Handberg, R., et al. 2015, PASP, 127, 1038 [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J. 2008, Ap&SS, 316, 113 [NASA ADS] [CrossRef] [Google Scholar]
 ChristensenDalsgaard, J. 2012, in Progress in Solar/Stellar Physics with Helio and Asteroseismology, eds. H. Shibahashi, M. Takata, & A. E. LynasGray, ASP Conf. Ser., 462, 503 [Google Scholar]
 Corsaro, E., & De Ridder, J. 2014, A&A, 571, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Corsaro, E., Stello, D., Huber, D., et al. 2012, ApJ, 757, 190 [NASA ADS] [CrossRef] [Google Scholar]
 Corsaro, E., De Ridder, J., & García, R. A. 2015a, A&A, 579, A83 [NASA ADS] [CrossRef] [EDP Sciences] [PubMed] [Google Scholar]
 Corsaro, E., De Ridder, J., & García, R. A. 2015b, A&A, 578, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., García, R. A., Chaplin, W. J., et al. 2012, ApJ, 756, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Deheuvels, S., Doğan, G., Goupil, M. J., et al. 2014, A&A, 564, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Deheuvels, S., Ballot, J., Beck, P. G., et al. 2015, A&A, 580, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Di Mauro, M. P., Ventura, R., Cardini, D., et al. 2016, ApJ, 817, 65 [NASA ADS] [CrossRef] [Google Scholar]
 Freytag, B., Ludwig, H.G., & Steffen, M. 1996, A&A, 313, 497 [Google Scholar]
 García, R. A., & Stello, D. 2015, in Extraterrestrial Seismology, eds. V. C. H. Tong & R. A. García (Cambridge University Press) [Google Scholar]
 Gough, D. O. 1990, in Progress of Seismology of the Sun and Stars, eds. Y. Osaki & H. Shibahashi (Berlin: Springer Verlag), Lect. Notes Phys., 367, 283 [Google Scholar]
 Grevesse, N., & Sauval, A. J. 1998, Space Sci. Rev., 85, 161 [NASA ADS] [CrossRef] [Google Scholar]
 Guillot, T., & Havel, M. 2011, A&A, 527, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Herwig, F. 2000, A&A, 360, 952 [NASA ADS] [Google Scholar]
 Houdek, G., & Gough, D. O. 2007, MNRAS, 375, 861 [NASA ADS] [CrossRef] [Google Scholar]
 Howell, S. B., Sobeck, C., Haas, M., et al. 2014, PASP, 126, 398 [NASA ADS] [CrossRef] [Google Scholar]
 Huber, D., Chaplin, W. J., ChristensenDalsgaard, J., et al. 2013, ApJ, 767, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
 Jiang, C., & ChristensenDalsgaard, J. 2014, MNRAS, 444, 3622 [NASA ADS] [CrossRef] [Google Scholar]
 Jiménez, A., García, R. A., Pérez Hernández, F., & Mathur, S. 2015, A&A, 583, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mathur, S., Metcalfe, T. S., Woitaszek, M., et al. 2012, ApJ, 749, 152 [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]
 Moravveji, E., Aerts, C., Pápics, P. I., Triana, S. A., & Vandoren, B. 2015, A&A, 580, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Goupil, M. J., Belkacem, K., et al. 2012, VizieR Online Data Catalog: J/A+A/540/A143 [Google Scholar]
 Mosser, B., Michel, E., Belkacem, K., et al. 2013, A&A, 550, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mosser, B., Benomar, O., Belkacem, K., et al. 2014, A&A, 572, L5 [NASA ADS] [CrossRef] [EDP Sciences] [PubMed] [Google Scholar]
 Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3 [Google Scholar]
 Pérez Hernández, F., & ChristensenDalsgaard, J. 1994, MNRAS, 269 [Google Scholar]
 Pinsonneault, M. H., An, D., MolendaŻakowicz, J., et al. 2012, ApJS, 199, 30 [NASA ADS] [CrossRef] [Google Scholar]
 Pinsonneault, M. H., Elsworth, Y., Epstein, C., et al. 2014, ApJS, 215, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Stello, D., Bruntt, H., Preston, H., & Buzasi, D. 2008, ApJ, 674, L53 [NASA ADS] [CrossRef] [Google Scholar]
 Stello, D., Huber, D., Bedding, T. R., et al. 2013, ApJ, 765, L41 [NASA ADS] [CrossRef] [Google Scholar]
 Stello, D., Huber, D., Sharma, S., et al. 2015, ApJ, 809, L3 [NASA ADS] [CrossRef] [Google Scholar]
 Vorontsov, S. V. 1988, in Advances in Helio and Asteroseismology, eds. J. ChristensenDalsgaard & S. Frandsen, IAU Symp., 123, 151 [Google Scholar]
All Tables
Correlation matrix for the grid of models that have no overshooting, corresponding to the first step in the minimization procedure.
All Figures
Fig. 1 Values of χ^{2} as function of age for an evolutionary sequence with stellar parameters that match those of KIC 003744043 within errors (letter A in Tables 1 and 2). The black open circles are computed values for models and frequencies, whereas the continuous black line corresponds to parameters and frequencies interpolated to intermediate ages. The filled black circle is the χ^{2} value derived from the bestfit model explicitly computed at the end of the procedure to exactly match the interpolated age. Also shown are values of (blue), (red), (green), and (magenta). Here the solid lines and the full and open circles have the same meaning as before. 

In the text 
Fig. 2 Frequency differences between the observations and the best models, [ ω(observations)−ω(best model) ] /ω (red points), and associated dimensionless energies in arbitrary units (blue squares) for radial oscillations. For a better comparison between modes in different stars, the horizontal axis gives the frequencies normalized to the frequency of maximum power. All nineteen stars are included. The black line is a secondorder polynomial fit. 

In the text 
Fig. 3 Frequency differences between the observations and the best models after subtracting the surface effects. Red points are for ℓ = 0, blue points for ℓ = 2 and green points for ℓ = 3. All nineteen stars are included. 

In the text 
Fig. 4 Period spacing obtained in this work compared to that in Corsaro et al. (2015a). Red points are values with errors lower than 0.1%. The continuous line is x = y and the dashed one y = x−2.5 s. 

In the text 
Fig. 5 Minimum values of χ^{2} as a function of M – changing Y_{ini}, Z_{ini} and α – for all the stars. Points with error bars correspond to χ^{2} values with 1σ uncertainties. The letters indicated in the legend correspond to those in Tables 1 and 2. 

In the text 
Fig. 6 Minimum values of χ^{2} corresponding to every Z_{ini} for all the stars. Points with error bars correspond to χ^{2} values with 1σ uncertainties. The letters indicated in the legend correspond to those in Tables 1 and 2. 

In the text 
Fig. 7 Minimum values of χ^{2} corresponding to every Y_{ini} for all the stars. Points with error bars correspond to χ^{2} values with 1σ uncertainties. The letters indicated in the legend correspond to those in Tables 1 and 2. 

In the text 
Fig. 8 Minimum values of χ^{2} corresponding to every α for all the stars. Points with error bars correspond to χ^{2} values with 1σ uncertainties. The letters indicated in the legend correspond to those in Tables 1 and 2. 

In the text 
Fig. 9 Values of χ^{2} against the overshooting parameter f_{ov}. The points with error bars correspond to χ^{2} values with 1σ uncertainties. The letters indicated in the legend correspond to those in Tables 1 and 2. The large circles correspond to models with M< 1 M_{⊙} and the large squares are for models with M> 1.39 M_{⊙}. 

In the text 
Fig. 10 Values of log g obtained in this work plotted against observed values. The blue points correspond to the spectroscopic values and were used in the minimization procedure while the vertical coordinate of the red points was obtained by using the observed values of ν_{max} and T_{eff}, assuming a relationship . The observed values of ν_{max} were not used in the minimization procedure. The continuous black line corresponds to x = y. 

In the text 
Fig. 11 Comparison between the masses obtained in the minimization procedure and those derived from simple scaling relations. The horizontal coordinate corresponds to the masses of the best models and the vertical coordinate gives the asteroseismic values derived from the scaling relationships. 

In the text 
Fig. 12 Comparison between the radii obtained in the minimization procedure and those derived from simple scaling relationships. The horizontal coordinate corresponds to the radii of the best models while the vertical coordinate gives the asteroseismic values derived from the scaling relationships. 

In the text 
Fig. 13 Acoustic depth of the HeII ionization zone versus log g. Blue points are for the best models and red points the values derived by Corsaro et al. (2015b) from a fit to the acoustic glitches. 

In the text 
Fig. 14 Overshooting parameter versus mass. Red and blue points follows the same criterion as in Fig. 4, the red points corresponding to stars where ΔΠ_{1} was obtained with a formal error lower than 0.1%. 

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.