We performed a theoretical analysis on the relevance of the small frequency separation $\delta \nu$ in determining stellar ages, masses, and radii. We adopted the SCEPtER pipeline for low-mass stars, [0.7, 1.05] Msun. Synthetic stars were generated and reconstructed assuming different relative precision in $\delta \nu$ (5% and 2%). The quantification of the systematic errors arising from a mismatch between synthetic stars and the recovery grid was also performed. The biases were negligible. The statistical error in age estimates was strongly dependent on the stellar evolutionary phase. It is at its maximum at ZAMS and decreases to about 11% and 6% ($\delta \nu$ at 5% and 2% level) when stars reach the 30% of their evolutionary MS lifetime, then vanishes in the last 20% of the MS. For mass and radius estimates we detected nearly identical behavior for an observational uncertainty of 5%. No benefit was detected for mass and radius determinations from a reduction $\delta \nu$ precision to 2%. The impact on the age of the initial helium abundance resulted in negligible results. On the other hand greater bias (2% and 1%) in mass and radius estimates were detected whenever $\delta \nu$ is in the observational pool. This result, together with the presence of further unexplored uncertainty sources, suggest that precision in the derived stellar quantities below these thresholds may possibly be overoptimistic. [...] Finally, we compared the estimates by the SCEPtER pipeline for 13 Kepler asteroseismic LEGACY sample stars with those given by six different pipelines. Overall, on the LEGACY sample data, we obtained a multi-pipeline precision of about 4.4%, 1.7%, and 11% on the estimated masses, radii, and ages, respectively.

We performed a theoretical investigation on the mixing-length parameter recovery from an eclipsing double-lined binary system. We focused on a syntetic system composed by a primary of mass M = 0.95 Msun and a secondary of M = 0.85 Msun. Monte Carlo simulations were conducted at three metallicities, and three evolutionary stages of the primary. For each configuration artificial data were sampled assuming an increasing difference between the mixing-length of the two stars. The mixing length values were reconstructed using three alternative set-ups. A first method, which assumes full independence between the two stars, showed a great difficulty to constrain the mixing-length values: the recovered values were nearly unconstrained with a standard deviation of 0.40. The second technique imposes the constraint of common age and initial chemical composition for the two stars in the fit. We found that $\alpha_{ml,1}$ values match the ones recovered under the previous configuration, but $\alpha_{ml,2}$ values are peaked around unbiased estimates. This occurs because the primary star provides a much more tight age constraint in the joint fit than the secondary. Within this second scenario we also explored, for systems sharing a common $\alpha_{ml}$, the difference in the mixing-length values of the two stars only due to random fluctuations owing to the observational errors. The posterior distribution of these differences was peaked around zero, with a large standard deviation of 0.3 (15\% of the solar-scaled value). The third technique also imposes the constraint of a common mixing-length value for the two stars, and served as a test for identification of wrong fitting assumptions. In this case the common mixing-length is mainly dictated by the value of $\alpha_{ml,2}$. [...] For $\Delta \alpha_{ml} > 0.4$ less than half of the systems can be recovered and only 20% at $\Delta \alpha_{ml} = 1.0$.

We analysed the theoretical foundation of the mixing-length calibration by means of classical and global asteroseismic observables of field stars. We discussed the soundness of inferring a metallicity dependence of the mixing-length parameter. We followed a theoretical approach based on mock datasets of artificial stars sampled from a grid of stellar models with a fixed $\alpha_{ml}$. [...] We verified that the $\alpha_{ml}$ estimates are affected by a huge spread, even in the case of perfect agreement between the mock data and the recovery grid. When the solar heavy-element mixture used to compute the models is different from that of the artificial stars, we also found a metallicity relationship in which $\alpha_{ml}$ increases by 0.4 for an increase of 1 dex in [Fe/H]. The origin of this trend was investigated considering the differences in the initial helium abundance Y -- [Fe/H] -- initial metallicity Z relation assumed in the models and data. A discrepancy between the adopted helium-to-metal enrichment ratio caused the appearance of relevant spurious trends in the estimated $\alpha_{ml}$. [...] A similar effect was caused by an offset in the [Fe/H] to Z conversion. An overestimation of [Fe/H] by 0.1 dex in the recovery grid forced an increasing trend of $\alpha_{ml}$ versus [Fe/H] of 0.2 per dex. We also explored the impact of some discrepancies between the adopted input physics in the recovery grid and mock data. We observed a trend with [Fe/H] of 0.3 per dex when the effect of the microscopic diffusion is neglected in the recovery grid [...]. Therefore, any attempt to calibrate the mixing-length parameter by means of classical and asteroseimic observables of field stars seems to be statistically poorly reliable. As such, any claim about the dependence of the mixing-length on the metallicity for field stars should be considered cautiously and critically.

