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

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