A stubbornly large mass of cold dust in the ejecta of Supernova 1987A

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

Spitzer Space Telescope spectra of post-AGB stars in the Large Magellanic Cloud ---polycyclic aromatic hydrocarbons at low metallicities

This paper reports variations of polycyclic aromatic hydrocarbons (PAHs) features that were found in Spitzer Space Telescope spectra of carbon-rich post-asymptotic giant branch (post-AGB) stars in the Large Magellanic Cloud (LMC). The paper consists of two parts. The first part describes our Spitzer spectral observing programme of 24 stars including post-AGB candidates. The latter half of this paper presents the analysis of PAH features in 20 carbon-rich post-AGB stars in the LMC, assembled from the Spitzer archive as well as from our own programme. We found that five post-AGB stars showed a broad feature with a peak at 7.7 micron, that had not been classified before. Further, the 10--13 micron PAH spectra were classified into four classes, one of which has three broad peaks at 11.3, 12.3 and 13.3 micron rather than two distinct sharp peaks at 11.3 and 12.7 micron, as commonly found in HII regions. Our studies suggest that PAHs are gradually processed while the central stars evolve from post-AGB phase to PNe, changing their composition before PAHs are incorporated into the interstellar medium. Although some metallicity dependence of PAH spectra exists, the evolutionary state of an object is more significant than its metallicity in determining the spectral characteristics of PAHs for LMC and Galactic post-AGB stars.

Herschel SPIRE and PACS observations of the red supergiant VY CMa: analysis of the molecular line spectra

We present an analysis of the far-infrared and submillimetre molecular emission line spectrum of the luminous M-supergiant VY CMa, observed with the SPIRE and PACS spectrometers aboard the Herschel Space Observatory. Over 260 emission lines were detected in the 190-650-micron SPIRE FTS spectra, with one-third of the observed lines being attributable to H2O. Other detected species include CO, 13CO, H2^18O, SiO, HCN, SO, SO2, CS, H2S, and NH3. Our model fits to the observed 12CO and 13CO line intensities yield a 12C/13C ratio of 5.6+-1.8, consistent with measurements of this ratio for other M supergiants, but significantly lower than previously estimated for VY CMa from observations of lower-J lines. The spectral line energy distribution for twenty SiO rotational lines shows two temperature components: a hot component at 1000 K, which we attribute to the stellar atmosphere and inner wind, plus a cooler ~200 K component, which we attribute to an origin in the outer circumstellar envelope. We fit the line fluxes of 12CO, 13CO, H2O and SiO, using the SMMOL non-LTE line transfer code, with a mass-loss rate of 1.85x10^-4 Msun yr^-1 between 9 R* and 350 R*. To fit the observed line fluxes of 12CO, 13CO, H2O and SiO with SMMOL non-LTE line radiative transfer code, along with a mass-loss rate of 1.85x10^-4 Msun yr^-1. To fit the high rotational lines of CO and H2O, the model required a rather flat temperature distribution inside the dust condensation radius, attributed to the high H2O opacity. Beyond the dust condensation radius the gas temperature is fitted best by an r^-0.5 radial dependence, consistent with the coolant lines becoming optically thin. Our H2O emission line fits are consistent with an ortho:para ratio of 3 in the outflow.

The global gas and dust budget of the Small Magellanic Cloud