We theoretically quantify the impact of observational errors on the recovered age for stars in MS and RGB phases. We adopted Teff, [Fe/H], Delta_nu and nu_max as observational constraints. Artificial stars were sampled from a reference isochrone and subjected to random perturbation to simulate observational errors. Their ages were then recovered by means of a MCMC approach. The differences between the recovered and true ages were modelled against the errors in the observables by means of linear models and projection pursuit regression models. From linear models we find that no age error source dominates. For MS the most important error source is Teff. An offset of 75 K accounts for an underestimation of the stellar age from 0.4 to 0.6 Gyr for initial and terminal MS. An error of 2.5% in nu_max accounted for about -0.3 Gyr. A 0.1 dex error in [Fe/H] resulted particularly important only at the end of the MS, with an age error of -0.4 Gyr. For the RGB phase the dominant source of uncertainty is nu_max, causing an underestimation of about 0.6 Gyr; the offset in the effective temperature and Delta_nu caused respectively an underestimation and overestimation of 0.3 Gyr. We find that the inference from the linear model is a good proxy for that from projection pursuit regression models, so that linear models can be safely used thanks to its broader generalizability. Finally, we explored the impact on age estimates of adding the luminosity to the other observational constraints. We assumed - for computational reasons - a 2.5% error in luminosity, much lower than the average error in the Gaia DR2 catalogue. Even in this optimistic case, the luminosity does not increase precision of age estimates. Moreover, the luminosity resulted as a major contributor to the variability in the estimated ages, accounting for an error of about -0.3 Gyr in the explored evolutionary phases.

We performed a theoretical investigation focused on the age estimate of RGB stars in OCs based on mixed classical surface and asteroseismic parameters. We evaluated the performances of three widely adopted methods (pure geometrical fit, maximum likelihood approach, and a single stars fit) in recovering stellar parameters. Artificial OCs were generated by means of a Monte Carlo procedure for two different ages (7.5 and 9.0 Gyr) and two different choices of the number of stars in the RGB evolutionary phase (35 and 80). The geometrical approach overestimated the age by about 0.3 and 0.2 Gyr for true ages of 7.5 and 9.0 Gyr, respectively. The ML approach provided similar biases (0.1 and 0.2 Gyr) but with a variance reduced by a factor of between two and four with respect to geometrical fit. The independent fit of single stars showed a very large variance. The most important difference between geometrical and ML approaches was the robustness against observational errors. For the geometrical method, we found that estimations starting from the same sample but with different Gaussian perturbations on the observables had about 0.3 Gyr random variability from one Monte Carlo run to another (45% of the intrinsic variability due to observational errors). On the other hand, for the ML method, this value was about 65% and 90% of the simulations failed to include the true parameter values in their estimated 1 sigma credible interval. Finally, we compared the performance of the three fitting methods for single RGB-star age estimation. The variability owing to the choice of the fitting method was minor, being about 15% of the variability caused by observational uncertainties.

