The chemical enrichment in the interstellar medium (ISM) of galaxies is regulated by several physical processes: stellar evolution, grain formation and destruction, galactic inflows and outflows. Understanding such processes is essential to follow the chemical enrichment of galaxies through the cosmic epochs, and to interpret the observations. Despite the importance of such topics, the efficiency of the different processes driving the evolution of baryons in galaxies, remain controversial. We revise the current description of metal and dust evolution in local low-metallicity dwarf galaxies and we develop a description for Lyman Break Galaxies. Our main goal is to reproduce i) the peak in the mass of dust over the mass of stars (sMdust) observed within few hundred Myrs; ii) the decrease of the sMdust at later time. The spectral energy distribution of the galaxies is fitted with the "Code Investigating GALaxies Emission" (CIGALE), through which the stellar and dust masses, and the star formation rate are estimated. For some of the dwarf galaxies, the metal and gas content are also available. We run different calculations of chemical evolution in galaxies, and we fit the observed properties through the model predictions. We show that i) a top-heavy initial mass function that favours massive stars and a dust condensation fraction for Type II Supernovae (SNe II) of 50% or more help to reproduce the peak of sMdust observed after 100 Myrs since the beginning of the cycle; ii) galactic outflows play a crucial role in reproducing the decline in sMdust with age, and they are more efficient than grain destruction from SNe II; iii) a star formation efficiency (mass of gas converted into stars) of few per cent is required to explain the metallicity of local dwarf galaxies; iv) dust growth in the ISM is not necessary to reproduce the sMdust and, if present, its effect is erased by galactic outflows.

(abridged)Variability is a key property of stars on the asymptotic giant branch (AGB). Their pulsation period is related to the luminosity and mass-loss rate of the star. The long-period variables (LPVs) and Mira variables are the most prominent of all types of variability of evolved stars. The reddest, most obscured AGB stars are too faint in the optical and have eluded large variability surveys. Selection criteria are derived based on colour-colour and colour-magnitude diagrams from the combination of VISTA Magellanic Cloud (VMC) survey, Spitzer IRAC and AllWISE data. After eliminating LPVs with known periods shorter than 450 days, a sample of 1299 candidate obscured AGB stars is selected. K-band LCs are constructed combining the epoch photometry available in the VMC survey with literature data, analysed for variability and fitted with a single period sine curve to derive mean magnitudes, amplitudes and periods. A subset of 254 stars are either new variables, or known variables where the period we find is better determined than the literature value, or longer than 1000 days. The spectral energy distributions (SEDs) of these stars are fitted to a large number of templates. For this purpose the SEDs and Spitzer IRS spectra of some non-AGB stars are also fitted to have templates of the most likely contaminants in the sample. A sample of 217 likely LPVs is found. Thirty-four stars have periods longer than 1000 days although some of them have alternative shorter periods. Mass-loss rates and luminosities are estimated from the template fitting. Period-luminosity relations are presented for C- and O-rich Miras that appear to be extensions of relations derived in the literature for shorter periods. The fit for the C-stars is particularly well defined (with 182 objects) and reads Mbol = (-2.27 \pm 0.20) log P + (1.45 \pm 0.54)mag with an rms of 0.41 mag.

The first generation of stars were born a few hundred million years after the big bang. These stars synthesized elements heavier than H and He, that are later expelled into the interstellar medium, initiating the rise of metals. Within this enriched medium, the first dust grains formed. This event is cosmological crucial for molecule formation as dust plays a major role by cooling low-metallicity star-forming clouds which can fragment to create lower mass stars. Collecting information on these first dust grains is difficult because of the negative alliance of large distances and low dust masses. We combine the observational information from galaxies at redshifts 5 < z < 10 to constrain their dust emission and theoretically understand the first evolutionary phases of the dust cycle. Spectral energy distributions (SEDs) are fitted with CIGALE and the physical parameters and their evolution are modelled. From this SED fitting, we build a dust emission template for this population of galaxies in the epoch of reionization. Our new models explain why some early galaxies are observed and others are not. We follow in time the formation of the first grains by supernovae later destroyed by other supernova blasts and expelled in the circumgalactic and intergalactic media. We have found evidence for the first dust grains formed in the universe. But, above all, this letter underlines the need to collect more data and to develop new facilities to further constrain the dust cycle in galaxies in the epoch of reionization.

