Researcher profile

M. A. T. Groenewegen

M. A. T. Groenewegen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
85works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

85 published item(s)

preprint2022arXiv

A new and homogeneous metallicity scale for Galactic classical Cepheids II. The abundance of iron and alpha elements

Classical Cepheids are the most popular distance indicators and tracers of young stellar populations. The key advantage is that they are bright and they can be easily identified in Local Group and Local Volume galaxies. Their evolutionary and pulsation properties depend on their chemical abundances. The main aim of this investigation is to perform a new and accurate abundance analysis of two tens of calibrating Galactic Cepheids using high spectral resolution (R$\sim$40,000-115,000) and high S/N spectra ($\sim$400) covering the entire pulsation cycle. We focus our attention on possible systematics affecting the estimate of atmospheric parameters and elemental abundances along the pulsation cycle. We cleaned the line list by using atomic transition parameters based on laboratory measurements and by removing lines that are either blended or display abundance variations along the pulsation cycle. The spectroscopic approach that we developed brings forward small dispersions in the variation of the atmospheric parameters ($σ$($T_{\rm eff}$)$\sim$50 K, $σ$($\log{g}$)$\sim$0.2 dex, and $σ$($ξ$)$\sim$0.2 km/s) and in the abundance of both iron ($\lesssim$ 0.05 dex) and alpha elements ($\lesssim$0.10 dex) over the entire pulsation cycle. We also provide new and accurate effective temperature templates by splitting the calibrating Cepheids into four different period bins, ranging from short to long periods. For each period bin, we performed an analytical fit with Fourier series providing $θ= 5040/{T_{\rm eff}}$ as a function of the pulsation phase. The current findings are a good viaticum to trace the chemical enrichment of the Galactic thin disk by using classical Cepheids and a fundamental stepping stone for further investigations into the more metal-poor regime typical of Magellanic Cepheids.

preprint2022arXiv

The VMC Survey -- XLVIII. Classical Cepheids unveil the 3D geometry of the LMC

We employed the {\it VISTA near-infrared $YJK_\mathrm{s}$ survey of the Magellanic System} (VMC), to analyse the $Y,\,J,\,K_\mathrm{s}$ light curves of $δ$ Cepheid stars (DCEPs) in the Large Magellanic Cloud (LMC). Our sample consists of 4408 objects accounting for 97 per cent of the combined list of OGLE\,IV and {\it Gaia}\,DR2 DCEPs. We determined a variety of period-luminosity ($PL$) and period-Wesenheit $PW$ relationships for Fundamental (F) and First Overtone (1O) pulsators. We discovered for the first time a break in these relationships for 1O DCEPs at $P$=0.58 d. We derived relative individual distances for DCEPs in the LMC with a precision of $\sim$1 kpc, calculating the position angle of the line of nodes and inclination of the galaxy: $θ$=145.6$\pm$1.0 deg and $i$=25.7$\pm$0.4 deg. The bar and the disc are seen under different viewing angles. We calculated the ages of the pulsators, finding two main episodes of DCEP formation lasting $\sim$40 Myr which happened 93 and 159 Myr ago. Likely as a result of its past interactions with the SMC, the LMC shows a non-planar distribution, with considerable structuring: the bar is divided into two distinct portions, the eastern and the western displaced by more than 1 kpc from each other. Similar behaviour is shown by the spiral arms. The LMC disc appears "flared" and thick, with a disc scale height of $h\sim 0.97$ kpc. This feature can be explained by strong tidal interactions with the Milky Way and/or the Small Magellanic Cloud or past merging events with now disrupted LMC satellites.

preprint2020arXiv

Analysing the spectral energy distributions of Galactic classical Cepheids

(abridged) Spectral energy distributions (SEDs) were constructed for a sample of 477 classical cepheids (CCs). The SEDs were fitted with a dust radiative transfer code. Four stars showed a large mid- or far-infrared excess and the fitting then included a dust component. These comprise the well-known case of RS Pup, and three stars that are (likely) Type-II cepheids (T2Cs), AU Peg, QQ Per, and FQ Lac. The remainder of the sample was fitted with a stellar photosphere to derive the best-fitting luminosity and effective temperature. Distance and reddening were taken from the literature. The stars were plotted in a Hertzsprung-Russell diagram and compared to evolutionary tracks for cepheids and theoretical instability strips. For the large majority of stars, the position in the HRD is consistent with the instability strip for a CC or T2C. About 5% of the stars are outliers in the sense that they are much hotter or cooler than expected. A comparison to effective temperatures derived from spectroscopy suggests in some cases that the photometrically derived temperature is not correct and that this is likely linked to an incorrectly adopted reddening. In this work the presence of a small NIR excess, as has been proposed in the literature for a few well-known cepheids, is investigated. Firstly, this was done by using a sample of about a dozen stars for which a mid-infrared spectrum is available. Secondly, the SEDs of all stars were fitted with a dust model to see if a statistically significant better fit could be obtained. The results were compared to recent work. Eight new candidates for exhibiting a NIR excess are proposed, solely based on the photometric SEDs. Obtaining mid-infrared spectra would be needed to confirm this excess. Finally, period-bolometric luminosity and period-radius relations are presented for samples of over 370 fundamental-mode CCs.

preprint2020arXiv

Are extreme AGB stars post-common envelope binaries?

Modelling dust formation in single stars evolving through the carbon-star stage of the asymptotic giant branch (AGB) reproduces well the mid-infrared colours and magnitudes of most of the C-rich sources in the Large Magellanic Cloud (LMC), apart from a small subset of extremely red objects (EROs). The analysis of EROs spectral energy distribution suggests the presence of large quantities of dust, which demand gas densities in the outflow significantly higher than expected from theoretical modelling. We propose that binary interaction mechanisms that involve common envelope (CE) evolution could be a possible explanation for these peculiar stars; the CE phase is favoured by the rapid growth of the stellar radius occurring after C$/$O overcomes unity. Our modelling of the dust provides results consistent with the observations for mass-loss rates $\dot M \sim 5\times 10^{-4}~\dot M/$yr, a lower limit to the rapid loss of the envelope experienced in the CE phase. We propose that EROs could possibly hide binaries of orbital periods $\sim$days and are likely to be responsible for a large fraction of the dust production rate in galaxies.

preprint2020arXiv

Characterization of M-stars in the LMC in the JWST era

We study the M-type asymptotic giant branch (AGB) population of the Large Magellanic Cloud (LMC) by characterizing the individual sources in terms of the main properties of the progenitors and of the dust present in the circumstellar envelope. To this aim, we compare the combination of the spectroscopic and photometric data collected by Spitzer, complemented by additional photometric results available in the literature, with results from AGB modelling that include the description of dust formation in the wind. To allow the interpretation of a paucity stars likely evolving through the post-AGB phase, we extended the available evolutionary sequences to reach the PN phase. The main motivation of the present analysis is to prepare the future observations of the evolved stellar populations of Local Group galaxies that will be done by the James Webb Space Telescope (JWST), by identifying the combination of filters that will maximize the possibilities of characterizing the observed sources. The present results show that for the M-star case the best planes to be used for this purpose are the colour-magnitude ([F770W]-[F2550W], [F770W]) and (Ks-[F770W], [F770W]) planes. In these observational diagrams the sequences of low-mass stars evolving in the AGB phases before the achievement of the C-star stage and of massive AGBs experiencing hot bottom burning are clearly separated and peculiar sources, such as post-AGB, dual-dust chemistry and iron-dust stars can be easily identified

preprint2020arXiv

DEATHSTAR: Nearby AGB stars with the Atacama Compact Array I. CO envelope sizes and asymmetries: A new hope for accurate mass-loss-rate estimates

This is the first publication of the DEATHSTAR project. The goal of the project is to reduce the uncertainties of observational estimates of mass-loss rates from Asymptotic Giant Branch (AGB) stars. Line emission from 12CO J=2-1 and 3-2 were mapped using the ACA. In this initial analysis, the emission distribution was fit to a Gaussian distribution in the uv-plane. Detailed radiative transfer analysis will be presented in the future. The axes of the best-fit Gaussian at the line center of the 12CO J=2-1 emission gives a first indication of the size of the emitting region. Furthermore, the fitting results, such as the major and minor axis, center position, and the goodness of fit across both lines, constrain the symmetry of the emission distribution. We find that the CO envelope sizes are, in general, larger for C-type than for M-type AGB stars, which is expected if the CO/H2 ratio is larger in C-type stars. Furthermore, a relation between the 12CO J=2-1 size and circumstellar density is shown that, while in broad agreement with photodissociation calculations, reveals large scatter and systematic differences between the stellar types. The majority of the sources have CO envelopes that are consistent with a spherically symmetric, smooth outflow. For about a third of the sources, indications of strong asymmetries are found. This is consistent with previous interferometric investigations of northern sources. Smaller scale asymmetries are found in a larger fraction of sources. These results for CO envelope radii and shapes can be used to constrain detailed radiative transfer modeling of the same stars so as to determine mass-loss rates that are independent of photodissociation models. For a large fraction of the sources, observations at higher spatial resolution will be necessary to further investigate the complex circumstellar dynamics revealed by our ACA observations.

preprint2020arXiv

Extended view on the dust shells around two carbon stars

Stars on the asymptotic giant branch (AGB) lose considerable amounts of matter through their dust-driven stellar winds. A number of such sources have been imaged by Herschel/PACS, revealing a diverse sample of different morphological types. Among them are a few examples which show geometrically thin, spherically symmetric shells which can be used to probe the mass loss history of their host stars. We aim to determine the physical properties of the dust envelope around the two carbon stars U Hya and W Ori. With the much-improved spatial constraints from the new far-infrared maps, our primary goal is to measure the dust masses contained in the shells and see how they fit the proposed scenarios of shell formation. We calculated the radiative transfer of the circumstellar dust envelope using the 1D code More of DUSTY (MoD). Adopting a parametrised density profile, we obtained a best-fit model in terms of the photometric and spectroscopic data, as well as a radial intensity profile based on Herschel/PACS data. For the case of U Hya, we also computed a grid of circumstellar envelopes by means of a stationary wind code and compare the results of the two modelling approaches. The Herschel/PACS maps show U Hya surrounded by a detached shell of $114''\ (0.12\,\mathrm{pc})$ in radius, confirming the observations from previous space missions. The dust masses calculated for the shell by the two approaches are consistent with respect to the adopted dust grain properties. In addition, around W Ori, we detect for the first time a weak spherically symmetric structure with a radius of $92''\ (0.17\,\mathrm{pc})$ and a dust mass of $(3.5\pm0.3)\times10^{-6}\,\mathrm{M_\odot}$.

preprint2020arXiv

Modelling the Milky Way. I -- Method and first results fitting the thick disk and halo with DES-Y3 data

We present MWFitting, a method to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) from TRILEGAL models to real data. We apply MWFitting to photometric data from the first three years of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud and the Sagittarius Stream, our main sample spans a total area of $\sim$2,300 deg$^2$ distributed across the DES footprint. We further explore a smaller subset ($\sim$ 1,300 deg$^2$) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3\%. Based on the best-fit models, we create simulated stellar catalogues covering the whole DES footprint down to $g = 24$ magnitude. Comparisons of data and simulations provide evidence for a break in the power law index describing the stellar density of the Milky Way (MW) halo. Several previously discovered stellar over-densities are recovered in the residual stellar density map, showing the reliability of MWFitting in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict MW star counts for surveys such as LSST and Euclid.

preprint2020arXiv

The flux-weighted gravity-luminosity relation of Galactic classical Cepheids

The flux-weighted gravity-luminosity relation (FWGLR) is investigated for a sample of 477 classical Cepheids (CCs), including stars that have been classified in the literature as such but are probably not. The luminosities are taken from the literature, based on the fitting of the spectral energy distributions (SEDs) assuming a certain distance and reddening. The flux-weighted gravity (FWG) is taken from gravity and effective temperature determinations in the literature based on high-resolution spectroscopy. There is a very good agreement between the theoretically predicted and observed FWG versus pulsation period relation that could serve in estimating the FWG (and $\log g$) in spectroscopic studies with a precision of 0.1~dex. As was known in the literature, the theoretically predicted FWGLR relation for CCs is very tight and is not very sensitive to metallicity (at least for LMC and solar values), rotation rate, and crossing of the instability strip. The observed relation has a slightly different slope and shows more scatter (0.54~dex). This is due both to uncertainties in the distances and to the pulsation phase averaged FWG values. Data from future Gaia data releases should reduce these errors, and then the FWGLR could serve as a powerful tool in Cepheid studies.

