We collected over 6000 high-resolution spectra of four dozen field RR Lyrae (RRL) variables pulsating either in the fundamental (39 RRab) or in the first overtone (9 RRc) mode. We measured radial velocities (RVs) of four strong metallic and four Balmer lines along the entire pulsational cycle and derived RV amplitudes with accuracies better than 1$-$2~\kmsec. The new amplitudes were combined with literature data for 23~RRab and 3~RRc stars (total sample 74 RRLs) which allowed us to investigate the variation of the Bailey diagram (photometric amplitude versus period) when moving from optical to mid-infrared bands and to re-cast the Bailey diagram in terms of RV amplitudes. We found that RV amplitudes for RRab are minimally affected by nonlinear phenomena (shocks) and multi-periodicity (Blazhko effect). The RV slope ($\log P$--A(V$_r$)) when compared with the visual slope ($\log P$--A($V$)) is shallower and the dispersion, at fixed period, decreases by a factor of two. We constructed homogeneous sets of Horizontal Branch evolutionary models and nonlinear, convective pulsation models of RRLs to constrain the impact of evolutionary effects on their pulsation properties. Evolution causes, on the Bailey diagram based on RV amplitudes, a modest variation in pulsation period and a large dispersion in amplitude. The broad dispersion in period of the Bailey diagram is mainly caused by variation in RRL intrinsic parameters (stellar mass, chemical composition). Empirical evidence indicates that RV amplitudes are an optimal diagnostic for tracing the mean effective temperature across the RRab instability strip.

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. The purpose of this paper is to analyze the effect of protostellar accretion on low metallicity low-mass stars with a focus on PMS lithium evolution. We computed the evolution from the protostar to the MS phase of accreting models with final masses of 0.7 and 0.8 M$_\odot$, and three metallicities Z=0.0001, Z=0.0010, and Z=0.0050. The effects of changing the main parameters affecting accreting models, that is the accretion energy (cold versus hot accretion), the initial seed mass $M_{seed}$ and radius $R_{seed}$, and the mass accretion rate $\dot{m}$, have been investigated in detail. As for the main stellar properties and the surface $^7 Li$ abundance, hot accretion models converge to standard non-accreting ones within 1 Myr, regardless of the actual value of $M_{seed}$, $R_{seed}$, and $\dot{m}$. Also, cold accretion models with a relatively large $M_{seed}$ ($\gtrsim 10~M_{jup}$) or $R_{seed}$ ($\gtrsim 1~R_\odot$) converge to standard non-accreting ones in less than about 10-20~Myr. A drastically different evolution occurs whenever a cold protostellar accretion process starts from small values of $M_{seed}$ and $R_{seed}$ ($M_{seed}\sim 1~M_{jup}$, $R_{seed} \lesssim 1~R_\odot$). These models almost entirely skip the standard Hayashi track evolution and deplete Li before the end of the accretion phase. The exact amount of depletion depends on the actual combination of the accretion parameters ($\dot{m}$, $M_{seed}$, and $R_{seed}$), achieving in some cases the complete exhaustion of Li in the whole star.

The cosmological lithium problem, i.e. 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. A possible astrophysical solution involves lithium burning due to protostellar mass accretion on Spite plateau stars. In present work, for the first time, we investigate with accurate evolutionary computations the impact of accretion on the lithium evolution in the metal-poor regime, that relevant for stars in the Spite plateau.

Studies of fundamental parameters of very low-mass objects are indispensable to provide tests of stellar evolution models that are used to derive theoretical masses of brown dwarfs and planets. However, only objects with dynamically determined masses and precise photometry can effectively evaluate the predictions of stellar models. AB Dor C (0.090 solar masses) has become a prime benchmark for calibration of theoretical evolutionary models of low-mass young stars. One of the ambiguities remaining in AB Dor C is the possible binary nature of this star. We observed AB Dor C with the VLTI/AMBER instrument in low-resolution mode at the J, H and K bands. The interferometric observables at the K-band are compatible with a binary brown dwarf system with tentative components AB Dor Ca/Cb with a K-band flux ratio of 5$\pm$1% and a separation of 38$\pm$1 mas. This implies theoretical masses of 0.072$\pm$0.013 M$_{\rm \odot}$ and 0.013$\pm$0.001 M$_{\rm \odot}$ for each component, near the hydrogen-burning limit for AB Dor Ca, and near the deuterium-burning limit, straddling the boundary between brown dwarfs and giant planets, for AB Dor Cb. The possible binarity of AB Dor C alleviates the disagreement between observed magnitudes and theoretical mass-luminosity relationships.

In the era of large high-resolution spectroscopic surveys, high-quality spectra can contribute to our understanding of the Galactic chemical evolution, providing chemical abundances belonging to the different nucleosynthesis channels, and also providing constraints to stellar age. Some abundance ratios have been proven to be excellent indicators of stellar ages. We aim at providing an empirical relationship between stellar ages and [C/N] using, as calibrators, open star clusters observed by both the Gaia-ESO and APOGEE surveys. We use stellar parameters and abundances from the Gaia-ESO and APOGEE of the Galactic field and open cluster stars. Ages of star clusters are retrieved from the literature sources and validated using a common set of isochrones. We use the same isochrones to determine, for each age and metallicity, the surface gravity at which the first dredge-up and red giant branch bump occur. We study the effect of extra-mixing processes in our sample of giant stars, and we derive the mean [C/N] in evolved stars, including only stars without evidence of extra-mixing. Combining the Gaia-ESO and APOGEE samples of open clusters, we derive a linear relationship between [C/N] and logarithmic cluster ages. We apply our relationship to selected giant field stars in both Gaia-ESO and APOGEE. We find an age separation between thin and thick disc stars and age trends within their populations, with an increasing age towards lower metallicity populations. With such empirical relationship, we are able to provide an age estimate for giant stars in which C and N abundances are measured. Isochrone fitting is less sensitive for giant than dwarf stars at the turn off. The present method can be thus considered as an additional tool to give an independent estimate of the age of giant stars, with uncertainties in their ages comparable to those obtained using isochrone fitting for dwarf stars.

