The initial-final mass relation (IFMR) links the birth mass of a star to the mass of the compact remnant left at its death. While the relevance of the IFMR across astrophysics is universally acknowledged, not all of its fine details have yet been resolved. A new analysis of a few carbon-oxygen white dwarfs in old open clusters of the Milky Way led us to identify a kink in the IFMR, located over a range of initial masses, $1.65 \lesssim M_{\rm i}/M_{\odot} \lesssim 2.10$. The kink's peak in WD mass of $\approx 0.70-0.75 \, M_{\odot}$ is produced by stars with $M_{\rm i} \simeq 1.8 - 1.9 \, M_{\odot}$, corresponding to ages of about $1.8 - 1.7 $ Gyr. Interestingly, this peak coincides with the initial mass limit between low-mass stars that develop a degenerate helium core after central hydrogen exhaustion, and intermediate-mass stars that avoid electron degeneracy. We interpret the IFMR kink as the signature of carbon star formation in the Milky Way. This finding is critical to constraining the evolution and chemical enrichment of low-mass stars, and their impact on the spectrophotometric properties of galaxies.

Among the reactions involved in the production and destruction of deuterium during Big Bang Nucleosynthesis, the deuterium-burning D(p,gamma)3He reaction has the largest uncertainty and limits the precision of theoretical estimates of primordial deuterium abundance. Here we report the results of a careful commissioning of the experimental setup used to measure the cross-section of the D(p,gamma)3He reaction at the Laboratory for Underground Nuclear Astrophysics of the Gran Sasso Laboratory (Italy). The commissioning was aimed at minimising all sources of systematic uncertainty in the measured cross sections. The overall systematic error achieved (< 3 %) will enable improved predictions of BBN deuterium abundance.