preprint2020arXiv

The VMC Survey -- XXXVII. Pulsation periods of dust enshrouded AGB stars in the Magellanic Clouds

(abridged)Variability is a key property of stars on the asymptotic giant branch (AGB). Their pulsation period is related to the luminosity and mass-loss rate of the star. The long-period variables (LPVs) and Mira variables are the most prominent of all types of variability of evolved stars. The reddest, most obscured AGB stars are too faint in the optical and have eluded large variability surveys. Selection criteria are derived based on colour-colour and colour-magnitude diagrams from the combination of VISTA Magellanic Cloud (VMC) survey, Spitzer IRAC and AllWISE data. After eliminating LPVs with known periods shorter than 450 days, a sample of 1299 candidate obscured AGB stars is selected. K-band LCs are constructed combining the epoch photometry available in the VMC survey with literature data, analysed for variability and fitted with a single period sine curve to derive mean magnitudes, amplitudes and periods. A subset of 254 stars are either new variables, or known variables where the period we find is better determined than the literature value, or longer than 1000 days. The spectral energy distributions (SEDs) of these stars are fitted to a large number of templates. For this purpose the SEDs and Spitzer IRS spectra of some non-AGB stars are also fitted to have templates of the most likely contaminants in the sample. A sample of 217 likely LPVs is found. Thirty-four stars have periods longer than 1000 days although some of them have alternative shorter periods. Mass-loss rates and luminosities are estimated from the template fitting. Period-luminosity relations are presented for C- and O-rich Miras that appear to be extensions of relations derived in the literature for shorter periods. The fit for the C-stars is particularly well defined (with 182 objects) and reads Mbol = (-2.27 \pm 0.20) log P + (1.45 \pm 0.54)mag with an rms of 0.41 mag.

preprint2016arXiv

Constraints on the H2O formation mechanism in the wind of carbon-rich AGB stars

Context. The recent detection of warm H$_2$O vapor emission from the outflows of carbon-rich asymptotic giant branch (AGB) stars challenges the current understanding of circumstellar chemistry. Two mechanisms have been invoked to explain warm H$_2$O vapor formation. In the first, periodic shocks passing through the medium immediately above the stellar surface lead to H$_2$O formation. In the second, penetration of ultraviolet interstellar radiation through a clumpy circumstellar medium leads to the formation of H$_2$O molecules in the intermediate wind. Aims. We aim to determine the properties of H$_2$O emission for a sample of 18 carbon-rich AGB stars and subsequently constrain which of the above mechanisms provides the most likely warm H$_2$O formation pathway. Methods, Results, and Conclusions. See paper.

preprint2016arXiv

Herschel observations of the nebula M1-67 around the Wolf-Rayet star WR 124

Infrared Herschel imaging and spectroscopic observations of the nebula M1-67 around the Wolf-Rayet star WR 124 have been obtained along with optical imaging observations. The infrared images reveal a clumpy dusty nebula that extends up to 1 pc. The comparison with the optical images shows that the ionized gas nebula coincides with the dust nebula, the dust and the gas being mixed together. A photodissociation region is revealed from the infrared spectroscopic analysis. The analysis of the infrared spectrum of the nebula, where forbidden emission lines of ionized elements were detected, showed that the nebula consists of mildly processed material with the calculated abundance number ratios being N/O = 1.0 +/- 0.5 and C/O = 0.46 +/- 0.27. Based on a radiative transfer model, the dust mass of the nebula was estimated to be 0.22 Msun with a population of large grains being necessary to reproduce the observations. The comparison of the mass-loss rate and the abundance ratios to theoretical models of stellar evolution led to the conclusion that the nebular ejection took place during a RSG/YSG evolutionary phase of a central star with an initial mass of 32 Msun.

preprint2016arXiv

The VMC Survey - XX. Identification of new Cepheids in the Small Magellanic Cloud

We present Ks -band light curves for 299 Cepheids in the Small Magellanic Cloud (SMC) of which 288 are new discoveries that we have identified using multi-epoch near-infrared photometry obtained by the VISTA survey of the Magellanic Clouds system (VMC). The new Cepheids have periods in the range from 0.34 to 9.1 days and cover the magnitude interval 12.9 <= Ks <= 17.6 mag. Our method was developed using variable stars previously identified by the optical microlensing survey OGLE. We focus on searching new Cepheids in external regions of the SMC for which complete VMC Ks-band observations are available and no comprehensive identification of different types of variable stars from other surveys exists yet.

preprint2016arXiv

The VMC Survey. XIX. Classical Cepheids in the Small Magellanic Cloud

The VISTA near infrared Y, J, Ks survey of the Magellanic System (VMC) is collecting deep Ks band time series photometry of pulsating variable stars hosted by the two Magellanic Clouds and their connecting Bridge. In this paper, we present Y, J, Ks light curves for a sample of 4172 Small Magellanic Cloud (SMC) Classical Cepheids (CCs). These data, complemented with literature V values, allowed us to construct a variety of period-luminosity (PL), period-luminosity-color (PLC), and period-wesenheit (PW) relationships, valid for Fundamental (F), First Overtone (FO) and Second Overtone (SO) pulsators. The relations involving V, J, Ks bands are in agreement with their counterparts in the literature. As for the Y band, to our knowledge we present the first CC PL, PW, and PLC relations ever derived using this filter. We also present the first near infrared PL, PW, and PLC relations for SO pulsators to date. We used PW(V,Ks) to estimate the relative SMC-LMC distance and, in turn, the absolute distance to the SMC. For the former quantity we find a value of Delta(mu)=0.55+-0.04 mag, in rather good agreement with other evaluations based on CCs, but significantly larger than the results obtained from older population II distance indicators. This discrepancy might be due to the different geometric distributions of young and old tracers in both Clouds. As for the absolute distance to the SMC, our best estimates are mu(SMC)=19.01+-0.05 mag and mu(SMC)=19.04+-0.06 mag, based on two distance measurements to the LMC, which rely on accurate CC and eclipsing Cepheid binary data, respectively.

preprint2015arXiv

Modelling the atmosphere of the carbon-rich Mira RU Vir

Context. We study the atmosphere of the carbon-rich Mira RU Vir using the mid-infrared high spatial resolution interferometric observations from VLTI/MIDI. Aims. The aim of this work is to analyse the atmosphere of the carbon-rich Mira RU Vir, with state of the art models, in this way deepening the knowledge of the dynamic processes at work in carbon-rich Miras. Methods. We compare spectro-photometric and interferometric measurements of this carbon-rich Mira AGB star, with the predictions of different kinds of modelling approaches (hydrostatic model atmospheres plus MOD-More Of Dusty, self-consistent dynamic model atmospheres). A geometric model fitting tool is used for a first interpretation of the interferometric data. Results. The results show that a joint use of different kind of observations (photometry, spectroscopy, interferometry) is essential to shed light on the structure of the atmosphere of a carbon-rich Mira. The dynamic model atmospheres fit well the ISO spectrum in the wavelength range λ = [2.9, 25.0] μm. Nevertheless, a discrepancy is noticeable both in the SED (visible), and in the visibilities (shape and level). A possible explanation are intra-/inter-cycle variations in the dynamic model atmospheres as well as in the observations. The presence of a companion star and/or a disk or a decrease of mass loss within the last few hundred years cannot be excluded but are considered unlikely.

preprint2015arXiv

New NIR light-curve templates for classical Cepheids

We present new near-infrared (NIR) light-curve templates for fundamental (FU, JHK) and first overtone (FO, J) Cepheids. The new templates together with PL and PW relations provide Cepheid distances from single-epoch observations with a precision only limited by the intrinsic accuracy of the method adopted. The templates rely on a very large set of Galactic and Magellanic Clouds (MCs) Cepheids (FU,~600; FO,~200) with well sampled NIR (IRSF data) and optical (V,I; OGLE data) light curves. To properly trace the change in the shape of the light curve as a function of period, we split the sample of calibrating Cepheids into 10 different period bins. The templates for the first time cover FO Cepheids and the FU short-period Cepheids (P<5 days). Moreover, the zero-point phase is anchored to the phase of the mean magnitude along the rising branch. The new approach has several advantages in sampling the light curve of bump Cepheids when compared with the phase of maximum light. We also provide new estimates of the NIR-to-optical amplitude ratios for FU and FO Cepheids. We perform detailed analytical fits using both 7th-order Fourier series and multi-Gaussian periodic functions. The latter are characterized by a smaller number of free parameters (9 vs 15). Mean NIR magnitudes based on the new templates are up to 80% more accurate than single-epoch measurements and up to 50% more accurate than mean magnitudes based on previous templates, with typical associated uncertainties ranging from 0.015 mag (J) to 0.019 mag (K). Moreover, the errors on individual distances of Small MC Cepheids derived from NIR PW relations, are essentially reduced to the intrinsic scatter of the adopted relations. Thus, the new templates are the ultimate tool to estimate precise Cepheid distances from NIR single-epoch observations, which can be adopted to derive the 3D structure of the MCs.

preprint2015arXiv

SAGE-Var: An Infrared Survey of Variability in the Magellanic Clouds

We present the first results from the SAGE-Var program, a follow on to the Spitzer legacy program Surveying the Agents of Galaxy Evolution (SAGE; Meixner, et al. 2006). We obtained 4 epochs of photometry at 3.6 & 4.5 microns covering the bar of the Large Magellanic Cloud (LMC) and the central region of the Small Magellanic Cloud (SMC) in order to probe the variability of extremely red sources missed by variability surveys conducted at shorter wavelengths, and to provide additional epochs of observation for known variables. Our 6 total epochs of observations allow us to probe infrared variability on 15 different timescales ranging from ~20 days to ~5 years. Out of a full catalog of 1,717,554 (LMC) and 457,760 (SMC) objects, we find 10 (LMC) and 6 (SMC) large amplitude AGB variables without optically measured variability owing to circumstellar dust obscuration. The catalog also contains multiple observations of known AGB variables, type I and II Cepheids, eclipsing variables, R CrB stars and young stellar objects which will be discussed in following papers. Here we present infrared Period-Luminosity (PL) relations for classical Cepheids in the Magellanic Clouds, as well as improved PL relationships for AGB stars pulsating in the fundamental mode using mean magnitudes constructed from 6 epochs of observations.

preprint2015arXiv

The Herschel view of the nebula around the luminous blue variable star AG Carinae

Far-infrared Herschel PACS imaging and spectroscopic observations of the nebula around the luminous blue variable (LBV) star AG Car have been obtained along with optical imaging in the Halpha+[NII] filter. In the infrared light, the nebula appears as a clumpy ring shell that extends up to 1.2 pc with an inner radius of 0.4 pc. It coincides with the Halpha nebula, but extends further out. Dust modeling of the nebula was performed and indicates the presence of large grains. The dust mass is estimated to be ~ 0.2 Msun. The infrared spectrum of the nebula consists of forbidden emission lines over a dust continuum. Apart from ionized gas, these lines also indicate the existence of neutral gas in a photodissociation region that surrounds the ionized region. The abundance ratios point towards enrichment by processed material. The total mass of the nebula ejected from the central star amounts to ~ 15 Msun, assuming a dust-to-gas ratio typical of LBVs. The abundances and the mass-loss rate were used to constrain the evolutionary path of the central star and the epoch at which the nebula was ejected, with the help of available evolutionary models. This suggests an ejection during a cool LBV phase for a star of ~ 55 Msun with little rotation.

preprint2014arXiv

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.

preprint2014arXiv

An Infrared Census of DUST in Nearby Galaxies with Spitzer (DUSTiNGS), I. Overview

Nearby resolved dwarf galaxies provide excellent opportunities for studying the dust-producing late stages of stellar evolution over a wide range of metallicity (-2.7 < [Fe/H] < -1.0). Here, we describe DUSTiNGS (DUST in Nearby Galaxies with Spitzer): a 3.6 and 4.5 micron post-cryogen Spitzer Space Telescope imaging survey of 50 dwarf galaxies within 1.5 Mpc that is designed to identify dust-producing Asymptotic Giant Branch (AGB) stars and massive stars. The survey includes 37 dwarf spheroidal, 8 dwarf irregular, and 5 transition-type galaxies. This near-complete sample allows for the building of statistics on these rare phases of stellar evolution over the full metallicity range. The photometry is >75% complete at the tip of the Red Giant Branch for all targeted galaxies, with the exception of the crowded inner regions of IC 10, NGC 185, and NGC 147. This photometric depth ensures that the majority of the dust-producing stars, including the thermally-pulsing AGB stars, are detected in each galaxy. The images map each galaxy to at least twice the half-light radius to ensure that the entire evolved star population is included and to facilitate the statistical subtraction of background and foreground contamination, which is severe at these wavelengths. In this overview, we describe the survey, the data products, and preliminary results. We show evidence for the presence of dust-producing AGB stars in 8 of the targeted galaxies, with metallicities as low as [Fe/H] = -1.9, suggesting that dust production occurs even at low metallicity.