Pismis 18 is a moderately populated, intermediate-age open cluster located within the solar circle at a Galactocentric distance of about 7 kpc. Few open clusters have been studied in detail in the inner disc region before the Gaia-ESO Survey. New data from the Gaia-ESO Survey allowed us to conduct an extended radial velocity membership study as well as spectroscopic metallicity and detailed chemical abundance measurements for this cluster. Gaia-ESO Survey data for 142 potential members, lying on the upper MS and on the red clump, yielded radial velocity measurements, which, together with proper motion measurements from the Gaia DR2, were used to determine the systemic velocity of the cluster and membership of individual stars. Photometry from Gaia DR2 was used to re-determine cluster parameters based on high confidence member stars only. Cluster abundance measurements of six radial-velocity member stars with UVES high-resolution spectroscopy are presented for 23 elements. According to the new estimates, based on high confidence members, Pismis 18 has an age of $700^{+40}_{-50}$ Myr, interstellar reddening of E(B-V) = $0.562^{+0.012}_{-0.026}$ mag and a de-reddened distance modulus of $DM_0 = 11.96^{+0.10}_{-0.24}$ mag. The median metallicity of the cluster (using the six UVES stars) is [Fe/H] = $+0.23 \pm 0.05$ dex, with [$\alpha$/Fe]= $0.07 \pm 0.13$ and a slight enhancement of s- and r- neutron-capture elements. With the present work, we fully characterized the open cluster Pismis 18, confirming its present location in the inner disc. We estimated a younger age than the previous literature values and gave, for the first time, its metallicity and its detailed abundances. Its [$\alpha$/Fe] and [s-process/Fe], both slightly super-solar, are in agreement with other inner-disc open clusters observed by the Gaia-ESO survey. [abridged]

The high-quality OmegaCAM photometry of the 3x3 deg around the Orion Nebula Cluster (ONC) in r, and i filters by Beccari et al.(2017) revealed three well-separated pre-main sequences in the color-magnitude diagram (CMD). The objects belonging to the individual sequences are concentrated towards the center of the ONC. The authors concluded that there are two competitive scenarios: a population of unresolved binaries and triples with an exotic mass ratio distribution, or three stellar populations with different ages. We use Gaia DR2 in combination with the photometric OmegaCAM catalog to test and confirm the presence of the putative three stellar populations. We also study multiple stellar systems in the ONC for the first time using Gaia DR2. We confirm that the second and third sequence members are more centrally concentrated towards the center of the ONC. In addition we find an indication that the parallax and proper motion distributions are different among the members of the stellar sequences. The age difference among stellar populations is estimated to be 1-2 Myr. We use Gaia measurements to identify and remove as many unresolved multiple system candidates as possible. Nevertheless we are still able to recover two well-separated sequences with evidence for the third one, supporting the existence of the three stellar populations. We were able to identify a substantial number of wide binary objects (separation between 1000-3000 au). This challenges previously inferred values that suggested no wide binary stars exist in the ONC. Our inferred wide-binary fraction is approx 5%. We confirm the three populations correspond to three separated episodes of star formation. Based on this result, we conclude that star formation is not happening in a single burst in this region. (abridged)

In very young clusters, stellar age distribution is the empirical proof of the duration of star formation (SF) and of the physical mechanisms involved in the process. We derived accurate stellar ages for the cluster NGC6530, associated with the Lagoon Nebula to infer its SF history. We use the Gaia-ESO survey observations and Gaia DR2 data, to derive cluster membership and fundamental stellar parameters. We identified 652 confirmed and 9 probable members. The reddening inferred for members and non-members allows us to distinguish MS stars and giants, in agreement with the distances inferred from Gaia DR2 data. The foreground and background stars show a spatial pattern that traces the 3D structure of the nebular dust component. We derive stellar ages for 382 confirmed cluster members and we find that the gravity-sensitive gamma index distribution for M stars is correlated with stellar age. For all members with Teff<5500 K, the mean logarithmic age is 5.84 (units of years) with a dispersion of 0.36 dex. The age distribution of stars with accretion and/or disk (CTTSe) is similar to that of stars without accretion and without disk (WTTSp). We interpret this dispersion as evidence of a real age spread since the total uncertainties on age determinations, derived from Monte Carlo simulations, are significantly smaller than the observed spread. This conclusion is supported by the evidence of a decreasing of the gravity-sensitive gamma index as a function of stellar ages. The presence of the age spread is also supported by the spatial distribution and the kinematics of old and young members. In particular, members with accretion and/or disk, formed in the last 1 Myr, show evidence of subclustering around the cluster center, in the Hourglass Nebula and in the M8-E region, suggesting a possible triggering of star formation events by the O-type star ionization fronts.

One of the aims of LSST is to perform a systematic survey of star clusters and star forming regions (SFRs) in our Galaxy. In particular, the observations obtained with LSST will make a big difference in Galactic regions that have been poorly studied in the past, such as the anticenter and the disk beyond the Galactic center, and they will have a strong impact in discovering new distant SFRs. These results can be achieved by exploiting the exquisite depth that will be attained if the wide-fast-deep (WFD) observing strategy of the main survey is also adopted for the Galactic plane, in the g, r, and i filters.

