Search

In this study we compute the equation of state and Rosseland mean opacity from temperatures of T~30000 K down to T~400 K, pushing the capabilities of the AESOPUS code (Marigo et al., 2022; Marigo & Aringer, 2009) into the regime where solid grains can form. The GGchem code (Woitke et al. 2018) is used to solve the chemistry for temperatures less than ~3000 K. Atoms, molecules, and dust grains in thermodynamic equilibrium are all included in the equation of state. To incorporate monochromatic atomic and molecular cross sections, an optimized opacity sampling technique is used. The Mie theory is employed to calculate the opacity of 43 grain species. Tables of Rosseland mean opacities for scaled-solar compositions are provided. Based on our computing resources, opacities for other chemical patterns, as well as various grain sizes, porosity, and shapes, can be easily computed upon user request to the corresponding author., Comment: 14 pages, 8 figures, resubmitted to ApJ following moderate revision

With an unprecedented astrometric and photometric data precision, Gaia EDR3 gives us, for the first time, the opportunity to systematically detect and map in the optical bands, the low mass populations of the star forming regions (SFRs) in the Milky Way. We provide a catalogue of the Gaia EDR3 data (photometry, proper motions and parallaxes) of the young stellar objects (YSOs) identified in the Galactic Plane (|b|<30 deg) within about 1.5 kpc. The catalogue of the SFRs to which they belong is also provided to study the properties of the very young clusters and put them in the context of the Galaxy structure. We applied the machine learning unsupervised clustering algorithm DBSCAN on a sample of Gaia EDR3 data photometrically selected on the region where very young stars (t<10 Myr) are expected to be found, with the aim to identify co-moving and spatially consistent stellar clusters. A subsample of 52 clusters, selected among the 7323 found with DBSCAN, has been used as template data set, to identify very young clusters from the pattern of the observed color-absolute magnitude diagrams through a pattern match process. We find 124440 candidate YSOs clustered in 354 SFRs and stellar clusters younger than 10 Myr and within about 1.5 Kpc. In addition, 65863 low mass members of 322 stellar clusters located within about 500 pc and with ages 10 Myr

Pre-main sequence models with inflated radii are needed to simultaneously reproduce the colour-magnitude diagram and the lithium depletion pattern in young open clusters. We tested a new set of PMS models including radius inflation due to starspots or magnetic inhibition of convection, using five clusters observed by the Gaia-ESO Survey, spanning the age range ~10-100 Myr where such effects could be important. Gaia-ESO radial velocities were combined with Gaia EDR3 astrometry to obtain clean lists of high-probability members for the five clusters. A Bayesian maximum likelihood method was adopted to derive the best model parameters and the cluster reddening and age. Models were calculated for different values of the mixing length parameter ($\alpha_{ML}=2.0$, 1.5 and 1.0), without spots or with effective spot coverage $\beta_{spot}=0.2$ and 0.4. To reproduce the CMD and the Li depletion pattern in Gamma Vel A and B and in 25 Ori we need both a reduced convection efficiency $\alpha_{ML}=1.0$ and an effective spot coverage of ~20%. We obtained ages of 18 Myr and 21 Myr for Gamma Vel A and B, respectively, and 19 Myr for 25 Ori. However, a single isochrone is not sufficient to account for the Li dispersion, and an increasing level of spot coverage as mass decreases seems to be required. The older clusters (NGC2451B at 30 Myr, NGC2547 at 35 Myr, and NGC2516 at 138 Myr) are consistent with standard models ($\alpha_{ML}=2.0$ and no spots) except at low masses: a 20% spot coverage seems to reproduce the sequence of M-type stars better and might explain the observed abundance spread. The quality of Gaia-ESO data combined with Gaia allows us to gain important insights on PMS evolution. Models including starspots can provide a consistent explanation of the cluster sequences and Li abundances of young clusters, although a range of starspot coverage is required to fully reproduce the data., Comment: Accepted for publication in A&A, 18 pages, 10 figures. Revised version after language editing, revised acknowledgements