preprint2014arXiv

An Infrared Census of DUST in Nearby Galaxies with Spitzer (DUSTiNGS), II. Discovery of Metal-poor Dusty AGB Stars

The DUSTiNGS survey (DUST in Nearby Galaxies with Spitzer) is a 3.6 and 4.5 micron imaging survey of 50 nearby dwarf galaxies designed to identify dust-producing Asymptotic Giant Branch (AGB) stars and massive stars. Using two epochs, spaced approximately six months apart, we identify a total of 526 dusty variable AGB stars (sometimes called &#34;extreme&#34; or x-AGB stars; [3.6]-[4.5]>0.1 mag). Of these, 111 are in galaxies with [Fe/H] < -1.5 and 12 are in galaxies with [Fe/H] < -2.0, making them the most metal-poor dust-producing AGB stars known. We compare these identifications to those in the literature and find that most are newly discovered large-amplitude variables, with the exception of approximately 30 stars in NGC 185 and NGC 147, one star in IC 1613, and one star in Phoenix. The chemical abundances of the x-AGB variables are unknown, but the low metallicities suggest that they are more likely to be carbon-rich than oxygen-rich and comparisons with existing optical and near-IR photometry confirms that 70 of the x-AGB variables are confirmed or likely carbon stars. We see an increase in the pulsation amplitude with increased dust production, supporting previous studies suggesting that dust production and pulsation are linked. We find no strong evidence linking dust production with metallicity, indicating that dust can form in very metal-poor environments.

preprint2014arXiv

Discovery of Time Variation of the Intensity of Molecular Lines in IRC+10216 in The Submillimeter and Far Infrared Domains

We report on the discovery of strong intensity variations in the high rotational lines of abundant molecular species towards the archetypical circumstellar envelope of IRC+10216. The observations have been carried out with the HIFI instrument on board \textit{Herschel}\thanks{\textit{Herschel} is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA} and with the IRAM\thanks{This work was based on observations carried out with the IRAM 30-meter telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain)} 30-m telescope. They cover several observing periods spreading over 3 years. The line intensity variations for molecules produced in the external layers of the envelope most probably result from time variations in the infrared pumping rates. We analyze the main implications this discovery has on the interpretation of molecular line emission in the envelopes of Mira-type stars. Radiative transfer calculations have to take into account both the time variability of infrared pumping and the possible variation of the dust and gas temperatures with stellar phase in order to reproduce the observation of molecular lines at different epochs. The effect of gas temperature variations with stellar phase could be particularly important for lines produced in the innermost regions of the envelope. Each layer of the circumstellar envelope sees the stellar light radiation with a different lag time (phase). Our results show that this effect must be included in the models. The sub-mm and FIR lines of AGB stars cannot anymore be considered as safe intensity calibrators.

preprint2014arXiv

Dusty shells surrounding the carbon variables S Scuti and RT Capricorni

For the Mass-loss of Evolved StarS (MESS) programme, the unprecedented spatial resolution of the PACS photometer on board the Herschel space observatory was employed to map the dusty environments of asymptotic giant branch (AGB) and red supergiant (RSG) stars. Among the morphologically heterogeneous sample, a small fraction of targets is enclosed by spherically symmetric detached envelopes. Based on observations in the 70 μm and 160 μm wavelength bands, we investigated the surroundings of the two carbon semiregular variables S Sct and RT Cap, which both show evidence for a history of highly variable mass-loss. S Sct exhibits a bright, spherically symmetric detached shell, 138&#34; in diameter and co-spatial with an already known CO structure. Moreover, weak emission is detected at the outskirts, where the morphology seems indicative of a mild shaping by interaction of the wind with the interstellar medium, which is also supported by the stellar space motion. Two shells are found around RT Cap that were not known so far in either dust emission or from molecular line observations. The inner shell with a diameter of 188&#34; shows an almost immaculate spherical symmetry, while the outer ~5&#39; structure is more irregularly shaped. MoD, a modification of the DUSTY radiative transfer code, was used to model the detached shells. Dust temperatures, shell dust masses, and mass-loss rates are derived for both targets.

preprint2014arXiv

Eclipsing binary stars in the Large Magellanic Cloud. Results from the EROS-2, OGLE and VMC surveys

We present a catalogue of 1768 eclipsing binary stars (EBs) detected in the Large Magellanic Cloud (LMC) by the second generation of the EROS survey (hereinafter EROS-2); 493 of them are new discoveries located in outer regions (out of the central bar) of the LMC. These sources were originally included in a list of candidate classical Cepheids (CCs) extracted from the EROS-2 catalogue on the basis of the period (0.89 $<P_{EROS}<$15.85 days) versus luminosity ($13.39 < \langle B_{EROS}\rangle<17.82$ mag) diagram. After visual inspection of the light curves we reclassified them as eclipsing binaries. They have blue colours ($B_{EROS} - R_{EROS} < $ 0.2 mag) hence we classed them as hot eclipsing binaries (HEBs) containing hot massive components: main sequence (MS) stars or blue giants. We present $K_{\rm s}$-band light curves for 999 binaries from our sample that have a counterpart in the VISTA near-infrared ESO public survey of the Magellanic Clouds system (VMC). We provide spectral classifications of 13 HEBs with existing spectroscopy. We divided our sample into contact-like binaries and detached/semi-detached systems based on both visual inspection and the parameters of the Fourier decomposition of the light curves and analysed the period-luminosity ($PL$) relations of the contact-like systems using the $R_{EROS}$ and $K_{\rm s}$ magnitudes at maximum light. The contact-like binaries in our sample do not follow $PL$ relations. We analysed the sample of contact binaries from the OGLE III catalogue and confirmed that $PL_I$ and $PL_{K_{\rm s}}$ sequences are defined only by eclipsing binaries containing a red giant component.

preprint2014arXiv

First detection of rotational CO line emission in an RGB star

Context: For stars with initial masses below about 1 solar mass, the mass loss during the first red giant branch (RGB) phase dominates mass loss in the later asymptotic giant branch (AGB) phase. Nevertheless, mass loss on the RGB is still often parameterised by a simple Reimers law in stellar evolution models. Aims: To try to detect CO thermal emission in a small sample of nearby RGB stars with reliable Hipparcos parallaxes that were shown to have infrared excess in an earlier paper. Methods: A sample of five stars was observed in the CO J=2-1 and J=3-2 lines with the IRAM and APEX telescopes. Results: One star, the one with the largest mass-loss rate based on the previous analysis of the spectral energy distribution, was detected. The expansion velocity is unexpectedly large at 12 km/s. The line profile and intensity are compared to the predictions from a molecular line emission code. The standard model predicts a double-peaked profile, while the observations indicate a flatter profile. A model that does fit the data has a much smaller CO envelope (by a factor of 3), and a CO abundance that is two times larger and/or a larger mass-loss rate than the standard model. This could indicate that the phase of large mass loss has only recently started. Conclusions: The detection of CO in an RGB star with a luminosity of only about 1300 solar luminosities and a mass-loss rate as low as a few 10(-9) solarmasses/yr is important and the results also raise new questions. However, ALMA observations are required in order to study the mass-loss process of RGB stars in more detail, both for reasons of sensitivity (6 hours of integration in superior weather at IRAM were needed to get a 4sigma detection in the object with the largest detection probability), and spatial resolution (to determine the size of the CO envelope).

preprint2014arXiv

Herschel imaging of the dust in the Helix Nebula (NGC 7293)

In our series of papers presenting the Herschel imaging of evolved planetary nebulae, we present images of the dust distribution in the Helix nebula (NGC 7293). Images at 70, 160, 250, 350, and 500 micron were obtained with the PACS and SPIRE instruments on board the Herschel satellite. The broadband maps show the dust distribution over the main Helix nebula to be clumpy and predominantly present in the barrel wall. We determined the spectral energy distribution of the main nebula in a consistent way using Herschel, IRAS, and Planck flux values. The emissivity index of 0.99 +/- 0.09, in combination with the carbon rich molecular chemistry of the nebula, indicates that the dust consists mainly of amorphous carbon. The dust excess emission from the central star disk is detected at 70 micron and the flux measurement agree with previous measurement. We present the temperature and dust column density maps. The total dust mass across the Helix nebula (without its halo) is determined to be 0.0035 solar mass at a distance of 216 pc. The temperature map shows dust temperatures between 22 and 42 K, which is similar to the kinetic temperature of the molecular gas, strengthening the fact that the dust and gas co-exist in high density clumps. Archived images are used to compare the location of the dust emission in the far infrared (Herschel) with the ionized (GALEX, Hbeta) and molecular hydrogen component. The different emission components are consistent with the Helix consisting of a thick walled barrel-like structure inclined to the line of sight. The radiation field decreases rapidly through the barrel wall.

preprint2014arXiv

Herschel spectral-mapping of the Helix Nebula (NGC 7293): Extended CO photodissociation and OH+ emission

The Helix Nebula (NGC 7293) is the closest planetary nebulae. Therefore, it is an ideal template for photochemical studies at small spatial scales in planetary nebulae. We aim to study the spatial distribution of the atomic and the molecular gas, and the structure of the photodissociation region along the western rims of the Helix Nebula as seen in the submillimeter range with Herschel. We use 5 SPIRE FTS pointing observations to make atomic and molecular spectral maps. We analyze the molecular gas by modeling the CO rotational lines using a non-local thermodynamic equilibrium (non-LTE) radiative transfer model. For the first time, we have detected extended OH+ emission in a planetary nebula. The spectra towards the Helix Nebula also show CO emission lines (from J= 4 to 8), [NII] at 1461 GHz from ionized gas, and [CI] (2-1), which together with the OH+ lines, trace extended CO photodissociation regions along the rims. The estimated OH+ column density is (1-10)x1e12 cm-2. The CH+ (1-0) line was not detected at the sensitivity of our observations. Non-LTE models of the CO excitation were used to constrain the average gas density (n(H2)=(1-5)x1e5 cm-3) and the gas temperature (Tk= 20-40 K). The SPIRE spectral-maps suggest that CO arises from dense and shielded clumps in the western rims of the Helix Nebula whereas OH+ and [CI] lines trace the diffuse gas and the UV and X-ray illuminated clumps surface where molecules reform after CO photodissociation. [NII] traces a more diffuse ionized gas component in the interclump medium.

preprint2014arXiv

Herschel/PACS observations of the 69 $μm$ band of crystalline olivine around evolved stars

We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $μ$m. This wavelength range covers the 69 $μ$m band of crystalline olivine ($\text{Mg}_{2-2x}\text{Fe}_{(2x)}\text{SiO}_{4}$). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, $\dot M \ge 10^{-5}$ M$_\odot$/yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. The goal of this study is to exploit the spectral properties of the 69 $μ$m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. We fit the 69 $μ$m band and use its width and position to probe the composition and temperature of the crystalline olivine. For 27 sources in the sample, we detected the 69 $μ$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $μ$m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 $μ$m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 $μ$m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. continued...

preprint2014arXiv

Large-scale environments of binary AGB stars probed by Herschel. II: Two companions interacting with the wind of pi1 Gruis