The main aim of this experiment is to provide a complete census of old (t > 10 Gyr, RR Lyrae, type II Cepheids, red horizontal branch), intermediate age (red clump, Miras) and young (classical Cepheids) stellar tracers across the Galactic Bulge. To fully exploit the unique photometric quality of LSST images, we plan to perform a Shallow minisurvey (ugrizy, -20 < l < 20 deg, -15 < b < 10 deg) and a Deep minisurvey (izy, -20 < l < 20 deg, -3 < b < 3 deg). The former one is aimed at constraining the 3D structure of the galactic Bulge across the four quadrants, and in particular, the transition between inner and outer Bulge. The u,g,r,i,z,y LSST bands provide fundamental diagnostics to constrain the evolutionary properties of low and intermediate-mass stars when moving from a metal-poor to a metal-rich regime. The deep minisurvey is aimed at tracing RR Lyrae, Red Clump stars, Miras and classical Cepheids in highly reddened regions of the Galactic center. These images will allow us to investigate the role that baryonic mass and dark matter played in the early formation and evolution of the MW.

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.

Most stars born in clusters and recent results suggest that star formation (SF) preferentially occurs in subclusters. Studying the morphology and SF history of young clusters is crucial to understanding early SF. We identify the embedded clusters of young stellar objects (YSOs) down to M stars, in the HII regions RCW33, RCW32 and RCW27 of the Vela Molecular Ridge. Our aim is to characterise their properties, such as morphology and extent of the clusters in the three HII regions, derive stellar ages and the connection of the SF history with the environment. Through public photometric surveys such as Gaia, VPHAS, 2MASS and Spitzer/GLIMPSE, we identify YSOs with IR, Halpha and UV excesses, as signature of circumstellar disks and accretion. In addition, we implement a method to distinguish M dwarfs and giants, by comparing the reddening derived in several optical/IR color-color diagrams, assuming suitable theoretical models. Since this diagnostic is sensitive to stellar gravity, the procedure allows us to identify pre-main sequence stars. We find a large population of YSOs showing signatures of circumstellar disks with or without accretion. In addition, with the new technique of M-type star selection, we find a rich population of young M stars with a spatial distribution strongly correlated to the more massive population. We find evidence of three young clusters, with different morphology. In addition, we identify field stars falling in the same region, by securely classifying them as giants and foreground MS stars. We identify the embedded population of YSOs, down to about 0.1 Msun, associated with the HII regions RCW33, RCW32 and RCW27 and the clusters Vela T2, Cr197 and Vela T1, respectively, showing very different morphologies. Our results suggest a decreasing SF rate in Vela T2 and triggered SF in Cr197 and Vela T1.

We performed an analysis of the main theoretical uncertainties that affect the radius of low- and very-low mass-stars predicted by current stellar models. We focused on stars in the mass range 0.1-1Msun, on both the zero-age main-sequence (ZAMS) and on 1, 2 and 5 Gyr isochrones. First, we quantified the impact on the radius of the uncertainty of several quantities, namely the equation of state, radiative opacity, atmospheric models, convection efficiency and initial chemical composition. Then, we computed the cumulative radius error stripe obtained by adding the radius variation due to all the analysed quantities. As a general trend, the radius uncertainty increases with the stellar mass. For ZAMS structures the cumulative error stripe of very-low mass stars is about $\pm 2$ and $\pm 3$ percent, while at larger masses it increases up to $\pm 4$ and $\pm 5$ percent. The radius uncertainty gets larger and age dependent if isochrones are considered, reaching for $M\sim 1$Msun about $+12(-15)$ percent at an age of 5 Gyr. We also investigated the radius uncertainty at a fixed luminosity. In this case, the cumulative error stripe is the same for both ZAMS and isochrone models and it ranges from about $\pm 4$ percent to $+7$ and $+9$($-5$) percent. We also showed that the sole uncertainty on the chemical composition plays an important role in determining the radius error stripe, producing a radius variation that ranges between about $\pm 1$ and $\pm 2$ percent on ZAMS models with fixed mass and about $\pm 3$ and $\pm 5$ percent at a fixed luminosity.

We present a new complete Near-Infrared (NIR, $JHK_s$) census of RR Lyrae stars (RRLs) in the globular $\omega$ Cen (NGC 5139). We collected 15,472 $JHK_s$ images with 4-8m class telescopes over 15 years (2000-2015) covering a sky area around the cluster center of 60x34 arcmin$^2$. These images provided calibrated photometry for 182 out of the 198 cluster RRL candidates with ten to sixty measurements per band. We also provide new homogeneous estimates of the photometric amplitude for 180 ($J$), 176 ($H$) and 174 ($K_s$) RRLs. These data were supplemented with single-epoch $JK_s$ magnitudes from VHS and with single-epoch $H$ magnitudes from 2MASS. Using proprietary optical and NIR data together with new optical light curves (ASAS-SN) we also updated pulsation periods for 59 candidate RRLs. As a whole, we provide $JHK_s$ magnitudes for 90 RRab (fundamentals), 103 RRc (first overtones) and one RRd (mixed--mode pulsator). We found that NIR/optical photometric amplitude ratios increase when moving from first overtone to fundamental and to long-period (P>0.7 days) fundamental RRLs. Using predicted Period-Luminosity-Metallicity relations, we derive a true distance modulus of 13.674$\pm$0.008$\pm$0.038 mag (statistical error and standard deviation of the median)---based on spectroscopic iron abundances---and of 13.698$\pm$0.004$\pm$0.048 mag---based on photometric iron abundances. We also found evidence of possible systematics at the 5-10% level in the zero-point of the PLs based on the five calibrating RRLs whose parallaxes had been determined with HST