We evaluated the bias and variability on the fitted age and convective core overshooting parameter for evolved binary stars accounting for observational and internal uncertainties. We considered a binary system composed of a 2.50 $M_{\sun}$ primary star coupled with a 2.38 $M_{\sun}$ secondary in three evolutionary stages (primary at the end of the central helium burning; at the bottom of the RGB; and in the helium core burning). The simulations have been carried out for two values of accuracy on the mass determination (1% and 0.1%). We found that the fitted age and overshooting efficiency are always biased towards low values. The underestimation is relevant for a primary in the central helium burning stage, reaching -8.5% in age and -0.04 (-25% relative error) in the overshooting parameter $\beta$. In the other scenarios, an undervaluation of the age by about 4% occurs. A large variability in the fitted values between simulations was found: for an individual system calibration, the value of $\beta$ can vary from 0.0 to 0.26. For an error of 0.1% on the masses the global variability is suppressed by a factor of two. We accounted for a systematic offset in the effective temperature of the stars by $\pm 150$ K. For a mass error of 1% $\beta$ is largely biased towards the edges of the explored range, while for the lower mass uncertainty it is basically unconstrained from 0.0 to 0.2. We evaluated the possibility of individually recovering the $\beta$ value for both stars. We found that this is impossible for a primary near to central hydrogen exhaustion, while in the other cases the fitted $\beta$ are consistent, but always biased. Finally, the possibility to distinguish between models computed with mild overshooting from models with no overshooting resulted in a reassuring power of 80%. However, the scenario with a primary in the central helium burning showed a power lower than 5%.

We aim to perform a theoretical evaluation of the impact of the mass loss indetermination on asteroseismic grid based estimates of masses, radii, and ages of stars in the red giant branch phase (RGB). We adopted the SCEPtER pipeline on a grid spanning the mass range [0.8; 1.8] Msun. As observational constraints, we adopted the star effective temperatures, the metallicity [Fe/H], the average large frequency spacing $\Delta \nu,$ and the frequency of maximum oscillation power $\nu_{\rm max}$. The mass loss was modelled following a Reimers parametrization with the two different efficiencies $\eta = 0.4$ and $\eta = 0.8$. In the RGB phase, the average error owing only to observational uncertainty on mass and age estimates is about 8% and 30% respectively. The bias in mass and age estimates caused by the adoption of a wrong mass loss parameter in the recovery is minor for the vast majority of the RGB evolution. The biases get larger only after the RGB bump. In the last 2.5% of the RGB lifetime the error on the mass determination reaches 6.5% becoming larger than the random error component in this evolutionary phase. The error on the age estimate amounts to 9%, that is, equal to the random error uncertainty. These results are independent of the stellar metallicity [Fe/H] in the explored range. Asteroseismic-based estimates of stellar mass, radius, and age in the RGB phase can be considered mass loss independent within the range ($\eta \in [0.0, 0.8]$) as long as the target is in an evolutionary phase preceding the RGB bump.

We attempt to constrain the initial helium abundance, the age and the efficiency of the convective core overshooting of the binary system TZ Fornacis. Our main aim is in pointing out the biases in the results due to not accounting for some sources of uncertainty. We adopt the SCEPtER pipeline, relying on stellar models computed with two stellar evolutionary codes (FRANEC and MESA). We found multiple independent groups of solutions. The best one suggests a system of age 1.10 $\pm$ 0.07 Gyr (primary star in the central helium burning stage, secondary in the sub-giant branch), with a helium-to-metal enrichment ratio of $\Delta Y/\Delta Z = 1$ and core overshooting parameter $\beta = 0.15 \pm 0.01$ (FRANEC) and $f_{\rm ov}= 0.013 \pm 0.001 $ (MESA). The second class of solutions, characterised by a worse goodness-of-fit, still suggest a primary star in the central helium-burning stage but a secondary in the overall contraction phase, at the end of the main sequence. The FRANEC grid provides an age of $1.16_{-0.02}^{+0.03}$ Gyr and $\beta = 0.25_{-0.01}^{+0.005}$, while the MESA grid gives $1.23 \pm 0.03$ Gyr and $f_{\rm ov} = 0.025 \pm 0.003$. We also show that very precise mass determinations with uncertainty of a few thousandths of solar mass are required to obtain reliable determinations of stellar parameters, as mass errors larger than approximately 1% lead to estimates that are not only less precise but also biased. Moreover, we show that a fit obtained with a grid of models computed at a fixed $\Delta Y/\Delta Z$ can provide biased age and overshooting estimates. The possibility of independent overshooting efficiencies for the two stars of the system is also explored. This work confirms that to constrain the core overshooting parameter by means of binary systems is a very difficult task that requires an observational precision still rarely achieved and robust statistical methods.

