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.

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.

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.

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.

We propose a combination of a modified Wide-Fast-Deep survey, a mini-survey of the South Celestial Pole, and a Deep Drilling-style survey to produce a 3-D map of the Magellanic System and to provide a detailed census of the transient and variable populations in the Clouds. We support modifying the Wide-Fast-Deep survey to cover the declination range $-72.25\deg<{\rm Dec}<12.4\deg$ and the Galactic latitude range $|b|>15\deg$, as proposed in a separate white paper. We additionally propose a mini-survey covering the 950$\deg^2$ with ${\rm Dec} < -72.25$ in $ugriz$ to the standard LSST single-exposure depth and with 40 visits per filter per field. Finally, we propose a mini-survey covering $\sim100 \deg^2$ of the main bodies of the Clouds with twelve total pointings, 2000 total visits per field, and shorter exposure time.

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.

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

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.

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