Spectroscopy from the final internal data release of the Gaia-ESO Survey (GES) has been combined with Gaia EDR3 to assign membership probabilities to targets observed towards 63 Galactic open clusters and 7 globular clusters. The membership probabilities are based chiefly on maximum likelihood modelling of the 3D kinematics of the targets, separating them into cluster and field populations. From 43211 observed targets, 13985 are identified as highly probable cluster members ($P>0.9$), with an average membership probability of 0.993. The addition of GES radial velocities successfully drives down the fraction of false positives and we achieve better levels of discrimination in most clusters over the use of astrometric data alone, especially those at larger distances. Since the membership selection is almost purely kinematic, the union of this catalogue with GES and Gaia is ideal for investigating the photometric and chemical properties of clusters as a function of stellar mass, age and Galactic position., Comment: Accepted version for publication in MNRAS. 16 pages + 38 pages of Appendices

The finding of multiple stellar populations (MP), defined by patterns in the stellar element abundances, is nowadays considered a distinctive feature of globular clusters. However, while data availability and quality improved in last decades, this is not always true for the techniques adopted to their analysis, rising problems of objectivity of the claims and reproducibility. Using NGC 2808 as test case we show the use of well established statistical clustering methods. We focus the analysis to the RGB phase, where two data sets are available from recent literature for low- and high-resolution spectroscopy. We adopt both hierarchical clustering and partition methods. We explicitly address the usually neglected problem of measurement errors. The results of the clustering algorithms were subjected to silhouette width analysis to compare the performance of the split into different number of MP. For both data sets the results are at odd with those reported in the literature. Two MP are detected for both data sets, while the literature reports five and four MP from high- and low-resolution spectroscopy respectively. The silhouette analysis suggests that the population sub-structure is reliable for high-resolution spectroscopy data, while the actual existence of MP is questionable for the low-resolution spectroscopy data. The discrepancy with literature claims is explainable due to the difference of methods adopted to MP characterisation. By means of Monte Carlo simulations and multimodality statistical tests we show that the often adopted study of the histogram of the differences in some key elements is prone to multiple false positive findings. The adoption of statistically grounded methods, which adopt all the available information to subset the data and explicitly address the problem of data uncertainty, is of paramount importance to present more robust and reproducible researches., Comment: Match the A&A accepted version. Acknowledgments fixed

Theoretical prediction of surface stellar abundances of light elements -- lithium, beryllium, and boron -- represents one of the most interesting open problems in astrophysics. As well known, several measurements of 7-Li abundances in stellar atmospheres point out a disagreement between predictions and observations in different stellar evolutionary phases, rising doubts about the capability of present stellar models to precisely reproduce stellar envelope characteristics. Light elements are burned at relatively low temperatures (from 2 to 5 MK) and thus in the evolutionary stages of a star they are gradually destroyed at different depths of stellar interior, in dependence on the stellar mass. Their surface abundances are strongly influenced by the nuclear cross sections, by the extension of the convective envelope and by the temperature at its bottom, which depend on the characteristics of the star (mass and chemical composition) as well as on the energy transport in the convective stellar envelope. In recent years, a great effort has been made to improve the precision of light element burning cross sections. However, theoretical predictions surface light element abundance are challenging because they are also influenced by the uncertainties in the input physics adopted in the calculations as well as the efficiency of several standard and non-standard physical processes active in young stars (i.e. diffusion, radiative levitation, magnetic fields, rotation). Moreover, it is still not completely clear how much the previous protostellar evolution affects the characteristics of a stellar model and thus the light element depletion. This paper presents the state-of-the-art of theoretical predictions for protostars and PMS stars, discussing the role of several input physics on the stellar evolution.

We aim to constrain the mixing processes in low-mass stars by investigating the behaviour of the Li surface abundance after the main sequence. We take advantage of the data from the sixth internal data release of Gaia-ESO, idr6, and from the Gaia Early Data Release 3, edr3. We select a sample of main sequence, sub-giant, and giant stars in which Li abundance is measured by the Gaia-ESO survey, belonging to 57 open clusters with ages from 120~Myr to about 7 Gyr and to Milky Way fields, covering a range in [Fe/H] between -1.0 and +0.5dex. We study the behaviour of the Li abundances as a function of stellar parameters. We compare the observed Li behaviour in field giant stars and in giant stars belonging to individual clusters with the predictions of a set of classical models and of models with mixing induced by rotation and thermohaline instability. The comparison with stellar evolution models confirms that classical models cannot reproduce the lithium abundances observed in the metallicity and mass regimes covered by the data. The models that include the effects of both rotation-induced mixing and thermohaline instability account for the Li abundance trends observed in our sample, in all metallicity and mass ranges. The differences between the results of the classical models and of the rotation models largely differ (up to ~2 dex), making lithium the best element to constrain stellar mixing processes in low-mass stars. For stars with well-determined masses, we find a better agreement between observed surface abundances and models with rotation-induced and thermohaline mixings, the former dominating during the main sequence and the first phases of the post-main sequence evolution and the latter after the bump in the luminosity function., Comment: Accepted for publication in A&A, 19 pages, 16 figures