We develop a statistical test on the expected difference in age estimates of two coeval stars in detached double-lined eclipsing binary systems that are only caused by observational uncertainties. We focus on stars in the mass range [0.8; 1.6] Msun, and on stars in the main-sequence phase. The ages were obtained by means of the maximum-likelihood SCEPtER technique. The observational constraints used in the recovery procedure are stellar mass, radius, effective temperature, and metallicity [Fe/H]. We defined the statistic W computed as the ratio of the absolute difference of estimated ages for the two stars over the age of the older one. We determined the critical values of this statistics above which coevality can be rejected. The median expected difference in the reconstructed age between the coeval stars of a binary system -- caused alone by the observational uncertainties -- shows a strong dependence on the evolutionary stage. This ranges from about 20% for an evolved primary star to about 75% for a near ZAMS primary. The median difference also shows an increase with the mass of the primary star from 20% for 0.8 Msun stars to about 50% for 1.6 Msun stars. The reliability of these results was checked by repeating the process with a grid of stellar models computed by a different evolutionary code. We show that the W test is much more sensible to age differences in the binary system components than the alternative approach of comparing the confidence interval of the age of the two stars. We also found that the distribution of W is, for almost all the examined cases, well approximated by beta distributions. The proposed method improves upon the techniques that are commonly adopted for judging the coevality of an observed system. It also provides a result founded on reliable statistics that simultaneously accounts for all the observational uncertainties.

In a robust statistical way, we quantify the uncertainty that affects the calibration of the overshooting efficiency parameter $\beta$ that is owing to the uncertainty on the observational data in double-lined eclipsing binary systems. We also quantify the bias that is caused by the lack of constraints on the initial helium content and on the efficiencies of the superadiabatic convection and microscopic diffusion. We adopted a modified grid-based SCEPtER pipeline using as observational constraints the effective temperatures, [Fe/H], masses, and radii of the two stars. In a reference scenario of mild overshooting $\beta = 0.2$ for the synthetic data, we found both large statistical uncertainties and biases on the estimated $\beta$. For the first 80% of the MS evolution, $\beta$ is biased and practically unconstrained in the whole explored range [0.0; 0.4]. In the last 5% of the MS the bias vanishes and the $1 \sigma$ error is about 0.05. For synthetic data computed with $\beta = 0.0$, the estimated $\beta$ is biased by about 0.12 in the first 80% of the MS evolution, and by 0.05 afterwards. Assuming an uncertainty of $\pm 1$ in the helium-to-metal enrichment ratio $\Delta Y/\Delta Z$, we found that in the terminal part of the MS evolution the error on the estimated $\beta$ values ranges from -0.05 to +0.10, while $\beta$ is basically unconstrained throughout the explored range at earlier evolutionary stages. A uniform variation of $\pm 0.24$ in the mixing-length parameter around the solar-calibrated value causes in last 5% of the MS an uncertainty from -0.09 to +0.15. A complete neglect of diffusion in the stellar evolution computations produces a $1 \sigma$ uncertainty of $\pm 0.08$ in the last 5% of the MS, while $\beta$ is practically unconstrained in the first 80% of the MS. Overall, the calibration appears poorly reliable.

We investigate the performance of grid-based techniques in estimating the age of stars in detached eclipsing binary systems. We evaluate the precision of the estimates due to the uncertainty in the observational constraints, and the systematic bias caused by the uncertainty in convective core overshooting, element diffusion, mixing-length value, and initial helium content. We adopted the SCEPtER grid, which includes stars with mass in the range [0.8; 1.6] Msun and evolutionary stages from the ZAMS to the central hydrogen depletion. Age estimates have been obtained by a generalisation of the technique described in our previous work. We showed that the typical 1 sigma random error in age estimates - due to the uncertainty on the observational constraints - is about +- 7%, which is nearly independent of the masses of the two stars. However, such an error strongly depends on the evolutionary phase and becomes larger and asymmetric for stars near the ZAMS where it ranges from about +90% to -25%. The systematic bias due to the including mild and strong convective core overshooting is about 50% and 120% of the error due to observational uncertainties. A variation of +- 1 in the helium-to-metal enrichment ratio accounts for about +- 150% of the random error. The neglect of microscopic diffusion accounts for a bias of about 60% of the random error. We also introduced a statistical test of the expected difference in the recovered age of two coeval stars in a binary system. We find that random fluctuations within the current observational uncertainties can lead genuine coeval binary components to appear to be non-coeval with a difference in age as high as 60%.

