Martin A. T. Groenewegen

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.5 × 10 −3 M⊙ yr −1 , up to 2.1 × 10 −2 M⊙ yr −1 . This number could be increased up to 2.7×10 −2 M⊙ yr −1 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–4×10 −2 M⊙ yr −1 . 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–800Myrs) is significant less than the dust mass in the ISM. Another dust source is required, possibly related to star-forming regions.

MESS (Mass-loss of Evolved StarS), a Herschel key program

MESS (Mass-loss of Evolved StarS) is a guaranteed time key program that uses the PACS and SPIRE instruments on board the Herschel space observatory to observe a representative sample of evolved stars, that include asymptotic giant branch (AGB) and post-AGB stars, planetary nebulae and red supergiants, as well as luminous blue variables, Wolf-Rayet stars and supernova remnants. In total, of order 150 objects are observed in imaging and about 50 objects inspectroscopy. This paper describes the target selection and target list, and the observing strategy. Key science projects are described, and illustrated using results obtained during Herschel’s science demonstration phase. Aperture photometry is given for the 70 AGB and post-AGB stars observed up to October 17, 2010, which constitutes the largest single uniform database of far-IR and sub-mm fluxes for late-type stars.

Post-AGB stars with hot circumstellar dust: binarity of the low-amplitude pulsators

Context. The influence of binarity on the late stages of stellar evolut ion. Aims. While the first binary post-AGB stars were serendipitously d iscovered, the distinct characteristics of their Spectral Energy Distribution (SED) allowed us to launch a more systematic search for binaries. We selected post-AGB objects which show a broad dust excess often starting already at H or K, pointing to the presence of a gravitationally bound dusty disc in the system. We started a very extensive multi-wavelength study of those systems and here we report on our radial velocity and photometric monitoring results for six stars of early F type, which are pulsators of small amplitude. Methods. To determine the radial velocity of low signal-to-noise time-series, we constructed dedicated auto-correlation masks based on high signal-to-noise spectra, used in our published chemical studies. The radial velocity variations were subjecte d to detailed analysis to differentiate between pulsational variability and variabilit y due to orbital motion. When available, the photometric monitoring data were used to complement the time series of radial velocity data and to establish the nature of the pulsation. Finally orbital minimalisation was performed to constrain the orbital elements. Results. All of the six objects are binaries, with orbital periods ran ging from 120 to 1800 days. Five systems have non-circular orbits. The mass functions range from 0.004 to 0.57 M⊙ and the companions are likely unevolved objects of (very) low initial mass. We argue that these binaries must have been subject to severe binary interaction when the primary was a cool supergiant. Although the origin of the circumstellar disc is not well understood, the disc is generally believed to be formed during this strong interaction phase. The eccentric orbits of these highly evolved objects remain poorly understood. In one object the line-of-sight is grazi ng the edge of the puffed-up inner rim of the disc. Conclusions. These results corroborate our earlier statement that evolved objects in binary stars create a Keplerian dusty circumbinary disc. With the measured orbits and mass functions we conclude that the circumbinary discs seem to have a major impact on the evolution of a significant fraction of binary systems.

On the metallicity gradient of the Galactic disk