Context. The Mass loss of Evolved StarS (MESS) sample observed with PACS on board the Herschel Space Observatory revealed that several asymptotic giant branch (AGB) stars are surrounded by an asymmetric circumstellar envelope (CSE) whose morphology is most likely caused by the interaction with a stellar companion. The evolution of AGB stars in binary systems plays a crucial role in understanding the formation of asymmetries in planetary nebulæ (PNe), but at present, only a handful of cases are known where the interaction of a companion with the stellar AGB wind is observed. Aims. We probe the environment of the very evolved AGB star $π^1$ Gruis on large and small scales to identify the triggers of the observed asymmetries. Methods. Observations made with Herschel/PACS at 70 $μ$m and 160 $μ$m picture the large-scale environment of $π^1$ Gru. The close surroundings of the star are probed by interferometric observations from the VLTI/AMBER archive. An analysis of the proper motion data of Hipparcos and Tycho-2 together with the Hipparcos Intermediate Astrometric Data help identify the possible cause for the observed asymmetry. Results. The Herschel/PACS images of $π^1$ Gru show an elliptical CSE whose properties agree with those derived from a CO map published in the literature. In addition, an arc east of the star is visible at a distance of $38^{\prime\prime}$ from the primary. This arc is most likely part of an Archimedean spiral caused by an already known G0V companion that is orbiting the primary at a projected distance of 460 au with a period of more than 6200 yr. However, the presence of the elliptical CSE, proper motion variations, and geometric modelling of the VLTI/AMBER observations point towards a third component in the system, with an orbital period shorter than 10 yr, orbiting much closer to the primary than the G0V star.

preprint2014arXiv

The VMC ESO Public Survey

The VISTA near-infrared YJKs survey of the Magellanic Clouds system (VMC) has entered its core phase: about 40% of the observations across the Large and Small Magellanic Clouds (LMC, SMC), the Magellanic Bridge and Stream have already been secured and the data are processed and analysed regularly. The initial analyses, concentrated in the first two completed tiles in the LMC (6_6 including 30 Doradus and 8_8 including the South Ecliptic Pole), show the superior quality of the data. The depth of the VMC survey allows the derivation of the star formation history (SFH) with unprecedented quality compared to previous wide-area surveys while reddening maps of high angular resolution are constructed using red clump stars. The multi-epoch Ks-band data reveal tight period-luminosity relations for variable stars and they permit the measurement of accurate proper motions of the stellar populations. The VMC survey continues to acquire data that will address many issues in the field of star and galaxy evolution.

preprint2014arXiv

The VMC Survey. XIII. Type II Cepheids in the Large Magellanic Cloud

The VISTA survey of the Magellanic Clouds System (VMC) is collecting deep $K_\mathrm{s}$--band time--series photometry of the pulsating variable stars hosted in the system formed by the two Magellanic Clouds and the Bridge connecting them. In this paper we have analysed a sample of 130 Large Magellanic Cloud (LMC) Type II Cepheids (T2CEPs) found in tiles with complete or near complete VMC observations for which identification and optical magnitudes were obtained from the OGLE III survey. We present $J$ and $K_\mathrm{s}$ light curves for all 130 pulsators, including 41 BL Her, 62 W Vir (12 pW Vir) and 27 RV Tau variables. We complement our near-infrared photometry with the $V$ magnitudes from the OGLE III survey, allowing us to build a variety of Period-Luminosity ($PL$), Period-Luminosity-Colour ($PLC$) and Period-Wesenheit ($PW$) relationships, including any combination of the $V, J, K_\mathrm{s}$ filters and valid for BL Her and W Vir classes. These relationships were calibrated in terms of the LMC distance modulus, while an independent absolute calibration of the $PL(K_\mathrm{s})$ and the $PW(K_\mathrm{s},V)$ was derived on the basis of distances obtained from $Hubble Space Telescope$ parallaxes and Baade-Wesselink technique. When applied to the LMC and to the Galactic Globular Clusters hosting T2CEPs, these relations seem to show that: 1) the two population II standard candles RR Lyrae and T2CEPs give results in excellent agreement with each other; 2) there is a discrepancy of $\sim$0.1 mag between population II standard candles and Classical Cepheids when the distances are gauged in a similar way for all the quoted pulsators. However, given the uncertainties, this discrepancy is within the formal 1$σ$ uncertainties.

preprint2014arXiv

The wind of W Hya as seen by Herschel - I. The CO envelope

Asymptotic giant branch (AGB) stars lose their envelopes by means of a stellar wind whose driving mechanism is not understood well. Characterizing the composition and thermal and dynamical structure of the outflow provides constraints that are essential for understanding AGB evolution, including the rate of mass loss and isotopic ratios. We modeled the CO emission from the wind of the low mass-loss rate oxygen-rich AGB star W Hya using data obtained by the HIFI, PACS, and SPIRE instruments onboard the Herschel Space Observatory and ground-based telescopes. $^{12}$CO and $^{13}$CO lines are used to constrain the intrinsic $^{12}$C/$^{13}$C ratio from resolved HIFI lines. The acceleration of the outflow up to about 5.5 km/s is quite slow and can be represented by a beta-type velocity law with index 5. Beyond this point, acceleration up the terminal velocity of 7 km/s is faster. Using the J=10-9, 9-8, and 6-5 transitions, we find an intrinsic $^{12}$C/$^{13}$C ratio of $18\pm10$ for W Hya, where the error bar is mostly due to uncertainties in the $^{12}$CO abundance and the stellar flux around 4.6 $μ$m. To match the low-excitation CO lines, these molecules need to be photo-dissociated at about 500 stellar radii. The radial dust emission intensity profile measured by PACS images at 70 $μ$m shows substantially stronger emission than our model predicts beyond 20 arcsec. The initial slow acceleration of the wind implies inefficient wind driving in the lower part of the envelope. The final injection of momentum in the wind might be the result of an increase in the opacity thanks to the late condensation of dust species. The derived intrinsic isotopologue ratio for W Hya is consistent with values set by the first dredge-up and suggestive of an initial mass of 2 M$_\odot$ or more. However, the uncertainty in the main-sequence mass derived based on this isotopologic ratio is large.

preprint2014arXiv

The wind of W Hya as seen by Herschel. II. The molecular envelope of W Hya

The evolution of low- and intermediate-mass stars on the asymptotic giant branch (AGB) is mainly controlled by the rate at which these stars lose mass in a stellar wind. Understanding the driving mechanism and strength of the stellar winds of AGB stars and the processes enriching their surfaces with products of nucleosynthesis are paramount to constraining AGB evolution and predicting the chemical evolution of galaxies. In a previous paper we have constrained the structure of the outflowing envelope of W Hya using spectral lines of the $^{12}$CO molecule. Here we broaden this study by modelling an extensive set of H$_{2}$O and $^{28}$SiO lines observed by the three instruments on board Herschel using a state-of-the-art molecular excitation and radiative transfer code. The oxygen isotopic ratios and the $^{28}$SiO abundance profile can be connected to the initial stellar mass and to crucial aspects of dust formation at the base of the stellar wind, respectively. The modelling of H$_{2}$O and $^{28}$SiO confirms the properties of the envelope model of W Hya derived from $^{12}$CO lines. We find an H$_2$O ortho-to-para ratio of 2.5\,$^{+2.5}_{-1.0}$, consistent with what is expected for an AGB wind. The O$^{16}$/O$^{17}$ ratio indicates that W Hya has an initial mass of about 1.5 M$_\odot$. Although the ortho- and para-H$_{2}$O lines observed by HIFI appear to trace gas of slightly different physical properties, a turbulence velocity of $0.7\pm0.1$ km s$^{-1}$ fits the HIFI lines of both spin isomers and those of $^{28}$SiO well. The ortho- and para-H$_2^{16}$O and $^{28}$SiO abundances relative to H$_{2}$ are $(6^{+3}_{-2}) \times 10^{-4}$, $(3^{+2}_{-1}) \times 10^{-4}$, and $(3.3\pm 0.8)\times 10^{-5}$, respectively. Assuming a solar silicon-to-carbon ratio, the $^{28}$SiO line emission model is consistent with about one-third of the silicon atoms being locked up in dust particles.

preprint2013arXiv

A Herschel study of NGC 650