Recent results by Bensby and collaborators on the ages of microlensed stars in the Galactic bulge have challenged the picture of an exclusively old stellar population. However, these age estimates have not been independently confirmed. In this paper we verify these results by means of a grid-based method and quantify the systematic biases that might be induced by some assumptions adopted to compute stellar models. We explore the impact of increasing the initial helium abundance, neglecting the element microscopic diffusion, and changing the mixing-length calibration in theoretical stellar track computations. We adopt the SCEPtER pipeline with a novel stellar model grid for metallicities [Fe/H] from -2.00 to 0.55 dex, and masses in the range [0.60; 1.60] Msun from the ZAMS to the helium flash at the red giant branch tip. We show for the considered evolutionary phases that our technique provides unbiased age estimates. Our age results are in good agreement with Bensby and collaborators findings and show 16 stars younger than 5 Gyr and 28 younger than 9 Gyr over a sample of 58. The effect of a helium enhancement as large as Delta Y/Delta Z = 5 is quite modest, resulting in a mean age increase of metal rich stars of 0.6 Gyr. Even simultaneously adopting a high helium content and the upper values of age estimates, there is evidence of 4 stars younger than 5 Gyr and 15 younger than 9 Gyr. For stars younger than 5 Gyr, the use of stellar models computed by neglecting microscopic diffusion or by assuming a super-solar mixing-length value leads to a mean increase in the age estimates of about 0.4 Gyr and 0.5 Gyr respectively. Even considering the upper values for the age estimates, there are four stars estimated younger than 5 Gyr is in both cases. Thus, the assessment of a sizeable fraction of young stars among the microlensed sample in the Galactic bulge appears robust.

We study the impact on stellar age determination by means of grid-based techniques adopting asteroseismic constraints of the uncertainty in the radiative opacity, in the initial helium abundance, in the mixing-length value, in the convective core overshooting, and in the microscopic diffusion efficiency adopted in stellar model computations. We extended our SCEPtER grid (Valle et al. 2014) to include stars with mass in the range [0.8; 1.6] Msun and evolutionary stages from the ZAMS to the central hydrogen depletion. The current typical uncertainty in the observations accounts for 1 sigma statistical relative error in age determination which in mean ranges from about -35% to +42%, depending on the mass. However, due to the strong dependence on the evolutionary phase, the age relative error can be higher than 120% for stars near the ZAMS, while it is typically of the order of 20% or lower in the advanced main-sequence phase. The systematic bias on age determination due to a variation of $\pm$ 1 in the helium-to-metal enrichment ratio Delta Y/Delta Z is about one-forth of the statistical error in the first 30% of the evolution while it is negligible for more evolved stages. The maximum bias due to the presence of the convective core overshooting is of -7% and -13% for mild and strong overshooting scenarios. For all the examined models the impact of a variation of $\pm$ 5 in the radiative opacity was found to be negligible. The most important source of bias are the uncertainty in the mixing-length value alpha_ml and the neglect of microscopic diffusion. Each of these effects accounts for a bias which is nearly equal to the random error uncertainty. Comparison of the results of our technique with other grid techniques on a set of common stars showed a general agreement. However, the adoption of a different grid can account for a variation in the mean estimated age up to 1 Gyr.