We present the first detection of diffuse dust in the intra-cluster medium of the Virgo cluster out to $\sim$0.4 virial radii, and study the radial variation of its properties on a radial scale of the virial radius. Analysing near-UV - $i$ colours for a sample of $\sim12000$ background galaxies with redshifts $0.02 < z < 0.8$, we find significant colour reddening and relate it to variation in $E(B-V)$ values. The $E(B-V)$ mean profile shows a dust component characterised by an average reddening $E(B-V)\sim0.042 \pm 0.004$ mag within 1.5 degrees ($\sim0.3\, r_{vir}$) from the cluster centre. Assuming a Large Magellanic Cloud extinction law, we derive an average visual extinction $A_{V} = 0.14\pm 0.01$ for a total dust mass, $M_{d} = 2.5\pm0.2\times10^{9}M_{\odot}$, hence a dust-to-gas mass ratio $M_{d}/M_{g} = 3.0\pm 0.3 \times 10^{-4}$. Based on the upper limits on the flux density $\mathrm{I_{250\mu m} = 0.1\, MJy sr^{-1}} $ derived from $Herschel$ data, we estimate an upper limit for the dust temperature of $T_{d} \sim 10\, K$. However, similar densities can be obtained with dust at higher temperatures with lower emissivities. The Virgo cluster has diffuse dust in its intra-cluster medium characterised by different physical properties as those characterising the Milky Way dust. The diffuse dust in Virgo is transported into the cluster space through similar phenomena (stripping) as those building up the optical intra-cluster light, and it constitutes an additional cooling agent of the cluster gas.

The properties of carbon stars in the Magellanic Clouds (MCs) and their total dust production rates are predicted by fitting their spectral energy distributions (SED) over pre-computed grids of spectra reprocessed by dust. The grids are calculated as a function of the stellar parameters by consistently following the growth for several dust species in their circumstellar envelopes, coupled with a stationary wind. Dust radiative transfer is computed taking as input the results of the dust growth calculations. The optical constants for amorphous carbon are selected in order to reproduce different observations in the infrared and optical bands of \textit{Gaia} Data Release 2. We find a tail of extreme mass-losing carbon stars in the Large Magellanic Cloud (LMC) with low gas-to-dust ratios that is not present in the Small Magellanic Cloud (SMC). Typical gas-to-dust ratios are around $700$ for the extreme stars, but they can be down to $\sim160$--$200$ and $\sim100$ for a few sources in the SMC and in the LMC, respectively. The total dust production rate for the carbon star population is $\sim 1.77\pm 0.45\times10^{-5}$~M$_\odot$~yr$^{-1}$, for the LMC, and $\sim 2.52\pm 0.96 \times 10^{-6}$~M$_\odot$~yr$^{-1}$, for the SMC. The extreme carbon stars observed with the Atacama Large Millimeter Array and their wind speed are studied in detail. For the most dust-obscured star in this sample the estimated mass-loss rate is $\sim 6.3 \times 10^{-5}$~M$_\odot$~yr$^{-1}$. The grids of spectra are available at: https://ambrananni085.wixsite.com/ambrananni/online-data-1 and included in the SED-fitting python package for fitting evolved stars https://github.com/s-goldman/Dusty-Evolved-Star-Kit .

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 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.