Aims. The iron abundance gradient in the Galactic stellar disk provides fundamental constraints on the chemical evolution of this important Galaxy component, however the spread around the mean slope is, at fixed Galactocentric distance, more than the estimated uncertainties. Methods. To provide quantitative constraints on these trends, we adopted iron abundances for 265 classical Cepheids (more than 50% of the currently known sample) based either on high-resolution spectra or on photometric metallicity indices. Homogeneous distances were estimated using near-infrared period-luminosity relations. The sample covers the four disk quadrants, and their Galactocentric distances range from similar to 5 to similar to 17 kpc. We provided a new theoretical calibration of the metallicity-index-color (MIC) relation based on Walraven and NIR photometric passbands. Results. We estimated the photometric metallicity of 124 Cepheids. Among them 66 Cepheids also have spectroscopic iron abundances and we found that the mean difference is -0.03 +/- 0.15 dex. We also provide new iron abundances, based on high-resolution spectra, for four metal-rich Cepheids located in the inner disk. The remaining iron abundances are based on high-resolution spectra collected by our group (73) or available in the literature (130). A linear regression over the entire sample provides an iron gradient of -0.051 +/- 0.004 dex kpc(-1). The above slope agrees quite well, within the errors, with previous estimates based either on Cepheids or on open clusters covering similar Galactocentric distances. However, Cepheids located in the inner disk systematically appear more metal-rich than the mean metallicity gradient. Once we split the sample into inner (R(G) <8 kpc) and outer disk Cepheids, the slope (-0.130 +/- 0.015 dex kpc(-1)) in the former region is approximate to 3 times steeper than the slope in the latter one (-0.042 +/- 0.004 dex kpc(-1)). In the outer disk the radial distribution of metal-poor (MP, [Fe/H] <-0.02 dex) and metal-rich (MR) Cepheids across the four disk quadrants does not show a clear trend when moving from the innermost to the external disk regions. The relative fractions of MP and MR Cepheids in the 1st and in the 3rd quadrants differ at the 8 sigma (MP) and 15 sigma (MR) levels. Finally, we found that iron abundances in two local overdensities of the 2nd and of the 4th quadrant cover individually a range in iron abundance of approximate to 0.5 dex. Conclusions. Current findings indicate that the recent chemical enrichment across the Galactic disk shows a clumpy distribution.

Unexpectedly large mass loss during the thermal pulse cycle of the red giant star R Sculptoris.

The asymptotic-giant-branch star R Sculptoris is surrounded by a detached shell of dust and gas. The shell originates from a thermal pulse during which the star underwent a brief period of increased mass loss. It has hitherto been impossible to constrain observationally the timescales and mass-loss properties during and after a thermal pulse—parameters that determine the lifetime of the asymptotic giant branch and the amount of elements returned by the star. Here we report observations of CO emission from the circumstellar envelope and shell around R Sculptoris with an angular resolution of 1.3″. What was previously thought to be only a thin, spherical shell with a clumpy structure is revealed to also contain a spiral structure. Spiral structures associated with circumstellar envelopes have been previously seen, leading to the conclusion that the systems must be binaries. Combining the observational data with hydrodynamic simulations, we conclude that R Sculptoris is a binary system that underwent a thermal pulse about 1,800 years ago, lasting approximately 200 years. About 3 × 10−3 solar masses of material were ejected at a velocity of 14.3 km s−1 and at a rate around 30 times higher than the pre-pulse mass-loss rate. This shows that about three times more mass was returned to the interstellar medium during and immediately after the pulse than previously thought.

The VMC survey - IV. The LMC star formation history and disk geometry from four VMC tiles

We derive the star formation history (SFH) for several regions of the Large Magellanic Cloud (LMC), using deep near-infrared data from the VISTA near-infrared Y JKs survey of the Magellanic system (VMC). The regions include three almost-complete 1.4 deg 2 tiles located ∼ 3.5 ◦ away from the LMC centre in distinct directions. They are split into 21.0 ′ × 21.5 ′ (0.12 deg 2 ) subregions, and each of these is analysed independently. To this dataset, we add two 11.3 ′ × 11.3 ′ (0.036 deg 2 ) subregions selected based on their small and uniform extinction inside the 30 Doradus tile. The SFH is derived from the simultaneous reconstruction of two different colour‐magnitude diagrams (CMDs), using the minimization code StarFISH together with a database of “partial models” representing the CMDs of LMC populations of various ages and metallicities, plus a partial model for the CMD of the Milky Way foreground. The distance modulus (m− M)0 and extinction AV is varied within intervals∼ 0.2 and∼ 0.5 mag wide, respectively, within which we identify the best-fitting star formation rate SFR( t) as a function of lookback time t, age‐metallicity relation (AMR), (m− M)0 and AV. Our results demonstrate that VMC data, due to the combination of depth and little sensitivity to differential reddening, allow the derivation of the space-reso lved SFH of the LMC with unprecedented quality compared to previous wide-area surveys. In particular, the data clea rly reveal the presence of peaks in the SFR(t) at ages log(t/yr)≃ 9.3 and 9.7, which appear in most of the subregions. The most recent SFR(t) is found to vary greatly from subregion to subregion, with the general trend of being more intense in the innermost LMC, except for the tile next to the N11 complex. In the bar region, the SFR(t) seems remarkably constant over the time interval from log(t/yr)≃ 8.4 to 9.7. The AMRs, instead, turn out to be remarkably similar across the LMC. Thanks to the accuracy in determining the distance modulus for every subregion ‐ with typical errors of just∼ 0.03 mag ‐ we make a first attempt to derive a spatial model of the LMC disk. The fields studied so far are fit extremel y well by a single disk of inclination i = 26.2± 2.0 ◦ , position angle of the line of nodesθ0 = 129.1± 13.0 ◦ , and distance modulus of (m− M)0 = 18.470± 0.006 mag (random errors only) up to the LMC centre. We show that once the (m− M)0 values or each subregion are assumed to be identical to those derived from this best-fitting plane, systematic errors in t he SFR(t) and AMR are reduced by a factor of about two.

