Modelling polarized light from dust shells surrounding asymptotic giant branch stars

Winds of asymptotic giant branch (AGB) stars are commonly assumed to be driven by radiative acceleration of dust grains. For M-type AGB stars, the nature of the wind-driving dust species has been a matter of intense debate. A proposed source of the radiation pressure triggering the outflows is photon scattering on Fe-free silicate grains. This wind-driving mechanism requires grain radii of about 0.1 - 1 micron in order to make the dust particles efficient at scattering radiation around the stellar flux maximum. Grain size is therefore an important parameter for understanding the physics behind the winds of M-type AGB stars. We seek to investigate the diagnostic potential of scattered polarized light for determining dust grain sizes. We have developed a new tool for computing synthetic images of scattered light in dust and gas shells around AGB stars, which can be applied to detailed models of dynamical atmospheres and dust-driven winds. We present maps of polarized light using dynamical models computed with the DARWIN code. The synthetic images clearly show that the intensity of the polarized light, the position of the inner edge of the dust shell, and the size of the dust grains near the inner edge are all changing with the luminosity phase. Non-spherical structures in the dust shells can also have an impact on the polarized light. We simulate this effect by combining different pulsation phases into a single 3D structure before computing synthetic images. An asymmetry of the circumstellar envelope can create a net polarization, which can be used as diagnostics for the grain size. The ratio between the size of the scattering particles and the observed wavelength determines at what wavelengths net polarization switches direction. If observed, this can be used to constrain average particle sizes.Comment: 9 page

Dynamic atmospheres and winds of cool luminous giants, I. Al2_2O3_3 and silicate dust in the close vicinity of M-type AGB stars

High spatial resolution techniques have given valuable insights into the mass loss mechanism of AGB stars, which presumably involves a combination of atmospheric levitation by pulsation-induced shock waves and radiation pressure on dust. Observations indicate that Al$_2$O$_3$ condenses at distances of about 2 stellar radii or less, prior to the formation of silicates. Al$_2$O$_3$ grains are therefore prime candidates for producing the scattered light observed in the close vicinity of several M-type AGB stars, and they may be seed particles for the condensation of silicates at lower temperatures. We have constructed a new generation of Dynamic Atmosphere & Radiation-driven Wind models based on Implicit Numerics (DARWIN), including a time-dependent treatment of grain growth & evaporation for both Al$_2$O$_3$ and Fe-free silicates (Mg$_2$SiO$_4$). The equations describing these dust species are solved in the framework of a frequency-dependent radiation-hydrodynamical model for the atmosphere & wind structure, taking pulsation-induced shock waves and periodic luminosity variations into account. Condensation of Al$_2$O$_3$ at the close distances and in the high concentrations implied by observations requires high transparency of the grains in the visual and near-IR region to avoid destruction by radiative heating. For solar abundances, radiation pressure due to Al$_2$O$_3$ is too low to drive a wind. Nevertheless, this dust species may have indirect effects on mass loss. The formation of composite grains with an Al$_2$O$_3$ core and a silicate mantle can give grain growth a head start, increasing both mass loss rates and wind velocities. Furthermore, our experimental core-mantle grain models lead to variations of visual and near-IR colors during a pulsation cycle which are in excellent agreement with observations.Comment: Accepted for publication in Astronomy & Astrophysics (18 pages, 9 figures

Constraining dust properties in circumstellar envelopes of C-stars in the Small Magellanic Cloud: optical constants and grain size of carbon dust

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 thermally pulsing asymptotic giant branch, for which we compute spectra and colours. Then, we compare our modelled colours 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 reproduce several colours 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 colours greater than 2 mag in some cases. We conclude that the complete set of observed NIR and MIR colours are best reproduced by small grains, with sizes between ~0.035 and ~0.12 μm, rather than by large grains between ~0.2 and 0.7 μm. The inability of large grains to reproduce NIR and MIR colours 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. © 2016 The Authors

Constraining the thermally pulsing asymptotic giant branch phase with resolved stellar populations in the Small Magellanic Cloud

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⊙. 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.This work is supported by the European Research Council (ERC) Consolidator Grant funding scheme (project STARKEY, G.A. n. 615604). We thank the entire VMC team for producing the spaceresolved SFH used in this work and J. Cummings and J. Kalirai for kindly providing us with their IFMR data. MRLC acknowledges funding from the ERC under the European Union’s Horizon 2020 research and innovation programme (G.A. n. 682115). Many thanks go to C. Maraston, S. Charlot, and G. Bruzual for providing us with their stellar population synthesis models

Constraining the thermally pulsing asymptotic giant branch phase with resolved stellar populations in the Small Magellanic Cloud

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

Constraining the thermally pulsing asymptotic giant branch phase with resolved stellar populations in the Large Magellanic Cloud

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 Ks-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 Mi 48 1.7 and 3 M 99 at Zi = 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

A new generation of Parsec-Colibri stellar isochrones including the TP-AGB phase

We introduce a new generation of PARSEC-COLIBRI stellar isochrones that includes a detailed treatment of the thermally pulsing asymptotic giant branch (TP-AGB) phase, 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 firstovertone modes; and 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

Dynamical atmospheres and winds of M-type AGB stars

Mass loss, in the form of slow stellar winds, is a decisive factor for the evolution of cool luminous giants, eventually turning them into white dwarfs. These dense outflows are also a key factor in the enrichment of the interstellar medium with newly produced elements from the interior of these stars. There are strong indications that these winds are accelerated by radiation pressure on dust grains, but the actual grain species responsible for driving the outflows in M-type Asymptotic Giant Branch stars are still a matter of debate. Observations of dust features in the circumstellar environment of these stars suggest that magnesium-iron silicates are possible wind-drivers. However, the optical properties of these silicate grains are strongly influenced by the Fe-content. Fe-bearing condensates heat up strongly when interacting with the radiation field and therefore cannot form close enough to the star to trigger outflows. Fe-free condensates, on the other hand, have a low absorption cross-section at near-IR wavelengths where AGB stars emit most of their flux.  To solve this conundrum, it has been suggested that winds of M-type AGB stars may be driven by photon scattering on Fe-free silicate grains with sizes comparable to the wavelength of the flux maximum, rather than by true absorption. In this thesis we investigate dynamical models of M-type AGB stars, using Fe-free silicates as the wind-driving dust species. According to our findings these models produce both dynamic and photometric properties consistent with observations. Especially noteworthy are the large photometric variations in the visual band during a pulsation cycle, seen both in the observed and synthetic fluxes. A closer examination of the models reveals that these variations are caused by changes in the molecular layers, and not by changes in the dust. This is a strong indication that stellar winds of M-type AGB stars are driven by dust materials that are very transparent in the visual and near-infrared wavelength regions, otherwise these molecular effects would not be visible