We present the \texttt{YBC} database of stellar bolometric corrections (BCs), available at \url{http://stev.oapd.inaf.it/YBC}. We homogenize widely-used theoretical stellar spectral libraries and provide BCs for many popular photometric systems, including the Gaia filters. The database can be easily extended to additional photometric systems and stellar spectral libraries. The web interface allows users to transform their catalogue of theoretical stellar parameters into magnitudes and colours of selected filter sets. The BC tables can be downloaded or also be implemented into large simulation projects using the interpolation code provided with the database. We compute extinction coefficients on a star-by-star basis, hence taking into account the effects of spectral type and non-linearity dependency on the total extinction. We illustrate the use of these BCs in \texttt{PARSEC} isochrones. We show that using spectral-type dependent extinction coefficients is quite necessary for Gaia filters whenever $\Av\gtrsim 0.5$\,mag. BC tables for rotating stars and tables of limb-darkening coefficients are also provided.

We perform a comprehensive study of the stellar populations in the young Large Magellanic Cloud cluster NGC 1866, combining the analysis of its best-studied Cepheids with that of a very accurate colour-magnitude diagram (CMD) obtained from the most recent Hubble Space Telescope photometry. We use a Bayesian method based on new PARSEC stellar evolutionary tracks with overshooting and rotation to obtain ages and initial rotation velocities of five well-studied Cepheids of the cluster. We find that four of the five Cepheids belong to an initially slowly rotating young population (of $ 176 \pm 5$ Myr), while the fifth one is significantly older, either $ 288 \pm 20$ Myr for models with high initial rotational velocity ($\omega_\mathrm{i} \sim 0.9$), or $ 202 \pm 5$ Myr for slowly rotating models. The complementary analysis of the CMD rules out the latter solution while strongly supporting the presence of two distinct populations of $\sim$176 Myr and $\sim$288 Myr, respectively. Moreover, the observed multiple main sequences and the turn-offs indicate that the younger population is mainly made of slowly rotating stars, as is the case of the four younger Cepheids, while the older population is made mainly of initially fast rotating stars, as is the case of the fifth Cepheid. Our study not only reinforces the notion that some young clusters like NGC 1866 harbor multiple populations, but gives also hints that the first population, the older, may inherit the angular momentum from the parent cloud while stars of the second one, the younger, do not.

Context: The period-luminosity diagram (PLD) has proven to be a powerful tool for studying populations of pulsating red giants. Gaia Data Release 2 (DR2) provides a large data set including many long-period variables (LPVs) on which this tool can be applied. Aims: We investigate the location of LPVs from the Large and Small Magellanic Clouds in the PLD using various optical and infrared luminosity indicators from Gaia and 2MASS, respectively. We thereby distinguish between stars of different masses and surface chemistry. Methods: The data set taken from the Gaia DR2 catalogue of LPVs allows for a homogeneous study from low- to high-mass LPVs. These sources are divided into sub-populations of asymptotic giant branch (AGB) stars according to their mass and their O- or C-rich nature using the Gaia-2MASS diagram developed by our group. This diagram uses a Wesenheit index Wrp based on Wesenheit functions in the Gaia and 2MASS photometric bands. Four different luminosity indicators are used to study the period-luminosity (P-L) relations. Results: We provide the first observational evidence of a P-L relation offset for both fundamental and 1O pulsators between low- and intermediate-mass O-rich stars, in agreement with published pulsation predictions. Among the luminosity indicators explored, sequence C' is the narrowest in the P-Wrp diagram, and is thus to be preferred over the other PLDs for the determination of distances using LPVs. The majority of massive asymptotic giant branch (AGB) stars and red supergiants form a smooth extension of sequence C of low- and intermediate-mass AGB stars in the P-Wrp diagram, suggesting that they pulsate in the fundamental mode. All results are similar in the two Magellanic Clouds.

Rapidly rotating stars have been recently recognized as having a major role in the interpretation of colour-magnitude diagrams of young and intermediate-age star clusters in the Magellanic Clouds and in the Milky Way. In this work, we evaluate the distinctive spectra and distributions in colour-colour space that follow from the presence of a substantial range in effective temperatures across the surface of fast rotators. The calculations are inserted in a formalism similar to the one usually adopted for non-rotating stars, which allows us to derive tables of bolometric corrections as a function not only of a reference effective temperature, surface gravity and metallicity, but also of the rotational speed with respect to the break-up value, $\omega$, and the inclination angle, $i$. We find that only very fast rotators ($\omega>0.95$) observed nearly equator-on ($i>45^\circ$) present sizable deviations from the colour-colour relations of non-rotating stars. In light of these results, we discuss the photometry of the $\sim$ 200-Myr-old cluster NGC 1866 and its split main sequence, which has been attributed to the simultaneous presence of slow and fast rotators. The small dispersion of its stars in colour-colour diagrams allow us to conclude that fast rotators in this cluster either have rotational velocities $\omega<0.95$, or are all observed nearly pole-on. Such geometric colour-colour effects, although small, might be potentially detectable in the huge, high-quality photometric samples in the post-Gaia era, in addition to the evolutionary effects caused by rotation-induced mixing.

We analyze resolved stellar near-infrared photometry of 21 HST fields in M31 to constrain the impact of metallicity on the formation of carbon stars. Observations of nearby galaxies show that the carbon stars are increasingly rare at higher metallicity. Models indicate that carbon star formation efficiency drops due to the decrease in dredge-up efficiency in metal-rich thermally-pulsing Asymptotic Giant Branch (TP-AGB) stars, coupled to a higher initial abundance of oxygen. However, while models predict a metallicity ceiling above which carbon stars cannot form, previous observations have not yet pinpointed this limit. Our new observations reliably separate carbon stars from M-type TP-AGB stars across 2.6-13.7 kpc of M31's metal-rich disk using HST WFC3/IR medium-band filters. We find that the ratio of C to M stars (C/M) decreases more rapidly than extrapolations of observations in more metal-poor galaxies, resulting in a C/M that is too low by more than a factor of 10 in the innermost fields and indicating a dramatic decline in C star formation efficiency at metallicities higher than [M/H] $\approx$ -0.1 dex. The metallicity ceiling remains undetected, but must occur at metallicities higher than what is measured in M31's inner disk ([M/H] $\gtrsim$ +0.06 dex).

We have computed a grid of hydrostatic spherical COMARCS models for C stars covering metallicities from [Z/H]=0 to -2 and values of the carbon excess [C-O] from 6.41 to 9.15, plus some temperature sequences, where the amount of oxygen and nitrogen is increased relative to a scaled solar element mixture. Such abundance variations may appear during the late stages of stellar evolution. Our study covers changes of [O/Z] and [N/Z] going up to +0.5. Based on the atmospheric structures we have calculated synthetic spectra and photometry for all of the models in a consistent way. The sequences with changed [O/Z] and [N/Z] can be used to derive correction terms, which are applied to the colours predicted for a certain combination of effective temperature, surface gravity, metallicity and carbon excess. If one neglects these shifts in case of a variable oxygen amount, taking [C-O] instead of C/O gives much better results, since the first quantity dominates the formation of many important molecular species. For the warmer C giants with weaker pulsation it is in principle possible to determine [C-O], [O/Z] or [N/Z] from high resolution spectra, when the opacities in the radiative transfer calculations for the models and observable properties are treated consistently. The corresponding changes due to the abundances often become significantly larger than the deviations caused by uncertainties of the stellar parameters or by a optically thin dust shell. Photometric data and low or medium resolution spectra are not sufficient to derive the mentioned quantities.

The purpose of this work is to present an extensive grid of dynamical atmosphere and wind models for M-type AGB stars, covering a wide range of relevant stellar parameters. We used the DARWIN code, which includes frequency-dependent radiation-hydrodynamics and a time-dependent description of dust condensation and evaporation, to simulate the dynamical atmosphere. The wind-driving mechanism is photon scattering on submicron-sized Mg$_2$SiO$_4$ grains. The grid consists of $\sim4000$ models, with luminosities from $L_\star=890\,{\mathrm{L}}_\odot$ to $L_\star=40000\,{\mathrm{L}}_\odot$ and effective temperatures from 2200K to 3400K. For the first time different current stellar masses are explored with M-type DARWIN models, ranging from 0.75M$_\odot$ to 3M$_\odot$. The modelling results are radial atmospheric structures, dynamical properties such as mass-loss rates and wind velocities, and dust properties (e.g. grain sizes, dust-to-gas ratios, and degree of condensed Si). We find that the mass-loss rates of the models correlate strongly with luminosity. They also correlate with the ratio $L_*/M_*$: increasing $L_*/M_*$ by an order of magnitude increases the mass-loss rates by about three orders of magnitude, which may naturally create a superwind regime in evolution models. There is, however, no discernible trend of mass-loss rate with effective temperature, in contrast to what is found for C-type AGB stars. We also find that the mass-loss rates level off at luminosities higher than $\sim14000\,{\mathrm{L}}_\odot$, and consequently at pulsation periods longer than $\sim800$ days. The final grain radii range from 0.25 micron to 0.6 micron. The amount of condensed Si is typically between 10% and 40%, with gas-to-dust mass ratios between 500 and 4000.

Models of the chemical evolution of the interstellar medium, galaxies, and the Universe rely on our understanding of the amounts and chemical composition of the material returned by stars and supernovae. Stellar yields are obtained from stellar-evolution models, which currently lack predictive prescriptions of stellar mass loss, although it significantly affects stellar lifetimes, nucleosynthesis, and chemical ejecta. Galaxy properties are derived from observations of the integrated light of bright member stars. Stars in the late stages of their evolution are among the infrared-brightest objects in galaxies. An unrealistic treatment of the mass-loss process introduces significant uncertainties in galaxy properties derived from their integrated light. We describe current efforts and future needs and opportunities to characterize AGB outflows: driving mechanisms, outflow rates, underlying fundamental physical and chemical processes such as dust grain formation, and dependency of these on metallicity.

The thermally-pulsing asymptotic giant branch (TP-AGB) experienced by low- and intermediate-mass stars is one of the most uncertain phases of stellar evolution and the models need to be calibrated with the aid of observations. To this purpose, we couple high-quality observations of resolved stars in the Small Magellanic Cloud (SMC) with detailed stellar population synthesis simulations computed with the TRILEGAL code. The strength of our approach relies on the detailed spatially-resolved star formation history of the SMC, derived from the deep near-infrared photometry of the VISTA survey of the Magellanic Clouds, as well as on the capability to quickly and accurately explore a wide variety of parameters and effects with the COLIBRI code for the TP-AGB evolution. Adopting a well-characterized set of observations -- star counts and luminosity functions -- we set up a calibration cycle along which we iteratively change a few key parameters of the TP-AGB models until we eventually reach a good fit to the observations. Our work leads to identify two best-fitting models that mainly differ in the efficiencies of the third dredge-up and mass loss in TP-AGB stars with initial masses larger than about 3 M$_{\odot}$. On the basis of these calibrated models we provide a full characterization of the TP-AGB stellar population in the SMC in terms of stellar parameters (initial masses, C/O ratios, carbon excess, mass-loss rates). Extensive tables of isochrones including these improved models are publicly available.

Double-line eclipsing binaries (DLEBs) have been recently used to constrain the amount of central mixing as a function of stellar mass, with contrasting results. In this work, we reanalyze the DLEB sample by Claret & Torres, using a Bayesian method and new PARSEC tracks that account for both convective core overshooting and rotational mixing. Using overshooting alone we obtain that, for masses larger than about 1.9 M$_{\odot}$, the distribution of the overshooting parameter, $\lambda_\mathrm{ov}$, has a wide dispersion between 0.3 and 0.8, with essentially no values below $\lambda_\mathrm{ov}=$ 0.3 - 0.4. While the lower limit supports a mild convective overshooting efficiency, the large dispersion derived is difficult to explain in the framework of current models of that process, which leave little room for large randomness. We suggest that a simple interpretation of our results can be rotational mixing: different initial rotational velocities, in addition to a fixed amount of overshooting, could reproduce the high dispersion derived for intermediate-mass stars. After a reanalysis of the data, we find good agreement with models computed with fixed overshooting parameter, $\lambda_\mathrm{ov}=0.4$, and initial rotational rates, $\omega$, uniformly distributed in a wide range between $0$ and $0.8$ times the break-up value, at varying initial mass. We also find that our best-fitting models for the components of $\alpha$ Aurigae and TZ Fornacis, agree with their observed rotational velocities, thus providing independent support to our hypothesis. We conclude that a constant efficiency of overshooting in concurrence with a star-to-star variation in the rotational mixing, might be crucial in the interpretation of such data.

The atmospheres and winds of C-type AGB stars are modeled with the 1D spherically symmetric radiation-hydrodynamical code DARWIN. To explore the metallicity-dependence of mass loss we calculate model grids at three different chemical abundances. Since carbon may be dredged up during the thermal pulses as AGB stars evolve, we keep the carbon abundance as a free parameter. The models in these three different grids all have a current mass of one solar mass; effective temperatures of 2600K, 2800K, 3000K, or 3200K; and stellar luminosities equal to log(L/Lsun)=3.70, 3.85, or 4.00. The models show that mass loss in carbon stars is facilitated by high luminosities, low effective temperatures, and a high carbon excess (C-O) at both solar and subsolar metallicities. Similar combinations of effective temperature, luminosity, and carbon excess produce outflows at both solar and subsolar metallicities. There are no large systematic differences in the mass-loss rates and wind velocities produced by these wind models with respect to metallicity, nor any systematic difference concerning the distribution of grain sizes or how much carbon is condensed into dust. DARWIN models at subsolar metallicity have approximately 15% lower mass-loss rates compared to DARWIN models at solar metallicity with the same stellar parameters and carbon excess. For both solar and subsolar environments typical grain sizes range between 0.1 and 0.5 micron, the degree of condensed carbon varies between 5% and 40%, and the gas-to-dust ratios between 500 and 10000. Furthermore, stellar evolution models can use the mass-loss rates calculated from DARWIN models at solar metallicity when modeling the AGB phase at subsolar metallicities if carbon excess is used as the critical abundance parameter instead of the C/O ratio.

We present a new grid of non-adiabatic, linear pulsation models of Long-Period Variables (LPVs), including periods and growth rates for radial modes from the fundamental to the fourth overtone. The models span a wide range in mass, luminosity, metallicity, C/O ratio and helium abundance, effectively covering the whole thermally-pulsing asymptotic giant branch (TP-AGB) evolution, and representing a significant update with respect to previous works. The main improvement is the inclusion of detailed atomic and molecular opacities, consistent with the models chemical mixture, that makes the present set of models the first to systematically account for variability in C-stars. We examine periods and growth rates in the models, and find that, while the fundamental mode is affected by the structure of the envelope, overtones are less sensitive to the interior and largely determined by the global properties. In the models, the frequency of the overtone with the largest degree of excitation is found to scale with the acoustic cut-off frequency at the stellar surface, a behaviour similar to that observed for the frequency of maximum oscillation power for solar-like oscillations in less evolved red giants. This allows us to provide a simple analytic prescription to predict the most-likely dominant mode as a function of stellar parameters. Best-fit relations for periods are also provided. By applying results of pulsation models to evolutionary tracks, we present a general picture of the evolution of long-period variability during the TP-AGB, that we find consistent with observations. Models are made public through a dedicated web interface.

Aims: We explore the wealth of high quality photometric data provided by data release 2 of the Gaia mission for long period variables (LPVs) in the Large Magellanic Cloud. Our goal is to identify stars of various types and masses along the Asymptotic Giant Branch. Methods: For this endeavour, we developed a new multi-band approach combining Wesenheit functions W_{RP,BP-RP} and W_{K_s,J-K_s} in the Gaia BP, RP and 2MASS J, K_s spectral ranges, respectively, and use a new diagram (W_{RP,BP-RP}-W_{K_s,J-K_s}) versus K_s to distinguish between different kinds of stars in our sample of LPVs. We used stellar population synthesis models to validate our approach. Results:We demonstrate the ability of the new diagram to discriminate between O-rich and C-rich objects, and to identify low-mass, intermediate-mass and massive O-rich red giants, as well as extreme C-rich stars. Stellar evolution and population synthesis models guide the interpretation of the results, highlighting the diagnostic power of the new tool to discriminate between stellar initial masses, chemical properties and evolutionary stages.

We recover the spatially resolved star formation history across the entire main body and Wing of the Small Magellanic Cloud (SMC), using fourteen deep tile images from the VISTA survey of the Magellanic Clouds (VMC), in the YJKs filters. The analysis is performed on 168 subregions of size 0.143 deg2, covering a total contiguous area of 23.57 deg2. We apply a colour-magnitude diagram (CMD) reconstruction method that returns the best-fitting star formation rate SFR(t), age--metallicity relation, distance and mean reddening, together with their confidence intervals, for each subregion. With respect to previous analyses, we use a far larger set of VMC data, updated stellar models, and fit the two available CMDs (Y-Ks versus Ks and J-Ks versus Ks) independently. The results allow us to derive a more complete and more reliable picture of how the mean distances, extinction values, star formation rate, and metallicities vary across the SMC, and provide a better description of the populations that form its Bar and Wing. We conclude that the SMC has formed a total mass of (5.31+-0.05)x10^8 Msun in stars over its lifetime. About two thirds of this mass is expected to be still locked in stars and stellar remnants. 50 per cent of the mass was formed prior to an age of 6.3 Gyr, and 80 per cent was formed between 8 and 3.5 Gyr ago. We also illustrate the likely distribution of stellar ages and metallicities in different parts of the CMD, to aid the interpretation of data from future astrometric and spectroscopic surveys of the SMC.

The experimental study of nuclear reactions of astrophysical interest is greatly facilitated by a low-background, high-luminosity setup. The Laboratory for Underground Nuclear Astrophysics (LUNA) 400 kV accelerator offers ultra-low cosmic-ray induced background due to its location deep underground in the Gran Sasso National Laboratory (INFN-LNGS), Italy, and high intensity, 250-500 $\mu$A, proton and $\alpha$ ion beams. In order to fully exploit these features, a high-purity, recirculating gas target system for isotopically enriched gases is coupled to a high-efficiency, six-fold optically segmented bismuth germanate (BGO) $\gamma$-ray detector. The beam intensity is measured with a beam calorimeter with constant temperature gradient. Pressure and temperature measurements have been carried out at several positions along the beam path, and the resultant gas density profile has been determined. Calibrated $\gamma$-intensity standards and the well-known $E_p$ = 278 keV $\mathrm{^{14}N(p,\gamma)^{15}O}$ resonance were used to determine the $\gamma$-ray detection efficiency and to validate the simulation of the target and detector setup. As an example, the recently measured resonance at $E_p$ = 189.5 keV in the $^{22}$Ne(p,$\gamma$)$^{23}$Na reaction has been investigated with high statistics, and the $\gamma$-decay branching ratios of the resonance have been determined.

Precise studies on the Galactic bulge, globular cluster, Galactic halo and Galactic thick disk require stellar models with alpha enhancement and various values of helium content. These models are also important for extra-Galactic population synthesis studies. For this purpose, we complement the existing PARSEC models, which are based on the solar partition of heavy elements, with alpha-enhanced partitions. We collect detailed measurements on the metal mixture and helium abundance for the two populations of 47 Tuc (NGC 104) from the literature, and calculate stellar tracks and isochrones with these alpha-enhanced compositions. By fitting the precise color-magnitude diagram with HST ACS/WFC data, from low main sequence till horizontal branch, we calibrate some free parameters that are important for the evolution of low mass stars like the mixing at the bottom of the convective envelope. This new calibration significantly improves the prediction of the RGB bump brightness. Comparison with the observed RGB and HB luminosity functions also shows that the evolutionary lifetimes are correctly predicted. As a further result of this calibration process, we derive the age, distance modulus, reddening, and the red giant branch mass loss for 47 Tuc. We apply the new calibration and alpha-enhanced mixtures of the two 47 Tuc populations ( [alpha/Fe] ~0.4 and 0.2) to other metallicities. The new models reproduce the RGB bump observations much better than previous models. This new PARSEC database, with the newly updated alpha-enhanced stellar evolutionary tracks and isochrones, will also be part of the new stellar products for Gaia.

We employ newly computed grids of spectra reprocessed by dust for estimating the total dust production rate (DPR) of carbon stars in the Small Magellanic Cloud (SMC). For the first time, the grids of spectra are computed as a function of the main stellar parameters, i.e. mass-loss rate, luminosity, effective temperature, current stellar mass and element abundances at the photosphere, following a consistent, physically grounded scheme of dust growth coupled with stationary wind outflow. The model accounts for the dust growth of various dust species formed in the circumstellar envelopes of carbon stars, such as carbon dust, silicon carbide and metallic iron. In particular, we employ some selected combinations of optical constants and grain sizes for carbon dust which have been shown to reproduce simultaneously the most relevant color-color diagrams in the SMC. By employing our grids of models, we fit the spectral energy distributions of $\approx$3100 carbon stars in the SMC, consistently deriving some important dust and stellar properties, i.e. luminosities, mass-loss rates, gas-to-dust ratios, expansion velocities and dust chemistry. We discuss these properties and we compare some of them with observations in the Galaxy and LMC. We compute the DPR of carbon stars in the SMC, finding that the estimates provided by our method can be significantly different, between a factor $\approx2-5$, than the ones available in the literature. Our grids of models, including the spectra and other relevant dust and stellar quantities, are publicly available at http://starkey.astro.unipd.it/web/guest/dustymodels

We report and discuss JHKs photometry for Sgr dIG, a very metal-deficient galaxy in the Local Group, obtained over 3.5 years with the Infrared Survey Facility in South Africa. Three large amplitude asymptotic giant branch variables are identified. One is an oxygen-rich star that has a pulsation period of 950 days, that was until recently undergoing hot bottom burning, with Mbol~-6.7. It is surprising to find a variable of this sort in Sgr dIG, given their rarity in other dwarf irregulars. Despite its long period the star is relatively blue and is fainter, at all wavelengths shorter than 4.5microns, than anticipated from period-luminosity relations that describe hot bottom burning stars. A comparison with models suggests it had a main sequence mass Mi~5 times solar and that it is now near the end of its AGB evolution. The other two periodic variables are carbon stars with periods of 670 and 503 days (Mbol~-5.7 and -5.3). They are very similar to other such stars found on the AGB of metal deficient Local Group Galaxies and a comparison with models suggests Mi~3 times solar. We compare the number of AGB variables in Sgr dIG to those in NGC6822 and IC1613, and suggest that the differences may be due to the high specific star formation rate and low metallicity of Sgr dIG.

Period-luminosity (PL) sequences of long period variables (LPVs) are commonly interpreted as different pulsation modes, but there is disagreement on the modal assignment. Here, we re-examine the observed PL sequences in the Large Magellanic Cloud, including the sequence of long secondary periods (LSPs), and their associated pulsation modes. Firstly, we theoretically model the sequences using linear, radial, non-adiabatic pulsation models and a population synthesis model of the LMC red giants. Then, we use a semi-empirical approach to assign modes to the pulsation sequences by exploiting observed multi-mode pulsators. As a result of the combined approaches, we consistently find that sequences B and C$^{\prime}$ both correspond to first overtone pulsation, although there are some fundamental mode pulsators at low luminosities on both sequences. The masses of these fundamental mode pulsators are larger at a given luminosity than the mass of the first overtone pulsators. These two sequences B and C$^{\prime}$ are separated by a small period interval in which large amplitude pulsation in a long secondary period (sequence D variability) occurs, meaning that the first overtone pulsation is not seen as the primary mode of pulsation. Observationally, this leads to the splitting of the first overtone pulsation sequence into the two observed sequences B and C$^{\prime}$. Our two independent examinations also show that sequences A$^{\prime}$, A and C correspond to third overtone, second overtone and fundamental mode pulsation, respectively.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

By means of the updated PARSEC database of evolutionary tracks of massive stars, we compute the integrated stellar light, the ionizing photon budget and the supernova rates of young simple stellar populations (SSPs), for different metallicities and IMF upper mass limits. Using CLOUDY we compute and include in the SSP spectra the neb- ular emission contribution. We also revisit the thermal and non-thermal radio emission contribution from young stars. Using GRASIL we can thus predict the panchromatic spectrum and the main recombination lines of any type of star-forming galaxy, including the effects of dust absorption and re-emission. We check the new models against the spectral energy distributions (SEDs) of selected well-observed nearby galaxies. From the best-fit models we obtain a consistent set of star formation rate (SFR) calibrations at wavelengths ranging from ultraviolet (UV) to radio. We also provide analytical calibrations that take into account the dependence on metallcity and IMF upper mass limit of the SSPs. We show that the latter limit can be well constrained by combining information from the observed far infrared, 24 {\mu}m, 33 GHz and H{\alpha} luminosities. Another interesting property derived from the fits is that, while in a normal galaxy the attenuation in the lines is significantly higher than that in the nearby continuum, in individual star bursting regions they are similar, supporting the notion that this effect is due to an age selective extinction. Since in these conditions the Balmer decrement method may not be accurate, we provide relations to estimate the attenuation from the observed 24 {\mu}m or 33 GHz fluxes. These relations can be useful for the analysis of young high redshift galaxies.

We introduce a new generation of PARSEC-COLIBRI stellar isochrones that include a detailed treatment of the thermally-pulsing asymptotic giant branch (TP-AGB) phase, and covering a wide range of initial metallicities (0.0001<Zi<0.06). Compared to previous releases, the main novelties and improvements are: use of new TP-AGB tracks and related atmosphere models and spectra for M and C-type stars; inclusion of the surface H+He+CNO abundances in the isochrone tables, accounting for the effects of diffusion, dredge-up episodes and hot-bottom burning; inclusion of complete thermal pulse cycles, with a complete description of the in-cycle changes in the stellar parameters; new pulsation models to describe the long-period variability in the fundamental and first overtone modes; new dust models that follow the growth of the grains during the AGB evolution, in combination with radiative transfer calculations for the reprocessing of the photospheric emission. Overall, these improvements are expected to lead to a more consistent and detailed description of properties of TP-AGB stars expected in resolved stellar populations, especially in regard to their mean photometric properties from optical to mid-infrared wavelengths. We illustrate the expected numbers of TP-AGB stars of different types in stellar populations covering a wide range of ages and initial metallicities, providing further details on the C-star island that appears at intermediate values of age and metallicity, and about the AGB-boosting effect that occurs at ages close to 1.6 Gyr for populations of all metallicities. The isochrones are available through a new dedicated web server.

Asteroseismic parameters allow us to measure the basic stellar properties of field giants observed far across the Galaxy. Most of such determinations are, up to now, based on simple scaling relations involving the large frequency separation, \Delta\nu, and the frequency of maximum power, \nu$_{max}$. In this work, we implement \Delta\nu\ and the period spacing, {\Delta}P, computed along detailed grids of stellar evolutionary tracks, into stellar isochrones and hence in a Bayesian method of parameter estimation. Tests with synthetic data reveal that masses and ages can be determined with typical precision of 5 and 19 per cent, respectively, provided precise seismic parameters are available. Adding independent information on the stellar luminosity, these values can decrease down to 3 and 10 per cent respectively. The application of these methods to NGC 6819 giants produces a mean age in agreement with those derived from isochrone fitting, and no evidence of systematic differences between RGB and RC stars. The age dispersion of NGC 6819 stars, however, is larger than expected, with at least part of the spread ascribable to stars that underwent mass-transfer events.

Modeling stellar atmospheres is a complex, intriguing task in modern astronomy. A systematic comparison of models with multi-technique observations is the only efficient way to constrain them. Aims: We performed a self-consistent modeling of the atmospheres of six C-rich AGB stars: R Lep, R Vol, Y Pav, AQ Sgr, U Hya and X TrA, with the aim of enlarging the knowledge of the dynamic processes occurring in their atmospheres. Methods. We used VLTI/MIDI interferometric observations, in combination with spectro-photometric data, and compared them with self-consistent dynamic models atmospheres. Results: The models can reproduce SED data well at wavelengths longwards of 1{\mu}m, and the interferometric observations between 8{\mu}m and 10{\mu}m. Discrepancies could be due to a combination of data- and model-related effects. The models best fitting the Miras are significantly extended, with a prominent shell-like structure, while the models best fitting the non-Miras are more compact, showing lower average mass-loss. The mass loss is of episodic or multi-periodic nature, but causes the visual amplitudes to be notably larger than the observed ones. Stellar parameters were derived from the model fitting: T_eff, L_bol, M, C/O, $\dot{\text{M}}$. Our findings agree well with literature values within the uncertainties. T_eff and L_bol are also in good agreement with the T derived from the angular diameter {\theta}_{(V-K)} and L_bol from the SED fitting, except for AQ Sgr. Finally, {\theta}_Ross and {\theta}_(V-K) agree with each other better for the Miras than for the non-Miras,which is probably connected to the episodic nature of the latter models. We also located the stars in the H-R diagram, comparing them with evolutionary tracks. We found that the main derived properties from the model fitting are in good agreement with TP-AGB evolutionary calculations for carbon stars (COLIBRI code).

We investigate the impact of the new LUNA rate for the nuclear reaction $^{22}$Ne$(p,\gamma)^{23}$Na on the chemical ejecta of intermediate-mass stars, with particular focus on the thermally-pulsing asymptotic giant branch (TP-AGB) stars that experience hot-bottom burning. To this aim we use the PARSEC and COLIBRI codes to compute the complete evolution, from the pre-main sequence up to the termination of the TP-AGB phase, of a set of stellar models with initial masses in the range $3.0\,M_{\odot} - 6.0\,M_{\odot}$, and metallicities $Z_{\rm i}=0.0005$, $Z_{\rm i}=0.006$, and $Z_{\rm i} = 0.014$. We find that the new LUNA measures have much reduced the nuclear uncertainties of the $^{22}$Ne and $^{23}$Na AGB ejecta, which drop from factors of $\simeq 10$ to only a factor of few for the lowest metallicity models. Relying on the most recent estimations for the destruction rate of $^{23}$Na, the uncertainties that still affect the $^{22}$Ne and $^{23}$Na AGB ejecta are mainly dominated by evolutionary aspects (efficiency of mass-loss, third dredge-up, convection). Finally, we discuss how the LUNA results impact on the hypothesis that invokes massive AGB stars as the main agents of the observed O-Na anti-correlation in Galactic globular clusters. We derive quantitative indications on the efficiencies of key physical processes (mass loss, third dredge-up, sodium destruction) in order to simultaneously reproduce both the Na-rich, O-poor extreme of the anti-correlation, and the observational constraints on the CNO abundance. Results for the corresponding chemical ejecta are made publicly available.

Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud. Methods: ALMA was used to observe 2 OH/IR stars and 4 carbon stars in the LMC in the CO J= 2-1 line. Results: We present the first measurement of expansion velocities in extragalactic carbon stars. All four C-stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. Gas-to-dust ratios are derived that make the predicted velocities agree with the observed ones. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with "extreme" C-stars that have properties more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than derived from the dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions: Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear similar. This is in agreement with dynamical dust-driven wind models.

We present a new approach aimed at constraining the typical size and optical properties of carbon dust grains in Circumstellar envelopes (CSEs) of carbon-rich stars (C-stars) in the Small Magellanic Cloud (SMC). To achieve this goal, we apply our recent dust growth description, coupled with a radiative transfer code to the CSEs of C-stars evolving along the TP-AGB, for which we compute spectra and colors. Then we compare our modeled colors in the near- and mid-infrared (NIR and MIR) bands with the observed ones, testing different assumptions in our dust scheme and employing several data sets of optical constants for carbon dust available in the literature. Different assumptions adopted in our dust scheme change the typical size of the carbon grains produced. We constrain carbon dust properties by selecting the combination of grain size and optical constants which best reproduces several colors in the NIR and MIR at the same time. The different choices of optical properties and grain size lead to differences in the NIR and MIR colors greater than two magnitudes in some cases. We conclude that the complete set of observed NIR and MIR colors are best reproduced by small grains, with sizes between $\sim$0.035 and $\sim$0.12~$\mu$m, rather than by large grains between $\sim0.2$ and $0.7$~$\mu$m. The inability of large grains to reproduce NIR and MIR colors seems independent of the adopted optical data set. We also find a possible trend of the grain size with mass-loss and/or carbon excess in the CSEs of these stars.

Thermally-Pulsing Asymptotic Giant Branch (TP-AGB) stars are relatively short lived (less than a few Myr), yet their cool effective temperatures, high luminosities, efficient mass-loss and dust production can dramatically effect the chemical enrichment histories and the spectral energy distributions of their host galaxies. The ability to accurately model TP-AGB stars is critical to the interpretation of the integrated light of distant galaxies, especially in redder wavelengths. We continue previous efforts to constrain the evolution and lifetimes of TP-AGB stars by modeling their underlying stellar populations. Using Hubble Space Telescope (HST) optical and near-infrared photometry taken of 12 fields of 10 nearby galaxies imaged via the ACS Nearby Galaxy Survey Treasury and the near-infrared HST/SNAP follow-up campaign, we compare the model and observed TP-AGB luminosity functions as well as the number ratio of TP-AGB to red giant branch stars. We confirm the best-fitting mass-loss prescription, introduced by Rosenfield et al. 2014, in which two different wind regimes are active during the TP-AGB, significantly improves models of many galaxies that show evidence of recent star formation. This study extends previous efforts to constrain TP-AGB lifetimes to metallicities ranging -1.59 < [Fe/H] < -0.56 and initial TP-AGB masses up to ~ 4 Msun, which include TP-AGB stars that undergo hot-bottom burning.

Based on a grid of hydrostatic spherical COMARCS models for cool stars we have calculated observable properties of these objects, which will be mainly used in combination with stellar evolution tracks and population synthesis tools. The high resolution opacity sampling and low resolution convolved spectra as well as bolometric corrections for a large number of filter systems are made electronically available. We exploit those data to study the effect of mass, C/O ratio and nitrogen abundance on the photometry of K and M giants. Depending on effective temperature, surface gravity and the chosen wavelength ranges variations of the investigated parameters cause very weak to moderate and, in the case of C/O values close to one, even strong shifts of the colours. For the usage with stellar evolution calculations they will be treated as correction factors applied to the results of an interpolation in the main quantities. When we compare the synthetic photometry to observed relations and to data from the Galactic Bulge, we find in general a good agreement. Deviations appear for the coolest giants showing pulsations, mass loss and dust shells, which cannot be described by hydrostatic models.

We check the performance of the {\sl\,PARSEC} tracks in reproducing the blue loops of intermediate age and young stellar populations at very low metallicity. We compute new evolutionary {\sl\,PARSEC} tracks of intermediate- and high-mass stars from 2\Msun to 350\Msun with enhanced envelope overshooting (EO), EO=2\HP and 4\HP, for very low metallicity, Z=0.0005. The input physics, including the mass-loss rate, has been described in {\sl\,PARSEC}~V1.2 version. By comparing the synthetic color-magnitude diagrams (CMDs) obtained from the different sets of models with envelope overshooting EO=0.7\HP (the standard {\sl\,PARSEC} tracks), 2\HP and 4\HP, with deep observations of the Sagittarius dwarf irregular galaxy (SagDIG), we find an overshooting scale EO=2\HP to best reproduce the observed loops. This result is consistent with that obtained by \citet{Tang_etal14} for Z in the range 0.001-0.004. We also discuss the dependence of the blue loop extension on the adopted instability criterion and find that, contrary to what stated in literature, the Schwarzschild criterion, instead of the Ledoux criterion, favours the development of blue loops. Other factors that could affect the CMD comparisons such as differential internal extinction or the presence of binary systems are found to have negligible effects on the results. We thus confirm that, in presence of core overshooting during the H-burning phase, a large envelope overshooting is needed to reproduce the main features of the central He-burning phase of intermediate- and high-mass stars.

We investigate the formation of hydrogen cyanide (HCN) in the inner circumstellar envelopes of thermally pulsing asymptotic giant branch (TP-AGB) stars. A dynamic model for periodically shocked atmospheres, which includes an extended chemo-kinetic network, is for the first time coupled to detailed evolutionary tracks for the TP-AGB phase computed with the COLIBRI code. We carried out a calibration of the main shock parameters (the shock formation radius and the effective adiabatic index) using the circumstellar HCN abundances recently measured for a populous sample of pulsating TP-AGB stars. Our models recover the range of the observed HCN concentrations as a function of the mass-loss rates, and successfully reproduce the systematic increase of HCN moving along the M-S-C chemical sequence of TP-AGB stars, that traces the increase of the surface C/O ratio. The chemical calibration brings along two important implications: i) the first shock should emerge very close to the photosphere, and ii) shocks are expected to have a dominant isothermal character in the denser region close to the star (within ~ 3-4 R), implying that radiative processes should be quite efficient. Our analysis also suggests that the HCN concentrations in the inner circumstellar envelopes are critically affected by the H-H2 chemistry during the post-shock relaxation stages.