A NEW GENERATION OF PARSEC-COLIBRI STELLAR ISOCHRONES INCLUDING THE TP-AGB PHASE

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.

Dust Formation in a Galaxy with Primitive Abundances

Interstellar dust plays a crucial role in the evolution of galaxies. It governs the chemistry and physics of the interstellar medium. In the local universe, dust forms primarily in the ejecta from stars, but its composition and origin in galaxies at very early times remain controversial. We report observational evidence of dust forming around a carbon star in a nearby galaxy with a low abundance of heavy elements, 25 times lower than the solar abundance. The production of dust by a carbon star in a galaxy with such primitive abundances raises the possibility that carbon stars contributed carbonaceous dust in the early universe.

The VMC survey - XIV : First results on the look-back time star formation rate tomography of the Small Magellanic Cloud

We analyse deep images from the VISTA survey of the Magellanic Clouds in the YJKs filters, covering 14 sqrdeg (10 tiles), split into 120 subregions, and comprising the main body and Wing of the Small Magellanic Cloud (SMC). We apply a colour--magnitude diagram reconstruction method that returns their best-fitting star formation rate SFR(t), age-metallicity relation (AMR), distance and mean reddening, together with 68% confidence intervals. The distance data can be approximated by a plane tilted in the East-West direction with a mean inclination of 39 deg, although deviations of up to 3 kpc suggest a distorted and warped disk. After assigning to every observed star a probability of belonging to a given age-metallicity interval, we build high-resolution population maps. These dramatically reveal the flocculent nature of the young star-forming regions and the nearly smooth features traced by older stellar generations. They document the formation of the SMC Wing at ages <0.2 Gyr and the peak of star formation in the SMC Bar at 40 Myr. We clearly detect periods of enhanced star formation at 1.5 Gyr and 5 Gyr. The former is possibly related to a new feature found in the AMR, which suggests ingestion of metal-poor gas at ages slightly larger than 1 Gyr. The latter constitutes a major period of stellar mass formation. We confirm that the SFR(t) was moderately low at even older ages.

ON THE EXCITATION AND FORMATION OF CIRCUMSTELLAR FULLERENES

We compare and analyze the Spitzer mid-infrared spectrum of three fullerene-rich planetary nebulae in the Milky Way and the Magellanic Clouds: Tc1, SMP SMC 16, and SMP LMC 56. The three planetary nebulae share many spectroscopic similarities. The strongest circumstellar emission bands correspond to the infrared active vibrational modes of the fullerene species C60 and little or no emission is present from polycyclic aromatic hydrocarbons. The strengths of the fullerene bands in the three planetary nebulae are very similar, while the ratios of the [Ne III]15.5 μm/[Ne II]12.8 μm fine structure lines, an indicator of the strength of the radiation field, are markedly different. This raises questions about their excitation mechanism and we compare the fullerene emission to fluorescent and thermal models. In addition, the spectra show other interesting and common features, most notably in the 6-9 μm region, where a broad plateau with substructure dominates the emission. These features have previously been associated with mixtures of aromatic/aliphatic hydrocarbon solids. We hypothesize on the origin of this band, which is likely related to the fullerene formation mechanism, and compare it with modeled hydrogenated amorphous carbon that present emission in this region.