As part of the Herschel Guaranteed Time Key Project MESS (Mass loss of Evolved StarS) we have imaged a sample of planetary nebulae. In this paper we present the PACS and SPIRE images of the classical bipolar planetary nebula NGC 650. We used these images to derive a temperature map of the dust. We also constructed a photoionization and dust radiative transfer model using the spectral synthesis code Cloudy. To constrain this model, we used the PACS and SPIRE fluxes and combined these with hitherto unpublished IUE and Spitzer IRS spectra as well as various other data from the literature. The temperature map combined with the photoionization model were used to study various aspects of the central star, the nebula, and in particular the dust grains in the nebula. The central star parameters are determined to be T_eff = 208 kK and L = 261 L_sol assuming a distance of 1200 pc. The stellar temperature is much higher than previously published values. We confirm that the nebula is carbon-rich with a C/O ratio of 2.1. The nebular abundances are typical for a type IIa planetary nebula. With the photoionization model we determined that the grains in the ionized nebula are large (assuming single-sized grains, they would have a radius of 0.15 micron. Most likely these large grains were inherited from the asymptotic giant branch phase. The PACS 70/160 micron temperature map shows evidence for two radiation components heating the grains. The first component is direct emission from the central star, while the second component is diffuse emission from the ionized gas (mainly Ly alpha). We show that previous suggestions that there is a photo-dissociation region surrounding the ionized region are incorrect. The neutral material resides in dense clumps inside the ionized region. These may also harbor stochastically heated very small grains in addition to the large grains.

preprint2013arXiv

ALMA observations of the variable 12CO/13CO ratio around the asymptotic giant branch star R Sculptoris

[abridged] The 12CO/13CO ratio is often used as a measure of the 12C/13C ratio in the circumstellar environment, carrying important information about the stellar nucleosynthesis. External processes can change the 12CO and 13CO abundances, and spatially resolved studies of the 12CO/13CO ratio are needed to quantify the effect of these processes on the globally determined values. Additionally, such studies provide important information on the conditions in the circumstellar environment. The detached-shell source R Scl, displaying CO emission from recent mass loss, in a binary-induced spiral structure as well as in a clumpy shell produced during a thermal pulse, provides a unique laboratory for studying the differences in CO isotope abundances throughout its recent evolution. We observed both the 12CO(J=3-2) and the 13CO(J=3-2) line using ALMA. We find significant variations in the 12CO/13CO intensity ratios and consequently in the abundance ratios. The average CO isotope abundance ratio is at least a factor three lower in the shell (~19) than that in the present-day (<300 years) mass loss (>60). Additionally, variations in the ratio of more than an order of magnitude are found in the shell itself. We attribute these variations to the competition between selective dissociation and isotope fractionation in the shell, of which large parts cannot be warmer than ~35 K. However, we also find that the 12CO/13CO ratio in the present-day mass loss is significantly higher than the 12C/13C ratio determined in the stellar photosphere from molecular tracers (~19). The origin of this discrepancy is still unclear, but we speculate that it is due to an embedded source of UV-radiation that is primarily photo-dissociating 13CO. This radiation source could be the hitherto hidden companion. Alternatively, the UV-radiation could originate from an active chromosphere of R Scl itself....

preprint2013arXiv

Carbon Monoxide in the Cold Debris of Supernova 1987A

We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1&#39;&#39; of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations.

preprint2013arXiv

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.

preprint2013arXiv

The VMC Survey - X. Cepheids, RR Lyrae stars and binaries as probes of the Magellanic System&#39;s structure

The VMC survey is obtaining multi-epoch photometry in the Ks band of the Magellanic System down to a limiting magnitude of Ks ~ 19.3 for individual epoch data. The observations are spaced in time such as to provide optimal sampling of the light curves for RR Lyrae stars and for Cepheids with periods up to 20-30 days. We present examples of the Ks-band light curves of Classical Cepheids and RR Lyrae stars we are obtaining from the VMC data and outline the strategy we put in place to measure distances and infer the System three-dimensional geometry from the variable stars. For this purpose the near-infrared Period-Luminosity, Period-Wesenheit, and Period-Luminosity-Colour relations of the system RR Lyrae stars and Cepheids are used. We extensively exploit the catalogues of the Magellanic Clouds&#39; variable stars provided by the EROS-2 and OGLE III/IV microlensing surveys. By combining these surveys we present the currently widest-area view of the Large Magellanic Cloud as captured by the galaxy Cepheids, RR Lyrae stars and binaries. This reveals the full extent of the main structures (bar/s - spiral arms) that have only been vaguely guessed before. Our work strengthens the case for a detailed study of the Large Magellanic Cloud three-dimensional geometry.

preprint2013arXiv

The VMC survey VII. Reddening map of the 30 Doradus field and the structure of the cold interstellar medium

Context: The details of how galaxies have evolved are imprinted in their star formation history, chemical enrichment and morpho-kinematic structure. Parameters behind these can be measured by combining photometric techniques with modelling. However, there are uncertainties from the ambiguity of colour and magnitude and the effects of interstellar reddening. Aims: In this paper we present a detailed reddening map of the central 30 Doradus region of the Large Magellanic Cloud; for community use and to test the methods used. The reddening, a measurement of dust extinction, acts as a tracer of the interstellar medium (ISM). Methods: Near infrared photometry of red clump stars are used to measure reddening as extinction is the main cause of their colour and magnitude variance. The star formation history is used to convert colour to reddening values which are subsequently converted to visual extinction. Results: Presented is a dust map for the 30 Doradus field. This map samples a region of 1x1.5 deg, containing ~1.5x10^5 red clump stars which probe reddening up to AV=6 mag. We compare our map with maps from the literature, including optical extinction maps and radio, mid- and far-infrared maps of atomic hydrogen and dust emission. Through estimation of column density we locate molecular clouds. Conclusions: This new reddening map shows correlation with equivalent maps in the literature, validating the method of red clump star selection. We make our reddening map available for community use. In terms of ISM the red clump stars appear to be more affected by the cooler dust measured by 70 micron emission because there is stronger correlation between increasing emission and extinction due to red clump stars not being located near hot stars that would heat the dust. The transition from atomic hydrogen to molecular hydrogen occurs between densities of NH=4x10^21 to 6x10^21 cm^-2.

preprint2013arXiv

The VMC Survey. VIII. First results for Anomalous Cepheids

The VISTA near-infrared $YJK_\mathrm{s}$ survey of the Magellanic Clouds System (VMC, PI M.-R. L. Cioni) is collecting deep $K_\mathrm{s}$-band time-series photometry of the pulsating variable stars hosted in the system formed by the two Magellanic Clouds and the Bridge connecting them. In this paper we present for the first time $K_\mathrm{s}$-band light curves for Anomalous Cepheid (AC) variables. In particular, we have analysed a sample of 48 Large Magellanic Cloud ACs, for which identification and optical magnitudes were obtained from the OGLE III and IV catalogues. The VMC $K_\mathrm{s}$-band light curves for ACs are well sampled, with the number of epochs ranging from 8 to 16, and allowing us to obtain very precise mean $K_\mathrm{s}$ magnitudes with errors on average of the order of 0.01 mag. The $\langle K_\mathrm{s} \rangle$ values were used to build the first Period-Luminosity and Period-Wesenheit relations in the near-infrared for fundamental-mode and first overtone ACs. At the same time we exploited the optical ($V,I$) OGLE data to build accurate Period-Luminosity, Period-Luminosity-Colour and Period-Wesenheit relations both for fundamental-mode and first overtone ACs. For the first time these relations were derived from a sample of pulsators which uniformly cover the whole AC instability strip. The application of the optical Period-Wesenheit relation to a sample of dwarf galaxies hosting a significant population of ACs revealed that this relation is a valuable tool for deriving distances within the Local Group. Due to its lower dispersion, we expect the $K_\mathrm{s}$ Period-Wesenheit relations first derived in this paper to represent a valuable tool for measuring accurate distances to galaxies hosting ACs when more data in near-infrared filters become available.

preprint2012arXiv

An independent distance estimate to CW Leo

CW Leo has been observed six times between October 2009 and June 2012 with the SPIRE instrument on board the Herschel satellite. Variability has been detected in the flux emitted by the central star with a period of 639 \pm 4 days, in good agreement with determinations in the literature. Variability is also detected in the bow shock around CW Leo that had previously been detected in the ultraviolet and Herschel PACS/SPIRE data. Although difficult to prove directly, our working hypothesis is that this variability is directly related to that of the central star. In this case, fitting a sine curve with the period fixed to 639 days results in a time-lag in the variability between bow shock and the central star of 402 \pm 37 days. The orientation of the bow shock relative to the plane of the sky is unknown (but see below). For an inclination angle of zero degrees, the observed time-lag translates into a distance to CW Leo of 130 \pm 13 pc, and for non-zero inclination angles the distance is smaller. Fitting the shape of the bow shock with an analytical model (Wilkin 1996), the effect of the inclination angle on the distance may be estimated. Making the additional assumption that the relative peculiar velocity between the interstellar medium (ISM) and CW Leo is determined entirely by the star space velocity with respect to the local standard of rest (i.e. a stationary ISM), the inclination angle is found to be (-33.3 \pm 0.8) degrees based on the observed proper motion and radial velocity. Using the Wilkin model, our current best estimate of the distance to CW Leo is 123 \pm 14 pc. For a distance of 123 pc, we derive a mean luminosity of 7790 \pm 150 Lsol (internal error).

preprint2012arXiv

Baade-Wesselink distances to Galactic and Magellanic Cloud Cepheids and the effect of metallicity

(abridged) The aim of this paper is to investigate the metallicity dependence of the PL relation in V and K, based on a sample of 128 Galactic, 36 LMC, and 6 SMC Cepheids with individual Baade-Wesselink (BW) distances and individually determined metallicities from high-resolution spectroscopy. The p-relation finally adopted is 1.50 -0.24log P. The slope of this relation is based on the condition that the distance to the LMC does not depend on period or (V-K) colour and that the slope of the PL relation based on the BW distances agrees with that based on apparent magnitude. The zero point of the relation is tight to the Cepheids with HST and revised Hipparcos parallaxes as well as to Cepheids in clusters. The slope of the Galactic and LMC K-band relation formally agrees within the errors, and combining all Cepheids (including the SMC) results in a negligible metallicity dependence. A similar conclusion is found for the reddening-free Wesenheit relation. In the V-band the situation is more complex. The slope of the LMC and the Galactic PL relation differ at the 3sigma level. Combining the sample nevertheless results in a metallicity term significant at the 2sigma level. The details of the comparison of BW-based distances and Cepheids with HST and revised Hipparcos parallaxes also play a role. The method used by Storm et al. would lead to larger DM of 18.37 and 18.81 for the LMC and SMC, respectively. The LMC DM is shorter than the currently accepted value, which is in the range 18.42 to 18.55 (Walker 2012), and it is speculated that the p-factor may depend on metallicity.

preprint2012arXiv

Carbon-rich dust production in metal-poor galaxies in the Local Group

We have observed a sample of 19 carbon stars in the Sculptor, Carina, Fornax, and Leo I dwarf spheroidal galaxies with the Infrared Spectrograph on the Spitzer Space Telescope. The spectra show significant quantities of dust around the carbon stars in Sculptor, Fornax, and Leo I, but little in Carina. Previous comparisons of carbon stars with similar pulsation properties in the Galaxy and the Magellanic Clouds revealed no evidence that metallicity affected the production of dust by carbon stars. However, the more metal-poor stars in the current sample appear to be generating less dust. These data extend two known trends to lower metallicities. In more metal-poor samples, the SiC dust emission weakens, while the acetylene absorption strengthens. The bolometric magnitudes and infrared spectral properties of the carbon stars in Fornax are consistent with metallicities more similar to carbon stars in the Magellanic Clouds than in the other dwarf spheroidals in our sample. A study of the carbon budget in these stars reinforces previous considerations that the dredge-up of sufficient quantities of carbon from the stellar cores may trigger the final superwind phase, ending a star&#39;s lifetime on the asymptotic giant branch.

preprint2012arXiv

Discovery of multiple dust shells beyond 1 arcmin in the circumstellar envelope of IRC+10216 using Herschel/PACS

We present new Herschel/PACS images at 70, 100, and 160 micron of the well-known, nearby, carbon-rich asymptotic giant branch star IRC+10216 revealing multiple dust shells in its circumstellar envelope. For the first time, dust shells (or arcs) are detected until 320 arcsec. The almost spherical shells are non-concentric and have an angular extent between 40 deg and 200 deg. The shells have a typical width of 5 arcsec - 8 arcsec, and the shell separation varies in the range of 10 arcsec - 35 arcsec, corresponding to 500-1700 yr. Local density variations within one arc are visible. The shell/intershell density contrast is typically 4, and the arcs contain some 50% more dust mass than the smooth envelope. The observed (nested) arcs record the mass-loss history over the past 16 000 yr, but Rayleigh-Taylor and Kelvin-Helmholtz instabilities in the turbulent astropause and astrosheath will erase any signature of the mass-loss history for at least the first 200 000 yr of mass loss. Accounting for the bowshock structure, the envelope mass around IRC+10216 contains >2Msun of gas and dust mass. It is argued that the origin of the shells is related to non-isotropic mass-loss events and clumpy dust formation.

preprint2012arXiv

Infrared excess around nearby RGB stars and Reimers law

(Abridged) The spectral energy distributions of a well-defined sample of 54 RGB stars are constructed, and fitted with the dust radiative transfer model DUSTY. The central stars are modeled by MARCS model atmospheres. In a first step, the best-fit MARCS model is derived, determining the effective temperature. In a second step, models with a finite dust optical depth are fitted and it is determined whether the reduction in chi2 in such models with one additional free parameter is statistically significant. 23 stars are found to have a significant infrared excess, which is interpreted as mass loss. The dust optical depths are translated into mass-loss rates assuming a typical expansion velocity of 10 km/s and a dust-to-gas ratio of 0.005. The mass-loss rates are compared to those derived for luminous stars in globular clusters, by fitting both the infrared excess, as in the present paper, and the chromospheric lines. There is excellent agreement between these values and the mass-loss rates derived from the chromospheric activity. There is a systematic difference with the literature mass-loss rates derived from modeling the infrared excess, and this has been traced to technical details on how the DUSTY radiative transfer model is run. If the present results are combined with those from modeling the chromospheric emission lines, we obtain the fits Log Mdot = (1.0 +- 0.3) Log L + (-12.0 +- 0.9) and Log Mdot = (0.6 +- 0.2) Log (LR/M) + (-11.9 +- 0.9). The predictions of these mass-loss rate formula are tested against the RGB mass loss determination in NGC 6791. Using a scaling factor of (8 +- 5), both relations can fit this value. That the scaling factor is larger than unity suggests that the expansion velocity and/or dust-to-gas ratio, or even the dust opacities, are different from the values adopted.

preprint2012arXiv

Large-scale environments of binary AGB stars probed by Herschel - I. Morphology statistics and case studies of R Aquarii and W Aquilae

The Mass loss of Evolved StarS (MESS) sample offers a selection of 78 Asymptotic Giant Branch (AGB) stars and Red Supergiants (RSGs) observed with the PACS photometer on-board Herschel at 70 and 160 μm. For most of these objects, the dusty AGB wind differs from spherical symmetry and the wind shape can be subdivided into four classes. In the present paper we concentrate on the influence of a companion on the morphology of the stellar wind. Literature was searched to find binaries in the MESS sample and these are subsequently linked to their wind-morphology class to assert that the binaries are not distributed equally among the classes. In the second part of the paper we concentrate on the circumstellar environment of the two prominent objects R Aqr and W Aql. Each shows a characteristic signature of a companion interaction with the stellar wind. For the symbiotic star R Aqr, PACS revealed two perfectly opposing arms which in part reflect the previously observed ring-shaped nebula in the optical. However, from the far-IR there is evidence that the emitting region is elliptical rather than circular. The outline of the wind of W Aql seems to follow a large Archimedean spiral formed by the orbit of the companion but also shows strong indications of an interaction with the interstellar medium. The nature of the companion of W Aql was investigated and the magnitude of the orbital period supports the size of the spiral outline.

preprint2012arXiv

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 SMC16, and SMP LMC56. 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 (PAHs). The strength of the fullerene bands in the three planetary nebulae is very similar, while the ratio of the [NeIII]15.5um/[NeII]12.8um fine structure lines, an indicator of the strength of the radiation field, is 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-9um 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.

preprint2012arXiv

On the physical structure of IRC+10216. Ground-based and Herschel observations of CO and CCH

The C-rich AGB star IRC+10216 undergoes strong mass loss, and quasi-periodic density enhancements in the circumstellar matter have been reported. CO is ubiquitous in the CSE, while CCH emission comes from a spatially confined shell. With the IRAM 30m telescope and Herschel/HIFI, we recently detected unexpectedly strong emission from the CCH N=4-3, 6-5, 7-6, 8-7, and 9-8 transitions, challenging the available chemical and physical models. We aim to constrain the physical properties of IRC+10216&#39;s CSE, including the effect of episodic mass loss on the observed emission. In particular, we aim to determine the excitation region and conditions of CCH and to reconcile these with interferometric maps of the N=1-0 transition. Via radiative-transfer modelling, we provide a physical description of the CSE, constrained by the SED and a sample of 20 high-resolution and 29 low-resolution CO lines. We further present detailed radiative-transfer analysis of CCH. Assuming a distance of 150pc, the SED is modelled with a stellar luminosity of 11300Lsun and a dust-mass-loss rate of 4.0\times10^{-8}Msun/yr. Based on the analysis of 20 high resolution CO observations, an average gas-mass-loss rate for the last 1000yrs of 1.5\times10^{-5}Msun/yr is derived. This gives a gas-to-dust-mass ratio of 375, typical for an AGB star. The gas kinetic temperature throughout the CSE is described by 3 powerlaws: it goes as r^{-0.58} for r<9R*, as r^{-0.40} for 9<=r<=65R*, and as r^{-1.20} for r>65R*. This model successfully describes all 49 CO lines. We show the effect of wind-density enhancements on the CCH-abundance profile, and the good agreement of the model with the CCH N=1-0 transition and with the lines observed with the 30m telescope and HIFI. We report on the importance of radiative pumping to the vibrationally excited levels of CCH and the significant effect this has on the excitation of all levels of the CCH-molecule.

preprint2012arXiv

Oxygen-rich dust production in IC 10

We report the detection of oxygen-rich circumstellar envelopes in stars of the nearby (700 kpc) starburst galaxy IC 10. The star formation history and the chemical environment of this galaxy makes it an ideal target to observe dust production by high-mass stars in a low-metallicity environment. The goal of this study is to identify oxygen-rich stars in IC 10 and to constrain their nature between asymptotic giant branch stars (AGBs), red supergiants (RSGs), and other infrared bright sources. We examine the mass-loss rate of the stars and compare to results obtained for the Magellanic Clouds. Our objectives are to (1) assess whether RSGs can be significant dust producers in IC 10, and (2), solve the discrepancy between the star formation history of IC 10 and the relatively low number of RSGs detected in the optical. We search for silicate dust in emission by using the spectral map observed with the Infrared Spectrograph on board the Spitzer Space Telescope. The optical (UBVRI) and infrared (JHK, Spitzer/IRAC and Spitzer/MIPS) photometry is used to assert the membership of the stars to IC 10 and disentangle between AGBs and RSGs. Radiative models are used to infer mass-loss rates and stellar luminosities. The luminosity and colors of at least 9 silicate emission sources are consistent with stars within IC 10. Furthermore, the photometry of 2 of these sources is consistent with RSGs. We derive dust mass-loss rates similar to the values found in the Magellanic Clouds. Accounting for the sample completeness, RSGs are not important contributors to the dust mass budget in IC 10.

preprint2012arXiv

The enigmatic nature of the circumstellar envelope and bow shock surrounding Betelgeuse as revealed by Herschel. I. Evidence of clumps, multiple arcs, and a linear bar-like structure

Context. The interaction between stellar winds and the interstellar medium (ISM) can create complex bow shocks. The photometers on board the Herschel Space Observatory are ideally suited to studying the morphologies of these bow shocks. Aims. We aim to study the circumstellar environment and wind-ISM interaction of the nearest red supergiant, Betelgeuse. Methods. Herschel PACS images at 70, 100, and 160 micron and SPIRE images at 250, 350, and 500 micron were obtained by scanning the region around Betelgeuse. These data were complemented with ultraviolet GALEX data, near-infrared WISE data, and radio 21 cm GALFA-HI data. The observational properties of the bow shock structure were deduced from the data and compared with hydrodynamical simulations. Results. The infrared Herschel images of the environment around Betelgeuse are spectacular, showing the occurrence of multiple arcs at 6-7 arcmin from the central target and the presence of a linear bar at 9 arcmin. Remarkably, no large-scale instabilities are seen in the outer arcs and linear bar. The dust temperature in the outer arcs varies between 40 and 140 K, with the linear bar having the same colour temperature as the arcs. The inner envelope shows clear evidence of a non-homogeneous clumpy structure (beyond 15 arcsec), probably related to the giant convection cells of the outer atmosphere. The non-homogeneous distribution of the material even persists until the collision with the ISM. A strong variation in brightness of the inner clumps at a radius of 2 arcmin suggests a drastic change in mean gas and dust density some 32 000 yr ago. Using hydrodynamical simulations, we try to explain the observed morphology of the bow shock around Betelgeuse. Conclusions: [abbreviated]

preprint2012arXiv

The Spitzer Spectroscopic Survey of S-type Stars

S-type AGB stars are thought to be in the transitional phase between M-type and C-type AGB stars. Because of their peculiar chemical composition, one may expect a strong influence of the stellar C/O ratio on the molecular chemistry and the mineralogy of the circumstellar dust. In this paper, we present a large sample of 87 intrinsic galactic S-type AGB stars, observed at infrared wavelengths with the Spitzer Space Telescope, and supplemented with ground-based optical data. On the one hand, we derive the stellar parameters from the optical spectroscopy and photometry, using a grid of model atmospheres. On the other, we decompose the infrared spectra to quantify the flux-contributions from the different dust species. Finally, we compare the independently determined stellar parameters and dust properties. For the stars without significant dust emission, we detect a strict relation between the presence of SiS absorption in the Spitzer spectra and the C/O ratio of the stellar atmosphere. These absorption bands can thus be used as an additional diagnostic for the C/O ratio. For stars with significant dust emission, we define three groups, based on the relative contribution of certain dust species to the infrared flux. We find a strong link between group-membership and C/O ratio. We show that these groups can be explained by assuming that the dust-condensation can be cut short before silicates are produced, while the remaining free atoms and molecules can then form the observed magnesium sulfides or the carriers of the unidentified 13 and 20 micron features. Finally, we present the detection of emission features attributed to molecules and dust characteristic to C-type stars, such as molecular SiS, hydrocarbons and magnesium sulfide grains. We show that we often detect magnesium sulfides together with molecular SiS and we propose that it is formed by a reaction of SiS molecules with Mg.

preprint2012arXiv

The VMC Survey - VI. Quasars behind the Magellanic system

The number and spatial distribution of confirmed quasi-stellar objects (QSOs) behind the Magellanic system is limited. This undermines their use as astrometric reference objects for different types of studies. We have searched for criteria to identify candidate QSOs using observations from the VISTA survey of the Magellanic Clouds system (VMC) that provides photometry in the YJKs bands and 12 epochs in the Ks band. The (Y-J) versus (J-Ks) diagram has been used to distinguish QSO candidates from Milky Way stars and stars of the Magellanic Clouds. Then, the slope of variation in the Ks band has been used to identify a sample of high confidence candidates. These criteria were developed based on the properties of 117 known QSOs presently observed by the VMC survey. VMC YJKs magnitudes and Ks light-curves of known QSOs behind the Magellanic system are presented. About 75% of them show a slope of variation in Ks>10^-4 mag/day and the shape of the light-curve is in general irregular and without any clear periodicity. The number of QSO candidates found in tiles including the South Ecliptic Pole and the 30 Doradus regions is 22 and 26, respectively, with a ~20% contamination by young stellar objects, planetary nebulae, stars and normal galaxies. By extrapolating the number of QSO candidates to the entire VMC survey area we expect to find about 1200 QSOs behind the LMC, 400 behind the SMC, 200 behind the Bridge and 30 behind the Stream areas, but not all will be suitable for astrometry. Further, the Ks band light-curves can help support investigations of the mechanism responsible for the variations.

preprint2012arXiv

The VMC Survey. V. First results for Classical Cepheids

The VISTA Magellanic Cloud (VMC, PI M.R. Cioni) survey is collecting deep Ks-band time-series photometry of the pulsating variable stars hosted by the system formed by the two Magellanic Clouds (MCs) and the &#34;bridge&#34; connecting them. In this paper we present the first results for Classical Cepheids, from the VMC observations of two fields in the Large Magellanic Cloud (LMC). The VMC Ks-band light curves of the Cepheids are well sampled (12-epochs) and of excellent precision. We were able to measure for the first time the Ks magnitude of the faintest Classical Cepheids in the LMC (Ks\sim17.5 mag), which are mostly pulsating in the First Overtone (FO) mode, and to obtain FO Period-Luminosity (PL), Period-Wesenheit (PW), and Period-Luminosity-Color (PLC) relations, spanning the full period range from 0.25 to 6 day. Since the longest period Cepheid in our dataset has a variability period of 23 day, we have complemented our sample with literature data for brighter F Cepheids. On this basis we have built a PL relation in the Ks band that, for the first time, includes short period pulsators, and spans the whole range from 1.6 to 100 days in period. We also provide the first ever empirical PW and PLC relations using the (V-Ks) color and time-series Ks photometry. The very small dispersion (\sim0.07 mag) of these relations makes them very well suited to study the three-dimensional (3D) geometry of the Magellanic system. The use of &#34;direct&#34; (parallax- and Baade-Wesselink- based) distance measurements to both Galactic and LMC Cepheids, allowed us to calibrate the zero points of the PL, PW, and PLC relations obtained in this paper, and in turn to estimate an absolute distance modulus of (m-M)0=18.46\pm0.03 for the LMC. This result is in agreement with most of the latest literature determinations based on Classical Cepheids.

preprint2012arXiv

Unexpectedly large mass loss during the thermal pulse cycle of the red giant 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 undergoes 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 on 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 arcsec. What was hitherto thought to be only a thin, spherical shell with a clumpy structure, is revealed to contain a spiral structure. Spiral structures associated with circumstellar envelopes have been seen previously, from which it was concluded that the systems must be binaries. Using the data, combined with hydrodynamic simulations, we conclude that R Sculptoris is a binary system that underwent a thermal pulse approximately 1800 years ago, lasting approximately 200 years. About 0.003 Msun of mass was ejected at a velocity of 14.3 km s-1 and at a rate approximately 30 times higher than the prepulse mass-loss rate. This shows that approximately 3 times more mass is returned to the interstellar medium during and immediately after a pulse than previously thought.

preprint2011arXiv

A far-infrared survey of bow shocks and detached shells around AGB stars and red supergiants

Far-infrared Herschel/PACS images at 70 and 160 micron of a sample of 78 Galactic evolved stars are used to study the (dust) emission structures, originating from stellar wind-ISM interaction. In addition, two-fluid hydrodynamical simulations of the coupled gas and dust in wind-ISM interactions are used to compare with the observations. Four distinct classes of wind-ISM interaction (i.e. &#34;fermata&#34;, &#34;eyes&#34;, &#34;irregular&#34;, and &#34;rings&#34;) are identified and basic parameters affecting the morphology are discussed. We detect bow shocks for ~40% of the sample and detached rings for ~20%. De-projected stand-off distances (R_0) -- defined as the distance between the central star and the nearest point of the interaction region -- of the detected bow shocks (&#34;fermata&#34; and &#34;eyes&#34;) are derived from the PACS images and compared to previous results, model predictions and the simulations. All observed bow shocks have stand-off distances smaller than 1 pc. Observed and theoretical stand-off distances are used together to independently derive the local ISM density. Both theoretical (analytical) models and hydrodynamical simulations give stand-off distances for adopted stellar properties that are in good agreement with the measured de-projected stand-off distance of wind-ISM bow shocks. The possible detection of a bow shock -- for the distance limited sample -- appears to be governed by its physical size as set roughly by the stand-off distance. In particular the star&#39;s peculiar space velocity and the density of the ISM appear decisive in detecting emission from bow shocks or detached rings. Tentatively, the &#34;eyes&#34; class objects are associated to (visual) binaries, while the &#34;rings&#34; generally appear not to occur for M-type stars, only for C or S-type objects that have experienced a thermal pulse.

preprint2011arXiv

Ground-Based BVRI Follow-Up Observations of the Cepheid V1154 Cyg in Kepler&#39;s Field

The almost un-interrupted observations from Kepler Space Telescope can be used to search for Earth-size and larger planets around other stars, as well as for stellar variability and asteroseismological study. However, the Kepler&#39;s observations are carried out with a single broad-band filter, and ground-based follow-up observations are needed to complement Kepler&#39;s light curves to fully characterize the properties of the target stars. Here we present ground-based optical (BVRI) followup observations of V1154 Cyg, the only Cepheid in the Kepler field of view, and deriving basic properties of this star.

preprint2011arXiv

Herschel&#39;s view into Mira&#39;s head

Herschel&#39;s PACS instrument observed the environment of the binary system Mira Ceti in the 70 and 160 micron bands. These images reveal bright structures shaped as five broken arcs and fainter filaments in the ejected material of Mira&#39;s primary star. The overall shape of the IR emission around Mira deviates significantly from the expected alignment with Mira&#39;s exceptionally high space velocity. The observed broken arcs are neither connected to each other nor are they of a circular shape; they stretch over angular ranges of 80 to 100 degrees. By comparing Herschel and GALEX data, we found evidence for the disruption of the IR arcs by the fast outflow visible in both Halpha and the far UV. Radial intensity profiles are derived, which place the arcs at distances of 6-85&#34; (550 - 8000 AU) from the binary. Mira&#39;s IR environment appears to be shaped by the complex interaction of Mira&#39;s wind with its companion, the bipolar jet, and the ISM.

preprint2011arXiv

Planetary Nebulae in the VISTA Magellanic Cloud (VMC) Survey

The multi-epoch YJKs sub-arcsecond photometry of the VMC survey provides a long anticipated deep near-infrared (NIR) window into further understanding the stellar populations of the Magellanic Clouds. The first year of observations consisted of six tiles covering ~9% of the Large Magellanic Cloud (LMC) survey region and contains 102 objects previously classified as planetary nebulae (PNe). A large proportion of the sample were found to be contaminated by non-PNe. These initial results underline the importance of establishing a clean catalogue of LMC PNe before they are applied in areas such as the planetary nebula luminosity function (PNLF) and searches for binary central stars. As the VMC survey progresses it will play a fundamental role in cleaning extant PN catalogues and a complementary role in the discovery of new PNe.

preprint2011arXiv

The VMC Survey - I. Strategy and First Data

The new VISual and Infrared Telescope for Astronomy (VISTA) has started operations. Over its first five years it will be collecting data for six public surveys, one of these is the near-infrared YJKsVISTA survey of the Magellanic Clouds system (VMC). This survey comprises the LMC, the SMC, the Bridge connecting the two galaxies and two fields in the Stream. This paper provides an overview of the VMC survey strategy and presents first science results. The main goals of the VMC survey are the determination of the spatially resolved SFH and 3D structure of the Magellanic system. Therefore, the VMC survey is designed to reach stars as faint as the oldest main sequence turn-off point and to constrain the mean magnitude of pulsating variable stars such as RR Lyrae stars and Cepheids. This paper focuses on observations of VMC fields in the LMC obtained between November 2009 and March 2010. These observations correspond to a 7% completeness of the LMC fields. The VMC observations consist of multi-epoch measurements organised following a specific structure. The data were reduced using the VISTA Data Flow System pipeline whose source catalogues were produced and made available via the VISTA Science Archive. The analysis of the data shows that the sensitivity in each wave band agrees with expectations. Uncertainties and completeness of the data are also derived. The first science results, aimed at assessing the science quality of the VMC data, include an overview of the distribution of stars in colour-magnitude and colour-colour diagrams, the detection of planetary nebulae and of stellar clusters, and the Ks band light-curve of variable stars. The VMC survey represents a tremendous improvement, in spatial resolution and sensitivity, on previous panoramic observations of the Magellanic system in the near-infrared and complements nicely the deep observations at other wavelengths. (Abridged)

preprint2011arXiv

The VMC survey III. Mass-loss rates and luminosities of LMC AGB stars

Dust radiative transfer models are presented for all 374 AGB stars candidates in one of the fields observed by the new VISTA survey of the Magellanic Clouds (VMC). Mass-loss rates, luminosities and a classification of C- and O-rich stars are derived by fitting the models to the spectral energy distribution obtained by combining VMC data with existing optical, near-, and mid-infrared photometry. This exploratory study shows that our method provides reliable mass-loss rates, luminosities and chemical classifications for all AGB stars. These results offer already important constraints to AGB evolutionary models. Most of our conclusions, especially for the rarer dust-enshrouded extreme AGB stars, are however strongly limited by the relatively small area covered by our study. Forthcoming VMC observations will easily remove this limitation. [abridged]

preprint2011arXiv

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

We derive the star formation history for several regions of the LMC, using deep near-infrared data from the VISTA near-infrared YJKs survey of the Magellanic system (VMC). The regions include three almost-complete 1.4 sqdeg tiles located 3.5 deg away from the LMC centre in distinct directions. To this dataset, we add two 0.036 sqdeg subregions inside the 30 Doradus tile. The SFH is derived from the simultaneous reconstruction of two different CMDs, using the minimization code StarFISH. 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), 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-resolved SFH of the LMC with unprecedented quality compared to previous wide-area surveys. In particular, the data clearly 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 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 seems remarkably constant over the time interval from 8.4 to 9.7. The AMRs, instead, turn out to be remarkably similar across the LMC. The fields studied so far are fit extremely well by a single disk of inclination 26.2+-2.0 deg, position angle of the line of nodes 129.1+-13.0 deg, and distance modulus of 18.470+-0.006 mag (random errors only) up to the LMC centre.

preprint2011arXiv

The VMC Survey. II. A multi-wavelength study of LMC planetary nebulae and their mimics

The VISTA Magellanic Cloud (VMC) survey is assembling a deep, multi-epoch atlas of YJKs photometry across the Magellanic Clouds. Prior to the VMC survey only the brightest Magellanic Cloud PNe (MCPNe) were accessible at near-infrared (NIR) wavelengths. It is now possible for the first time to assemble the NIR properties of MCPNe and to identify contaminating non-PNe mimics which are best revealed at NIR wavelengths (e.g. HII regions and symbiotic stars). To maintain the unique scientific niche that MCPNe occupy these contaminants must be removed. Here we conduct a VMC-led, multi-wavelength study of 102 objects previously classified as PNe that are located within the first six VMC tiles observed. We present images, photometry, lightcurves, diagnostic colour-colour diagrams and spectral energy distributions used to analyse the entire sample. At least five PNe have newly resolved nebula morphologies, a task previously only possible with the HST. A total 45/67 (67%) of Reid & Parker (RP) catalogued objects were reclassified as non-PNe, most of which were located in the vicinity of 30 Doradus. This sample included 16 field stars, 5 emission line stars, 19 HII regions, 4 symbiotic star candidates and 1 young stellar object. We discuss possible selection effects responsible for their inclusion in the RP catalogue and the implications for binary central star surveys targeting LMC PNe. A total of five new LMC symbiotic star candidates identified, compared to eight previously known, underlines the important role the VMC survey will have in advancing Magellanic symbiotic star studies.

preprint2011arXiv

Warm water vapour in the sooty outflow from a luminous carbon star

In 2001, the discovery of circumstellar water vapour around the ageing carbon star IRC+10216 was announced. This detection challenged the current understanding of chemistry in old stars, since water vapour was predicted to be absent in carbon-rich stars. Several explanations for the occurrence of water vapour were postulated, including the vaporization of icy bodies (comets or dwarf planets) in orbit around the star, grain surface reactions, and photochemistry in the outer circumstellar envelope. However, the only water line detected so far from one carbon-rich evolved star can not discriminate, by itself, between the different mechanisms proposed. Here we report on the detection by the Herschel satellite of dozens of water vapour lines in the far-infrared and sub-millimetre spectrum of IRC+10216, including some high-excitation lines with energies corresponding to ~1000 K. The emission of these high-excitation water lines can only be explained if water vapour is present in the warm inner region of the envelope. A plausible explanation for the formation of warm water vapour appears to be the penetration of ultraviolet (UV) photons deep into a clumpy circumstellar envelope. This mechanism triggers also the formation of other molecules such as ammonia, whose observed abundances are much higher than hitherto predicted.

preprint2011arXiv

X Her and TX Psc: Two cases of ISM interaction with stellar winds observed by Herschel

The asymptotic giant branch (AGB) stars X Her and TX Psc have been imaged at 70 and 160 microns with the PACS instrument onboard the Herschel satellite, as part of the large MESS (Mass loss of Evolved StarS) Guaranteed Time Key Program. The images reveal an axisymmetric extended structure with its axis oriented along the space motion of the stars. This extended structure is very likely to be shaped by the interaction of the wind ejected by the AGB star with the surrounding interstellar medium (ISM). As predicted by numerical simulations, the detailed structure of the wind-ISM interface depends upon the relative velocity between star+wind and the ISM, which is large for these two stars (108 and 55 km/s for X Her and TX Psc, respectively). In both cases, there is a compact blob upstream whose origin is not fully elucidated, but that could be the signature of some instability in the wind-ISM shock. Deconvolved images of X Her and TX Psc reveal several discrete structures along the outermost filaments, which could be Kelvin-Helmholtz vortices. Finally, TX Psc is surrounded by an almost circular ring (the signature of the termination shock?) that contrasts with the outer, more structured filaments. A similar inner circular structure seems to be present in X Her as well, albeit less clearly.

preprint2010arXiv

870 um observations of evolved stars with LABOCA

During their evolution, Asymptotic Giant Branch (AGB) stars experience a high mass-loss which leads to the formation of a Circumstellar Envelope (CSE) of dust and gas. The mass-loss process is the most important phenomenon during this evolutionary stage. In order to understand it, it is important to study the physical parameters of the CSE. The emission of the CSE in the (sub)millimetre range is dominated by the dust continuum. This means that (sub)millimetre observations are a key tool in tracing the dust and improving our knowledge of the mass-loss process. In this paper we analyse new sub-millimetre observations of 9 evolved stars in order to constrain the CSE physical parameters. The data were taken by the new APEX bolometer LABoCa. The fluxes at 870 um are derived and the extended emission is investigated. We compute the spectral energy distribution (SEDs) using a 1D radiative transfer code, DUSTY which we compared to literature data. Grain properties are calculated using both spherical grains distribution and a Continuous Distribution of Ellipsoids (CDE) and a comparison between the two is drawn. Synthetic surface brightness maps have been derived from the modelling and were compared to the LABoCa brightness maps. We detected the presence of extended emission around four stars. Physical parameters of the circumstellar envelope are derived from SED modelling, such as the dust chemical composition, the dust condensation temperature and the total mass of the envelope. It proves difficult however to fit the SED and the intensity profile simultaneously.

preprint2010arXiv

Detection of anhydrous hydrochloric acid, HCl, in IRC+10216 with the Herschel SPIRE and PACS spectrometers

We report on the detection of anhydrous hydrochloric acid (hydrogen chlorine, HCl) in the carbon-rich star IRC+10216 using the spectroscopic facilities onboard the Herschel satellite. Lines from J=1-0 up to J=7-6 have been detected. From the observed intensities, we conclude that HCl is produced in the innermost layers of the circumstellar envelope with an abundance relative to H2 of 5x10^-8 and extends until the molecules reach its photodissociation zone. Upper limits to the column densities of AlH, MgH, CaH, CuH, KH, NaH, FeH, and other diatomic hydrides have also been obtained.

preprint2010arXiv

Herschel images of NGC 6720: H2 formation on dust grains

Herschel PACS and SPIRE images have been obtained of NGC 6720 (the Ring Nebula). This is an evolved planetary nebula with a central star that is currently on the cooling track, due to which the outer parts of the nebula are recombining. From the PACS and SPIRE images we conclude that there is a striking resemblance between the dust distribution and the H2 emission, which appears to be observational evidence that H2 forms on grain surfaces. We have developed a photoionization model of the nebula with the Cloudy code which we used to determine the physical conditions of the dust and investigate possible formation scenarios for the H2. We conclude that the most plausible scenario is that the H2 resides in high density knots which were formed after the recombination of the gas started when the central star entered the cooling track. Hydrodynamical instabilities due to the unusually low temperature of the recombining gas are proposed as a mechanism for forming the knots. H2 formation in the knots is expected to be substantial after the central star underwent a strong drop in luminosity about one to two thousand years ago, and may still be ongoing at this moment, depending on the density of the knots and the properties of the grains in the knots.

preprint2010arXiv

Herschel PACS and SPIRE imaging of CW Leo

Herschel PACS and SPIRE images have been obtained over a 30&#39;x30&#39; area around the well-known carbon star CW Leo (IRC +10 216). An extended structure is found in an incomplete arc of ~22&#39; diameter, which is cospatial with the termination shock due to interaction with the interstellar medium (ISM) as defined by Sahai & Chronopoulos from ultraviolet GALEX images. Fluxes are derived in the 70, 160, 250, 350, and 550 um bands in the region where the interaction with the ISM takes place, and this can be fitted with a modified black body with a temperature of 25+-3 K. Using the published proper motion and radial velocity for the star, we derive a heliocentric space motion of 25.1 km/s. Using the PACS and SPIRE data and the analytical formula of the bow shock structure, we infer a de-projected standoff distance of the bow shock of R0 = (8.0+-0.3)x10^17 cm. We also derive a relative velocity of the star with respect to the ISM of (106.6+-8.7)/sqrt(n_ISM) km/s, where n_ISM is the number density of the local ISM.

preprint2010arXiv

Herschel-SPIRE FTS spectroscopy of the carbon-rich objects AFGL 2688, AFGL 618 and NGC 7027

We present far-infrared and submillimetre spectra of three carbon-rich evolved objects, AFGL 2688, AFGL 618 and NGC 7027. The spectra were obtained with the SPIRE Fourier transform spectrometer on board the Herschel Space Observatory, and cover wavelengths from 195-670 um, a region of the electromagnetic spectrum hitherto difficult to study in detail. The far infrared spectra of these objects are rich and complex, and we measure over 150 lines in each object. Lines due to 18 different species are detected. We determine physical conditions from observations of the rotational lines of several molecules, and present initial large velocity gradient models for AFGL 618. We detect water in AFGL 2688 for the first time, and confirm its presence in AFGL 618 in both ortho and para forms. In addition, we report the detection of the J=1-0 line of CH+ in NGC 7027.

preprint2010arXiv

PACS and SPIRE Spectroscopy of the Red Supergiant VY CMa

With a luminosity > 10^5 Lsun and a mass-loss rate of about 2.10-4 Msun/yr, the red supergiant VY CMa truly is a spectacular object. Because of its extreme evolutionary state, it could explode as supernova any time. Studying its circumstellar material, into which the supernova blast will run, provides interesting constraints on supernova explosions and on the rich chemistry taking place in such complex circumstellar envelopes. We have obtained spectroscopy of VYCMa over the full wavelength range offered by the PACS and SPIRE instruments of Herschel, i.e. 55 to 672 micron. The observations show the spectral fingerprints of more than 900 spectral lines, of which more than half belong to water. In total, we have identified 13 different molecules and some of their isotopologues. A first analysis shows that water is abundantly present, with an ortho-to-para ratio as low as 1.3:1, and that chemical non-equilibrium processes determine the abundance fractions in the inner envelope.

preprint2010arXiv

Silicon in the dust formation zone of IRC +10216 as observed with PACS and SPIRE on board Herschel

The interstellar medium is enriched primarily by matter ejected from evolved low and intermediate mass stars. The outflows from these stars create a circumstellar envelope in which a rich gas-phase and dust-nucleation chemistry takes place. We observed the nearest carbon-rich evolved star, IRC+10216, using the PACS (55-210 μm) and SPIRE (194-672 μm) spectrometers on board Herschel. We find several tens of lines from SiS and SiO, including lines from the v=1 vibrational level. For SiS these transitions range up to J=124-123, corresponding to energies around 6700K, while the highest detectable transition is J=90-89 for SiO, which corresponds to an energy around 8400K. Both species trace the dust formation zone of IRC+10216, and the broad energy ranges involved in their detected transitions permit us to derive the physical properties of the gas and the particular zone in which each species has been formed. This allows us to check the accuracy of chemical thermodynamical equilibrium models and the suggested depletion of SiS and SiO due to accretion onto dust grains.

preprint2010arXiv

The detached dust shells of AQ And, U Ant, and TT Cyg

Detached circumstellar dust shells are detected around three carbon variables using Herschel-PACS. Two of them are already known on the basis of their thermal CO emission and two are visible as extensions in IRAS imaging data. By model fits to the new data sets, physical sizes, expansion timescales, dust temperatures, and more are deduced. A comparison with existing molecular CO material shows a high degree of correlation for TT Cyg and U Ant but a few distinct differences with other observables are also found.

preprint2009arXiv

Luminosities and mass-loss rates of SMC and LMC AGB stars and Red Supergiants

(Abridged) Dust radiative transfer models are presented for 101 carbon stars and 86 oxygen-rich evolved stars in the Magellanic Clouds for which 5-35 \mum\ {\it Spitzer} IRS spectra are available. The spectra are complemented with available optical and infrared photometry to construct the spectral energy distribution. A minimisation procedure is used to fit luminosity, mass-loss rate and dust temperature at the inner radius. Different effective temperatures and dust content are also considered. Periods from the literature and from new OGLE-III data are compiled and derived. The O-rich stars are classified in foreground objects, AGB stars and Red Super Giants. For the O-rich stars silicates based on laboratory optical constants are compared to &#34;astronomical silicates&#34;. Overall, the grain type by Volk & Kwok (1988) fit the data best. However, the fit based on laboratory optical constants for the grains can be improved by abandoning the small-particle limit. The influence of grain size, core-mantle grains and porosity are explored. Relations between mass-loss rates and luminosity and pulsation period are presented and compared to the predictions of evolutionary models, those by Vassiliadis & Wood (1993) and their adopted mass-loss recipe, and those based on a Reimers mass-loss law with a scaling of a factor of five. The Vassiliadis & Wood models describe the data better, although there are also some deficiencies, in particular to the maximum adopted mass-loss rate. The OGLE-III data reveal an O-rich star in the SMC with a period of 1749 days. Its absolute magnitude of $M_{\rm bol}= -8.0$ makes it a good candidate for a super-AGB star.

preprint2009arXiv

Near-infrared spectroscopy of AGB star candidates in Fornax, Sculptor and NGC 6822

Context: The Asymptotic Giant Branch (AGB) phase is characterised by substantial mass loss that is accompanied by the formation of dust. In extreme cases this will make the star no longer visible in the optical. For a better understanding of AGB evolution it is important to identify and characterise these very red AGB stars. Aims: The first aim of this article is to improve the census of red AGB stars in three Local Group galaxies, based on near-IR spectroscopic observations of new candidates with red 2MASS (J-K) colours. The opportunity is taken to compare the near-IR spectra with those of Milky Way stars. Methods: We used ISAAC on the ESO VLT to take J and H-band spectra of 36 targets in Fornax, Sculptor and NGC 6822. Results: Twelve new C-stars are found in Fornax, and one is confirmed in Sculptor. All C-stars have (J-K) > 1.6, and are brighter than -3.55 in bolometric magnitude. Ten new oxygen-rich late-type giant stars are identified in Fornax, but none is extremely red or very luminous. Five luminous O-rich AGB stars are identified in NGC 6822, of which 3 show water absorption, indicative of spectral type M. Again, none is as red as Milky Way OH/IR stars, but in this galaxy the list of candidate AGB stars is biased against very red objects. In some C-stars with (J-K)>2 an extremely strong 1.53 $μ$m absorption band is found. These stars are probably all Mira variables and the feature is related to the low temperature, high density chemistry that is a first step towards dust formation and mass loss.

preprint2009arXiv

The Dependency of the Cepheid Period-Luminosity Relation on Chemical Composition

The dependency of the Cepheid Period-Luminosity Relation on chemical composition at different wavelengths is assessed via direct detailed abundance analysis of Galactic and Magellanic Cepheids, as derived from high resolution, high signal-to-noise spectra. Our measurements span one order of magnitude in iron content and allow to rule out at the ~ 9 sigma level the universality of the Period-Luminosity Relation in the V band, with metal rich stars being fainter than metal poor ones by ~0.3 mag. The dependency is less pronounced in the K band. Its magnitude and statistical significance decisively depend on detailed distance measurements to individual stars, as inferred via the Infrared Surface Brightness Method.

preprint2009arXiv

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&#39; 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.

preprint2008arXiv

The Magellanic Clouds as a template for the study of stellar populations and galaxy interactions

The Magellanic System represents one of the best places to study the formation and evolution of galaxies. Photometric surveys of various depths, areas and wavelengths have had a significant impact on our understanding of the system; however, a complete picture is still lacking. VMC (the VISTA near-infrared YJKs survey of the Magellanic System) will provide new data to derive the spatially resolved star formation history and to construct a three-dimensional map of the system. These data combined with those from other ongoing and planned surveys will give us an absolutely unique view of the system opening up the doors to truly new science!

preprint2007arXiv

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.

preprint2005arXiv

ESO-VLT and Spitzer spectroscopy of IRAS05328-6827: a massive young stellar object in the Large Magellanic Cloud

We present the first thermal-infrared spectra of an extra-galactic Young Stellar Object (YSO), IRAS05328-6827 in the HII region LHA 120-N148 in the Large Magellanic Cloud. The observed and modelled spectral energy distribution reveals a massive YSO, M=20 Msun, which is heavily-embedded and probably still accreting. The reduced dust content as a consequence of the lower metallicity of the LMC allows a unique view into this object, and together with a high C/O ratio may be responsible for the observed low abundance of water ice and relatively high abundances of methanol and CO2 ices.

preprint2005arXiv

The Chemical Composition of Cepheids in the Milky Way and the Magellanic Clouds

We have measured the elemental abundances of 68 Galactic and Magellanic Cepheids from FEROS and UVES high-resolution and high signal-to-noise spectra in order to establish the influence of the chemical composition on the properties of these stars (see Romaniello et al. 2005). Here we describe the robust analytical procedure we have developed to accurately determine them. The resulting iron abundances span a range between ~ -0.80 dex for stars in the Small Magellanic Cloud and ~ +0.20 dex for the most metal-rich ones in the Galaxy. While the average values for each galaxy are in good agreement with non-pulsating stars of similar age, Cepheids display a significant spread. Thus it is fundamental to measure the metallicity of individual stars.

preprint2001arXiv

ESO Imaging Survey. Deep Public Survey: Multi-Color Optical Data for the Chandra Deep Field South

This paper presents multi-passband optical data obtained from observations of the Chandra Deep Field South (CDF-S), located at alpha ~ 3h 32m, delta ~ -27d 48m. The observations were conducted at the ESO/MPG 2.2m telescope at La Silla using the 8kx8k Wide-Field Imager (WFI). This data set, taken over a period of one year, represents the first field to be completed by the ongoing Deep Public Survey (DPS) being carried out by the ESO Imaging Survey (EIS) project. This paper describes the optical observations, the techniques employed for un-supervised pipeline processing and the general characteristics of the final data set. The paper includes data taken in six different filters U&#39;UBVRI. The data cover an area of about 0.25 square degrees reaching 5 sigma limiting magnitudes of U&#39;_AB=26.0, U_AB=25.7, B_AB=26.4$, V_AB=25.4, R_AB=25.5 and I_AB= 24.7 mag, as measured within a 2xFWHM aperture. The optical data covers the area of ~ 0.1

preprint1999arXiv

The absolute magnitudes of RR Lyrae stars from Hipparcos parallaxes

Using the method of ``reduced parallaxes&#39;&#39; for the Halo RR Lyrae stars in the Hipparcos catalogue we derive a zero point of 0.77 $\pm$ 0.26 mag for an assumed slope of 0.18 in the $M_{\rm V}$-[Fe/H] relation. This is 0.28 magnitude brighter than the value Fernley et al. (1998a) derived by employing the method of statistical parallax for the {\it identical} sample and using the same slope. We point out that a similar difference exists between the ``reduced parallaxes&#39;&#39; method and the statistical parallax method for the Cepheids in the Hipparcos catalogue. We also determine the zero point for the $M_{\rm K}$-$\log P_{0}$ relation, and obtain a value of -1.16 $\pm$ 0.27 mag (for a slope of -2.33). The distance moduli to the Hipparcos RR Lyrae stars derived from the two relations agree well. The derived distance scale is in good agreement with the results from the Main Sequence fitting distances of Galactic globular clusters and with the results of theoretical Horizontal Branch models, and implies a distance modulus to the LMC of 18.61 $\pm$ 0.28 mag.

preprint1998arXiv

The absolute magnitude of K0V stars from HIPPARCOS parallaxes

We investigate the properties of K0V stars with Hipparcos parallaxes and spectral types taken from the Michigan Spectral Survey. The sample of 200 objects allows the empirical investigation of the magnitude selection (Malmquist) bias, which appears clearly present. By selecting those objects that are not affected by bias, we find a mean absolute magnitude of Mv~5.7, a downward revision from 5.9 mag. listed in Schmidt-Kaler (1982). Some objects have absolute magnitudes far brighter than Mv~5.7, and it is suggested that these objects (~20% of the total sample) are K0IV stars which may have been mis-classified as a K0V star. The presence of the Malmquist bias in even this high quality sample suggests that no sample can be expected to be bias-free.