..., Comment: Accepted for publication in A&A. Abstract abridged

We tested the capability of a Bayesian procedure to calibrate both the helium abundance and the mixing length parameter ($\alpha_{ML}$), using precise photometric data for main-sequence (MS) stars in a cluster with negligible reddening and well-determined distance. The method has been applied first to a mock data set generated to mimic Hyades MS stars and then to the real Hyades cluster. We tested the impact on the results of varying the number of stars in the sample, the photometric errors, and the estimated [Fe/H]. The analysis of the synthetic data set shows that $\alpha_{ML}$ is recovered with a very good precision in all the analysed cases (with an error of few percent), while [Fe/H] and the helium-to-metal enrichment ratio $\Delta Y/\Delta Z$ are more problematic. If spectroscopic determinations of [Fe/H] are not available and thus [Fe/H] has to be recovered alongside with $\Delta Y/\Delta Z$ and $\alpha_{ML}$, the well-known degeneracy between [Fe/H]-$\Delta Y/\Delta Z$-$\alpha_{ML}$ could result in a large uncertainty on the recovered parameters, depending on the portion of the MS used for the analysis. On the other hand, the prior knowledge of an accurate [Fe/H] value puts a strong constraint on the models, leading to a more precise parameters recovery. Using the current set of PISA models, the most recent [Fe/H] value and the Gaia photometry and parallaxes for the Hyades cluster, we obtained the average values $<\alpha_{ML}>=2.01\pm0.05$ and $<\Delta Y/\Delta Z>=2.03\pm0.33$, sensitively reducing the uncertainty in these important parameters., Comment: Accepted for publication in MNRAS

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., Comment: 14 pages, 4 figures, 1 table, Accepted on ApJ Letter

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., Comment: 4 pages, 5 figures. Proceedings of the conference "Lithium in the Universe: to Be or not to Be", 18-22 November 2019, Rome (Italy)

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., Comment: 8 pages, 5 figures, accepted for publication in ApJ Letters

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., Comment: 28 pages, 24 figures. Accepted by A&A

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], Comment: 15 pages, 10 figures, 8 tables (Table 2 available at the CDS)

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., Comment: Accepted for publication in A&A; 23 pages, 17 figures, 3 tables

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., Comment: Contribution in response of the Call for White Papers on LSST Cadence Optimization

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., Comment: Submitted White Paper for the cadence optimization call of LSST, November 30 2018

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

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., Comment: Accepted for publication in A&A; 20 pages, 22 figures, 6 tables

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., Comment: 18 pages, 20 figures, 1 table; accepted for publication in MNRAS

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, Comment: Accepted for publication in AJ (2018), 36 pages, 22 figures

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., Comment: Accepted for publication by A&A. 22 pages main text plus 8 pages appendices, 12 figures main text, 7 figures appendices

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., Comment: 15 pages, 19 figures, accepted for publication in A&A

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., Comment: Accepted for publication on Astronomy and Astrophysics

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., Comment: AJ, accepted

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., Comment: Accepted for publication in MNRAS. 13 pages, 3 tables, 9 figures

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., Comment: Accepted for pubblication in MNRAS. 19 pages, 15 figures, 2 tables

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., Comment: 9 pages, 4 tables, 5 figures. Accepted for publication by Astronomy & Astrophysics

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., Comment: accepted

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., Comment: 8 pages, 6 fugures, full Table 1 online. Accepted for publication in MNRAS on 2013 November 30

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., Comment: 15 pages, 13 figures, accepted by A&A

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., Comment: 13 pages, 10 figures, 3 tables. Accepted for publication in A&A. This version includes referee comments. Planet's name changed from KOI-127b to Kepler-77b

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., Comment: 10 pages, 5 figures. Accepted for publication in A&A

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., Comment: The database is available at http://astro.df.unipi.it/stellar-models/