We complement the PARSEC data base of stellar evolutionary tracks with new models of massive stars, from the pre-main sequence phase to the central carbon ignition. We consider a broad range of metallicities, 0.0001$\leq Z \leq$0.04 and initial masses up to $M_{\rm ini}=350\,M_\odot$. The main difference with respect to our previous models of massive stars is the adoption of a recent formalism accounting for the mass-loss enhancement when the ratio of the stellar to the Eddington luminosity, $\Gamma_e$, approaches unity. With this new formalism, the models are able to reproduce the Humphreys-Davidson limit observed in the Galactic and Large Magellanic Cloud colour-magnitude diagrams, without an ad hoc mass-loss enhancement. We also follow the predictions of recent wind models indicating that the metallicity dependence of the mass-loss rates becomes shallower when $\Gamma_e$ approaches unity. We thus find that the more massive stars may suffer from substantial mass-loss even at low metallicity. We also predict that the Humphreys-Davidson limit should become brighter at decreasing metallicity. We supplement the evolutionary tracks with new tables of theoretical bolometric corrections, useful to compare tracks and isochrones with the observations. For this purpose, we homogenize existing stellar atmosphere libraries of hot and cool stars (PoWR, ATLAS9 and Phoenix) and we add, where needed, new atmosphere models computed with WM-basic. The mass, age and metallicity grids are fully adequate to perform detailed investigations of the properties of very young stellar systems, both in local and distant galaxies. The new tracks supersede the previous old Padova models of massive stars.

By examining the diffusion of young white dwarfs through the core of the globular cluster 47 Tucanae, we estimate the time when the progenitor star lost the bulk of its mass to become a white dwarf. According to stellar evolution models of the white-dwarf progenitors in 47 Tucanae, we find this epoch to coincide approximately with the star ascending the asymptotic giant branch ($3.0 \pm 8.1$ Myr before the tip of the AGB) and more than ninety million years after the helium flash (with ninety-percent confidence). From the diffusion of the young white dwarfs we can exclude the hypothesis that the bulk of the mass loss occurs on the red-giant branch at the four-sigma level. Furthermore, we find that the radial distribution of horizontal branch stars is consistent with that of the red-giant stars and upper-main-sequence stars and inconsistent with the loss of more than 0.2 solar masses on the red-giant branch at the six-sigma level.

Low-mass stars in the He-core-burning phase (HeCB) play a major role in stellar, galactic, and extragalactic astrophysics. The ability to predict accurately the properties of these stars, however, depends on our understanding of convection, which remains one of the key open questions in stellar modelling. We argue that the combination of the luminosity of the AGB bump (AGBb) and the period spacing of gravity modes (DP) during the HeCB phase, provides us with a decisive test to discriminate between competing models of these stars. We use the MESA, BaSTI, and PARSEC stellar evolution codes to model a typical giant star observed by Kepler. We explore how various near-core-mixing scenarios affect the predictions of the above-mentioned constraints, and we find that DP depends strongly on the prescription adopted. Moreover we show that the detailed behaviour of DP shows the signature of sharp variations in the Brunt-Vaisala frequency, which could potentially give additional information about near-core features. We find evidence for the AGBb among Kepler targets, and a first comparison with observations shows that, even if standard models are able to reproduce the luminosity distribution, no standard model can account for satisfactorily the period spacing of HeCB stars. Our analysis allows us to outline a candidate model to describe simultaneously the two observed distributions: a model with a moderate overshooting region characterized by an adiabatic thermal stratification. This prescription will be tested in the future on cluster stars, to limit possible observational biases.

Lithium abundance derived in metal-poor main sequence stars is about three times lower than the value of primordial Li predicted by the standard Big Bang nucleosynthesis when the baryon density is taken from the CMB or the deuterium measurements. This disagreement is generally referred as the lithium problem. We here reconsider the stellar Li evolution from the pre-main sequence to the end of the main sequence phase by introducing the effects of convective overshooting and residual mass accretion. We show that $^7$Li could be significantly depleted by convective overshooting in the pre-main sequence phase and then partially restored in the stellar atmosphere by a tail of matter accretion which follows the Li depletion phase and that could be regulated by EUV photo-evaporation. By considering the conventional nuclear burning and microscopic diffusion along the main sequence we can reproduce the Spite plateau for stars with initial mass $m_0=0.62 - 0.80 M_{\odot}$, and the Li declining branch for lower mass dwarfs, e.g, $m_0=0.57 - 0.60 M_{\odot}$, for a wide range of metallicities (Z=0.00001 to Z=0.0005), starting from an initial Li abundance $A({\rm Li}) =2.72$. This environmental Li evolution model also offers the possibility to interpret the decrease of Li abundance in extremely metal-poor stars, the Li disparities in spectroscopic binaries and the low Li abundance in planet hosting stars.

Recent high-quality photometry of many star clusters in the Magellanic Clouds with ages of 1$\,-\,$2 Gyr revealed main sequence turnoffs (MSTOs) that are significantly wider than can be accounted for by a simple stellar population (SSP). Such extended MSTOs (eMSTOs) are often interpreted in terms of an age spread of several $10^8$ yr, challenging the traditional view of star clusters as being formed in a single star formation episode. Li et al. and Bastian & Niederhofer recently investigated the sub-giant branches (SGBs) of NGC 1651, NGC 1806, and NGC 1846, three star clusters in the Large Magellanic Cloud (LMC) that exhibit an eMSTO. They argued that the SGB of these star clusters can be explained only by a SSP. We study these and two other similar star clusters in the LMC, using extensive simulations of SSPs including unresolved binaries. We find that the shapes of the cross-SGB profiles of all star clusters in our sample are in fact consistent with their cross-MSTO profiles when the latter are interpreted as age distributions. Conversely, SGB morphologies of star clusters with eMSTOs are found to be inconsistent with those of simulated SSPs. Finally, we create PARSEC isochrones from tracks featuring a grid of convective overshoot levels and a very fine grid of stellar masses. A comparison of the observed photometry with these isochrones shows that the morphology of the red clump (RC) of such star clusters is also consistent with that implied by their MSTO in the age spread scenario. We conclude that the SGB and RC morphologies of star clusters featuring eMSTOs are consistent with the scenario in which the eMSTOs are caused by a distribution of stellar ages.

The efficiency of convection in stars affects many aspects of their evolution and remains one of the key-open questions in stellar modelling. In particular, the size of the mixed core in core-He-burning low-mass stars is still uncertain and impacts the lifetime of this evolutionary phase and, e.g., the C/O profile in white dwarfs. One of the known observables related to the Horizontal Branch (HB) and Asymptotic Giant Branch (AGB) evolution is the AGB bump. Its luminosity depends on the position in mass of the helium-burning shell at its first ignition, that is affected by the extension of the central mixed region. In this preliminary work we show how various assumptions on near-core mixing and on the thermal stratification in the overshooting region affect the luminosity of the AGB bump, as well as the period spacing of gravity modes in core-He-burning models.

We analyse deep images from the VISTA survey of the Magellanic Clouds in the YJKs filters, covering 14 sqrdeg (10 tiles), split into 120 subregions, and comprising the main body and Wing of the Small Magellanic Cloud (SMC). We apply a colour--magnitude diagram reconstruction method that returns their best-fitting star formation rate SFR(t), age-metallicity relation (AMR), distance and mean reddening, together with 68% confidence intervals. The distance data can be approximated by a plane tilted in the East-West direction with a mean inclination of 39 deg, although deviations of up to 3 kpc suggest a distorted and warped disk. After assigning to every observed star a probability of belonging to a given age-metallicity interval, we build high-resolution population maps. These dramatically reveal the flocculent nature of the young star-forming regions and the nearly smooth features traced by older stellar generations. They document the formation of the SMC Wing at ages <0.2 Gyr and the peak of star formation in the SMC Bar at 40 Myr. We clearly detect periods of enhanced star formation at 1.5 Gyr and 5 Gyr. The former is possibly related to a new feature found in the AMR, which suggests ingestion of metal-poor gas at ages slightly larger than 1 Gyr. The latter constitutes a major period of stellar mass formation. We confirm that the SFR(t) was moderately low at even older ages.

