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 Opacity Project At Livermore (OPAL) data at high temperatures to the molecular opacities computed with our æsopus code. 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., 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 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 < Z i < 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 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 selfconsistently: (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) the conditions under which sphericity effects may significantly affect the atmospheres of giant stars; (iii) the core mass-luminosity relation and its possible break-down due to the occurrence of hot bottom burning (HBB) in the most massive AGB stars, by taking properly into account the nuclear energy generation in the H-burning shell and in the deepest layers of the convective envelope; (iv) the HBB nucleosynthesis via the solution of a complete nuclear network (including the pp chains, and the CNO, NeNa, MgAl cycles) coupled to a diffusive description of mixing, suitable to follow also the synthesis of 7 Li via the Cameron-Fowler beryllium transport mechanism; (v) the intershell abundances left by each thermal pulse via the solution of a complete nuclear network applied to a simple model of the pulsedriven convective zone; (vi) the onset and quenching of the third dredge-up, with a temperature criterion that is applied, at each thermal pulse, to the result of envelope integrations at the stage of the post-flash luminosity peak.At the same time COLIBRI pioneers new techniques in the treatment of the physics of stellar interiors, not yet adopted in full TP-AGB models. 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 atmosphere and the deep envelope. This ensures a complete consistency, step by step, of both EoS and opacity with the evolution of the chemical abundances caused by the third dredge-up and HBB. Another distinguishing aspect of COLIBRI is its high computational speed, that allows to generate complete grids of TP-AGB models in just a few hours. This feature is absolutely necessary for calibrating the many uncertain parameters and processes that characterize the TP-AGB phase.We illustrate the many unique features of COLIBRI by means of detailed evolutionary tracks computed for several choices of model parameters, including initial star masses, chemical abundances, nuclear reaction rates, efficiency of the third dredge-up, overshooting at the base of the pulse-driven convection zone, etc. Future papers in this series will deal with the calibration of all these and other parameters using observational data of AGB stars in the Galaxy and in nearby systems, a step that is of paramount importance for producing reliable stellar population synthesis models...

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 H 2 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 T cond ∼1400 K) than those found when chemisputtering is included (T cond ∼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. If confirmed by more observational data, it would indicate that either the amount of the carbon excess, determined by the complex interplay between mass loss, third dredge-up and hot bottom burning, or the homogeneous growth scheme should be revised. Finally, we analyze the differences in the total dust production of M-stars that arise from the use of the two approaches (i.e. with or without chemisputtering). We find that, in spite of the differences in the expected dust stratification, for a given set of TP-AGB models, the ejecta are only weakly sensitive to the specific assumption. This work also shows that the properties of TP-AGB circumstellar envelopes are important diagnostic tools that may be profitably added to the traditional calibrators for setting further constraints on this complex phase of stellar evolution. our galaxy and in the nearby ones (Knapp 1985;Matsuura et al. 2009Matsuura et al. , 2012.Comparing the dust mass loss rates derived from far-infrared observations with the gas mass loss rates obtained from CO observations for Galactic AGB stars, Knapp (1985) have found typical dust-to-gas ratios of ∼6 × 10 −3 for oxygen-rich (M) star...