In order to understand the evolution of the interstellar medium (ISM) of a galaxy, we have analysed the gas and dust budget of the Small Magellanic Cloud (SMC). Using the Spitzer Space Telescope, we measured the integrated gas mass-loss rate across asymptotic giant branch (AGB) stars and red supergiants (RSGs) in the SMC, and obtained a rate of 1.4x10^-3 Msun yr-1. This is much smaller than the estimated gas ejection rate from type II supernovae (SNe) (2-4x10^-2 Msun yr-1). The SMC underwent a an increase in starformation rate in the last 12 Myrs, and consequently the galaxy has a relatively high SN rate at present. Thus, SNe are more important gas sources than AGB stars in the SMC. The total gas input from stellar sources into the ISM is 2-4x10^-2 Msun yr-1. This is slightly smaller than the ISM gas consumed by starformation (~8x10^-2 Msun yr-1). Starformation in the SMC relies on a gas reservoir in the ISM, but eventually the starformation rate will decline in this galaxy, unless gas infalls into the ISM from an external source. The dust injection rate from AGB and RSG candidates is 1x10^-5 Msun yr-1. Dust injection from SNe is in the range of 0.2--11x10^-4 Msun yr-1, although the SN contribution is rather uncertain. Stellar sources could be important for ISM dust (3x10^5 Msun yr-1) in the SMC, if the dust lifetime is about 1.4 Gyrs. We found that the presence of poly-aromatic hydrocarbons (PAHs) in the ISM cannot be explained entirely by carbon-rich AGB stars. Carbon-rich AGB stars could inject only 7x10^-9 Msun yr-1 of PAHs at most, which could contribute up to 100 Msun of PAHs in the lifetime of a PAH. The estimated PAH mass of 1800 Msun in the SMC can not be explained. Additional PAH sources, or ISM reprocessing should be needed.

Radiation Transfer in the Cavity and Shell of Planetary Nebulae

We develop an approximate analytical solution for the transfer of line-averaged radiation in the hydrogen recombination lines for the ionized cavity and molecular shell of a spherically symmetric planetary nebula. The scattering problem is treated as a perturbation, using a mean intensity derived from a scattering-free solution. The analytical function was fitted to Halpha and Hbeta data from the planetary nebula NGC6537. The position of the maximum in the intensity profile produced consistent values for the radius of the cavity as a fraction of the radius of the dusty nebula: 0.21 for Halpha and 0.20 for Hbeta. Recovered optical depths were broadly consistent with observed optical extinction in the nebula, but the range of fit parameters in this case is evidence for a clumpy distribution of dust.

Observational studies of mass loss from AGB stars

It is important to properly describe the mass-loss rate of AGB stars, in order to understand their evolution from the AGB to PN phase. The primary goal of this study is to investigate the influence of metallicity on the mass-loss rate, under well determined luminosities. The luminosity of the star is a crucial parameter for the radiative driven stellar wind. Many efforts have been invested to constrain the AGB mass-loss rate, but most of the previous studies use Galactic objects, which have poorly known distances, thus their luminosities. To overcome this problem, we have studied mass loss from AGB stars in the Galaxies of the Local Group. The distance to the stars have been independently measured, thus AGB stars in these galaxies are ideal for understanding the mass-loss rate. Moreover, these galaxies have a lower metallicity than the Milky Way, providing an ideal target to study the influence of metallicity on the mass-loss rate. We report our analysis of mass loss, using the Spitzer Space Telescope and the Herschel Space Observatory. We will discuss the influence of AGB mass-loss on stellar evolution, and explore AGB and PN contribution to the lifecycle of matter in galaxies.

Herschel Detects a Massive Dust Reservoir in Supernova 1987A

We report far-infrared and submillimeter observations of Supernova 1987A, the star that exploded on February 23, 1987 in the Large Magellanic Cloud, a galaxy located 160,000 lightyears away. The observations reveal the presence of a population of cold dust grains radiating with a temperature of ~17-23 K at a rate of about 220 solar luminosity. The intensity and spectral energy distribution of the emission suggests a dust mass of ~0.4-0.7 solar mass. The radiation must originate from the SN ejecta and requires the efficient precipitation of all refractory material into dust. Our observations imply that supernovae can produce the large dust masses detected in young galaxies at very high redshifts.

The global gas and dust budget of the Large Magellanic Cloud: AGB stars and supernovae, and the impact on the ISM evolution