We describe the methodologies that, taking advantage of Gaia-DR1 and the Gaia-ESO Survey data, enable the comparison of observed open star cluster sequences with stellar evolutionary models. The final, long-term goal is the exploitation of open clusters as age calibrators. We perform a homogeneous analysis of eight open clusters using the Gaia-DR1 TGAS catalogue for bright members, and information from the Gaia-ESO Survey for fainter stars. Cluster membership probabilities for the Gaia-ESO Survey targets are derived based on several spectroscopic tracers. The Gaia-ESO Survey also provides the cluster chemical composition. We obtain cluster parallaxes using two methods. The first one relies on the astrometric selection of a sample of bona fide members, while the other one fits the parallax distribution of a larger sample of TGAS sources. Ages and reddening values are recovered through a Bayesian analysis using the 2MASS magnitudes and three sets of standard models. Lithium depletion boundary (LDB) ages are also determined using literature observations and the same models employed for the Bayesian analysis. For all but one cluster, parallaxes derived by us agree with those presented in Gaia Collaboration et al. (2017), while a discrepancy is found for NGC 2516; we provide evidence supporting our own determination. Inferred cluster ages are robust against models and are generally consistent with literature values. The systematic parallax errors inherent in the Gaia DR1 data presently limit the precision of our results. Nevertheless, we have been able to place these eight clusters onto the same age scale for the first time, with good agreement between isochronal and LDB ages where there is overlap. Our approach appears promising and demonstrates the potential of combining Gaia and ground-based spectroscopic datasets.

We present new theoretical period-luminosity-metallicity (PLZ) relations for RR Lyrae stars (RRL) at Spitzer and WISE wavelengths. The PLZ relations were derived using nonlinear, time-dependent convective hydrodynamical models for a broad range in metal abundances (Z=0.0001 to 0.0198). In deriving the light curves, we tested two sets of atmospheric models (Brott & Hauschildt 2005, Castelli & Kurucz 2003) and found no significant difference between the resulting mean magnitudes. We also compare our theoretical relations to empirical relations derived from RRL in both the field and in the globular cluster M4. Our theoretical PLZ relations were combined with multi-wavelength observations to simultaneously fit the distance modulus, mu_0, and extinction, Av, of both the individual Galactic RRL and of the cluster M4. The results for the Galactic RRL are consistent with trigonometric parallax measurements from Gaia's first data release. For M4, we find a distance modulus of $\mu_0=11.257 \pm 0.035$ mag with $A_V = 1.45 \pm 0.12$ mag, which is consistent with measurements from other distance indicators. This analysis has shown that when considering a sample covering a range of iron abundances, the metallicity spread introduces a dispersion in the PL relation on the order of 0.13 mag. However, if this metallicity component is accounted for in a PLZ relation, the dispersion is reduced to ~0.02 mag at MIR wavelengths.

Context. Precise determination of stellar masses is necessary to test the validity of pre-main-sequence (PMS) stellar evolutionary models, whose predictions are in disagreement with measurements for masses below 1.2 Msun. To improve such a test, and based on our previous studies, we selected the AB Doradus moving group (AB Dor-MG) as the best-suited association on which to apply radio-based high-precision astrometric techniques to study binary systems. Aims. We seek to determine precise estimates of the masses of a set of stars belonging to the AB Dor-MG using radio and infrared observations. Methods. We observed in phase-reference mode with the Very Large Array (VLA) at 5 GHz and with the European VLBI Network (EVN) at 8.4 GHz the stars HD 160934, EK Dra, PW And, and LO Peg. We also observed some of these stars with the near-infrared CCD AstraLux camera at the Calar Alto observatory to complement the radio observations. Results. We determine model-independent dynamical masses of both components of the star HD 160934, A and c, which are 0.70+/-0.07 Msun and 0.45+/-0.04 Msun , respectively. We revised the orbital parameters of EK Dra and we determine a sum of the masses of the system of 1.38+/-0.08 Msun. We also explored the binarity of the stars LO Peg and PW And. Conclusions. We found observational evidence that PMS evolutionary models underpredict the mass of PMS stars by 10%-40%, as previously reported by other authors. We also inferred that the origin of the radio emission must be similar in all observed stars, that is, extreme magnetic activity of the stellar corona that triggers gyrosynchrotron emission from non-thermal, accelerated electrons.

The young (~2 Myr) cluster Chamaeleon I is one of the closest laboratories to study the early stages of star cluster dynamics in a low-density environment. We studied its structural and kinematical properties combining parameters from the high-resolution spectroscopic survey Gaia-ESO with data from the literature. Our main result is the evidence of a large discrepancy between the velocity dispersion (sigma = 1.14 \pm 0.35 km s^{-1}) of the stellar population and the dispersion of the pre-stellar cores (~0.3 km s^{-1}) derived from submillimeter observations. The origin of this discrepancy, which has been observed in other young star clusters is not clear. It may be due to either the effect of the magnetic field on the protostars and the filaments, or to the dynamical evolution of stars driven by two-body interactions. Furthermore, the analysis of the kinematic properties of the stellar population put in evidence a significant velocity shift (~1 km s^{-1}) between the two sub-clusters located around the North and South main clouds. This result further supports a scenario, where clusters form from the evolution of multiple substructures rather than from a monolithic collapse. Using three independent spectroscopic indicators (the gravity indicator $\gamma$, the equivalent width of the Li line, and the H_alpha 10\% width), we performed a new membership selection. We found six new cluster members located in the outer region of the cluster. Starting from the positions and masses of the cluster members, we derived the level of substructure Q, the surface density \Sigma and the level of mass segregation $\Lambda_{MSR}$ of the cluster. The comparison between these structural properties and the results of N-body simulations suggests that the cluster formed in a low density environment, in virial equilibrium or supervirial, and highly substructured.