We present color-magnitude diagram analysis of deep Hubble Space Telescope imaging of a mass-limited sample of 18 intermediate-age (1 - 2 Gyr old) star clusters in the Magellanic Clouds, including 8 clusters for which new data was obtained. We find that ${\it all}$ star clusters in our sample feature extended main sequence turnoff (eMSTO) regions that are wider than can be accounted for by a simple stellar population (including unresolved binary stars). FWHM widths of the MSTOs indicate age spreads of 200-550 Myr. We evaluate dynamical evolution of clusters with and without initial mass segregation. Our main results are: (1) the fraction of red clump (RC) stars in secondary RCs in eMSTO clusters scales with the fraction of MSTO stars having pseudo-ages $\leq 1.35$ Gyr; (2) the width of the pseudo-age distributions of eMSTO clusters is correlated with their central escape velocity $v_{\rm esc}$, both currently and at an age of 10 Myr. We find that these two results are unlikely to be reproduced by the effects of interactive binary stars or a range of stellar rotation velocities. We therefore argue that the eMSTO phenomenon is mainly caused by extended star formation within the clusters; (3) we find that $v_{\rm esc} \geq 15$ km/s out to ages of at least 100 Myr for ${\it all}$ clusters featuring eMSTOs, while $v_{\rm esc} \leq 12$ km/s at all ages for two lower-mass clusters in the same age range that do ${\it not}$ show eMSTOs. We argue that eMSTOs only occur for clusters whose early escape velocities are higher than the wind velocities of stars that provide material from which second-generation stars can form. The threshold of 12-15 km/s is consistent with wind velocities of intermediate-mass AGB stars in the literature.

We discuss the dust chemistry and growth in the circumstellar envelopes (CSEs) of Thermally Pulsing Asymptotic Giant Branch (TP-AGB) star models computed with the COLIBRI code, at varying initial mass and metallicity (Z=0.001, 0.008, 0.02, 0.04, 0.06). A relevant result of our analysis deals with the silicate production in M-stars. We show that, in order to reproduce the observed trend between terminal velocities and mass-loss rates in Galactic M-giants, one has to significantly reduce the efficiency of chemisputtering by H2 molecules, usually considered as the most effective dust destruction mechanism. This indication is also in agreement with the most recent laboratory results, which show that silicates may condense already at T=1400 K, instead than at Tcond=1000 K, as obtained by models that include chemisputtering. From the analysis of the total dust ejecta, we find that the total dust-to-gas ejecta of intermediate-mass stars are much less dependent on metallicity than usually assumed. In a broader context, our results are suitable to study the dust enrichment of the interstellar medium provided by TP-AGB stars in both nearby and high redshift galaxies.

In the last ten years three main facts about the thermally pulsing asymptotic giant branch (TP-AGB) have become evident: 1) the modelling of the TP-AGB phase is critical for the derivation of basic galaxy properties (e.g. mass and age) up to high redshift, with consequent cosmological implications; 2) current TP-AGB calibrations based on Magellanic Cloud (MC) clusters come out not to work properly for other external galaxies, yielding a likely TP-AGB overestimation; 3) the significance of the TP-AGB contribution in galaxies, hence their derived properties, are strongly debated, with conflicting claims in favour of either a heavy or a light TP-AGB. The only way out of this condition of persisting uncertainty is to perform a reliable calibration of the TP-AGB phase as a function of the star's initial mass (hence age) over a wide range of metallicity, from very low to super-solar values. In this context, I will review recent advancements and ongoing efforts towards a physically-sound TP-AGB calibration that, moving beyond the classical use of the MC clusters, combines increasingly refined TP-AGB stellar models with exceptionally high-quality data for resolved TP-AGB stars in nearby galaxies. Preliminary results indicate that a sort of "TP-AGB island" emerges in the age-metallicity plane, where the contribution of these stars is especially developed, embracing preferentially solar- and MC-like metallicities, and intermediate ages (~ few Gyr).

The evolution and lifetimes of thermally pulsating asymptotic giant branch (TP-AGB) stars suffer from significant uncertainties. We present a detailed framework for constraining model luminosity functions of TP-AGB stars using resolved stellar populations. We show an example of this method that compares various TP-AGB mass-loss prescriptions that differ in their treatments of mass loss before the onset of dust-driven winds (pre-dust). We find that models with more efficient pre-dust driven mass loss produce results consistent with observations, as opposed to more canonical mass-loss models. Efficient pre-dust driven mass-loss predicts for [Fe/H] < -1.2, lower mass TP-AGB stars (M < 1 Msun) must have lifetimes less than about 1.2 Myr.

We extend the {\sl\,PARSEC} library of stellar evolutionary tracks by computing new models of massive stars, from 14\Msun to 350\Msun. The input physics is the same used in the {\sl\,PARSEC}~V1.1 version, but for the mass-loss rate which is included by considering the most recent updates in literature. We focus on low metallicity, $Z$=0.001 and $Z$=0.004, for which the metal poor dwarf irregular star forming galaxies, Sextans A, WLM and NCG6822, provide simple but powerful workbenches. The models reproduce fairly well the observed CMDs but the stellar colour distributions indicate that the predicted blue loop is not hot enough in models with canonical extent of overshooting. In the framework of a mild extended mixing during central hydrogen burning, the only way to reconcile the discrepancy is to enhance the overshooting at the base of the convective envelope (EO) during the first dredge-UP. The mixing scales required to reproduce the observed loops, EO=2\HP or EO=4\HP, are definitely larger than those derived from, e.g., the observed location of the RGB bump in low mass stars. This effect, if confirmed, would imply a strong dependence of the mixing scale below the formal Schwarzschild border, on the stellar mass or luminosity. Reproducing the features of the observed CMDs with standard values of envelope overshooting would require a metallicity significantly lower than the values measured in these galaxies. Other quantities, such as the star formation rate and the initial mass function, are only slightly sensitive to this effect. Future investigations will consider other metallicities and different mixing schemes.

Using the stellar population synthesis tool TRILEGAL, we discuss the expected stellar populations in the Kepler and CoRoT fields

In spite of the great effort made in the last decades to improve our understanding of stellar evolution, significant uncertainties remain due to our poor knowledge of some complex physical processes that require an empirical calibration, such as the efficiency of the interior mixing related to convective overshoot. Here we review the impact of convective overshoot on the evolution of stars during the main Hydrogen and Helium burning phases.

Many stellar models present difficulties in reproducing basic observational relations of very low mass stars (VLMS), including the mass--radius relation and the optical colour--magnitudes of cool dwarfs. Here, we improve PARSEC models on these points. We implement the T--tau relations from PHOENIX BT-Settl model atmospheres as the outer boundary conditions in the PARSEC code, finding that this change alone reduces the discrepancy in the mass--radius relation from 8 to 5 per cent. We compare the models with multi--band photometry of clusters Praesepe and M67, showing that the use of T--tau relations clearly improves the description of the optical colours and magnitudes. But anyway, using both Kurucz and PHOENIX model spectra, model colours are still systematically fainter and bluer than the observations. We then apply a shift to the above T--tau relations, increasing from 0 at T_eff = 4730 K to ~14% at T_eff = 3160 K, to reproduce the observed mass--radius radius relation of dwarf stars. Taking this experiment as a calibration of the T--tau relations, we can reproduce the optical and near infrared CMDs of low mass stars in the old metal--poor globular clusters NGC6397 and 47Tuc, and in the intermediate--age and young solar--metallicity open clusters M67 and Praesepe. Thus, we extend PARSEC models using this calibration, providing VLMS models more suitable for the lower main sequence stars over a wide range of metallicities and wavelengths. Both sets of models are available on PARSEC webpage.

The evolution and lifetimes of thermally pulsating asymptotic giant branch (TP-AGB) stars suffer from significant uncertainties. In this work, we analyze the numbers and luminosity functions of TP-AGB stars in six quiescent, low metallicity ([Fe/H] $\lesssim -0.86$) galaxies taken from the ANGST sample, using HST photometry in both optical and near-infrared filters. The galaxies contain over 1000 TP-AGB stars (at least 60 per field). We compare the observed TP-AGB luminosity functions and relative numbers of TP-AGB and RGB stars, to models generated from different suites of TP-AGB evolutionary tracks after adopting star formation histories (SFH) derived from the HST deep optical observations. We test various mass-loss prescriptions that differ in their treatments of mass-loss before the onset of dust-driven winds (pre-dust). These comparisons confirm that pre-dust mass-loss is important, since models that neglect pre-dust mass-loss fail to explain the observed TP-AGB/RGB ratio or the luminosity functions. In contrast, models with more efficient pre-dust mass-loss produce results consistent with observations. We find that for [Fe/H]$\lesssim-0.86$, lower mass TP-AGB stars ($M\lesssim1M_\odot$) must have lifetimes of $\sim0.5$ Myr and higher masses ($M\lesssim 3M_\odot$) must have lifetimes $\lesssim 1.2$ Myr. In addition, assuming our best-fitting mass-loss prescription, we show that the third dredge up has no significant effect on TP-AGB lifetimes in this mass and metallicity range.

We present HI spectral-line imaging of the extremely metal-poor galaxy DDO 68. This system has a nebular oxygen abundance of only 3% Z$_{\odot}$, making it one of the most metal-deficient galaxies known in the local volume. Surprisingly, DDO 68 is a relatively massive and luminous galaxy for its metal content, making it a significant outlier in the mass-metallicity and luminosity-metallicity relationships. The origin of such a low oxygen abundance in DDO 68 presents a challenge for models of the chemical evolution of galaxies. One possible solution to this problem is the infall of pristine neutral gas, potentially initiated during a gravitational interaction. Using archival HI spectral-line imaging obtained with the Karl G. Jansky Very Large Array, we have discovered a previously unknown companion of DDO 68. This low-mass (M$_{\rm HI}$ $=$ 2.8$\times$10$^{7}$ M$_{\odot}$), recently star-forming (SFR$_{\rm FUV}$ $=$ 1.4$\times$10$^{-3}$ M$_{\odot}$ yr$^{-1}$, SFR$_{\rm H\alpha}$ $<$ 7$\times$10$^{-5}$ M$_{\odot}$ yr$^{-1}$) companion has the same systemic velocity as DDO 68 (V$_{\rm sys}$ $=$ 506 km s$^{-1}$; D $=$ 12.74$\pm$0.27 Mpc) and is located at a projected distance of 42 kpc. New HI maps obtained with the 100m Robert C. Byrd Green Bank Telescope provide evidence that DDO 68 and this companion are gravitationally interacting at the present time. Low surface brightness HI gas forms a bridge between these objects.

We establish new constraints on the intermediate-mass range of the initial-final mass relation by studying white dwarfs in four young star clusters, and apply the results to study the evolution of stars on the thermally pulsing asymptotic giant branch (TP-AGB). We show that the stellar core mass on the AGB grows rapidly from 10% to 30% for stars with $M_{\rm initial}$ = 1.6 to 2.0 $M_\odot$. At larger masses, the core-mass growth decreases steadily to $\sim$10% at $M_{\rm initial}$ = 3.4 $M_\odot$. These observations are in excellent agreement with predictions from the latest TP-AGB evolutionary models in Marigo et al. (2013). We also compare to models with varying efficiencies of the third dredge-up and mass loss, and demonstrate that the process governing the growth of the core is largely the stellar wind, while the third dredge-up plays a secondary, but non-negligible role. Based on the new white dwarf measurements, we perform an exploratory calibration of the most popular mass-loss prescriptions in the literature. Finally, we estimate the lifetime and the integrated luminosity of stars on the TP-AGB to peak at $t$ $\sim$ 3 Myr and $E$ = 1.2 $\times$ 10$^{10}$ $L_\odot$ yr for $M_{\rm initial}$ $\sim$ 2 $M_\odot$ ($t$ $\sim$ 2 Myr for luminosities brighter than the RGB tip at $\log(L/L_{\odot})$ $>$ 3.4), decreasing to $t$ = 0.4 Myr and $E$ = 6.1 $\times$ 10$^{9}$ $L_\odot$ yr for stars with $M_{\rm initial}$ $\sim$ 3.5 $M_\odot$. The implications of these results are discussed with respect to general population synthesis studies that require correct modeling of the TP-AGB phase of stellar evolution.

We extend the formalism presented in our recent calculations of dust ejecta from the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase, to the case of super-solar metallicity stars. The TP-AGB evolutionary models are computed with the COLIBRI code. We adopt our preferred scheme for dust growth. For M-giants, we neglect chemisputtering by H$_2$ molecules and, for C-stars we assume a homogeneous growth scheme which is primarily controlled by the carbon over oxygen excess. At super-solar metallicities, dust forms more efficiently and silicates tend to condense significantly closer to the photosphere (r~1.5 R$_*$) - and thus at higher temperatures and densities - than at solar and sub-solar metallicities (r~2-3 R$_*$). In such conditions, the hypothesis of thermal decoupling between gas and dust becomes questionable, while dust heating due to collisions plays an important role. The heating mechanism delays dust condensation to slightly outer regions in the circumstellar envelope. We find that the same mechanism is not significant at solar and sub-solar metallicities. The main dust products at super-solar metallicities are silicates. We calculate the total dust ejecta and dust-to-gas ejecta, for various values of the stellar initial masses and initial metallicities Z=0.04, 0.06. Merging these new calculations with those for lower metallicities it turns out that, contrary to what often assumed, the total dust-to-gas ejecta of intermediate-mass stars exhibit only a weak dependence on the initial metal content.

(Abridged) In the recent controversy about the role of TP-AGB stars in evolutionary population synthesis (EPS) models of galaxies, one particular aspect is puzzling: TP-AGB models aimed at reproducing the lifetimes and integrated fluxes of the TP-AGB phase in Magellanic Cloud (MC) clusters, when incorporated into EPS models, are found to overestimate the TP-AGB contribution in resolved star counts and integrated spectra of galaxies. In this paper, we call attention to a particular evolutionary aspect that in all probability is the main cause of this conundrum. As soon as stellar populations intercept the ages at which RGB stars first appear, a sudden change in the lifetime of the core He-burning phase causes a temporary boost in the production rate of subsequent evolutionary phases, including the TP-AGB. For a timespan of about 0.1 Gyr, triple TP-AGB branches develop at slightly different initial masses, causing their frequency and contribution to the integrated luminosity of the stellar population to increase by a factor of 2. The boost occurs just in the proximity of the expected peak in the TP-AGB lifetimes, and for ages of 1.6 Gyr. Coincidently, this relatively narrow age interval happens to contain the few very massive MC clusters that host most of the TP-AGB stars used to constrain stellar evolution and EPS models. This concomitance makes the AGB-boosting particularly insidious in the context of present EPS models. The effect brings about three main consequences. (1) Present estimates of the TP-AGB contribution to the integrated light of galaxies derived from MC clusters, are biased towards too large values. (2) The relative TP-AGB contribution of single-burst populations falling in this critical age range cannot be accurately derived by the fuel consumption theorem. (3) A careful revision of AGB star populations in intermediate-age MC clusters is urgently demanded.