We report on an analysis of the gas and dust budget in the the interstellar medium (ISM) of the Large Magellanic Cloud (LMC). Recent observations from the Spitzer Space Telescope enable us to study the mid-infrared dust excess of asymptotic giant branch (AGB) stars in the LMC. This is the first time we can quantitatively assess the gas and dust input from AGB stars over a complete galaxy, fully based on observations. The integrated mass-loss rate over all intermediate and high mass-loss rate carbon-rich AGB candidates in the LMC is 8.5x10^-3 solar mass per year, up to 2.1x10^-2 solar mass per year. This number could be increased up to 2.7x10^-2 solar mass per year, if oxygen-rich stars are included. This is overall consistent with theoretical expectations, considering the star formation rate when these low- and intermediate-mass stars where formed, and the initial mass functions. AGB stars are one of the most important gas sources in the LMC, with supernovae (SNe), which produces about 2-4x10^-2 solar mass per year. At the moment, the star formation rate exceeds the gas feedback from AGB stars and SNe in the LMC, and the current star formation depends on gas already present in the ISM. This suggests that as the gas in the ISM is exhausted, the star formation rate will eventually decline in the LMC, unless gas is supplied externally. Our estimates suggest `a missing dust-mass problem in the LMC, which is similarly found in high-z galaxies: the accumulated dust mass from AGB stars and possibly SNe over the dust life time (400--800 Myrs) is significant less than the dust mass in the ISM. Another dust source is required, possibly related to star-forming regions.

Spitzer Space Telescope spectral observations of AGB stars in the Fornax dwarf spheroidal galaxy

We have observed five carbon-rich AGB stars in the Fornax dwarf spheroidal (dSph) galaxy, using the Infrared Spectrometer on board the Spitzer Space Telescope. The stars were selected from a near-infrared survey of Fornax and include the three reddest stars, with presumably the highest mass-loss rates, in that galaxy. Such carbon stars probably belong to the intermediate-age population (2-8 Gyr old and metallicity of [Fe/H] -1) of Fornax. The primary aim of this paper is to investigate mass-loss rate, as a function of luminosity and metallicity, by comparing AGB stars in several galaxies with different metallicities. The spectra of three stars are fitted with a radiative transfer model. We find that mass-loss rates of these three stars are 4-7x10^-6 Msun yr-1. The other two stars have mass-loss rates below 1.3x10^-6 Msun yr-1. We find no evidence that these rates depend on metallicity, although we do suggest that the gas-to-dust ratio could be higher than at solar metallicity, in the range 240 to 800. The C2H2 bands are stronger at lower metallicity because of the higher C/O ratio. In contrast, the SiC fraction is reduced at low metallicity, due to low silicon abundance. The total mass-loss rate from all known carbon-rich AGB stars into the interstellar medium of this galaxy is of the order of 2x10^-5 Msun yr-1. This is much lower than that of the dwarf irregular galaxy WLM, which has a similar visual luminosity and metallicity. The difference is attributed to the younger stellar population of WLM. The suppressed gas-return rate to the ISM accentuates the difference between the relatively gas-rich dwarf irregular and the gas-poor dwarf spheroidal galaxies. Our study will be useful to constrain gas and dust recycling processes in low metallicity galaxies.

VLT / Infrared Integral Field Spectrometer Observations of Molecular Hydrogen Lines in the Knots in the Planetary Nebula NGC 7293 (the Helix Nebula)

Knots are commonly found in nearby planetary nebulae (PNe) and star forming regions. Within PNe, knots are often found to be associated with the brightest parts of the nebulae and understanding the physics involved in knots may reveal the processes dominating in PNe. As one of the closest PNe, the Helix Nebula (NGC 7293) is an ideal target to study such small-scale (~300 AU) structures. We have obtained infrared integral spectroscopy of a comet-shaped knot in the Helix Nebula using SINFONI on the Very Large Telescope at high spatial resolution (50-125 mas). With spatially resolved 2 micron spectra, we find that the H2 rotational temperature within the cometary knots is uniform. The rotational-vibrational temperature of the cometary knot (situated in the innermost region of the nebula, 2.5 arcmin away from the central star), is 1800 K, higher than the temperature seen in the outer regions (5-6 arcmin from the central star) of the nebula (900 K), showing that the excitation temperature varies across the nebula. The obtained intensities are reasonably well fitted with 27 km s-1 C-type shock model. This ambient gas velocity is slightly higher than the observed [HeII] wind velocity of 13 km s-1. The gas excitation can also be reproduced with a PDR (photo dominant region) model, but this requires an order of magnitude higher UV radiation. Both models have limitations, highlighting the need for models that treats both hydrodynamical physics and the PDR.