We study the temporal evolution of the habitable zone (HZ) of low-mass stars - only due to stellar evolution - and evaluate the related uncertainties. These uncertainties are then compared with those due to the adoption of different climate models. We computed stellar evolutionary tracks from the pre-main sequence phase to the helium flash at the red-giant branch tip for stars with masses in the range [0.70 - 1.10] Msun, metallicity Z in the range [0.005 - 0.04], and various initial helium contents. We evaluated several characteristics of the HZ, such as the distance from the host star at which the habitability is longest, the duration of this habitability, the width of the zone for which the habitability lasts one half of the maximum, and the boundaries of the continuously habitable zone (CHZ) for which the habitability lasts at least 4 Gyr. We developed analytical models, accurate to the percent level or lower, which allowed to obtain these characteristics in dependence on the mass and the chemical composition of the host star. The metallicity of the host star plays a relevant role in determining the HZ. The initial helium content accounts for a variation of the CHZ boundaries as large as 30% and 10% in the inner and outer border. The computed analytical models allow the first systematic study of the variability of the CHZ boundaries that is caused by the uncertainty in the estimated values of mass and metallicity of the host star. An uncertainty range of about 30% in the inner boundary and 15% in the outer one were found. We also verified that these uncertainties are larger than that due to relying on recently revised climatic models, which leads to a CHZ boundaries shift within 5% with respect to those of our reference scenario. We made an on-line tool available that provides both HZ characteristics and interpolated stellar tracks.

Some aspects of the systematic and statistical errors affecting grid-based estimation of stellar masses and radii have still not been investigated well. We study the impact on mass and radius determination of the uncertainty in the input physics, in the mixing-length value, in the initial helium abundance, and in the microscopic diffusion efficiency adopted in stellar model computations. We consider stars with mass in the range [0.8 - 1.1] Msun and evolutionary stages from the zero-age main sequence to the central hydrogen exhaustion. Stellar parameters were recovered by a maximum-likelihood technique, comparing the observations constraints to a grid of stellar models. Synthetic grids with perturbed input were adopted to estimate the systematic errors arising from the current uncertainty in model computations. We found that the statistical error components, owing to the current typical uncertainty in the observations, are nearly constant in all cases at about 4.5% and 2.2% on mass and radius determination, respectively. The systematic bias on mass and radius determination due to a variation of $\pm$ 1 in Delta Y/Delta Z is $\pm$ 2.3% and $\pm$ 1.1%; the one due to a change of $\pm$ 0.24 in the value of the mixing-length is $\pm$ 2.1% and $\pm$ 1.0%; the one due to a variation of $\pm$ 5% in the radiative opacity is $\mp$ 1.0% and $\mp$ 0.45%. An important bias source is to neglect microscopic diffusion, which accounts for errors of about 3.7% and 1.5% on mass and radius. The cumulative effects of the considered uncertainty sources can produce biased estimates of stellar characteristics. Comparison of the results of our technique with other grid techniques shows that the systematic biases induced by the differences in the estimation grids are generally greater than the statistical errors involved.

We extend our work on the effects of the uncertainties on the main input physics for the evolution of low-mass stars. We analyse the dependence of the cumulative physical uncertainty affecting stellar tracks on the chemical composition. We calculated more than 6000 stellar tracks and isochrones, with metallicity ranging from Z = 0.0001 to 0.02, by changing the following physical inputs within their current range of uncertainty: 1H(p,nu e+)2H, 14N(p,gamma)15O and triple-alpha reaction rates, radiative and conductive opacities, neutrino energy losses, and microscopic diffusion velocities. The analysis was performed using a latin hypercube sampling design. We examine in a statistical way the dependence on the variation of the physical inputs of the turn-off (TO) luminosity, the central hydrogen exhaustion time (t_H), the luminosity and the helium core mass at the red-giant branch (RGB) tip, and the zero age horizontal branch (ZAHB) luminosity in the RR Lyrae region. For the stellar tracks, an increase from Z = 0.0001 to Z = 0.02 produces a cumulative physical uncertainty in TO luminosity from 0.028 dex to 0.017 dex, while the global uncertainty on t_H increases from 0.42 Gyr to 1.08 Gyr. For the RGB tip, the cumulative uncertainty on the luminosity is almost constant at 0.03 dex, whereas the one the helium core mass decreases from 0.0055 M_sun to 0.0035 M_sun. The dependence of the ZAHB luminosity error is not monotonic with Z, and it varies from a minimum of 0.036 dex at Z = 0.0005 to a maximum of 0.047 dex at Z = 0.0001. Regarding stellar isochrones of 12 Gyr, the cumulative physical uncertainty on the predicted TO luminosity and mass increases respectively from 0.012 dex to 0.014 dex and from 0.0136 M_sun to 0.0186 M_sun. Consequently, for ages typical of galactic globular clusters, the uncertainty on the age inferred from the TO luminosity increases from 325 Myr to 415 Myr.