We use medium-band near-infrared (NIR) Hubble Space Telescope WFC3 photometry with model NIR spectra of Asymptotic Giant Branch (AGB) stars to develop a new tool for efficiently distinguishing carbon-rich (C-type) AGB stars from oxygen-rich (M-type) AGB stars in galaxies at the edge of and outside the Local Group. We present the results of a test of this method on a region of the inner disk of M31, where we find a surprising lack of C stars, contrary to the findings of previous C star searches in other regions of M31. We find only 1 candidate C star (plus up to 6 additional, less certain C stars candidates), resulting in an extremely low ratio of C to M stars (C/M = (3.3(+20,-0.1))x10^-4) that is 1-2 orders of magnitude lower than other C/M estimates in M31. The low C/M ratio is likely due to the high metallicity in this region which impedes stars from achieving C/O > 1 in their atmospheres. These observations provide stringent constraints evolutionary models of metal-rich AGB stars and suggest that there is a metallicity threshold above which M stars are unable to make the transition to C stars, dramatically affecting AGB mass loss and dust production and, consequently, the observed global properties of metal-rich galaxies.

We present the dust ejecta of the new stellar models for the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase computed with the COLIBRI code. We use a formalism of dust growth coupled with a stationary wind for both M and C-stars. In the original version of this formalism, the most efficient destruction process of silicate dust in M-giants is chemisputtering by H2 molecules. For these stars we find that dust grains can only form at relatively large radial distances (r~5 R*), where they cannot be efficiently accelerated, in agreement with other investigations. In the light of recent laboratory results, we also consider the alternative case that the condensation temperature of silicates is determined only by the competition between growth and free evaporation processes (i.e. no chemisputtering). With this latter approach we obtain dust condensation temperatures that are significantly higher (up to Tcond~1400 K) than those found when chemisputtering is included (Tcond~900 K), and in better agreement with condensation experiments. As a consequence, silicate grains can remain stable in inner regions of the circumstellar envelopes (r~2 R*), where they can rapidly grow and can be efficiently accelerated. With this modification, our models nicely reproduce the observed trend between terminal velocities and mass loss rates of Galactic M-giants. For C-stars the formalism is based on the homogeneous growth scheme where the key role is played by the carbon over oxygen excess. The models reproduce fairly well the terminal velocities of Galactic stars and there is no need to invoke changes in the standard assumptions. At decreasing metallicity the carbon excess becomes more pronounced and the efficiency of dust formation increases. This trend could be in tension with recent observational evidence in favour of a decreasing efficiency, at decreasing metallicity.

We present the COLIBRI code for computing the evolution of stars along the TP-AGB phase. Compared to purely synthetic TP-AGB codes, COLIBRI relaxes a significant part of their analytic formalism in favour of a detailed physics applied to a complete envelope model, in which the stellar structure equations are integrated from the atmosphere down to the bottom of the hydrogen-burning shell. This allows to predict self-consistently: (i) the effective temperature, and more generally the convective envelope and atmosphere structures, correctly coupled to the changes in the surface chemical abundances and gas opacities; (ii) sphericity effects in the atmospheres; (iii) the core mass-luminosity relation and its break-down due to hot bottom burning (HBB) in the most massive AGB stars, (iv) the HBB nucleosynthesis via the solution of a complete nuclear network (pp chains, and the CNO, NeNa, MgAl cycles), including also the production of 7Li via the Cameron-Fowler beryllium transport mechanism; (v) the chemical composition of the pulse-driven convective zone; (vi) the onset and quenching of the third dredge-up, with a suitable temperature criterion. At the same time COLIBRI pioneers new techniques in the treatment of the physics of stellar interiors. It is the first evolutionary code ever to use accurate on-the-fly computation of the equation of state for roughly 800 atoms, ions, molecules, and of the Rosseland mean opacities throughout the deep envelope. Another distinguishing aspect of COLIBRI is its high computational speed. This feature is necessary for calibrating the uncertain parameters and processes that characterize the TP-AGB phase, a step of paramount importance for producing reliable stellar population synthesis models of galaxies up to high redshift. (abridged)

Based on new observations with the Wide Field Camera 3 onboard the Hubble Space Telescope, we report the discovery of an extended main sequence turn-off (eMSTO) in the intermediate-age star cluster NGC411. This is the second case of an eMSTO being identified in a star cluster belonging to the Small Magellanic Cloud (SMC), after NGC419. Despite the present masses of these two SMC clusters differ by a factor of 4, the comparison between their colour--magnitude diagrams (CMD) shows striking similarities, especially regarding the shape of their eMSTOs. The loci of main CMD features are so similar that they can be well described, in a first approximation, by the same mean metallicity, distance and extinction. NGC411, however, presents merely a trace of secondary red clump as opposed to its prominent manifestation in NGC419. This could be due either to the small number statistics in NGC411, or by the star formation in NGC419 having continued for 60 Myr longer than in NGC411. Under the assumption that the eMSTOs are caused by different generations of stars at increasing age, both clusters are nearly coeval in their first episodes of star formation. The initial period of star formation, however, is slightly more marked in NGC419 than in NGC411. We discuss these findings in the context of possible scenarios for the origin of eMSTOs.

In spite of the great effort made in the last decades to improve our understanding of stellar evolution, significant uncertainties still remain due to our poor knowledge of some complex physical processes that still require an empirical calibration, such as the efficiency of convective heat transport and interior mixing. Here we will review the impact of these uncertainties on the evolution of red giant stars.

NGC1846 and NGC1783 are two massive star clusters in the Large Magellanic Cloud, hosting both an extended main sequence turn-off and a dual clump of red giants. They present similar masses but differ mainly in angular size. Starting from their high-quality ACS data in the F435W, F555W and F814W filters, and updated sets of stellar evolutionary tracks, we derive their star formation rates as a function of age, SFR(t), by means of the classical method of CMD reconstruction which is usually applied to nearby galaxies. The method confirms the extended periods of star formation derived from previous analysis of the same data. When the analysis is performed for a finer resolution in age, we find clear evidence for a 50-Myr long hiatus between the oldest peak in the SFR(t), and a second prolonged period of star formation, in both clusters. For the more compact cluster NGC1846, there seems to be no significant difference between the SFR(t) in the cluster centre and in an annulus with radii between 20 and 60 arcsec (from 4.8 to 15.4 pc). The same does not occur in the more extended NGC1783 cluster, where the outer ring (between 33 and 107 arcsec, from 8.0 to 25.9 pc) is found to be slightly younger than the centre. We also explore the best-fitting slope of the present-day mass function and binary fraction for the different cluster regions, finding hints of a varying mass function between centre and outer ring in NGC1783. These findings are discussed within the present scenarios for the formation of clusters with multiple turn-offs.

Stellar parameter determination is a challenging task when dealing with galactic giant stars. The combination of different investigation techniques has proven to be a promising approach. We analyse archive spectra obtained with the Short-Wavelength-Spectrometer (SWS) onboard of ISO, and new interferometric observations from the Very Large Telescope MID-infrared Interferometric instrument (VLTI/MIDI) of a very well studied carbon-rich giant: TX Psc. The aim of this work is to determine stellar parameters using spectroscopy and interferometry. The observations are used to constrain the model atmosphere, and eventually the stellar evolutionary model in the region where the tracks map the beginning of the carbon star sequence. Two different approaches are used to determine stellar parameters: (i) the 'classic' interferometric approach where the effective temperature is fixed by using the angular diameter in the N-band (from interferometry) and the apparent bolometric magnitude; (ii) parameters are obtained by fitting a grid of state-of-the-art hydrostatic models to spectroscopic and interferometric observations. We find a good agreement between the parameters of the two methods. The effective temperature and luminosity clearly place TX Psc in the carbon-rich AGB star domain in the H-R-diagram. Current evolutionary tracks suggest that TX Psc became a C-star just recently, which means that the star is still in a 'quiet' phase compared to the subsequent strong-wind regime. This is in agreement with the C/O ratio being only slightly larger than 1.

We present the updated version of the code used to compute stellar evolutionary tracks in Padova. It is the result of a thorough revision of the major input physics, together with the inclusion of the pre-main sequence phase, not present in our previous releases of stellar models. Another innovative aspect is the possibility of promptly generating accurate opacity tables fully consistent with any selected initial chemical composition, by coupling the OPAL opacity data at high temperatures to the molecular opacities computed with our AESOPUS code (Marigo & Aringer 2009). In this work we present extended sets of stellar evolutionary models for various initial chemical compositions, while other sets with different metallicities and/or different distributions of heavy elements are being computed. For the present release of models we adopt the solar distribution of heavy elements from the recent revision by Caffau et al. (2011), corresponding to a Sun's metallicity Z=0.0152. From all computed sets of stellar tracks, we also derive isochrones in several photometric systems. The aim is to provide the community with the basic tools to model star clusters and galaxies by means of population synthesis techniques.

We reconsider the role of pre-main sequence (pre-MS) Li depletion on the basis of new observational and theoretical evidence: i) new observations of Halpha emissions in young clusters show that mass accretion could be continuing till the first stages of the MS, ii) theoretical implications from helioseismology suggest large overshooting values below the bottom of the convective envelopes. We argue here that a significant pre-MS 7Li destruction, caused by efficient overshoot mixing, could be followed by a matter accretion after 7Li depletion has ceased on MS thus restoring Li almost to the pristine value. As a test case we show that a halo dwarf of 0.85 Msun with an extended overshooting envelope starting with an initial abundance of A(Li) = 2.74 would burn Li completely, but an accretion rate of the type 1e-8xe^{-t/3e6} Msun yr$^{-1}$ would restore Li to end with an A(Li) = 2.31. A self-regulating process is required to produce similar final values in a range of different stellar masses to explain the PopII Spite plateau. However, this framework could explain why open cluster stars have lower Li abundances than the pre-solar nebula, the absence of Li in the most metal poor dwarfs and a number of other features which lack of a satisfactory explanation.

Using high spatial resolution HST WFC3 and ACS imaging of resolved stellar populations, we constrain the contribution of thermally-pulsing asymptotic giant branch (TP-AGB) stars and red helium burning (RHeB) stars to the 1.6 um near-infrared (NIR) luminosities of 23 nearby galaxies. The TP-AGB phase contributes as much as 17% of the integrated F160W flux, even when the red giant branch is well populated. The RHeB population contribution can match or even exceed the TP-AGB contribution, providing as much as 21% of the integrated F160W light. The NIR mass-to-light (M/L) ratio should therefore be expected to vary significantly due to fluctuations in the star formation rate over timescales from 25 Myr to several Gyr. We compare our observational results to predictions based on optically derived star formation histories and stellar population synthesis (SPS) models, including models based on the Padova isochrones (used in popular SPS programs). The SPS models generally reproduce the expected numbers of TP-AGB stars in the sample. The same SPS models, however, give a larger discrepancy in the F160W flux contribution from the TP-AGB stars, over-predicting the flux by a weighted mean factor of 2.3 +/-0.8. This larger offset is driven by the prediction of modest numbers of high luminosity TP-AGB stars at young (<300 Myrs) ages. The best-fit SPS models simultaneously tend to under-predict the numbers and fluxes of stars on the RHeB sequence, typically by a factor of 2.0+/-0.6 for galaxies with significant numbers of RHeBs. Coincidentally, over-prediction of the TP-AGB and under-prediction of the RHeBs result in a NIR M/L ratio largely unchanged for a rapid star formation rate. However, the NIR-to-optical flux ratio of galaxies could be significantly smaller than AGB-rich models would predict, an outcome that has been observed in some intermediate redshift post-starburst galaxies. (Abridged)

We derive the star formation history for several regions of the LMC, using deep near-infrared data from the VISTA near-infrared YJKs survey of the Magellanic system (VMC). The regions include three almost-complete 1.4 sqdeg tiles located 3.5 deg away from the LMC centre in distinct directions. To this dataset, we add two 0.036 sqdeg subregions inside the 30 Doradus tile. The SFH is derived from the simultaneous reconstruction of two different CMDs, using the minimization code StarFISH. The distance modulus (m-M)_0 and extinction Av is varied within intervals 0.2 and 0.5 mag wide, respectively, within which we identify the best-fitting star formation rate SFR(t), age-metallicity relation (AMR), (m-M)_0 and Av. Our results demonstrate that VMC data, due to the combination of depth and little sensitivity to differential reddening, allow the derivation of the space-resolved SFH of the LMC with unprecedented quality compared to previous wide-area surveys. In particular, the data clearly reveal the presence of peaks in the SFR(t) at ages log(t/yr)=9.3 and 9.7, which appear in most of the subregions. The most recent SFR is found to vary greatly from subregion to subregion, with the general trend of being more intense in the innermost LMC, except for the tile next to the N11 complex. In the bar region, the SFR seems remarkably constant over the time interval from 8.4 to 9.7. The AMRs, instead, turn out to be remarkably similar across the LMC. The fields studied so far are fit extremely well by a single disk of inclination 26.2+-2.0 deg, position angle of the line of nodes 129.1+-13.0 deg, and distance modulus of 18.470+-0.006 mag (random errors only) up to the LMC centre.

We present near-infrared (NIR) color-magnitude diagrams (CMDs) for the resolved stellar populations within 26 fields of 23 nearby galaxies (<4 Mpc), based on F110W and F160W images from Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). The CMDs sample both old dormant and young star-forming populations. We match key NIR CMD features with their counterparts in optical CMDs, and identify the red core Helium burning (RHeB) sequence as a significant contributor to the NIR flux in stellar populations younger than a few 100 Myrs old, suggesting that star formation can drive surprisingly rapid variations in the NIR mass-to-light ratio. The NIR luminosity of star forming galaxies is therefore not necessarily proportional to the stellar mass. We note that these individual bright RHeB stars may be misidentified as old stellar clusters in low resolution imaging. We also discuss the CMD location of asymptotic giant branch (AGB) stars, and the separation of AGB sub-populations using a combination of optical and NIR colors. We empirically calibrate the NIR magnitude of the tip of the red giant branch (TRGB) as a function of color, allowing this widely adopted filter to be used for distance measurements. We find a clear trend between NIR RGB color and metallicity. However, it appears unlikely that the slope of the NIR RGB can be used as a metallicity indicator in extragalactic systems with comparable data. Finally, we discuss scattered light in the WFC3, which becomes significant for exposures taken close to a bright earth limb.

We include an energy term based on Dark Matter (DM) self-annihilation during the cooling and subsequent collapse of the metal-free gas, in halos hosting the formation of the first stars in the Universe. We have found that the feedback induced on the chemistry of the cloud does modify the properties of the gas throughout the collapse. However, the modifications are not dramatic, and the typical Jeans mass within the halo is conserved throughout the collapse, for all the DM parameters we have considered. This result implies that the presence of Dark Matter annihilations does not substantially modify the Initial Mass Function of the First Stars, with respect to the standard case in which such additional energy term is not taken into account. We have also found that when the rate of energy produced by the DM annihilations and absorbed by the gas equals the chemical cooling (at densities yet far from the actual formation of a proto-stellar core) the structure does not halt its collapse, although that proceeds more slowly by a factor smaller than few per cent of the total collapse time.

