Search

Context. The cosmological lithium problem, that is, the discrepancy between the lithium abundance predicted by the Big Bang nucleosynthesis and the one observed for the stars of the “Spite plateau”, is one of the long standing problems of modern astrophysics. Recent hints for a possible solution involve lithium burning induced by protostellar mass accretion on Spite plateau stars. However, to date, most of the protostellar and pre-main sequence stellar models that take mass accretion into account have been computed at solar metallicity, and a detailed analysis on the impact of protostellar accretion on the lithium evolution in the metal-poor regime, which is relevant for stars in the Spite plateau, is completely missing. Aims. The purpose of this paper is to fill this gap, analysing, in detail, for the first time the effect of protostellar accretion on low metallicity low-mass stars with a focus on pre-main sequence lithium evolution. Methods. We computed the evolution from the protostar to the main-sequence phase of accreting models with final masses equal to 0.7 and 0.8 M⊙, and three metallicities Z = 0.0001, Z = 0.0010, and Z = 0.0050, corresponding to [Fe/H] ∼ −2.1, −1.1 (typical of Spite plateau stars), and [Fe/H] ∼ −0.42, respectively. We followed the temporal evolution of the chemical composition by considering nuclear burning, convective mixing, and diffusion. The effects of changing some of the main parameters affecting accreting models, that is the accretion energy (i.e. cold versus hot accretion), the initial seed mass Mseed and radius Rseed, and the mass accretion rate ṁ (also considering episodic accretion), have been investigated in detail. Results. As for the main stellar properties and in particular the surface 7Li abundance, hot accretion models converge to standard non-accreting ones within 1 Myr, regardless of the actual value of Mseed, Rseed, and ṁ. Also, cold accretion models with a relatively large Mseed (≳10 MJ) or Rseed (≳1 R⊙) converge to standard non-accreting ones in less than about 10−20 Myr. However, a drastically different evolution occurs whenever a cold protostellar accretion process starts from small values of Mseed and Rseed (Mseed ∼ 1 MJ, Rseed ≲ 1 R⊙). These models almost entirely skip the standard Hayashi track evolution and deplete lithium before the end of the accretion phase. The exact amount of depletion depends on the actual combination of the accretion parameters (ṁ, Mseed, and Rseed), achieving in some cases the complete exhaustion of lithium in the whole star. Finally, the lithium evolution in models accounting for burst accretion episodes or for an initial hot accretion followed by a cold accretion phase closely resemble that of standard non-accreting ones. Conclusions. To significantly deplete lithium in low-mass metal poor stars by means of protostellar accretion, a cold accretion scenario starting from small initial Mseed and Rseed is required. Even in this extreme configuration leading to a non-standard evolution that misses almost entirely the standard Hayashi track, an unsatisfactory fine tuning of the parameters governing the accretion phase is required to deplete lithium in stars of different mass and metallicity – starting from the Big Bang nucleosynthesis abundance – in such a way as to produce the observed Spite plateau.

The direct measurements of nuclear reaction cross sections at astrophysical energies are usually difficult to be performed in terrestrial laboratories, mainly because of the Coulomb barrier. An alternative way of studying these processes is via the indirect methods among which the Trojan Horse (THM) represents a powerful tool for studying charged-particles induced reactions at astrophysically energies without the need of extrapolation procedures. The recent THM studies about the lithium and boron destroying reactions will be here shown and discussed.

Aims. The capability of grid-based techniques to estimate the age together with the convective core overshooting efficiency of stars in detached eclipsing binary systems for main sequence stars has previously been investigated. We have extended this investigation to later evolutionary stages and have evaluated the bias and variability on the recovered age and convective core overshooting parameter accounting for both observational and internal uncertainties. Methods. We considered synthetic binary systems, whose age and overshooting efficiency should be recovered by applying the SCEPtER pipeline to the same grid of models used to build the mock stars. We focus our attention on a binary system composed of a 2.50 M⊙ primary star coupled with a 2.38 M⊙ secondary. To explore different evolutionary scenarios, we performed the estimation at three different times: when the primary is at the end of the central helium burning, when it is at the bottom of the RGB, and when it is in the helium core burning phase. The Monte Carlo simulations have been carried out for two typical values of accuracy on the mass determination, that is, 1% and 0.1%. Results. Adopting typical observational uncertainties, we found that the recovered age and overshooting efficiency are biased towards low values in all three scenarios. For an uncertainty on the masses of 1%, the underestimation is particularly 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 β. In the other scenarios, an undervaluation of the age by about 4% occurs. A large variability in the fitted values between Monte Carlo simulations was found: for an individual system calibration, the value of the overshooting parameter can vary from β = 0.0 to β = 0.26. When adopting a 0.1% error on the masses, the biases remain nearly unchanged but the global variability is suppressed by a factor of about two. We also explored the effect of a systematic discrepancy between the artificial systems and the model grid by accounting for an offset in the effective temperature of the stars by ± 150 K. For a mass error of 1% the overshooting parameter 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 also evaluate the possibility of individually recovering the β value for both binary stars. We found that this is impossible for a primary near to central hydrogen exhaustion owing to huge biases for the primary star of + 0.14 (90% relative error), while in the other cases the fitted β are consistent, but always biased by about − 0.04 (− 25% relative error). Finally, the possibility to distinguish between models computed with mild overshooting from models with no overshooting was evaluated, resulting in a reassuring power of distinction greater than 80%. However, the scenario with a primary in the central helium burning was a notable exception, showing a power of distinction lower than 5%.