New accurate and homogeneous optical UBVRI photometry has been obtained for variable stars in the Galactic globular $\omega$ Cen (NGC 5139). We secured 8202 CCD images covering a time interval of 24 years and a sky area of 84x48 arcmin. The current data were complemented with data available in the literature and provided new, homogeneous pulsation parameters (mean magnitudes, luminosity amplitudes,periods) for 187 candidate $\omega$ Cen RR Lyrae (RRLs). Among them we have 101 RRc (first overtone), 85 RRab (fundamental) and a single candidate RRd (double-mode) variables. Candidate Blazhko RRLs show periods and colors that are intermediate between RRc and RRab variables, suggesting that they are transitional objects. The comparison of the period distribution and of the Bailey diagram indicates that RRLs in $\omega$ Cen show a long-period tail not present in typical Oosterhoff II (OoII) globulars. The RRLs in dwarf spheroidals and in ultra faint dwarfs have properties between Oosterhoff intermediate and OoII clusters. Metallicity plays a key role in shaping the above evidence. These findings do not support the hypothesis that $\omega$ Cen is the core remnant of a spoiled dwarf galaxy. Using optical Period-Wesenheit relations that are reddening-free and minimally dependent on metallicity we find a mean distance to $\omega$ Cen of 13.71$\pm$0.08$\pm$0.01 mag (semi-empirical and theoretical calibrations). Finally, we invert the I-band Period-Luminosity-Metallicity relation to estimate individual RRLs metal abundances. The metallicity distribution agrees quite well with spectroscopic and photometric metallicity estimates available in the literature.

