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 .
Reliable models of the thermally pulsing asymptotic giant branch (TP-AGB) phase are of critical importance across astrophysics, including our interpretation of the spectral energy distribution of ...galaxies, cosmic dust production, and enrichment of the interstellar medium. With the aim of improving sets of stellar isochrones that include a detailed description of the TP-AGB phase, we extend our recent calibration of the AGB population in the Small Magellanic Cloud (SMC) to the more metal rich Large Magellanic Cloud (LMC). We model the LMC stellar populations with the TRILEGAL code, using the spatially-resolved star formation history derived from the VISTA survey. We characterize the efficiency of the third dredge-up by matching the star counts and the \(K_{\rm s}\)-band luminosity functions of the AGB stars identified in the LMC. In line with previous findings, we confirm that, compared to the SMC, the third dredge-up in AGB stars of the LMC is somewhat less efficient, as a consequence of the higher metallicity. The predicted range of initial mass of C-rich stars is between \(M_{\rm i} \approx 1.7 - 3~\mathrm{M}_{\odot}\) at \(Z_{\rm i} = 0.008\). We show how the inclusion of new opacity data in the carbon star spectra will improve the performance of our models. We discuss the predicted lifetimes, integrated luminosities and mass-loss rate distributions of the calibrated models. The results of our calibration are included in updated stellar isochrones publicly available.
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 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.
In this work we study how to employ the upcoming Legacy Survey of Space and Time (LSST) data to constrain physical properties of normal, star forming galaxies. We use simulated LSST data and existing ...real observations to test the estimations of the physical properties of galaxies, such as star formation rate (SFR), stellar mass (\(M_{star}\)), and dust luminosity (\(L_{dust}\)). We focus on normal star-forming galaxies, as they form the majority of the galaxy population in the universe and therefore are more likely to be observed by the LSST. We perform a simulation of LSST observations and uncertainties of 50,385 real galaxies within redshift range \(0<z<2.5\). In order to achieve this goal, we used the unique multi-wavelength data from the Herschel Extragalactic Legacy Project (HELP) survey. Our analysis focus on two fields: ELAIS-N1 and COSMOS. To obtain galaxy physical properties we fit their Spectral Energy Distributions (SEDs) using the Code Investigating GALaxy Emission (CIGALE). We compare the main galaxy physical properties obtained from the fit of the observed multi-wavelength photometry of galaxies (from UV to FIR) to the ones obtained from the simulated LSST optical measurements only. The stellar masses estimated based on the LSST measurements are in agreement with the full UV-FIR SED estimations, as they depend mainly on the UV and optical emission, well covered by LSST in the considered redshift range. We obtain a clear overestimation of SFR, \(L_{dust}\), \(M_{dust}\) estimated with LSST only, highly correlated with redshift. We investigate the cause of this overestimation and we conclude that it is related to an overestimation of the dust attenuation, both UV and NIR. We find that it is necessary to employ auxiliary rest-frame mid-infrared observations, simulated UV observations, or FUV attenuation (AFUV)- Mstar relation, to correct the overestimation.
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 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.
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 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 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.