No Arabic abstract
In order to determine the composition of the dust in the circumstellar envelopes of oxygen-rich asymptotic giant branch (AGB) stars we have computed a grid of modust radiative-transfer models for a range of dust compositions, mass-loss rates, dust shell inner radii and stellar parameters. We compare the resulting colours with the observed oxygen-rich AGB stars from the SAGE-Spec Large Magellanic Cloud (LMC) sample, finding good overall agreement for stars with a mid-infrared excess. We use these models to fit a sample of 37 O-rich AGB stars in the LMC with optically thin circumstellar envelopes, for which 5$-$35-$mu$m Spitzer infrared spectrograph (IRS) spectra and broadband photometry from the optical to the mid-infrared are available. From the modelling, we find mass-loss rates in the range $sim 8times10^{-8}$ to $5times10^{-6}$ M$_{odot} mathrm{yr}^{-1}$, and we show that a grain mixture consisting primarily of amorphous silicates, with contributions from amorphous alumina and metallic iron provides a good fit to the observed spectra. Furthermore, we show from dust models that the AKARI [11]$-$[15] versus [3.2]$-$[7] colour-colour diagram, is able to determine the fractional abundance of alumina in O-rich AGB stars.
The circumstellar ammonia (NH$_3$) chemistry in evolved stars is poorly understood. Previous observations and modelling showed that NH$_3$ abundance in oxygen-rich stars is several orders of magnitude above that predicted by equilibrium chemistry. In this article, we characterise the spatial distribution and excitation of NH$_3$ in the O-rich circumstellar envelopes (CSEs) of four diverse targets: IK Tau, VY CMa, OH 231.8+4.2, and IRC +10420 with multi-wavelength observations. We observed the 1.3-cm inversion line emission with the Very Large Array (VLA) and submillimetre rotational line emission with the Heterodyne Instrument for the Far-Infrared (HIFI) aboard Herschel from all four targets. For IK Tau and VY CMa, we observed the rovibrational absorption lines in the $ u_2$ band near 10.5 $mu$m with the Texas Echelon Cross Echelle Spectrograph (TEXES) at the NASA Infrared Telescope Facility (IRTF). We also attempted to search for the rotational transition within the $v_2=1$ state near 2 mm with the IRAM 30m Telescope towards IK Tau. Non-LTE radiative transfer modelling, including radiative pumping to the vibrational state, was carried out to derive the radial distribution of NH$_3$ in these CSEs. Our modelling shows that the NH$_3$ abundance relative to molecular hydrogen is generally of the order of $10^{-7}$, which is a few times lower than previous estimates that were made without considering radiative pumping and is at least 10 times higher than that in the C-rich CSE of IRC +10216. Incidentally, we also derived a new period of IK Tau from its $V$-band light curve. NH$_3$ is again detected in very high abundance in O-rich CSEs. Its emission mainly arises from localised spatial-kinematic structures that are probably denser than the ambient gas. Circumstellar shocks in the accelerated wind may contribute to the production of NH$_3$. (Abridged abstract)
We combine variability information from the MAssive Compact Halo Objects (MACHO) survey of the Large Magellanic Cloud (LMC) with infrared photometry from the Spitzer Space Telescope Surveying the Agents of a Galaxys Evolution (SAGE) survey to create a dataset of ~30 000 variable red sources. We photometrically classify these sources as being on the first ascent of the Red Giant Branch (RGB), or as being in one of three stages along the Asymptotic Giant Branch (AGB): oxygen-rich, carbon-rich, or highly reddened with indeterminate chemistry (extreme AGB candidates). We present linear period-luminosity relationships for these sources using 8 separate infrared bands (J, H, K, 3.6, 4.5, 5.8, 8.0, and 24 micron) as proxies for the luminosity. We find that the wavelength dependence of the slope of the period-luminosity relationship is different for different photometrically determined classes of AGB stars. Stars photometrically classified as O-rich show the least variation of slope with wavelength, while dust enshrouded extreme AGB stars show a pronounced trend toward steeper slopes with increasing wavelength. We find that O-rich AGB stars pulsating in the fundamental mode obey a period-magnitude relation with a slope of -3.41 +/- 0.04 when magnitude is measured in the 3.6 micron band, in contrast to C-rich AGB stars, which obey a relation of slope -3.77 +/- 0.05.
We present a 2Dust model for the dust shell around a LMC long-period variable (LPV) previously studied as part of the OGLE survey. OGLE LMC LPV 28579 (SAGE J051306.40-690946.3) is a carbon-rich asymptotic giant branch (AGB) star for which we have photometry and spectra from the Spitzer SAGE and SAGE-Spec programs along with UBVIJHK_s photometry. By modeling this source, we obtain a baseline set of dust properties to be used in the construction of a grid of models for carbon stars. We reproduce its spectral energy distribution using a mixture of AmC and SiC (15% by mass). The grain sizes are distributed according to the KMH model. The best-fit model has an optical depth of 0.28 for the shell at the peak of the SiC feature, with R_in~1430 R_sun or 4.4 R_star. The temperature at this inner radius is 1310 K. Assuming an expansion velocity of 10 km s^-1, we obtain a dust mass-loss rate of 2.5x10^-9 M_sun yr-1. We calculate a 15% variation in this rate by testing the fit sensitivity against variation in input parameters. We also present a simple model for the molecular gas in the extended atmosphere that could give rise to the 13.7 mu m feature seen in the spectrum. We find that a combination of CO and C_2H_2 gas at an excitation temperature of about 1000 K and column densities of 3x10^21 cm^-2 and 10^19 cm^-2 respectively are able to reproduce the observations. Given that the excitation temperature is close to T_dust(R_in), most of the molecular contribution probably arises from the inner shell region. The luminosity corresponding to the first epoch of SAGE observations is 6580 L_sun. For an effective temperature of about 3000 K, this implies a stellar mass of 1.5-2 M_sun and an age of 1-2.5 Gyr. For a gas:dust ratio of 200, we obtain a gas mass-loss rate of 5.0x10^-7 M_sun yr^-1, consistent with the gas mass-loss rates estimated from the period, color and 8 mu m flux of the source.
For decades ever since the early detection in the 1990s of the emission spectral features of crystalline silicates in oxygen-rich evolved stars, there is a long-standing debate on whether the crystallinity of the silicate dust correlates with the stellar mass loss rate. To investigate the relation between the silicate crystallinities and the mass loss rates of evolved stars, we carry out a detailed analysis of 28 nearby oxygen-rich stars. We derive the mass loss rates of these sources by modeling their spectral energy distributions from the optical to the far infrared. Unlike previous studies in which the silicate crystallinity was often measured in terms of the crystalline-to-amorphous silicate mass ratio, we characterize the silicate crystallinities of these sources with the flux ratios of the emission features of crystalline silicates to that of amorphous silicates. This does not require the knowledge of the silicate dust temperatures which are the major source of uncertainties in estimating the crystalline-to-amorphous silicate mass ratio. With a Pearson correlation coefficient of ~0.24, we find that the silicate crystallinities and the mass loss rates of these sources are not correlated. This supports the earlier findings that the dust shells of low mass-loss rate stars can contain a significant fraction of crystalline silicates without showing the characteristic features in their emission spectra.
A sample of 28 oxygen-rich evolved stars is selected based on the presence of crystalline silicate emission features in their ISO/SWS spectra. The crystallinity, measured as the flux fraction of crystalline silicate features, is found not related to mass loss rate that is derived from fitting the spectral energy distribution.