We analysed the effects of planet ingestion on the characteristics of a pre-MS star similar to the Gamma Velorum cluster member 2MASS J08095427--4721419 (#52). We discussed the effects of changing the age $t_0$ at which the accretion episode occurs, the mass of the ingested planet and its chemical composition. We showed that the mass of the ingested planet required to explain the current [Fe/H]^#52 increases by decreasing the age $t_0$ and/or by decreasing the Iron content of the accreted matter. We compared the predictions of a simplified accretion method -- where only the variation of the surface chemical composition is considered -- with that of a full accretion model that properly accounts for the modification of the stellar structure. We showed that the two approaches result in different convective envelope extension which can vary up to 10 percent. We discussed the impact of the planet ingestion on a stellar model in the colour-magnitude diagram, showing that a maximum shift of about 0.06 dex in the colour and 0.07 dex in magnitude are expected and that such variations persist even much later the accretion episode. We also analysed the systematic bias in the stellar mass and age inferred by using a grid of standard non accreting models to recover the characteristics of an accreting star. We found that standard non accreting models can safely be adopted for mass estimate, as the bias is <= 6 percent, while much more caution should be used for age estimate where the differences can reach about 60 percent.

The Goldreich-Schubert-Fricke (GSF) instability is a rotational instability that is thought to contribute to the transfer of angular momentum in differentially rotating stars. It has been included in recent codes of stellar evolution in a diffusion-like approximation, under the assumption that the kinematic viscosity $\nu$ is unimportant for the development of the instability. As noted previously by other authors, for most stellar applications this may not be a valid approximation. We discuss this issue in detail, solving the dispersion relation of the perturbed modes for realistic values of $\nu$ in the bulk of the radiative zone of the Sun and of three red giant stars at various evolutionary stages. We find that the GSF instability is triggered only in layers of extremely strong shear. In a simple case study, we also investigate the effect of a small deviation from axisymmetry or a small background magnetic field. We find that, like the viscosity, these have a stabilising effect. We conclude that this instability is probably far less efficient in transporting angular momentum than is often assumed, and may not even be present.

There is a growing body of evidences for the presence of crystalline material in comets. These crystals are believed to have been annealed in the inner part of the proto-solar nebula, while comets should have been formed in the outer regions. Several transport processes have been proposed to reconcile these two facts; among them a migration driven by photophoresis. The primarily goal of this work is to assess whether disk irradiation by a Pre-Main Sequence star would influence the photophoretic transport. To do so, we have implemented an evolving 1+1D model of an accretion disk, including advanced numerical techniques, undergoing a time-dependent irradiation, consistent with the evolution of the proto-Sun along the Pre-Main Sequence. The photophoresis is described using a formalism introduced in several previous works. Adopting the opacity prescription used in these former studies, we find that the disk irradiation enhances the photophoretic transport: the assumption of a disk central hole of several astronomical units in radius is no longer strictly required, whereas the need for an ad hoc introduction of photoevaporation is reduced. However, we show that a residual trail of small particles could annihilate the photophoretic driven transport via their effect on the opacity. We have also confirmed that the thermal conductivity of transported aggregates is a crucial parameter which could limit or even suppress the photophoretic migration and generate several segregation effects.

We investigated the possibility of producing helium enhanced stars in globular clusters by accreting polluted matter during the pre-main sequence phase. We followed the evolution of two different classes of pre-main sequence accreting models, one which neglects and the other that takes into account the protostellar evolution. We analysed the dependence of the final central helium abundance, of the tracks position in the HR diagram and of the surface lithium abundance evolution on the age at which the accretion of polluted material begins and on the main physical parameters that govern the protostellar evolution. The later is the beginning of the late accretion and the lower are both the central helium and the surface lithium abundances at the end of the accretion phase and in ZAMS (Zero Age Main Sequence). In order to produce a relevant increase of the central helium content the accretion of polluted matter should start at ages lower than 1 Myr. The inclusion of the protostellar evolution has a strong impact on the ZAMS models too. The adoption of a very low seed mass (i.e. 0.001 M$_{\odot}$) results in models with the lowest central helium and surface lithium abundances. The higher is the accretion rate and the lower is the final helium content in the core and the residual surface lithium. In the worst case -- i.e. seed mass 0.001 M$_\odot$ and accretion rate $\ge 10^{-5}$ M$_\odot$ yr$^{-1}$ -- the central helium is not increased at all and the surface lithium is fully depleted in the first few million years.

Age and mass determinations for isolated stellar objects remain model-dependent. While stellar interior and atmospheric theoretical models are rapidly evolving, we need a powerful tool to test them. Open clusters are good candidates for this role. We complement previous studies on the Hyades multiplicity by Lucky Imaging observations with the AstraLux Norte camera. This allows us to exclude possible binary and multiple systems with companions outside 2--7 AU separation and to create a "single-star sequence" for the Hyades. The sequence encompasses 250 main-sequence stars ranging from A5V to M6V. Using the "Tool for Astrophysical Data Analysis" (TA-DA), we create various theoretical isochrones applying different combinations of interior and atmospheric models. We compare the isochrones with the observed Hyades single-star sequence on J vs. J - K_s, J vs. J - H and K_s vs. H - K_s color-magnitude diagrams. As a reference we also compute absolute fluxes and magnitudes for all stars from X-ray to mid-infrared based on photometric measurements available in the literature(ROSAT X-ray, GALEX UV, APASS gri, 2MASS JHK_s, and WISE W1 to W).We find that combinations of both PISA and DARTMOUTH stellar interior models with BT-Settl 2010 atmospheric models describe the observed sequence well. The full sequence covers the mass range 0.13 to 2.3 Msun, and effective temperatures between 3060 K and 8200 K. Within the measurement uncertainties, the current generation of models agree well with the single-star sequence. The primary limitations are the uncertainties in the measurement of the distance to individual Hyades members, and uncertainties in the photometry. Additionally, a small (~0.05 mag) systematic offset can be noted on J vs. J - K and K vs. H - K diagrams - the observed sequence is shifted to redder colors from the theoretical predictions.

We investigated the absolute age of the Galactic globular cluster M71 (NGC 6838) by using optical ground-based images (u',g',r',i',z') collected with the MegaCam camera at the Canada-France-Hawaii-Telescope (CFHT). We performed a robust selection of field and cluster stars by applying a new method based on the 3D (r',u'-g',g'-r') Color-Color-Magnitude-Diagram. The comparison between the Color-Magnitude-Diagram of the candidate cluster stars and a new set of isochrones, at the locus of the Main Sequence Turn Off (MSTO), suggests an absolute age of 12+/-2 Gyr. The absolute age was also estimated using the difference in magnitude between the MSTO and the so-called main sequence knee, a well defined bending occurring in the lower main sequence. This feature was originally detected in the near-infrared (NIR) bands and explained as a consequence of an opacity mechanism (collisionally induced absorption of molecular hydrogen) in the atmosphere of cool low-mass stars (Bono et al. 2010). The same feature was also detected in the r',u'-g' and in the r',g'-r' CMD, thus supporting previous theoretical predictions by Borysow et al. (1997). The key advantage in using the Delta^(Knee)_(TO) as an age diagnostic is that it is independent of uncertainties affecting the distance, the reddening and the photometric zero-point. We found an absolute age of 12+/-1 Gyr that agrees, within the errors, with similar age estimates, but the uncertainty is on average a factor of two smaller. We also found that the Delta^(Knee)_(TO) is more sensitive to the metallicity than the MSTO, but the dependence becomes vanishing using the difference in color between the MSK and the MSTO.

Context. ABDoradus is the main system of the ABDoradus moving group. It is a quadruple system formed by two widely separated binaries of pre-main-sequence (PMS) stars: ABDor A/C and ABDor Ba/Bb. The pair ABDor A/C has been extensively studied and its dynamical masses have been determined with high precision, thus making of ABDor C a benchmark for calibrating PMS stellar models. If the orbit and dynamical masses of the pair ABDor Ba/Bb can be determined, they could not only play a similar role to that of ABDor C in calibrating PMS models, but would also help to better understand the dynamics of the whole ABDoradus system. Aims. We aim to determine the individual masses of the pair ABDor Ba/Bb using VLBI observations and archive infrared data, as part of a larger program directed to monitor binary systems in the ABDoradus moving group. Methods. We observed the system ABDor B between 2007 and 2013 with the Australian Long Baseline Array (LBA), at a frequency of 8.4 GHz in phase-reference mode. Results. We detected, for the first time, compact radio emission from both stars in the binary, ABDor Ba and ABDor Bb. This result allowed us to determine the orbital parameters of both the relative and absolute orbits and, consequently, their individual dynamical masses: 0.28+/-0.05 Msun and 0.25+/-0.05 Msun, respectively. Conclusions. Comparisons of the dynamical masses with the prediction of PMS evolutionary models show that the models underpredict the dynamical masses of the binary components Ba and Bb by 10-30% and 10-40%, respectively, although they all still agree at the 2-sigma level. Some of stellar models considered favour an age between 50 and 100 Myr for this system, meanwhile others predict older ages. We also discuss the evolutionary status of ABDor Ba/Bb in terms of an earlier double-double star scenario that might explain the strong radio emission detected in both components.

We performed a detailed analysis of the main theoretical uncertainties affecting the age at the lithium depletion boundary (LDB). To do that we computed almost 12000 pre-main sequence models with mass in the range [0.06, 0.4] M_sun by varying input physics (nuclear reaction cross-sections, plasma electron screening, outer boundary conditions, equation of state, and radiative opacity), initial chemical elements abundances (total metallicity, helium and deuterium abundances, and heavy elements mixture), and convection efficiency (mixing length parameter, alpha_ML). As a first step, we studied the effect of varying these quantities individually within their extreme values. Then, we analysed the impact of simultaneously perturbing the main input/parameters without an a priori assumption of independence. Such an approach allowed us to build for the first time the cumulative error stripe, which defines the edges of the maximum uncertainty region in the theoretical LDB age. We found that the cumulative error stripe is asymmetric and dependent on the adopted mixing length value. For alpha_ML = 1.00, the positive relative age error ranges from 5 to 15 per cent, while for solar-calibrated mixing length, the uncertainty reduces to 5-10 per cent. A large fraction of such an error (about 40 per cent) is due to the uncertainty in the adopted initial chemical elements abundances.

The p(p,e^+ \nu_e)^2H reaction rate is an essential ingredient for theoretical computations of stellar models. In the past several values of the corresponding S-factor have been made available by different authors. Prompted by a recent evaluation of S(E), we analysed the effect of the adoption of different proton-proton reaction rates on stellar models, focusing, in particular, on the age of mid and old stellar clusters (1-12 Gyr) and on standard solar model predictions. By comparing different widely adopted p(p,e^+ \nu_e)^2H reaction rates, we found a maximum difference in the temperature regimes typical of main sequence hydrogen-burning stars (5x10^6 - 3x10^7 K) of about 3%. Such a variation translates into a change of cluster age determination lower than 1%. A slightly larger effect is observed in the predicted solar neutrino fluxes with a maximum difference, in the worst case, of about 8%. Finally we also notice that the uncertainty evaluation of the present proton-proton rate is at the level of few \permil, thus the p(p,e^+ \nu_e)^2H reaction rate does not constitute anymore a significant uncertainty source in stellar models.

The dust cloud TGU H645 P2 and embedded in it young open cluster NGC 7129 are investigated using the results of medium-band photometry of 159 stars in the Vilnius seven-colour system down to V = 18.8 mag. The photometric data were used to classify about 50 percent of the measured stars in spectral and luminosity classes. The extinction A_V vs. distance diagram for the 20x20 arcmin area is plotted for 155 stars with two-dimensional classification from the present and the previous catalogues. The extinction values found range between 0.6 and 3.4 mag. However, some red giants, located in the direction of the dense parts of the cloud, exhibit the infrared extinction equivalent up to A_V = 13 mag. The distance to the cloud (and the cluster) is found to be 1.15 kpc (the true distance modulus 10.30 mag). For determining the age of NGC 7129, a luminosity vs. temperature diagram for six cluster members of spectral classes B3 to A1 was compared with the Pisa pre-main-sequence evolution tracks and the Palla birthlines. The cluster can be as old as about 3 Myr, but star forming continues till now as witnessed by the presence in the cloud of many younger pre-main-sequence objects identified with photometry from 2MASS, Spitzer and WISE infrared surveys.

Embedded clusters like W3 Main are complex and dynamically evolving systems that represent an important phase of the star formation process. We aim at the characterization of the entire stellar content of W3 Main in a statistical sense to identify possible differences in evolutionary phase of the stellar populations and find clues about the formation mechanism of this massive embedded cluster. Methods. Deep JHKs imaging is used to derive the disk fraction, Ks-band luminosity functions and mass functions for several subregions in W3 Main. A two dimensional completeness analysis using artificial star experiments is applied as a crucial ingredient to assess realistic completeness limits for our photometry. We find an overall disk fraction of 7.7 $\pm$ 2.3%, radially varying from 9.4 $\pm$ 3.0 % in the central 1 pc to 5.6 $\pm$ 2.2 % in the outer parts of W3 Main. The mass functions derived for three subregions are consistent with a Kroupa and Chabrier mass function. The mass function of IRSN3 is complete down to 0.14 Msun and shows a break at M $\sim$ 0.5 Msun. We interpret the higher disk fraction in the center as evidence for a younger age of the cluster center. We find that the evolutionary sequence observed in the low-mass stellar population is consistent with the observed age spread among the massive stars. An analysis of the mass function variations does not show evidence for mass segregation. W3 Main is currently still actively forming stars, showing that the ionizing feedback of OB stars is confined to small areas ($\sim$ 0.5 pc). The FUV feedback might be influencing large regions of the cluster as suggested by the low overall disk fraction.

We report the discovery of Kepler-77b (alias KOI-127.01), a Saturn-mass transiting planet in a 3.6-day orbit around a metal-rich solar-like star. We combined the publicly available Kepler photometry (quarters 1-13) with high-resolution spectroscopy from the Sandiford@McDonald and FIES@NOT spectrographs. We derived the system parameters via a simultaneous joint fit to the photometric and radial velocity measurements. Our analysis is based on the Bayesian approach and is carried out by sampling the parameter posterior distributions using a Markov chain Monte Carlo simulation. Kepler-77b is a moderately inflated planet with a mass of Mp=0.430+/-0.032 Mjup, a radius of Rp=0.960+/-0.016 Rjup, and a bulk density of 0.603+/-0.055 g/cm^3. It orbits a slowly rotating (P=36+/-6 days) G5V star with M*=0.95+/-0.04 Msun, R*=0.99+/-0.02 Rsun, Teff=5520+/-60 K, [M/H]=0.20+/-0.05, that has an age of 7.5+/-2.0 Gyr. The lack of detectable planetary occultation with a depth higher than about 10 ppm implies a planet geometric and Bond albedo of Ag<0.087+/-0.008 and Ab<0.058+/-0.006, respectively, placing Kepler-77b among the gas-giant planets with the lowest albedo known so far. We found neither additional planetary transit signals nor transit-timing variations at a level of about 0.5 minutes, in accordance with the trend that close-in gas giant planets seem to belong to single-planet systems. The 106 transits observed in short-cadence mode by Kepler for nearly 1.2 years show no detectable signatures of the planet's passage in front of starspots. We explored the implications of the absence of detectable spot-crossing events for the inclination of the stellar spin-axis, the sky-projected spin-orbit obliquity, and the latitude of magnetically active regions.

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.

We present a new method to determine the age spread of resolved stellar populations in a starburst cluster. The method relies on a two-step process. In the first step, kinematic members of the cluster are identified based on multi-epoch astrometric monitoring. In the second step, a Bayesian analysis is carried out, comparing the observed photometric sequence of cluster members with sets of theoretical isochrones. When applying this methodology to optical and near-infrared high angular resolution Hubble Space Telescope (HST) and adaptive optics observations of the ~5 Myr old starburst cluster Westerlund 1 and ~2 Myr old starburst cluster NGC 3603 YC, we derive upper limits for the age spreads of 0.4 and 0.1 Myr, respectively. The results strongly suggest that star formation in these starburst clusters happened almost instantaneously.

Located at the tip of the wing of the Small Magellanic Cloud (SMC), the star-forming region NGC602/N90 is characterized by the HII nebular ring N90 and the young cluster of pre--main-sequence (PMS) and early-type main sequence stars NGC602. We present a thorough cluster analysis of the stellar sample identified with HST/ACS camera in the region. We show that apart from the central cluster, low-mass PMS stars are congregated in thirteen additional small compact sub-clusters at the periphery of NGC602. We find that the spatial distribution of the PMS stars is bimodal, with an unusually large fraction (~60%) of the total population being clustered, while the remaining is diffusely distributed in the inter-cluster area. From the corresponding color-magnitude diagrams we disentangle an age-difference of ~2.5Myr between NGC602 and the compact sub-clusters which appear younger. The diffuse PMS population appears to host stars as old as those in NGC602. Almost all detected PMS sub-clusters appear to be centrally concentrated. When the complete PMS stellar sample, including both clustered and diffused stars, is considered in our cluster analysis, it appears as a single centrally concentrated stellar agglomeration, covering the whole central area of the region. Considering also the hot massive stars of the system, we find evidence that this agglomeration is hierarchically structured. Based on our findings we propose a scenario, according to which the region NGC602/N90 experiences an active clustered star formation for the last ~5Myr. The central cluster NGC602 was formed first and rapidly started dissolving into its immediate ambient environment, possibly ejecting also massive stars found away from its center. Star formation continued in sub-clusters of a larger stellar agglomeration, introducing an age-spread of the order of 2.5Myr among the PMS populations.

Pre-main sequence (PMS) models provide invaluable tools for the study of star forming regions as they allow to assign masses and ages to young stars. Thus it is of primary importance to test the models against observations of PMS stars with dynamically determined mass. We developed a Bayesian method for testing the present generation of PMS models which allows for a quantitative comparison with observations, largely superseding the widely used isochrones and tracks qualitative superposition. Using the available PMS data we tested the newest PISA PMS models establishing their good agreement with the observations. The data cover a mass range from ~0.3 to ~3.1 Msun, temperatures from ~3x10^3 to ~1.2x10^4 K and luminosities from ~3x10^-2 to ~60 Lsun. Masses are correctly predicted within 20% of the observed values in most of the cases and for some of them the difference is as small as 5%. Nevertheless some discrepancies are also observed and critically discussed. By means of simulations, using typical observational errors, we evaluated the spread of log \tau_sim - log \tau_rec, i.e. simulated minus recovered ages distribution of the single objects. We also found that stars in binary systems simulated as coeval might be recovered as non coeval, due to observational errors. The actual fraction of fake non coevality is a complex function of the simulated ages, masses and mass ratios. We demonstrated that it is possible to recover the systems' ages with better precision than for single stars using the composite age-probability distribution, i.e. the product of the components' age distributions. Using this valuable tool we estimated the ages of the presently observed PMS binary systems.

In recent years new observations of pre-main sequence stars (pre-MS) with Z <= Zsun have been made available. To take full advantage of the continuously growing amount of data of pre-MS stars in different environments, we need to develop updated pre-MS models for a wide range of metallicity to assign reliable ages and masses to the observed stars. We present updated evolutionary pre-MS models and isochrones for a fine grid of mass, age, metallicity, and helium values. We use a standard and well-tested stellar evolutionary code (i.e. FRANEC), that adopts outer boundary conditions from detailed and realistic atmosphere models. In this code, we incorporate additional improvements to the physical inputs related to the equation of state and the low temperature radiative opacities essential to computing low-mass stellar models. We make available via internet a large database of pre-MS tracks and isochrones for a wide range of chemical compositions (Z = 0.0002 - 0.03), masses (M = 0.2 - 7.0 Msun), and ages (1 - 100 Myr) for a solar-calibrated mixing length parameter \alpha (i.e. 1.68). For each chemical composition, additional models were computed with two different mixing length values, namely \alpha = 1.2 and 1.9. Moreover, for Z >= 0.008, we also provided models with two different initial deuterium abundances. The characteristics of the models have been discussed in detail and compared with other work in the literature. The main uncertainties affecting theoretical predictions have been critically discussed. Comparisons with selected data indicate that there is close agreement between theory and observation.