Context. A promising technique to derive the physical parameters and the distance of pulsating stars is the fit of the observed light curves by nonlinear pulsation models. Aims. We apply this technique to a subsample of the RR Lyrae belonging to the Galactic globular cluster M3. The application of the method to cluster pulsators has the advantage of dealing with a homogeneous sample at the same distance and with the same chemical composition allowing to be checked the internal consistency of pulsational calculations. Methods. We selected seven pulsators (three RR c and four RR ab ) which cover a significant period range and show detailed light curves in the B, V and in some cases I bands. For four of them, with different periods, pulsation modes and light curve properties, we analyze the dependence of the theoretical light curve variations on the model input parameters. Results. For all selected objects, except the reddest one, we are able to theoretically reproduce the observed light curve morphology for self-consistent ranges of intrinsic stellar parameters, in agreement with the evolutionary predictions for the corresponding metal abundance. It is worth noting that, even if the theoretical reproduction of individual light curves gives for each variable only a range of stellar parameters and distances, the analysis of several variables belonging to the same cluster provides a mean distance modulus, namely μ v = 15.10 ± 0.1 mag, and also checks the self-consistency of the adopted theoretical scenario. Taking also into account the evolutionary constraints, the range of the accepted distance modulus is significantly reduced giving a weighted mean value of 15.05 ± 0.02. Our estimates are in agreement with available results in the literature obtained from independent methods.

We present recent work undertaken by the Evolution and Seismic Tools Activity (ESTA) team of the CoRoT Seismology Working Group. The new ESTA-Task 3 aims at testing, comparing and optimising stellar evolution codes which include microscopic diffusion of the chemical elements resulting from pressure, temperature and concentration gradients. The results already obtained are globally satisfactory, but some differences between the different numerical tools appear that require further investigations., to appear in "Stellar Evolution and Seismic Tools for Asteroseismology", (Eds) C. W. Straka, Y. Lebreton and M. J. P. F. G. Monteiro, EAS Publications Series, 2007

We supplement current evolutionary computations concerning Magellanic Cloud stars by exploring the evolutionary behavior of canonical stellar models (i.e.,with inefficient core overshooting) with metallicities suitable for stars in the Clouds. After discussing the adequacy of the adopted evolutionary scenario, we present evolutionary sequences as computed following a selected sample of stellar models in the mass range 0.8-8 Mo from the Main Sequence till the C ignition or the onset of thermal pulses in the advanced Asymptotic Giant Branch phase. On this basis, cluster isochrones covering the range of ages from 100 Myr to 15 Gyr are presented and discussed. To allow a comparison with evolutionary investigations appeared in the recent literature, we computed additional sets of models which take into account moderate core overshooting during the H burning phase, discussing the comparison in terms of current uncertainties in the stellar evolutionary models. Selected predictions constraining the cluster ages are finally discussed, presenting a calibration of the difference in magnitude between the luminous MS termination and the He burning giants in terms of cluster age. Both evolutionary tracks and isochrones have been made available at the node http://gipsy.cjb.net as a first step of a planned ``Pisa Evolutionary Library''., 11 pages, 9 eps figures, A&A accepted, evolutionary tracks and isochrones available at http://gipsy.cjb.net at the link ``Pisa Evolutionary Library''

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., Comment: Abstract shortened. Accepted for publication in A&A. One reference fixed

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., Comment: Abstract abridged; accepted for publication in A&A

The Trojan Horse Method (THM) has been largely adopted for investigating astrophysically relevant charged-particle induced reactions at Gamow energies. Indeed, THM allows one to by pass extrapolation procedures, thus overcoming this source of uncertainty. Here, the recent THM results and their impact in astrophysics are going to be discussed.

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%., A&A accepted

We present a new approach to estimate the ages of very young clusters based on the Pre-Main Sequence (PMS) Turn-On (TOn), point in the colour magnitude diagram where the PMS joins the main sequence (MS). In order to detect the TOn magnitude, we propose to monitor the spatial distribution of MS stars: any drop in the MS spatial density is a necessary signature for a successful TOn identification and, accordingly, a robust cluster dating. We use Hubble Space Telescope photometry of the star forming region NGC 346 as a benchmark for the method.

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., Accepteted for pubblication in A&A. The abstract is shortened to fill in the arxiv abstract field

The disagreement between theoretical predictions and observations for surface lithium abundance in stars is a long-standing problem, which indicates that the adopted physical treatment is still lacking in some points. However, thanks to the recent improvements in both models and observations, it is interesting to analyse the situation to evaluate present uncertainties. We present a consistent and quantitative analysis of the theoretical uncertainties affecting surface lithium abundance in the current generation of models. By means of an up-to-date and well tested evolutionary code, FRANEC, theoretical errors on surface 7Li abundance predictions, during the pre-main sequence (pre-MS) and main sequence (MS) phases, are discussed in detail. Then, the predicted surface 7Li abundance was tested against observational data for five open clusters, namely Ic 2602, \alpha Per, Blanco1, Pleiades, and Ngc 2516, and for four detached double-lined eclipsing binary systems. Stellar models for the aforementioned clusters were computed by adopting suitable chemical composition, age, and mixing length parameter for MS stars determined from the analysis of the colour-magnitude diagram of each cluster. We restricted our analysis to young clusters, to avoid additional uncertainty sources such as diffusion and/or radiative levitation efficiency. We confirm the disagreement, within present uncertainties, between theoretical predictions and 7Li observations for standard models. However, we notice that a satisfactory agreement with observations for 7Li abundance in both young open clusters and binary systems can be achieved if a lower convection efficiency is adopted during the pre-MS phase with respect to the MS one.