In this work the optical properties of amorphous carbon (amC) dust condensed around carbon(C)-stars are constrained by comparing the observations for the Large Magellanic Cloud C-stars from the Two Micron All Sky Survey (2MASS) and from the Gaia data release 2 (DR2) with the synthetic photometry obtained by computing dust growth and radiative transfer in their circumstellar envelopes. The set of optical constants of amC dust considered have been pre-selected according to their ability to reproduce the infrared colour-colour diagrams in the Small Magellanic Cloud. Only two combinations of the optical data set and grain size are able to reproduce the infrared photometry and the Gaia observations simultaneously. The analysis presented provides information about the properties of amC dust grains that might be characterized by a diamond-like structure, rather than a graphite-like one, at least around the most dust-enshrouded C-stars, or be composed of "small" grains of size less than 0.04 $\mu$m. The selected data sets will be adopted to compute grids of spectra as a function of the stellar parameters that will be employed to estimate the dust return and mass-loss rates of C-stars by fitting their spectral energy distribution, and to study the resolved stellar populations of nearby objects.

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.

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 present the results of targeted observations and a survey of 1612-, 1665-, and 1667-MHz circumstellar OH maser emission from asymptotic giant branch (AGB) stars and red supergiants (RSGs) in the Small Magellanic Cloud (SMC), using the Parkes and Australia Telescope Compact Array radio telescopes. No clear OH maser emission has been detected in any of our observations targeting luminous, long-period, large-amplitude variable stars, which have been confirmed spectroscopically and photometrically to be mid- to late-M spectral type. These observations have probed 3 - 4 times deeper than any OH maser survey in the SMC. Using a bootstrapping method with LMC and Galactic OH/IR star samples and our SMC observation upper limits, we have calculated the likelihood of not detecting maser emission in any of the two sources considered to be the top maser candidates to be less than 0.05%, assuming a similar pumping mechanism as the LMC and Galactic OH/IR sources. We have performed a population comparison of the Magellanic Clouds and used Spitzer IRAC and MIPS photometry to confirm that we have observed all high luminosity SMC sources that are expected to exhibit maser emission. We suspect that, compared to the OH/IR stars in the Galaxy and LMC, the reduction in metallicity may curtail the dusty wind phase at the end of the evolution of the most massive cool stars. We also suspect that the conditions in the circumstellar envelope change beyond a simple scaling of abundances and wind speed with metallicity.

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.

We present the results of our survey of 1612 MHz circumstellar OH maser emission from asymptotic giant branch (AGB) stars and red supergiants (RSGs) in the Large Magellanic Cloud. We have discovered four new circumstellar maser sources in the LMC, and increased the number of reliable wind speeds from IR stars in the LMC from 5 to 13. Using our new wind speeds, as well as those from Galactic sources, we have derived an updated relation for dust driven winds: $v_{exp} \propto Z L^{0.4}$. We compare the sub-solar metallicity LMC OH/IR stars with carefully selected samples of more metal-rich OH/IR stars, also at known distances, in the Galactic Centre and Galactic Bulge. For 8 of the Bulge stars we derive pulsation periods for the first time, using near-IR photometry from the VVV survey. We have modeled our LMC OH/IR stars and developed an empirical method of deriving gas-to-dust ratios and mass loss rates by scaling the models to the results from maser profiles. We have done this also for samples in the Galactic Centre and Bulge and derived a new mass loss prescription that includes luminosity, pulsation period, and gas-to-dust ratio $\dot{M} = 1.06^{+3.5}_{-0.8} \rm{ \cdot }10^{-5}\,(L/10^4\,\rm{L}_\odot)^{0.9\pm0.1}(P/500\,\rm{d})^{0.75\pm0.3} (r_{gd}/200)^{-0.03\pm0.07}\,\rm{M_{\odot}}\, yr^{-1}$. The tightest correlation is found between mass loss rate and luminosity. We find that the gas-to-dust ratio has little effect on the mass loss of oxygen-rich AGB stars and RSGs within the Galaxy and the LMC. This suggests that mass loss of oxygen-rich AGB stars and RSGs is (nearly) independent of metallicity between a half and twice solar.

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.

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 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.

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.

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)

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.

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.