We present comprehensive catalogs of hot star candidates in the Milky Way, selected from GALEX far-UV (FUV, 1344-1786 AA) and near-UV (NUV, 1771-2831 AA) imaging. The FUV and NUV photometry allows us to extract the hottest stellar objects, in particular hot white dwarfs (WD), which are elusive at other wavelengths because of their high temperatures and faint optical luminosities. We generated catalogs of UV sources from two GALEX's surveys: AIS (All-Sky Imaging Survey, depth ABmag~19.9/20.8 in FUV/NUV) and MIS (Medium-depth Imaging Survey, depth ~22.6/22.7mag). The two catalogs (from GALEX fifth data release) contain 65.3/12.6 million (AIS/MIS) unique UV sources with error(NUV)<0.5mag, over 21,435/1,579 square degrees. We also constructed subcatalogs of the UV sources with matched optical photometry from SDSS (7th data release): these contain 0.6/0.9million (AIS/MIS) sources with errors <0.3mag in both FUV and NUV, excluding sources with multiple optical counterparts, over an area of 7,325/1,103 square degrees. All catalogs are available online. We then selected 28,319(AIS)/9,028(MIS) matched sources with FUV-NUV<-0.13; this color cut corresponds to stellar Teff hotter than ~18,000 K. An additional color cut of NUV-r>0.1 isolates binaries with largely differing Teff's, and some intruding QSOs. Available spectroscopy for a subsample indicates that hot-star candidates with NUV-r<0.1 have negligible contamination by non-stellar objects. We discuss the distribution of sources in the catalogs, and the effects of error and color cuts on the samples. The density of hot-star candidates increases from high to low Galactic latitudes, but drops on the MW plane due to dust extinction. Our hot-star counts at all latitudes are better matched by Milky Way models computed with an initial-final mass relation that favours lower final masses. (ABRIDGED)

In an attempt to constrain evolutionary models of the asymptotic giant branch (AGB) phase at the limit of low masses and low metallicities, we have examined the luminosity functions and number ratio between AGB and red giant branch (RGB) stars from a sample of resolved galaxies from the ACS Nearby Galaxy Survey Treasury (ANGST). This database provides HST optical photometry together with maps of completeness, photometric errors, and star formation histories for dozens of galaxies within 4 Mpc. We select 12 galaxies characterized by predominantly metal-poor populations as indicated by a very steep and blue RGB, and which do not present any indication of recent star formation in their color--magnitude diagrams. Thousands of AGB stars brighter than the tip of the RGB (TRGB) are present in the sample (between 60 and 400 per galaxy), hence the Poisson noise has little impact in our measurements of the AGB/RGB ratio. We model the photometric data with a few sets of thermally pulsing AGB (TP-AGB) evolutionary models with different prescriptions for the mass loss. This technique allows us to set stringent constraints to the TP-AGB models of low-mass metal-poor stars (with M<1.5 Msun, [Fe/H]<~-1.0). Indeed, those which satisfactorily reproduce the observed AGB/RGB ratios have TP-AGB lifetimes between 1.2 and 1.8 Myr, and finish their nuclear burning lives with masses between 0.51 and 0.55 Msun. This is also in good agreement with recent observations of white dwarf masses in the M4 old globular cluster. These constraints can be added to those already derived from Magellanic Cloud star clusters as important mileposts in the arduous process of calibrating AGB evolutionary models.

We investigate the main physical properties of low-metallicity Asymptotic Giant Branch stars, with the aim of quantifying the uncertainties that presently affect the predicted chemical yields of these stars, associated to mass loss and description of molecular opacities. We find that above a threshold mass, M ~ 3.5Msun for Z=0.001, the results are little dependent on the opacity treatment, as long as hot-bottom burning prevents the surface C/O ratio from exceeding unity; the yields of these massive AGB stars are expected to be mostly determined by the efficiency of convection, with a relatively mild dependence on the mass-loss description. A much higher degree of uncertainty is associated to the yields of less massive models, which critically depend on the adopted molecular opacities. An interval of masses exists, say 2.0-3.0Msun, (the exact range depends on mass loss), in which HBB may be even extinguished following the cooling produced by the opacity of C-bearing molecules. The yields of these stars are the most uncertain, the variation range being the largest (up to ~ 2dex) for the nitrogen and sodium yields. For very low-mass models, not experiencing hot-bottom burning (M< 1.5Msun),the description of mass loss and the treatment of the convective boundaries are crucial for the occurrence of the third dredge-up, with sizable consequences on the CNO yields.

A near-infrared imaging study of the evolved stellar populations in the dwarf spheroidal galaxy Leo I is presented. Based on JHK observations obtained with the WFCAM wide-field array at the UKIRT telescope, we build a near-infrared photometric catalogue of red giant branch (RGB) and asymptotic giant branch (AGB) stars in Leo I over a 13.5 arcmin square area. The V-K colours of RGB stars, obtained by combining the new data with existing optical observations, allow us to derive a distribution of global metallicity [M/H] with average [M/H] = -1.51 (uncorrected) or [M/H] = -1.24 +/- 0.05 (int) +/- 0.15 (syst) after correction for the mean age of Leo I stars. This is consistent with the results from spectroscopy once stellar ages are taken into account. Using a near-infrared two-colour diagram, we discriminate between carbon- and oxygen-rich AGB stars and obtain a clean separation from Milky Way foreground stars. We reveal a concentration of C-type AGB stars relative to the red giant stars in the inner region of the galaxy, which implies a radial gradient in the intermediate-age (1-3 Gyr) stellar populations. The numbers and luminosities of the observed carbon- and oxygen-rich AGB stars are compared with those predicted by evolutionary models including the thermally-pulsing AGB phase, to provide new constraints to the models for low-metallicity stars. We find an excess in the predicted number of C stars fainter than the RGB tip, associated to a paucity of brighter ones. The number of O-rich AGB stars is roughly consistent with the models, yet their predicted luminosity function is extended to brighter luminosity. It appears likely that the adopted evolutionary models overestimate the C star lifetime and underestimate their K-band luminosity.