We present the R package stellaR, which is designed to access and manipulate publicly available stellar evolutionary tracks and isochrones from the Pisa low-mass database. The procedures of the extraction of important stages in the evolution of a star from the database, of the isochrones construction from stellar tracks and of the interpolation among tracks are discussed and demonstrated.

Using our updated stellar evolutionary code, we quantitatively evaluate the effects of the uncertainties in the main physical inputs on the evolutionary characteristics of low mass stars from the main sequence to the zero age horizontal branch (ZAHB). We calculated more than 3000 stellar tracks and isochrones, with updated solar mixture, by changing the following physical inputs within their current range of uncertainty: 1H(p,nu e+)2H, 14N(p,gamma)15O, and triple-alpha reaction rates, radiative and conductive opacities, neutrino energy losses, and microscopic diffusion velocities. We performed a systematic variation on a fixed grid, in a way to obtain a full crossing of the perturbed input values. The effect of the variations of the chosen physical inputs on relevant stellar evolutionary features, such as the turn-off luminosity, the central hydrogen exhaustion time, the red-giant branch (RGB) tip luminosity, the helium core mass, and the ZAHB luminosity in the RR Lyrae region are statistically analyzed. For a 0.9 Msun model, the cumulative uncertainty on the turn-off, the RGB tip, and the ZAHB luminosities accounts for $\pm$ 0.02 dex, $\pm$ 0.03 dex, and $\pm$ 0.045 dex respectively, while the central hydrogen exhaustion time varies of about $\pm$ 0.7 Gyr. The most relevant effect is due to the radiative opacities uncertainty; for the later evolutionary stages the second most important effect is due to the triple-alpha reaction rate uncertainty. For an isochrone of 12 Gyr, we find that the isochrone turn-off log luminosity varies of $\pm$ 0.013 dex, the mass at the isochrone turn-off varies of $\pm$ 0.015 Msun, and the difference between ZAHB and turn-off log-luminosity varies of $\pm$ 0.05 dex. The effect of the physical uncertainty affecting the age inferred from turn-off luminosity and from the vertical method are of $\pm$ 0.375 Gyr and $\pm$ 1.25 Gyr respectively.

The last decade showed an impressive observational effort from the photometric and spectroscopic point of view for ancient stellar clusters in our Galaxy and beyond. The theoretical interpretation of these new observational results requires updated evolutionary models and isochrones spanning a wide range of chemical composition. With this aim we built the new "Pisa Stellar Evolution Database" of stellar models and isochrones by adopting a well-tested evolutionary code (FRANEC) implemented with updated physical and chemical inputs. In particular, our code adopts realistic atmosphere models and an updated equation of state, nuclear reaction rates and opacities calculated with recent solar elements mixture. A total of 32646 models have been computed in the range of initial masses 0.30 - 1.10 Msun for a grid of 216 chemical compositions with the fractional metal abundance in mass, Z, ranging from 0.0001 to 0.01, and the original helium content, Y, from 0.25 to 0.42. Models were computed for both solar-scaled and alpha-enhanced abundances with different external convection efficiencies. Correspondingly, 9720 isochrones were computed in the age range 8 - 15 Gyr, in time steps of 0.5 Gyr. The whole database is available to the scientific community on the web. Models and isochrones were compared with recent calculations available in the literature and with the color-magnitude diagram of selected Galactic globular clusters. The dependence of relevant evolutionary quantities on the chemical composition and convection efficiency were analyzed in a quantitative statistical way and analytical formulations were made available for reader's convenience.