We introduce a new tool - AESOPUS: Accurate Equation of State and OPacity Utility Software - for computing the equation of state and the Rosseland mean (RM) opacities of matter in the ideal gas phase. Results are given as a function of one pair of state variables, (i.e. temperature T in the range 3.2 <= log(T) <= 4.5, and parameter R= rho/(T/10^6 K)^3 in the range -8 <= log(R) <= 1), and arbitrary chemical mixture. The chemistry is presently solved for about 800 species, consisting of almost 300 atomic and 500 molecular species. The gas opacities account for many continuum and discrete sources, including atomic opacities, molecular absorption bands, and collision-induced absorption. Several tests made on AESOPUS have proved that the new opacity tool is accurate in the results,flexible in the management of the input prescriptions, and agile in terms of computational time requirement. We set up a web-interface (http://stev.oapd.inaf.it/aesopus) which enables the user to compute and shortly retrieve RM opacity tables according to his/her specific needs, allowing a full degree of freedom in specifying the chemical composition of the gas. Useful applications may regard RM opacities of gas mixtures with i) scaled-solar abundances of metals, choosing among various solar mixture compilations available in the literature; ii) varying CNO abundances, suitable for evolutionary models of red and asymptotic giant branch stars and massive stars in the Wolf-Rayet stages; iii) various degrees of enhancement in alpha-elements, and C-N, Na-O and Mg-Al abundance anti-correlations, necessary to properly describe the properties of stars in early-type galaxies and Galactic globular clusters; iv) zero-metal abundances appropriate for studies of gas opacity in primordial conditions.

We extend our theoretical computations for low-mass stars to intermediate-mass and massive stars, for which few databases exist in the literature. Evolutionary tracks and isochrones are computed from 2.50 to 20 solar masses for agrid of 37 chemical compositions with metal content Z between 0.0001 and 0.070 and helium content Y between 0.23 and 0.40. Synthetic TP-AGB models allow stellar tracks and isochrones to be extended until the end of the thermal pulses along the AGB. We provide software tools for the bidimensional interpolation (in Y and Z) of the isochrones. We present tracks for scaled-solar abundances and the corresponding isochrones from very old ages down to about 10 million years. The extension of the blue loops and the instability strip of Cepheid stars are compared and the Cepheid mass-discrepancy is discussed. The location of red supergiants in the H-R diagram is in good agreement with the evolutionary tracks for masses from 10 to 20 solar masses. Tracks and isochrones are available in tabular form for the adopted grid of chemical compositions in the extended plane Z-Y in three photometric systems. An interactive web interface allows users to obtain isochrones of any chemical composition inside the provided Z-Y range and also to simulate stellar populations with different Y(Z) helium-to-metal enrichment laws.

The studies focused on the Thermally-Pulsing Asymptotic Giant Branch phase experienced by low- and intermediate-mass stars are extremely important in many astrophysical contexts. In particular, a detailed computation of their chemical yields is essential for several issues, ranging from the chemical evolution of galaxies, to the mechanisms behind the formation of globular clusters. Among all the uncertainties affecting the theoretical modelling of this phase, and described in the literature, it remains to be fully clarified which results are severely affected by the use of inadequate low-temperature opacities, that are in most cases calculated on the basis of the original chemical composition of the stars, and do not consider the changes in the surface chemistry due to the occurrence of the third dredge-up and hot-bottom burning. Our investigation is aimed at investigating this point. By means of full evolutionary models including new set of molecular opacities computed specifically with the AESOPUS tool, we highlight which stellar models, among those present in the literature, need a substantial revision, mainly in relation to the predicted chemical yields. The interplay among convection, hot bottom burning and the low-temperature opacity treatment is also discussed

We study the spectroscopic and photometric properties of carbon stars. In the first paper of this series we focus on objects that can be described by hydrostatic models neglecting dynamical phenomena like pulsation and mass loss. As a consequence, the reddening due to circumstellar dust is not included. Our results are collected in a database, which can be used in conjunction with stellar evolution and population synthesis calculations involving the AGB. We have computed a grid of 746 spherically symmetric COMARCS atmospheres covering effective temperatures between 2400 and 4000K, surface gravities from log(g) = 0.0 to -1.0, metallicities ranging from the solar value down to one tenth of it and C/O ratios in the interval between 1.05 and 5.0. Subsequently, we used these models to create synthetic low resolution spectra and photometric data for a large number of filter systems. The tables including the results are electronically available. We have selected some of the most commonly used colours in order to discuss their behaviour as a function of the stellar parameters. A comparison with measured data shows that down to 2800K the agreement between predictions and observations of carbon stars is good. Below this limit the synthetic colours are much too blue. The obvious reason for these problems is the neglect of circumstellar reddening and structural changes due to pulsation and mass loss.

MASYS is the AKARI spectroscopic survey of Symbiotic Stars in the Magellanic Clouds, and one of the European Open Time Observing Programmes approved for the AKARI (Post-Helium) Phase-3. It is providing the first ever near-IR spectra of extragalactic symbiotic stars. The observations are scheduled to be completed in July 2009.

We study the effects of WIMP dark matter (DM) annihilations on the thermal and chemical evolution of the gaseous clouds where the first generation of stars in the Universe is formed. We follow the collapse of the gas inside a typical halo virializing at very high redshift, from well before virialization until a stage where the heating from DM annihilations exceeds the gas cooling rate. The DM energy input is estimated by inserting the energy released by DM annihilations (as predicted by an adiabatic contraction of the original DM profile) in a spherically symmetric radiative transfer scheme. In addition to the heating effects of the energy absorbed, we include its feedback upon the chemical properties of the gas, which is critical to determine the cooling rate in the halo, and hence the fragmentation scale and Jeans mass of the first stars. We find that DM annihilation does alter the free electron and especially the H2 fraction when the gas density is n>~ 10^4 cm^-3, for our fiducial parameter values. However, even if the change in the H2 abundance and the cooling efficiency of the gas is large (sometimes exceeding a factor 100), the effects on the temperature of the collapsing gas are far smaller (a reduction by a factor <~1.5), since the gas cooling rate depends very strongly on temperature: then, the fragmentation mass scale is reduced only slightly, hinting towards no dramatic change in the initial mass function of the first stars.

We study the evolution of the first stars in the universe (Population III) from the early pre-Main Sequence until the end of helium burning in the presence of WIMP dark matter annihilation inside the stellar structure. The two different mechanisms that can provide this energy source are the contemporary contraction of baryons and dark matter, and the capture of WIMPs by scattering off the gas with subsequent accumulation inside the star. We find that the first mechanism can generate an equilibrium phase, previously known as a "dark star", which is transient and present in the very early stages of pre-MS evolution. The mechanism of scattering and capture acts later, and can support the star virtually forever, depending on environmental characteristic of the dark matter halo and on the specific WIMP model.

We study the effects of WIMP dark matter (DM) on the collapse and evolution of the first stars in the Universe. Using a stellar evolution code, we follow the pre-Main Sequence (MS) phase of a grid of metal-free stars with masses in the range 5-600 solar mass forming in the centre of a 1e6 solar mass halo at redhisft z=20. DM particles of the parent halo are accreted in the proto-stellar interior by adiabatic contraction and scattering/capture processes, reaching central densities of order 1e12 GeV/cm3 at radii of the order of 10 AU. Energy release from annihilation reactions can effectively counteract the gravitational collapse, in agreement with results from other groups. We find this stalling phase (known as "dark" star) is transients and lasts from 2.1e3 yr (M=600 solar mass) to 1.8e4 yr (M=9 solar mass). Later in the evolution, DM scattering/capture rate becomes high enough that energy deposition from annihilations significantly alters the pre-MS evolution of the star in a way that depends on DM (i) velocity dispersion, (ii) density, (iii) elastic scattering cross section with baryons. For our fiducial set of parameters (10 km/s, 1e11 GeV/cm3, 1e-38 cm2) we find that the evolution of stars of mass lower than 40 solar masses "freezes" on the HR diagram before reaching the ZAMS. Stars with bigger masses manage to ignite nuclear reactions; however, DM "burning" prolonges their lifetimes by a factor 2 (5) for a 600 (40) solar mass star.

We present a study of the evolved stellar populations in the dwarf spheroidal galaxy LeoII, based on JHKs observations obtained with the near-infrared array WFCAM at the UKIRT telescope. Combining the new data with optical data, we derived photometric estimates of the distribution of global metallicity [M/H] of individual red giant stars from their V-Ks colours. We used the new data to derive the properties of a nearly complete sample of asymptotic giant branch (AGB) stars in LeoII. Using a near-infrared two-colour diagram, we were able to obtain a clean separation from Milky Way foreground stars and discriminate between carbon- and oxygen-rich AGB stars. We simulate the JHKs data with the TRILEGAL population synthesis code together with the most updated thermally pulsing AGB models, and using the star formation histories derived from independent work based on deep HST photometry. After scaling the mass of LeoII models to the observed number of upper RGB stars, we find that present models predict too many O-rich TP-AGB stars of higher luminosity due to a likely under-estimation of either their mass-loss rates at low metallicity, and/or their degree of obscuration by circumstellar dust. On the other hand, the TP-AGB models are able to reproduce the observed number and luminosities of carbon stars satisfactorily well, indicating that in this galaxy the least massive stars that became carbon stars should have masses as low as ~1 M_sun.

Tracks and isochrones have been computed in the range of initial masses 0.15 - 20 Mo for a grid of 39 chemical compositions with the metal content Z between 0.0001 and 0.070, and helium content Y between 0.23 and 0.46. The Padova stellar evolution code has been implemented with updated physics. New synthetic TP-AGB models allow the extension of stellar models and isochrones until the end of the thermal pulses along the AGB. Software tools for the bidimensional interpolation (in Y and Z) of the tracks have been tuned. This first paper presents tracks for low mass stars (from 0.15 to 2.5 Mo) with scaled-solar abundances and the corresponding isochrones from very old ages down to about 1 Gyr. Tracks and isochrones are made available in tabular form for the adopted grid of chemical compositions in the plane Z-Y. An interactive web interface will allow users to obtain isochrones of whatever chemical composition and also simulated stellar populations with different Y(Z) helium-to-metal enrichment laws.

We present a large set of theoretical isochrones, whose distinctive features mostly reside on the greatly improved treatment of the thermally pulsing asymptotic giant branch (TP-AGB) phase. Essentially, we have coupled the TP-AGB tracks described in Paper I, at their stages of pre-flash quiescent H-shell burning, with the evolutionary tracks for the previous evolutionary phases from Girardi et al. (2000). Theoretical isochrones for any intermediate value of age and metallicity are then derived by interpolation in the grids. We take care that the isochrones keep, to a good level of detail, the several peculiarities present in these TP-AGB tracks. Theoretical isochrones are then converted to about 20 different photometric systems -- including traditional ground-based systems, and those of recent major wide-field surveys such as SDSS, OGLE, DENIS, 2MASS, UKIDSS, etc., -- by means of synthetic photometry applied to an updated library of stellar spectra, suitably extended to include C-type stars. Finally, we correct the predicted photometry by the effect of circumstellar dust during the mass-losing stages of the AGB evolution, which allows us to improve the results for the optical-to-infrared systems, and to simulate mid- and far-IR systems such as those of Spitzer and AKARI. Access to the data is provided both via a web repository of static tables (http://stev.oapd.inaf.it/dustyAGB07 and CDS), and via an interactive web interface (http://stev.oapd.inaf.it/cmd) that provides tables for any intermediate value of age and metallicity, for several photometric systems, and for different choices of dust properties.

We present new synthetic models of the TP-AGB evolution. They are computed for 7 values of initial metal content (Z from 0.0001 to 0.03) and for initial masses between 0.5 and 5.0 Msun, thus extending the low- and intermediate-mass tracks of Girardi et al. (2000) until the beginning of the post-AGB phase. The calculations are performed by means of a synthetic code that incorporates many recent improvements, among which we mention: (1) the use of detailed and revised analytical relations to describe the evolution of quiescent luminosity, inter-pulse period, third dredge-up, hot bottom burning, pulse cycle luminosity variations, etc.; (2) the use of variable molecular opacities -- i.e. opacities consistent with the changing photospheric chemical composition -- in the integration of a complete envelope model, instead of the standard choice of scaled-solar opacities; (3) the use of formalisms for the mass-loss rates derived from pulsating dust-driven wind models of C- and O-rich AGB stars; and (4) the switching of pulsation modes between the first overtone and the fundamental one along the evolution, which has consequences in terms of the history of mass loss. It follows that, in addition to the time evolution on the HR diagram, the new models predict in a consistent fashion also variations in surface chemical compositions, pulsation modes and periods, and mass-loss rates. The onset and efficiency of the third dredge-up process are calibrated in order to reproduce basic observables like the carbon star luminosity functions in the Magellanic Clouds, and TP-AGB lifetimes (star counts) in Magellanic Cloud clusters. Forthcoming papers will present the theoretical isochrones and chemical yields derived from these tracks, and additional tests performed with the aid of a complete population synthesis code.

Marigo (2002) has highlighted the crucial importance of molecular opacities in modelling the evolution of AGB stars at varying surface C/O ratio. In particular, it has been shown the large inadequacy of solar-scaled opacities when applied to models of carbon stars, and hence the need for correctly coupling the molecular opacities to the current surface chemical composition of AGB stars. The aim of the present follow-up study is to investigate the effects of variable molecular opacities on the evolutionary properties of luminous AGB stars with massive envelopes, i.e. with initial masses from ~3.5 Msun up to 5-8 Msun, which are predicted to experience both the third dredge-up and hot-bottom burning. It is found that if the dredge-up of carbon is efficient enough to lead to an early transition from C/O<1 to C/O>1, then hot-bottom burning may be weakened, extinguished, or even prevented. The physical conditions for this occurrence are analysed and a few theoretical and observational implications are discussed. Importantly, it is found that the inclusion of variable molecular opacities could significantly change the current predictions for the chemical yields contributed by intermediate-mass AGB stars, with M~3.5 - 4.0 Msun that make as much as ~ 30-50 % of all stars expected to undergo hot-bottom burning.

I will briefly review the state of the art of evolutionary population synthesis (EPS) models that include the contribution from AGB stars.

Extended and updated grids of TP-AGB tracks have been implemented in the TRILEGAL population synthesis code, which generates mock stellar catalogues for a galaxy given its mass, distance, star formation history and age-metallicity relation, including also the Milky Way foreground population. Among the stellar parameters that are simulated, we now include the surface chemistry, mass-loss rates, pulsation modes and periods of LPVs. This allows us to perform a series of consistency checks between AGB model predictions and observations, that we are just starting to explore. We present a few examples of model--data comparisons, mostly regarding the near-infrared and variability data for AGB stars in the Magellanic Clouds.

Using available data for C and M giants with Mbol<-3.6 in Magellanic Cloud clusters, we derive limits to the lifetimes of the corresponding evolutionary phases, as a function of stellar mass. The C-star phase is found to have a duration between 2 and 3 Myr for stars in the mass range from 1.5 to 2.8 Msun. There is also an indication that the peak of C-star lifetime shifts to lower masses (from slightly above to slightly below 2 Msun) as we move from LMC to SMC metallicities. The M-giant lifetimes also peak at 2 Msun in the LMC, with a maximum value of about 4 Myr, whereas in the SMC their lifetimes appear much shorter but, actually, they are poorly constrained by the data. These numbers constitute useful constraints to theoretical models of the TP-AGB phase. We show that several models in the literature underestimate the duration of the C-star phase at LMC metallicities.

This article compares the distribution of Ks magnitude of Large Magellanic Cloud (LMC) asymptotic giant branch (AGB) stars obtained from the DENIS and 2MASS data with theoretical distributions. These have been constructed using up-to-date stellar evolution calculations for low and intermediate-mass stars, and in particular for thermally pulsing AGB stars. A fit of both the magnitude distribution of carbon- and oxygen-rich AGB stars allowed us to constrain the metallicity distribution across the LMC and its star formation rate (SFR). The LMC stellar population is found to be on average 5-6 Gyr old and is consistent with a mean metallicity corresponding to Z=0.006. These values may however be affected by systematic errors in the underlying stellar models, and by the limited exploration of the possible SFR histories. Our method should be, instead, particularly useful for detecting variations in the mean metallicity and SFR across the LMC disk. In fact, there are well defined regions where both the metallicity and the mean-age of the underlying stellar population span the whole range of grid parameters. The C/M ratio discussed in paper I is a tracer of the metallicity distribution if the underlying stellar population is older than about a few Gyr. A similar study across the Small Magellanic Cloud is given in paper III of this series.

This article compares the Ks magnitude distribution of Small Magellanic Cloud asymptotic giant branch stars obtained from the DENIS and 2MASS data with theoretical distributions. Theoretical Ks magnitude distributions have been constructed using up-to-date stellar evolution calculations for low and intermediate-mass stars, and in particular for thermally pulsing asymptotic giant branch stars. Separate fits of the magnitude distributions of carbon- and oxygen-rich stars allowed us to constrain the metallicity distribution across the galaxy and its star formation rate. The Small Magellanic Cloud stellar population is found to be on average 7-9 Gyr old but older stars are present at its periphery and younger stars are present in the direction of the companion galaxy the Large Magellanic Cloud. The metallicity distribution traces a ring-like structure that is more metal rich than the inner region of the galaxy. The C/M ratio discussed in Paper I is a tracer of the metallicity distribution only if the underlying stellar population is of intermediate-age.

We investigate the bright cut-off of the [OIII]l5007 planetary nebula luminosity function (PNLF), that has been suggested as a powerful extragalactic distance indicator. Theoretical PNLFs are constructed via Monte-Carlo simulations of populations of PNe, whose individual properties are described with the aid of recent PN synthetic models (Marigo et al. 2001), coupled to a detailed photoionisation code (CLOUDY). The basic dependences of the cut-off magnitude M* are then discussed. We find that: (i) In galaxies with recent or ongoing star formation, the modelled PNLF present M* values between -4 and -5, very close to the observationally-calibrated value for the LMC. (ii) In these galaxies, the PNLF cut-off is produced by PNe with progenitor masses of about 2.5 Msun, while less massive stars give origin to fainter PNe. As a consequence M* is expected to depend strongly on the age of the last burst of star formation, dimming by as much as 5 mag as we go from young to 10-Gyr old populations. (iii) Rather than on the initial metallicity of a stellar population, M* depends on the actual [O/H] of the observed PNe, a quantity that may differ significantly from the initial value (due to dredge-up episodes), especially in young and intermediate-age PN populations. (iv) Also the transition time from the end of AGB to the PN phase, and the nuclear-burning properties (i.e. H- or He-burning) of the central stars introduce non-negligible effects on M*. The strongest indication derived from the present calculations is a serious difficulty to explain the age-invariance of the cut-off brightness over an extended interval, say from 1 to 13 Gyr, that observations of PNLFs in galaxies of late-to-early type seem to suggest.

The elemental abundances of ten planetary nebulae, derived with high accuracy including ISO and IUE spectra, are analysed with the aid of synthetic evolutionary models for the TP-AGB phase. Model prescriptions are varied until we achieve the simultaneous reproduction of all elemental features, which allows placing important constraints on the characteristic masses and nucleosynthetic processes experienced by the stellar progenitors. First of all, it is possible to separate the sample into two groups of PNe, one indicating the occurrence of only the third dredge-up during the TP-AGB phase, and the other showing also the chemical signature of hot-bottom burning. The former group is reproduced by stellar models with variable molecular opacities (see Marigo 2002), adopting initial solar metallicity, and typical efficiency of the third dredge-up 0.3-0.4. The latter group of PNe, with extremely high He content 0.15<=He/H<=0.20 and marked oxygen deficiency, is consistent with original sub-solar metallicity (i.e. LMC composition). Moreover, we are able to explain quantitatively both the N/H-He/H correlation and the N/H-C/H anti-correlation, thus solving the discrepancy pointed out long ago by Becker & Iben (1980). This is obtained only under the hypothesis that intermediate-mass TP-AGB progenitors (M > 4.5-5.0 Mo) with LMC composition have suffered a number of very efficient, carbon-poor, dredge-up events. Finally, the neon abundances of the He-rich PNe can be recovered by invoking a significant production of 22Ne during thermal pulses, which would imply a reduced role of the 22Ne(alpha,n)25Mg reaction as neutron source to the s-process nucleosynthesis in these stars.

Carbon stars are known to exhibit systematically redder near-infrared colours with respect to M-type stars. In the near-infrared colour-magnitude diagrams provided by the 2MASS and DENIS surveys, the LMC C-type stars draw a striking red tail, well separated from the sequences of O-rich giants. So far, this conspicuous feature has been absent from any set of available isochrones, even the few existing ones that include the TP-AGB evolution of low- and intermediate-mass stars. To investigate such issue we simulate the complete 2MASS Ks vs.(J-Ks) data towards the LMC by means of a population synthesis approach, that relies on extended libraries of published stellar evolutionary tracks, including the TP-AGB phase. The simulations provide quite a detailed description of the several vertical fingers and inclined sequences seen in 2MASS data, due to both Galactic foreground and LMC O-rich stars. Instead, as mentioned, the red tail of C-stars sets a major difficulty: we find that TP-AGB models with solar-scaled molecular opacities, the usual assumption of existing AGB calculations, do not succeed in reproducing this feature. Our tests indicate that the main reason for this failure should not be ascribed to empirical Teff - (J-K) transformations for C-type stars. Instead, the discrepancy is simply removed by adopting new evolutionary models that account for the changes in molecular opacities as AGB stars get enriched in carbon via the third dredge-up (Marigo 2002). In fact, simulations that adopt these models are able to reproduce, for the first time, the red tail of C-stars in near-infrared CMDs. Finally, we point out that these simulations also provide useful indications about the efficiency of the third dredge-up process, and the pulsation modes of long-period variables.

We use the available data for Magellanic Cloud carbon stars to constrain the efficiency of the third dredge-up process in TP-AGB models. We show that star counts in LMC clusters provide quite stringent limits to the lifetime of the C-star phase, with a duration between 2 and 3 Myr for stars in the mass range from 1.5 to 2.8 Msun. Together with the luminosity functions of field C stars, this information allows us to re-calibrate the third dredge-up parameters log T_b_dred and lambda in TP-AGB models that include variable molecular opacities (Marigo 2002). Preliminary results are presented here.

This study calls attention to the importance of properly coupling the molecular opacities to the actual surface abundances of TP-AGB stars that experience the third dredge-up and/or hot-bottom burning, i.e. with surface abundances of carbon and oxygen varying with time. New TP-AGB calculations with variable opacities -- replacing the usually adopted solar-scaled opacity tables -- have proven to reproduce, for the first time, basic observables of carbon stars, like their effective temperatures, C/O ratios, and near-infrared colours. Moreover, it turns out that the effect of envelope cooling -- due to the increase in molecular opacities -- may cause other important effects, namely: i) shortening of the C-star phase; ii) possible reduction or shut-down of the third dredge-up in low-mass carbon stars; and iii) weakening or even extinction of hot-bottom burning in intermediate-mass stars.

We present evolutionary models of zero-metallicity very massive objects, with initial masses in the range 120 Msun -- 1000 Msun, covering their quiescent evolution up to central carbon ignition. In the attempt of exploring the possible occurrence of mass loss by stellar winds, calculations are carried out with recently-developed formalisms for the mass-loss rates driven by radiation pressure (Kudritzki 2002) and stellar rotation (Maeder & Meynet 2000).The study completes the previous analysis by Marigo et al. (2001) on the constant-mass evolution of primordial stars. Our results indicate that radiation pressure (assuming a minimum metallicity Z = 10^{-4} Zsun)is not an efficient driving force of mass loss, except for very massive stars with M >= 750 Msun. On the other hand, stellar rotation might play a crucial role in triggering powerful stellar winds, once the (Omega-Gamma)-limit is approached. However, this critical condition of intense mass loss can be maintained just for short, as the loss of angular momentum due to mass ejection quickly leads to the spinning down of the star. As by-product to the present work, the wind chemical yields from massive zero-metallicity stars are presented. The helium and metal enrichments, and the resulting Delta(Y)/Delta(Z) ratio are briefly discussed.

We discuss the evolutionary properties of primordial massive and very massive stars, supposed to have formed from metal-free gas. Stellar models are presented over a large range of initial masses (8 Msun <= Mi <= 1000 Msun), covering the hydrogen- and helium-burning phases up to the onset of carbon burning. In most cases the evolution is followed at constant mass. To estimate the possible effect of mass loss via stellar winds, recent analytic formalisms for the mass-loss rates are applied to the very massive models (Mi >= 120 Msun).

We present results of a new investigation (Carraro et al. 2002) aimed at clarifying the mutual relationship between the most prominent young open clusters close to eta Carinae, namely Trumpler 16, Trumpler 14 and Collinder 232.

We provide tables of theoretical isochrones in several photometric systems. To this aim, the following steps are followed: (1) First, we re-write the formalism for converting synthetic stellar spectra into tables of bolometric corrections. The resulting formulas can be applied to any photometric system, provided that the zero-points are specified by means of either ABmag, STmag, VEGAmag, or a standard star system that includes well-known spectrophotometric standards. Interstellar absorption can be considered in a self-consistent way. (2) We assemble an extended and updated library of stellar intrinsic spectra. It is mostly based on non-overshooting ATLAS9 models, suitably extended to both low and high effective temperatures. This offers an excellent coverage of the parameter space of Teff, logg, and [M/H]. We briefly discuss the main uncertainties and points still deserving more improvement. (3) From the spectral library, we derive tables of bolometric corrections for Johnson-Cousins-Glass, HST/WFPC2, HST/NICMOS, Washington, and ESO Imaging Survey systems (this latter consisting on the WFI, EMMI, and SOFI filter sets). (4) These tables are used to convert several sets of Padova isochrones into the corresponding absolute magnitudes and colours, thus providing a useful database for several astrophysical applications. All data files are made available in electronic form.

We investigate the effects of molecular opacities on the evolution of TP-AGB stars that experience the third dredge-up, i.e. with surface abundances of carbon and oxygen varying with time. To this aim, a routine is constructed to derive the molecular concentrations through dissociation equilibrium calculations, and estimate the opacities due to H2, H2O, OH, C2, CN, and CO for any given density, temperature, and chemical composition of the gas. Then, synthetic TP-AGB models with dredge-up are calculated by either 1) adopting the newly developed routine, or 2) interpolating between fixed opacity tables for solar chemical composition, which is the usual prescription in AGB evolution calculations. The comparison between the two cases shows that the change in the dominant opacity sources, as the C/O ratio grows from below to above unity, crucially affects the evolution of the effective temperature, i.e. causing a notable cooling of the carbon-rich models (with C/O>1). From the comparison with observational data, it turns out that TP-AGB models with variable molecular opacities are able to reproduce the observed range of effective temperatures, mass-loss rates, and wind expansion velocities of C-type giants in the solar neighbourhood, otherwise failed if assuming fixed molecular opacities for solar-scaled mixtures. Finally, we mention other important evolutionary and observational effects that result from the adoption of the variable opacities, such as: i) significant shortening of the C-star phase due to the earlier onset of the super-wind; ii) consequent reduction of the carbon yields; iii) reproduction of the observed range of near-infrared colours of C-stars.

We report on $UBVRI$ CCD photometry of two overlapping fields in the region of the intermediate-age open cluster NGC 2158 down to V=21. By analyzing Colour-Colour (CC) and Colour-Magnitude Diagrams (CMD) we infer a reddening $E_{B-V}= 0.55\pm0.10$, a distance of $3600 \pm 400$ pc, and an age of about 2 Gyr. Synthetic CMDs performed with these parameters (but fixing $E_{B-V}=0.60$ and $[{\rm Fe/H}]=-0.60$), and including binaries, field contamination, and photometric errors, allow a good description of the observed CMD. The elongated shape of the clump of red giants in the CMD is interpreted as resulting from a differential reddening of about $\Delta E_{B-V}=0.06$ across the cluster, in the direction perpendicular to the Galactic plane. NGC 2158 turns out to be an intermediate-age open cluster with an anomalously low metal content. The combination of these parameters together with the analysis of the cluster orbit, suggests that the cluster belongs to the old thin disk population.

We present a new synthetic model to follow the evolution of a planetary nebula (PN) and its central star, starting from the onset of AGB phase up to the white dwarf cooling sequence. The model suitably combines various analytical prescriptions to account for different (but inter-related) aspects of planetary nebulae, such as: the dynamical evolution of the primary shell and surrounding ejecta, the photoionisation of H and He by the central star, the nebular emission of a few relevant optical lines (e.g. Hbeta; HeII4686; [OIII]5007). Predictions of the synthetic model are tested by comparison with both findings of hydrodynamical calculations, and observations of Galactic PNe. The sensitiveness of the results to the models parameters (e.g. transition time, mass of the central star, H-/He-burning tracks, etc.) is also discussed. We briefly illustrate the systematic differences that are expected in the luminosities and lifetimes of PNe with either H- or He-burning central stars, which result in different ``detection probabilities'' across the H-R diagram, in both Hbeta and [OIII]5007 lines. Adopting reasonable values of the model parameters, we are able to reproduce, in a satisfactory way, many general properties of PNe, like the ionised mass--nebular radius relationship, the trends of a few main nebular line ratios, and the observed ranges of nebular shell thicknesses, electron densities, and expansion velocities. The models naturally predict also the possible transitions from optically-thick to optically-thin configurations (and vice versa). In this context, the origin of the Zanstra discrepancy is also analysed. (abridged)

In this paper we emphasize the close connection between the chemical and spectrophotometric evolution of stellar systems: Chemical yields from stars correspond to a precise fraction of their emitted light. We translate this concept quantitatively. Starting from simple stellar populations, we derive useful analytical relations to calculate the stellar fuel consumption (emitted light) as a function of basic quantities predicted by stellar models, i.e. the mass of the core and the chemical composition of the envelope. The final formulas explicate the relation between integrated light contribution (total or limited to particular evolutionary phases), chemical yields and stellar remnants. We test their accuracy in the case of low- and intermediate-mass stars, and indicate the way to extend the analysis to massive stars. This formalism provides an easy tool to check the internal consistency between the different stellar inputs adopted in galaxy models: The fuel computed by means of the analytical formulas (corresponding to a given set of chemical yields) should be compared to the exact values given by the luminosity integration along the stellar evolutionary tracks or isochrones (corresponding to a given set of spectrophotometric models). Only if both estimates of the fuel are similar, the stellar inputs can be considered self-consistent in terms of their energetics. This sets an important requirement to galaxy models, also in consideration of the fact that different sources of input stellar data are frequently used to model their spectro-photometric and chemical evolution.

We present extensive evolutionary models of stars with initial zero-metallicity, covering a large range of initial masses (i.e. 0.7 <= M <= 100 Msun). Calculations are carried out at constant mass, with updated input physics, and applying an overshooting scheme to convective boundaries. The nuclear network includes all the important reactions of the p-p chain, CNO-cycle and alpha-captures, and is solved by means of a suitable semi-implicit method. The evolution is followed up to the thermally pulsing AGB in the case of low- and intermediate-mass stars, or to the onset of carbon burning in massive stars. The main evolutionary features of these models are discussed, also in comparison with models of non-zero metallicity. Among several interesting aspects, particular attention has been paid to describe: i) the first synthesis of 12C inside the stars, that may suddenly trigger the CNO-cycle causing particular evolutionary features; ii) the pollution of the stellar surface by the dredge-up events, that are effective only within particular mass ranges; iii) the mass limits which conventionally define the classes of low-, intermediate-, and high-mass stars on the basis of common evolutionary properties, including the upper mass limit for the achievement of super-Eddington luminosities before C-ignition in the high-mass regime; and iv) the expected pulsational properties of zero-metallicity stars. All relevant information referring to the evolutionary tracks and isochrones is made available in computer-readable format at http://pleiadi.pd.astro.it .

In this work we analyse the role of low- and intermediate-mass stars in contributing to the chemical enrichment of the interstellar medium. First we present new sets of stellar yields basing on the results of updated evolutionary calculations, which extend from the ZAMS up to the end of the AGB phase (Girardi et al. 2000; Marigo et al. 1999a). These new yields, that present a significant dependence on metallicity, are then compared to those of other available sets (Renzini & Voli 1981; van de Hoek & Groenewegen 1997). The resulting differences are explained in terms of different model assumptions -- i.e. treatment of convective boundaries, mass loss, dredge-up, hot-bottom burning --, and further discussed on the basis of important empirical constraints which should be reproduced by theory -- i.e. the initial-final mass relation, white dwarf mass distribution, carbon star luminosity function, and chemical abundances of planetary nebulae. We show that present models are able to reproduce such constraints in a satisfactory way.

The luminosity evolution of stars with highly condensed cores surrounded by nuclear-burning shell(s) is analytically investigated with the aid of homology relations. With respect to earlier works using a similar approach (e.g. Refsdal & Weigert 1970; Kippenhahn 1981), the major improvement is that we derive all the basic dependences (i.e. on core mass, core radius, and chemical composition) in a completely generalised fashion, then accounting for a large range of possible physical properties characterising the burning shell(s). Parameterised formulas for the luminosity are given as a function of the (i) relative contribution of the gas to the total pressure (gas plus radiation), (ii) opacity source, and (iii) dominant nuclear reaction rates. In this way, the same formalism can be applied to shell-burning stars of various metallicities and in different evolutionary phases. In particular, we present some applications concerning the luminosity evolution of RGB and AGB stars with different chemical compositions, including the case of initial zero metallicity. It turns out that homology predictions provide a good approximation to the results of stellar model calculations. Therefore, the proposed formalism is useful to understand the possible differences in the luminosity evolution of shell-burning stars within a unified interpretative framework, and can be as well adopted to improve the analytical description of stellar properties in synthetic models.

We present new sets of chemical yields from low- and intermediate-mass stars with 0.8 Msun <= M <= Mup ~ 5 Msun, and three choices of the metallicity, Z=0.02, Z=0.008, and Z=0.004 (Marigo 2000, in preparation). These are then compared with the yields calculated by other authors on the basis of different model prescriptions, and basic observational constraints which should be reproduced.

We present new evolutionary models for zero-metallicity stars, covering a large range of initial masses (from 0.8 to 100 M_sun). Models are computed with overshooting from stellar cores and convective envelopes, and assuming mass-loss from the most massive stars. We discuss the main evolutionary features of these stars, and provide estimates of the amount of newly-synthesized elements dredged-up to the stellar surface, and possibly lost by stellar winds from the most massive stars. Full details of these models will be given in Marigo et al. (2000, in preparation).

The core mass-luminosity (Mc-L) relation of TP-AGB stars is a key ingredient in synthetic calculations of their evolution. Recently, Herwig et al. (1998) have presented full calculations of TP-AGB models with strong dredge-up occurring already during early thermal pulses. The resulting luminosity evolution differs appreciably from the simple linear Mc-L relation. In this paper, we show that at least part of the luminosity evolution can be understood as being the result of two well-known effects: the gradual approach to the asymptotic behaviour that characterises the first thermal pulses, and the chemical composition changes of the envelope. Both effects are already implemented in the Mc-L relations used in synthetic models. Consequently, these models are able to reproduce the behaviour of full calculations. Whether additional effects, not yet taken into account, are present, can be decided only through additional calculations and data. We also comment on the validity of a linear Mc-L relation and its possible violation, as mentioned by Herwig et al.

We investigate the formation of carbon stars as a function of the stellar mass and parent metallicity. Theoretical modelling is based on an improved scheme for treating the third dredge-up in synthetic calculations of thermally pulsing asymptotic giant branch (TP-AGB) stars. In this approach, the usual criterion (based on a constant minimum core mass for the occurrence of dredge-up, M_c^min) is replaced by one on the minimum temperature at the base of the convective envelope, T_b^dred, at the stage of the post-flash luminosity maximum. Envelope integrations then allow determination of M_c^min as a function of stellar mass, metallicity, and pulse strength (see Wood 1981), thus inferring if and when dredge-up first occurs. Moreover, the final possible shut down of the process is predicted. Extensive grids of TP-AGB models were computed using this scheme. We present and discuss the calibration of the two dredge-up parameters (lambda and T_b^dred) aimed at reproducing the carbon star luminosity function (CSLF) in the LMC. It turns out that the faint tail is almost insensitive to the history of star formation rate (SFR) in the parent galaxy (it is essentially determined by T_b^dred), in contrast to the bright wing which may be more affected by the details of the recent SFR. Once the faint end is reproduced, the peak location is a stringent calibrator of lambda. The best fit to the observed CSLF in the LMC is obtained with Z=0.008, lambda=0.50, log(T_b^dred)=6.4, and a constant SFR up to 5x10^8 yr ago. A good fit to the CSLF in the SMC is then easily derived from the Z=0.004 models, with a single choice of parameters, and a constant SFR over the entire significant age interval. The results are consistent with the theoretical expectation that the third dredge-up is more efficient at lower Zs.

Until recently synthetic AGB models had not taken into account the break-down of the core mass-luminosity (Mc-L) relation due to the occurrence of envelope burning in the most massive (M > 3.5 Msun) and luminous (Mbol > -6) stars. Marigo et al. (1998) made the first attempt to consistently include the related over-luminosity effect (i.e. above the Mc-L relation) in synthetic TP-AGB calculations. In this paper the reliability of the solution scheme is tested by comparison with the results of complete evolutionary calculations for a 7 Msun AGB star undergoing envelope burning (e.g. Bloecker & Schoenberner 1991). Indeed, the method proves to be valid as it is able to reproduce with remarkable accuracy several evolutionary features of the 7 Msun star. We present extensive synthetic TP-AGB calculations for stars with initial masses of 3.5, 4.0, 4.5, and 5.0 Msun, and three choices of the initial metallicity, i.e. Z=0.019, Z=0.008, and Z=0.004. Three values of the mixing-length parameter are used, i.e. alpha=1.68, 2.0, 2.5. We investigate the dependence of envelope burning on such stellar parameters (M, Z, and alpha). The comparison between different cases gives hints on the interplay between envelope burning over-luminosity and mass loss, and related effects on TP-AGB lifetimes.

In this study we present a simple model of elliptical galaxies aimed at interpreting the gradients in colours and narrow band indices observed across these systems. Salient features of the model are the gradients in mass density and star formation and infall of primordial gas aimed at simulating the collapse of a galaxy into the potential well of dark matter. Adopting a multi-zone model we follow in detail the history of star formation, gas consumption, and chemical enrichment of the galaxy and also allow for the occurrence of galactic winds according to the classical supernova (and stellar winds) energy deposit. The outline of the model, the time scale of gas accretion and rate of star formation as a function of the galacto-centric distance in particular, seek to closely mimic the results from Tree-SPH dynamical models. Although same specific ingredients of the model can be questioned from many points of view (of which we are well aware) the model predictions have to be considered as a gross tool for exploring the consequences of different receipts of gas accretion and star formation in which the simple one-zone scheme is abandoned. With the aid of this model we discuss the observational data on the gradients in metallicity, colour, and narrow band indices across elliptical galaxies.

In this paper we focus on the TP-AGB evolution of intermediate-mass stars experiencing envelope burning (M=4-5Mo). The approach we have adopted is a semi-analytical one as it combines analytical relationships derived from complete models of TP-AGB stars with sole envelope models in which the physical structure is calculated from the photosphere down to the core. Our efforts are directed to analyse the effects produced by envelope burning, such as: i) the energy contribution which may drive significant deviations from the standard core mass-luminosity relationship; and ii) the changes in the surface chemical composition due to nuclear burning via the CNO cycle. Evolutionary models for stars with initial mass of 4.0, 4.5, 5.0 Mo and two choices of the initial chemical composition ([Y=0.28,Z=0.02] and [Y=0.25,Z=0.008]) are calculated from the first thermal pulse till the complete ejection of the envelope. We find that massive TP-AGB stars can rapidly reach high luminosities (-6>Mbol>-7), without exceeding, however, the classical limit to the AGB luminosity of Mbol~-7.1 corresponding to the Chandrasekhar value of the core mass. No carbon stars brighter than Mbol~-6.5 are predicted to form (the alternative of a possible transition from M-star to C-star during the final pulses is also explored), in agreement with observations which indicate that most of the very luminous AGB stars are oxygen-rich. Finally, new chemical yields from stars in the mass range 4-5 Mo are presented, so as to extend the sets of stellar yields from low-mass stars already calculated by Marigo et al. (1996). For each CNO element we give both the secondary and the primary components.