Researcher profile

M. Matsuura

M. Matsuura contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
74works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

74 published item(s)

preprint2025arXiv

The Nearby Evolved Stars Survey III: First data release of JCMT CO-line observations

Low- to intermediate-mass ($\sim$0.8$-$8 M$_\odot$) evolved stars contribute significantly to the chemical enrichment of the interstellar medium in the local Universe, making accurate mass-return estimates in their final stages crucial. The Nearby Evolved Stars Survey (NESS) is a large multi-telescope project targeting a volume-limited sample of $\sim$850 stars within 3 kpc in order to derive the dust and gas return rates in the Solar Neighbourhood, and to constrain the physics underlying these processes. We present an initial analysis of the CO-line observations, including detection statistics, carbon isotopic ratios, initial mass-loss rates, and gas-to-dust ratios. We describe a new data reduction pipeline to analyse the available NESS CO data from the JCMT, measuring line parameters and calculating empirical gas mass-loss rates. We present the first release of the available data on 485 sources, one of the largest homogeneous samples of CO data to date. Comparison with a large literature sample finds that high mass-loss rate and especially carbon-rich sources are over-represented in literature, while NESS is probing significantly more sources at low mass-loss rates, detecting 59 sources in CO for the first time and providing useful upper limits. CO line detection rates are 81% for the CO (2--1) line and 75% for CO (3--2). The majority (82%) of detected lines conform to the expected soft parabola shape, while eleven sources show a double wind. Calculated mass-loss rates show power-law relations with both the dust-production rates and expansion velocities up to $\sim 5 \times 10^{-6}$~\msunyr. Median gas-to-dust ratios of 250 and 680 are found for oxygen-rich and carbon-rich sources, respectively. Our analysis of CO observations in this first data release highlights the importance of our volume-limited approach in characterizing the local AGB population as a whole.

preprint2023arXiv

Quantifying the dust in SN 2012aw and iPTF14hls with ORBYTS

Core-collapse supernovae (CCSNe) are potentially capable of producing large quantities of dust, with strong evidence that ejecta dust masses can grow significantly over extended periods of time. Red-blue asymmetries in the broad emission lines of CCSNe can be modelled using the Monte Carlo radiative transfer code DAMOCLES, to determine ejecta dust masses. To facilitate easier use of DAMOCLES, we present a Tkinter graphical user interface (GUI) running DAMOCLES. The GUI was tested by high school students as part of the Original Research By Young Twinkle Students (ORBYTS) programme, who used it to measure the dust masses formed at two epochs in two Type IIP CCSNe: SN 2012aw and iPTF14hls, demonstrating that a wide range of people can contribute significantly to scientific advancement. Bayesian methods were used to quantify uncertainties on our model parameters. From the presence of a red scattering wing in the day 1863 H$α$ profile of SN 2012aw, we were able to constrain the dust composition to large (radius $>0.1 μ$m) silicate grains, with a dust mass of $6.0^{+21.9}_{-3.6}\times10^{-4} M_\odot$. From the day 1158 H$α$ profile of SN 2012aw, we found a dust mass of $3.0^{+14}_{-2.5}\times10^{-4}$ M$_\odot$. For iPTF14hls, we found a day 1170 dust mass of 8.1 $^{+81}_{-7.6}\times10^{-5}$ M$_{\odot}$ for a dust composition consisting of 50% amorphous carbon and 50% astronomical silicate. At 1000 days post explosion, SN 2012aw and iPTF14hls have formed less dust than SN 1987A, suggesting that SN 1987A could form larger dust masses than other Type IIP's.

preprint2022arXiv

Observation of topological $Z_2$ vortex fluctuations in the frustrated Heisenberg magnet NaCrO$_2$

Spin fluctuations in the triangular-lattice Heisenberg antiferromagnet NaCrO$_2$ are investigated by means of quasi-elastic neutron scattering with high energy resolution and wide energy band. Two components with the following features are captured separately. They are pronounced at intermediate temperatures of 20-50 K. One with quite an extended lifetime corresponding to $\sim$$0.001E_{\rm ex}$ ($E_{\rm ex}$ the exchange energy) nearly disappears at low temperature 10 K, and the other with an extended lifetime $\sim$$0.01E_{\rm ex}$ survives there, identified as free $Z_2$ vortex and $Z_2$-vortex pair, respectively, in harmony with the $Z_2$ vortex theory.

preprint2021arXiv

The Nearby Evolved Stars Survey II: Constructing a volume-limited sample and first results from the James Clerk Maxwell Telescope

The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $μ$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust.

preprint2020arXiv

Atacama Compact Array Observations of the Pulsar-Wind Nebula of SNR 0540-69.3

We present observations of the pulsar-wind nebula (PWN) region ofSNR 0540-69.3. The observations were made with the Atacama Compact Array (ACA) in Bands 4 and 6. We also add radio observations from the Australia Compact Array (ATCA) at 3 cm. For 1.449 - 233.50 GHz we obtain a synchrotron spectrum $F_ν \propto ν^{-α_ν}$, with the spectral index $α_ν = 0.17\pm{0.02}$. To conclude how this joins the synchrotron spectrum at higher frequencies we include hitherto unpublished AKARI mid-infrared data, and evaluate published data in the ultraviolet (UV), optical and infrared (IR). In particular, some broad-band filter data in the optical must be discarded from our analysis due to contamination by spectral line emission. For the UV/IR part of the synchrotron spectrum, we arrive at $α_ν = 0.87^{+0.08}_{-0.10}$. There is room for $2.5\times10^{-3}$ solar masses of dust with temperature $\sim 55$ K if there are dual breaks in the synchrotron spectrum, one around $\sim 9\times10^{10}$ Hz, and another at $\sim 2\times10^{13}$ Hz. The spectral index then changes at $\sim 9\times10^{10}$ Hz from $α_ν = 0.14\pm0.07$ in the radio, to $α_ν = 0.35^{-0.07}_{+0.05}$ in the millimetre to far-IR range. The ACA Band 6 data marginally resolves the PWN. In particular, the strong emission 1.5" south-west of the pulsar, seen at other wavelengths, and resolved in the 3-cm data with its 0.8" spatial resolution, is also strong in the millimeter range. The ACA data clearly reveal the supernova remnant shell 20-35 arcsec west of the pulsar, and for the shell we derive $α_ν = 0.64\pm{0.05}$ for the range $8.6-145$~GHz.

preprint2016arXiv

On the properties of dust and gas in the environs of V838 Monocerotis

Herschel FIR imaging and spectroscopy were taken at several epochs to probe the central point source and the extended environment of V838 Mon. PACS and SPIRE maps were used to obtain photometry of the near and far dust around V838 Mon. Fitting reveals 0.5-0.6 solar masses of ~19K dust in the environs (~2.7pc) surrounding the star. The surface-integrated infrared flux (signifying the thermal light echo) and derived dust properties do not vary significantly between the epochs. We also fit the SED of the point source. As the peak of the SED lies outside the Herschel spectral range, it is only by incorporating data from other observatories and epochs that we can perform useful fitting; with this we explicitly assume no evolution of the point source between the epochs. We find warm dust with a temperature of ~300K distributed over a radius of 150-200AU. PACS and SPIRE spectra were also used to detect emission lines from the extended environment around the star. We fit the far-infrared lines of CO arising from the point source, from an extended environment around V838 Mon. Assuming a model of a spherical shell for this gas, we find that the CO appears to arise from two temperature zones: a cold zone (Tkin ~18K) that could be associated with the ISM or possibly with a cold layer in the outermost part of the shell, and a warm (Tkin ~400K) zone that is associated with the extended environment of V838 Mon within a region of radius of ~210AU. The SiO lines arise from a warm/hot zone. We did not fit the lines of H2O as they are far more dependent on the model assumed.

preprint2016arXiv

The ALMA detection of CO rotational line emission in AGB stars in the Large Magellanic Cloud

Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud. Methods: ALMA was used to observe 2 OH/IR stars and 4 carbon stars in the LMC in the CO J= 2-1 line. Results: We present the first measurement of expansion velocities in extragalactic carbon stars. All four C-stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. Gas-to-dust ratios are derived that make the predicted velocities agree with the observed ones. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with "extreme" C-stars that have properties more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than derived from the dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions: Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear similar. This is in agreement with dynamical dust-driven wind models.

preprint2016arXiv

The infrared spectral properties of Magellanic carbon stars

The Infrared Spectrograph on the Spitzer Space Telescope observed 184 carbon stars in the Magellanic Clouds. This sample reveals that the dust-production rate (DPR) from carbon stars generally increases with the pulsation period of the star. The composition of the dust grains follows two condensation sequences, with more SiC condensing before amorphous carbon in metal-rich stars, and the order reversed in metal-poor stars. MgS dust condenses in optically thicker dust shells, and its condensation is delayed in more metal-poor stars. Metal-poor carbon stars also tend to have stronger absorption from C2H2 at 7.5 um. The relation between DPR and pulsation period shows significant apparent scatter, which results from the initial mass of the star, with more massive stars occupying a sequence parallel to lower-mass stars, but shifted to longer periods. Accounting for differences in the mass distribution between the carbon stars observed in the Small and Large Magellanic Clouds reveals a hint of a subtle decrease in the DPR at lower metallicities, but it is not statistically significant. The most deeply embedded carbon stars have lower variability amplitudes and show SiC in absorption. In some cases they have bluer colors at shorter wavelengths, suggesting that the central star is becoming visible. These deeply embedded stars may be evolving off of the AGB and/or they may have non-spherical dust geometries.

preprint2016arXiv

Three-dimensional distribution of ejecta in Supernova 1987A at 10 000 days

Due to its proximity, SN 1987A offers a unique opportunity to directly observe the geometry of a stellar explosion as it unfolds. Here we present spectral and imaging observations of SN 1987A obtained ~10,000 days after the explosion with HST/STIS and VLT/SINFONI at optical and near-infrared wavelengths. These observations allow us to produce the most detailed 3D map of H-alpha to date, the first 3D maps for [Ca II] λλ7292, 7324, [O I] λλ6300, 6364 and Mg II λλ9218, 9244, as well as new maps for [Si I]+[Fe II] 1.644 μm and He I 2.058 μm. A comparison with previous observations shows that the [Si I]+[Fe II] flux and morphology have not changed significantly during the past ten years, providing evidence that it is powered by 44Ti. The time-evolution of H-alpha shows that it is predominantly powered by X-rays from the ring, in agreement with previous findings. All lines that have sufficient signal show a similar large-scale 3D structure, with a north-south asymmetry that resembles a broken dipole. This structure correlates with early observations of asymmetries, showing that there is a global asymmetry that extends from the inner core to the outer envelope. On smaller scales, the two brightest lines, H-alpha and [Si I]+[Fe II] 1.644 μm, show substructures at the level of ~ 200 - 1000 km/s and clear differences in their 3D geometries. We discuss these results in the context of explosion models and the properties of dust in the ejecta.

preprint2015arXiv

Early Dust Formation and a Massive Progenitor for SN 2011ja?

SN 2011ja was a bright (I = -18.3) Type II supernova occurring in the nearby edge on spiral galaxy NGC 4945. Flat-topped and multi-peaked H-alpha and H-beta spectral emission lines appear between 64 - 84 days post-explosion, indicating interaction with a disc-like circumstellar medium inclined 30-45 degrees from edge-on. After day 84 an increase in the H- and K-band flux along with heavy attenuation of the red wing of the emission lines are strong indications of early dust formation, likely located in the cool dense shell created between the forward shock of the SN ejecta and the reverse shock created as the ejecta plows into the existing CSM. Radiative transfer modeling reveals both ~1.5 x 10^-4 Msun of pre-existing dust located ~ 10^16.7 cm away and ~ 5 x 10^-5 Msun of newly formed dust. Spectral observations after 1.5 years reveal the possibility that the fading SN is located within a young (3-6 Myr) massive stellar cluster, which when combined with tentative 56Ni mass estimates of 0.2 Msun may indicate a massive (> 25 Msun) progenitor for SN 2011ja.

preprint2015arXiv

Herschel observations of extreme OH/IR stars - the isotopic ratios of oxygen as a sign-post for the stellar mass

Aim: The late stages of stellar evolution are mainly governed by the mass of the stars. Low- and intermediate-mass stars lose copious amounts of mass during the asymptotic giant branch (AGB) which obscure the central star making it difficult to study the stellar spectra and determine the stellar mass. In this study, we present observational data that can be used to determine lower limits to the stellar mass. Method: Spectra of nine heavily reddened AGB stars taken by the Herschel Space Observatory display numerous molecular emission lines. The strongest emission lines are due to H2O. We search for the presence of isotopologues of H2O in these objects. Result: We detected the 16O and 17O isotopologues of water in these stars, but lines due to H2^{18}O are absent. The lack of 18O is predicted by a scenario where the star has undergone hot-bottom burning which preferentially destroys 18O relative to 16O and 17O. From stellar evolution calculations, this process is thought to occur when the stellar mass is above 5 Msun for solar metallicity. Hence, observations of different isotopologues of H2O can be used to help determine the lower limit to the initial stellar mass. Conclusion: From our observations, we deduce that these extreme OH/IR stars are intermediate-mass stars with masses of >= 5 Msun. Their high mass-loss rates of ~ 1.0e-4 Msun/yr may affect the enrichment of the interstellar medium and the overall chemical evolution of our Galaxy.

preprint2015arXiv

Spitzer Infrared Spectrographic point source classification in the Small Magellanic Cloud

The Magellanic clouds are uniquely placed to study the stellar contribution to dust emission. Individual stars can be resolved in these systems even in the mid-infrared, and they are close enough to allow detection of infrared excess caused by dust.We have searched the Spitzer Space Telescope data archive for all Infrared Spectrograph (IRS) staring-mode observations of the Small Magellanic Cloud (SMC) and found that 209 Infrared Array Camera (IRAC) point sources within the footprint of the Surveying the Agents of Galaxy Evolution in the Small Magellanic Cloud (SAGE-SMC) Spitzer Legacy programme were targeted, within a total of 311 staring mode observations. We classify these point sources using a decision tree method of object classification, based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership and variability information. We find 58 asymptotic giant branch (AGB) stars, 51 young stellar objects (YSOs), 4 post-AGB objects, 22 Red Supergiants (RSGs), 27 stars (of which 23 are dusty OB stars), 24 planetary nebulae (PNe), 10Wolf-Rayet (WR) stars, 3 Hii regions, 3 R Coronae Borealis (R CrB) stars, 1 Blue Supergiant and 6 other objects, including 2 foreground AGB stars. We use these classifications to evaluate the success of photometric classification methods reported in the literature.

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

Carbon-rich dust past the asymptotic giant branch: aliphatics, aromatics, and fullerenes in the Magellanic Clouds

Infrared spectra of carbon-rich objects which have evolved off the asymptotic giant branch reveal a range of dust properties, including fullerenes, polycyclic aromatic hydrocarbons (PAHs), aliphatic hydrocarbons, and several unidentified features, including the 21 um emission feature. To test for the presence of fullerenes, we used the position and width of the feature at 18.7-18.9 um and examined other features at 17.4 and 6-9 um. This method adds three new fullerene sources to the known sample, but it also calls into question three previous identifications. We confirm that the strong 11 um features seen in some sources arise primarily from SiC, which may exist as a coating around carbonaceous cores and result from photo-processing. Spectra showing the 21 um feature usually show the newly defined Class D PAH profile at 7-9 um. These spectra exhibit unusual PAH profiles at 11-14 um, with weak contributions at 12.7 um, which we define as Class D1, or show features shifted to ~11.4, 12.4, and 13.2 um, which we define as Class D2. Alkyne hydrocarbons match the 15.8 um feature associated with 21 um emission. Sources showing fullerene emission but no PAHs have blue colors in the optical, suggesting a clear line of sight to the central source. Spectra with 21 um features and Class D2 PAH emission also show photometric evidence for a relatively clear line of sight to the central source. The multiple associations of the 21 um feature to aliphatic hydrocarbons suggest that the carrier is related to this material in some way.

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

The Mass-Loss Return from Evolved Stars to the Large Magellanic Cloud II: Dust Properties for Oxygen-Rich Asymptotic Giant Branch Stars

We model multi-wavelength broadband UBVIJHKs and Spitzer IRAC and MIPS photometry and IRS spectra from the SAGE and SAGE-Spec observing programs of two oxygen-rich asymptotic giant branch (O-rich AGB) stars in the Large Magellanic Cloud (LMC) using radiative transfer models of dust shells around stars. We chose a star from each of the bright and faint O-rich AGB populations found by earlier studies of the SAGE sample in order to derive a baseline set of dust properties to be used in the construction of an extensive grid of radiative transfer models of the O-rich AGB stars found in the SAGE surveys. From the bright O-rich AGB population we chose HV 5715, and from the faint O-rich AGB population we chose SSTISAGE1C J052206.92-715017.6 (SSTSAGE052206). We found the complex indices of refraction of oxygen-deficient silicates from Ossenkopf et al (1992) and a "KMH"-like grain size distribution with gamma of -3.5, a_min of 0.01 microns, and a_0 of 0.1 microns to be reasonable dust properties for these models. There is a slight indication that the dust around the faint O-rich AGB may be more silica-rich than that around the bright O-rich AGB. Simple models of gas emission suggest a relatively extended gas envelope for the faint O-rich AGB star modeled, consistent with the relatively large dust shell inner radius for the same model. The models compute the dust mass-loss rates for the two stars to be 2.0x10^-9 M_Sun/yr and 2.3x10^-9 M_Sun/yr, respectively. When a dust-to-gas mass ratio of 0.002 is assumed for SSTSAGE052206 and HV 5715, the dust mass-loss rates imply total mass-loss rates of 1.0x10^-6 M_Sun/yr and 1.2x10^-6 M_Sun/yr, respectively. These properties of the dust shells and stars, as inferred from our models of the two stars, are found to be consistent with properties observed or assumed by detailed studies of other O-rich AGB stars in the LMC and elsewhere.

preprint2014arXiv

The timing and location of dust formation in the remnant of SN 1987A

The discovery with the {\it Herschel Space Observatory} of bright far infrared and submm emission from the ejecta of the core collapse supernova SN\,1987A has been interpreted as indicating the presence of some 0.4--0.7\,M$_\odot$ of dust. We have constructed radiative transfer models of the ejecta to fit optical to far-infrared observations from the literature at epochs between 615 days and 24 years after the explosion, to determine when and where this unexpectedly large amount of dust formed. We find that the observations by day 1153 are consistent with the presence of 3$\times$10$^{-3}$M$_\odot$ of dust. Although this is a larger amount than has previously been considered possible at this epoch, it is still very small compared to the amount present in the remnant after 24 years, and significantly higher dust masses at the earlier epochs are firmly ruled out by the observations, indicating that the majority of the dust must have formed at very late times. By 8515-9200 days after the explosion, 0.6--0.8\,M$_\odot$ of dust is present, and dust grains with radii greater than 2\,$μ$m are required to obtain a fit to the observed SED. This suggests that the dust mass increase at late times was caused by accretion onto and coagulation of the dust grains formed at earlier epochs. These findings provide further confirmation that core collapse supernovae can create large quantities of dust, and indicate that the reason for small dust masses being estimated in many cases is that the vast majority of the dust forms long after most supernovae have been detectable at mid-infrared wavelengths.

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.

preprint2014arXiv

VLT observations of the asymmetric Etched Hourglass Nebula, MyCn 18

Context. The mechanisms that form extreme bipolar planetary nebulae remain unclear. Aims. The physical properties, structure, and dynamics of the bipolar planetary nebula, MyCn 18, are investigated in detail with the aim of understanding the shaping mechanism and evolutionary history of this object. Methods. VLT infrared images, VLT ISAAC infrared spectra, and long-slit optical Echelle spectra are used to investigate MyCn 18. Morpho-kinematic modelling was used to firmly constrain the structure and kinematics of the source. A timescale analysis was used to determine the kinematical age of the nebula and its main components. Results. A spectroscopic study of MyCn 18's central and offset region reveals the detailed make-up of its nebular composition. Molecular hydrogen, atomic helium, and Bracket gamma emission are detected from the central regions of MyCn 18. ISAAC spectra from a slit position along the narrow waist of the nebula demonstrate that the ionised gas resides closer to the centre of the nebula than the molecular emission. A kinematical age of the nebula and its components were obtained by the P-V arrays and timescale analysis. Conclusions. The structure and kinematics of MyCn 18 are better understood using an interactive 3-D modelling tool called shape. A dimensional and timescale analysis of MyCn 18's major components provides a possible mechanism for the nebula's asymmetry. The putative central star is somewhat offset from the geometric centre of the nebula, which is thought to be the result of a binary system. We speculate that the engulfing and destruction of an exoplanet during the AGB phase may have been a key event in shaping MyCn 18 and generating of its hypersonic knotty outflow.

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'' 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

Detection of a Noble Gas Molecular Ion, {36}ArH^+, in the Crab Nebula

Noble gas molecules have not hitherto been detected in space. From spectra obtained with the Herschel Space Observatory, we report the detection of emission in the 617.5 GHz and 1234.6 GHz J = 1-0 and 2-1 rotational lines of {36}ArH^+ at several positions in the Crab Nebula, a supernova remnant known to contain both H2 molecules and regions of enhanced ionized argon emission. {36}Ar is believed to have originated from explosive nucleosynthesis in massive stars during core-collapse supernova events. Its detection in the Crab Nebula, the product of such a supernova event, confirms this expectation. The likely excitation mechanism for the observed {36}ArH^+ emission lines is electron collisions in partially ionized regions with electron densities of a few hundred per centimeter cubed.

preprint2013arXiv

Dust Production and Particle Acceleration in Supernova 1987A Revealed with ALMA

Supernova (SN) explosions are crucial engines driving the evolution of galaxies by shock heating gas, increasing the metallicity, creating dust, and accelerating energetic particles. In 2012 we used the Atacama Large Millimeter/Submillimeter Array to observe SN 1987A, one of the best-observed supernovae since the invention of the telescope. We present spatially resolved images at 450um, 870um, 1.4mm, and 2.8mm, an important transition wavelength range. Longer wavelength emission is dominated by synchrotron radiation from shock-accelerated particles, shorter wavelengths by emission from the largest mass of dust measured in a supernova remnant (>0.2Msun). For the first time we show unambiguously that this dust has formed in the inner ejecta (the cold remnants of the exploded star's core). The dust emission is concentrated to the center of the remnant, so the dust has not yet been affected by the shocks. If a significant fraction survives, and if SN 1987A is typical, supernovae are important cosmological dust producers.

preprint2013arXiv

OH/IR stars and their superwinds as observed by the Herschel Space Observatory

Aim : In order to study the history of mass loss in extreme OH/IR stars, we observed a number of these objects using CO as a tracer of the density and temperature structure of their circumstellar envelopes. Method : Combining CO observations from the Herschel Space Observatory with those from the ground, we trace mass loss rates as a function of radius in five extreme OH/IR stars. Using radiative transfer modelling, we modelled the dusty envelope as well as the CO emission. The high-rotational transitions of CO indicate that they originate in a dense superwind region close to the star while the lower transitions tend to come from a more tenuous outer wind which is a result of the mass loss since the early AGB phase. Result : The models of the circumstellar envelopes around these stars suggest that they have entered a superwind phase in the past 200 - 500 years. The low 18O/17O (~ 0.1 compared to the solar abundance ratio of ~ 5) and 12C/13C (3-30 cf. the solar value of 89) ratios derived from our study support the idea that these objects have undergone hot-bottom burning and hence that they are massive M >= 5 solar-mass AGB stars.

preprint2012arXiv

A Cool Dust Factory in the Crab Nebula: A Herschel Study of the Filaments

Whether supernovae are major sources of dust in galaxies is a long-standing debate. We present infrared and submillimeter photometry and spectroscopy from the Herschel Space Observatory of the Crab Nebula between 51 and 670 micron as part of the Mass Loss from Evolved StarS program. We compare the emission detected with Herschel with multiwavelength data including millimeter, radio, mid-infrared and archive optical images. We carefully remove the synchrotron component using the Herschel and Planck fluxes measured in the same epoch. The contribution from line emission is removed using Herschel spectroscopy combined with Infrared Space Observatory archive data. Several forbidden lines of carbon, oxygen and nitrogen are detected where multiple velocity components are resolved, deduced to be from the nitrogen-depleted, carbon-rich ejecta. No spectral lines are detected in the SPIRE wavebands; in the PACS bands, the line contribution is 5% and 10% at 70 and 100 micron and negligible at 160 micron. After subtracting the synchrotron and line emission, the remaining far-infrared continuum can be fit with two dust components. Assuming standard interstellar silicates, the mass of the cooler component is 0.24(+0.32)(-0.08) Msolar for T = 28.1(+5.5)(-3.2)K. Amorphous carbon grains require 0.11 +/- 0.01 Msolar of dust with T = 33.8(+2.3)(-1.8) K. A single temperature modified-blackbody with 0.14Msolar and 0.08Msolar for silicate and carbon dust respectively, provides an adequate fit to the far- infrared region of the spectral energy distribution but is a poor fit at 24-500 micron. The Crab Nebula has condensed most of the relevant refractory elements into dust, suggesting the formation of dust in core-collapse supernova ejecta is efficient.

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

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's lifetime on the asymptotic giant branch.

preprint2012arXiv

Dust in Historical Galactic Type Ia Supernova Remnants with Herschel

The origin of interstellar dust in galaxies is poorly understood, particularly the relative contributions from supernovae and the cool stellar winds of low-intermediate mass stars. Here, we present Herschel PACS and SPIRE photometry at 70-500um of the historical young supernova remnants: Kepler and Tycho; both thought to be the remnants of Type Ia explosion events. We detect a warm dust component in Kepler's remnant with T = 82K and mass 0.0031Msun; this is spatially coincident with thermal X-ray emission optical knots and filaments, consistent with the warm dust originating in the circumstellar material swept up by the primary blast wave of the remnant. Similarly for Tycho's remnant, we detect warm dust at 90K with mass 0.0086Msun. Comparing the spatial distribution of the warm dust with X-rays from the ejecta and swept-up medium, and Ha emission arising from the post-shock edge, we show that the warm dust is swept up interstellar material. We find no evidence of a cool (25-50 K) component of dust with mass >0.07Msun as observed in core-collapse remnants of massive stars. Neither the warm or cold dust components detected here are spatially coincident with supernova ejecta material. We compare the lack of observed supernova dust with a theoretical model of dust formation in Type Ia remnants which predicts dust masses of 0.088(0.017)Msun for ejecta expanding into surrounding densities of 1(5)cm-3. The model predicts that silicon- and carbon-rich dust grains will encounter the interior edge of the observed dust emission at 400 years confirming that the majority of the warm dust originates from swept up circumstellar or interstellar grains (for Kepler and Tycho respectively). The lack of cold dust grains in the ejecta suggests that Type Ia remnants do not produce substantial quantities of iron-rich dust grains and has important consequences for the 'missing' iron mass observed in ejecta.

preprint2012arXiv

Modified Cross-Correlation for Efficient White-Beam Inelastic Neutron Scattering Spectroscopy

We describe a method of white-beam inelastic neutron scattering for improved measurement efficiency. The method consists of matrix inversion and selective extraction. The former is to resolve each incident energy component from the white-beam data, and the latter eliminates contamination by elastic components, which produce strong backgrounds that otherwise obfuscate the inelastic scattering components. In this method, the optimal experimental condition to obtain high efficiency will strongly depend on the specific aim of the individual experiments.

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

Radiation Transfer in the Cavity and Shell of Planetary Nebulae

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

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.

preprint2011arXiv

Herschel Detects a Massive Dust Reservoir in Supernova 1987A

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

preprint2011arXiv

Temperature independent diffuse scattering and elastic lattice deformations in relaxor PbMg1/3Nb2/3O3

The results of diffuse neutron scattering experiment on PbMg1/3Nb2/3O3 single crystal above the Burns temperature are reported. It is shown that the high temperature elastic diffuse component is highly anisotropic in low-symmetry Brillouin zones and this anisotropy can be described using Huang scattering formalism assuming that the scattering originates from mesoscopic lattice deformations due to elastic defects. The qualitative agreement between this model and the experimental data is achieved with simple isotropic defects. It is demonstrated that weak satellite maxima near the Bragg reflections can be interpreted as the finite resolution effect.

preprint2011arXiv

The Circumstellar Environment of R Coronae Borealis: White Dwarf Merger or Final Helium Shell Flash?

In 2007, R Coronae Borealis (R CrB) went into an historically deep and long decline. In this state, the dust acts like a natural coronagraph at visible wavelengths, allowing faint nebulosity around the star to be seen. Imaging has been obtained from 0.5 to 500 micron with Gemini/GMOS, HST/WFPC2, Spitzer/MIPS, and Herschel/SPIRE. Several of the structures around R CrB are cometary globules caused by wind from the star streaming past dense blobs. The estimated dust mass of the knots is consistent with their being responsible for the R CrB declines if they form along the line of sight to the star. In addition, there is a large diffuse shell extending up to 4 pc away from the star containing cool 25 K dust that is detected all the way out to 500 micron. The SED of R CrB can be well fit by a 150 AU disk surrounded by a very large diffuse envelope which corresponds to the size of the observed nebulosity. The total masses of the disk and envelope are 10^-4 and 2 M(Sun), respectively, assuming a gas-to-dust ratio of 100. The evidence pointing toward a white-dwarf merger or a final-helium-shell flash origin for R CrB is contradictory. The shell and the cometary knots are consistent with a fossil planetary nebula. Along with the fact that R CrB shows significant Lithium in its atmosphere, this supports the final-helium-shell flash. However, the relatively high inferred mass of R CrB and its high fluorine abundance support a white-dwarf merger.

preprint2011arXiv

The Dust Properties of Two Hot R Coronae Borealis Stars and a Wolf-Rayet Central Star of a Planetary Nebula: in Search of a Possible Link

We present new Spitzer/IRS spectra of two hot R Coronae Borealis (RCB) stars, one in the Galaxy,V348 Sgr, and one lying in the LMC, HV 2671. These two objects may constitute a link between the RCB stars and the late Wolf-Rayet ([WCL]) class of central stars of planetary nebula (CSPNe) such as CPD -56 8032 that has little or no hydrogen in their atmospheres. HV 2671 and V348 Sgr are members of a rare subclass that has significantly higher effective temperatures than most RCB stars, but sharing the traits of hydrogen deficiency and dust formation that define the cooler RCB stars. The [WC] CSPNe star, CPD -56 8032, displays evidence for dual-dust chemistry showing both PAHs and crystalline silicates in its mid-IR spectrum. HV 2671 shows strong PAH emission but shows no sign of having crystalline silicates. The spectrum of V348 Sgr is very different from those of CPD -56 8032 and HV 2671. The PAH emission seen strongly in the other two stars is not present. Instead, the spectrum is dominated by a broad emission centered at about 8.2 micron. The mid-IR spectrum of CPD -56 8032 shows emission features that may be associated with C60. The other two stars do not show evidence for C60. HV 2671 has also been detected by Herschel/PACS and SPIRE. V348 Sgr and CPD -56 8032 have been detected by AKARI/FIS. These data were combined with Spitzer, IRAS, 2MASS and other photometry to produce their spectral energy distributions from the visible to the far-IR. Monte Carlo radiative transfer modeling was used to study the circumstellar dust around these stars. HV 2671 and CPD -56 8032 require both a flared inner disk with warm dust and an extended diffuse envelope with cold dust to to fit their SEDs. The SED of V348 Sgr can be fit with a much smaller disk and envelope.

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.

preprint2010arXiv

AKARI Near- to Mid-Infrared Imaging and Spectroscopic Observations of the Small Magellanic Cloud. I. Bright Point Source List

We carried out a near- to mid-infrared imaging and spectroscopic observations of the patchy areas in the Small Magellanic Cloud using the Infrared Camera on board AKARI. Two 100 arcmin2 areas were imaged in 3.2, 4.1, 7, 11, 15, and 24 um and also spectroscopically observed in the wavelength range continuously from 2.5 to 13.4 um. The spectral resolving power (lambda/Delta lambda) is about 20, 50, and 50 at 3.5, 6.6 and 10.6 um, respectively. Other than the two 100 arcmin2 areas, some patchy areas were imaged and/or spectroscopically observed as well. In this paper, we overview the observations and present a list of near- to mid-infrared photometric results, which lists ~ 12,000 near-infrared and ~ 1,800 mid-infrared bright point sources detected in the observed areas. The 10 sigma limits are 16.50, 16.12, 13.28, 11.26, 9.62, and 8.76 in Vega magnitudes at 3.2, 4.1, 7, 11, 15, and 24 um bands, respectively.

preprint2010arXiv

AKARI's infrared view on nearby stars : Using AKARI Infrared Camera All-Sky Survey, 2MASS, and Hipparcos catalog

--Results-- We found that the (B-V) v.s. (V-S9W) color-color diagram is useful to identify the stars with infrared excess emerged from circumstellar envelopes/disks. Be stars with infrared excess are well separated from other types of stars in this diagram. Whereas (J-L18W) v.s. (S9W-L18W) diagram is a powerful tool to classify several object-types. Carbon-rich asymptotic giant branch (AGB) stars and OH/IR stars form distinct sequences in this color-color diagram. Young stellar objects (YSOs), pre-main sequence (PMS) stars, post-AGB stars and planetary nebulae (PNe) have largest mid-infrared color-excess, and can be identified in infrared catalog. Finally, we plot L18W v.s. (S9W-L18W) color-magnitude diagram, using the AKARI data together with Hipparcos parallaxes. This diagram can be used to identify low-mass YSOs, as well as AGB stars. We found that this diagram is comparable to the [24] vs ([8.0]-[24]) diagram of Large Magellanic Cloud sources using the Spitzer Space Telescope data. Our understanding of Galactic objects will be used to interpret color-magnitude diagram of stellar populations in nearby galaxies which Spitzer Space Telescope has observed. --Conclusions-- Our study of the AKARI color-color and color-magnitude will be used to explore properties of unknown objects in future. In addition, our analysis highlights a future key project to understand stellar evolution with circumstellar envelope, once the forthcoming astronometrical data with GAIA are available.

preprint2010arXiv

AKARI's infrared view on nearby stars : Using AKARI Infrared Camera All-Sky Survey, 2MASS, and Hipparcos catalog

Results: We found that the (B-V) v.s. (V-S9W) color-color diagram is useful to identify the stars with infrared excess emerged from circumstellar envelopes/disks. Be stars with infrared excess are well separated from other types of stars in this diagram. Whereas (J-L18W) v.s. (S9W-L18W) diagram is a powerful tool to classify several object-types. Carbon-rich asymptotic giant branch (AGB) stars and OH/IR stars form distinct sequences in this color-color diagram. Young stellar objects (YSOs), pre-main sequence (PMS) stars, post-AGB stars and planetary nebulae (PNe) have largest mid-infrared color-excess, and can be identified in infrared catalog. Finally, we plot L18W v.s. (S9W-L18W) color-magnitude diagram, using the AKARI data together with Hipparcos parallaxes. This diagram can be used to identify low-mass YSOs, as well as AGB stars. We found that this diagram is comparable to the [24] vs ([8.0]-[24]) diagram of Large Magellanic Cloud sources using the Spitzer Space Telescope data. Our understanding of Galactic objects will be used to interpret color-magnitude diagram of stellar populations in nearby galaxies which Spitzer Space Telescope has observed. Conclusions: Our study of the AKARI color-color and color-magnitude will be used to explore properties of unknown objects in future. In addition, our analysis highlights a future key project to understand stellar evolution with circumstellar envelope, once the forthcoming astronometrical data with GAIA are available.

preprint2010arXiv

Clouds, filaments and protostars: the Herschel Hi-GAL Milky Way

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2° x 2° tiles approximately centered at l=30° and l=59°. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59° field; a A_V value between 5 and 10 is found for the l=30° field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

preprint2010arXiv

Cold dust clumps in dynamically hot gas

We present clumps of dust emission from Herschel observations of the Large Magellanic Cloud (LMC) and their physical and statistical properties. We catalog cloud features seen in the dust emission from Herschel observations of the LMC, the Magellanic type irregular galaxy closest to the Milky Way, and compare these features with HI catalogs from the ATCA+Parkes HI survey. Using an automated cloud-finding algorithm, we identify clouds and clumps of dust emission and examine the cumulative mass distribution of the detected dust clouds. The mass of cold dust is determined from physical parameters that we derive by performing spectral energy distribution fits to 250, 350, and 500 micronm emission from SPIRE observations using DUSTY and GRASIL radiative transfer calculation with dust grain size distributions for graphite/silicate in low-metallicity extragalactic environments. The dust cloud mass spectrum follows a power law distribution with an exponent of gamma=-1.8 for clumps larger than 400 solar mass and is similar to the HI mass distribution. This is expected from the theory of ISM structure in the vicinity of star formation.

preprint2010arXiv

Cold Dust in Three Massive Evolved Stars in the LMC

Massive evolved stars can produce large amounts of dust, and far-infrared (IR) data are essential for determining the contribution of cold dust to the total dust mass. Using Herschel, we search for cold dust in three very dusty massive evolved stars in the Large Magellanic Cloud: R71 is a Luminous Blue Variable, HD36402 is a Wolf-Rayet triple system, and IRAS05280-6910 is a red supergiant. We model the spectral energy distributions using radiative transfer codes and find that these three stars have mass-loss rates up to 10^-3 solar masses/year, suggesting that high-mass stars are important contributors to the life-cycle of dust. We found far-IR excesses in two objects, but these excesses appear to be associated with ISM and star-forming regions. Cold dust (T < 100 K) may thus not be an important contributor to the dust masses of evolved stars.

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

Determining dust temperatures and masses in the Herschel era: the importance of observations longward of 200 micron

The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (>100 micron). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500 micron that is beyond the peak (~160 micron) of dust emission in most galaxies. We investigate the differences in the fitted dust temperatures and masses determined using only <200 micron data and then also including >200 micron data (new SPIRE observations) to determine how important having >200 micron data is for deriving these dust properties. We fit the 100 to 350 micron observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single temperature and fixed emissivity law. The data used are existing observations at 100 and 160 micron (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE Key Project as part of the Herschel Science Demonstration phase. The dust temperatures and masses computed using only 100 and 160 micron data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350 micron data. We find that an emissivity law proportional to lambda^-1.5 minimizes the 100-350 micron fractional residuals. We find that the emission at 500 micron is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500 micron excess is weakly anti-correlated with MIPS 24 micron flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500 micron excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent emissivity law or a population of very cold dust only detectable at lambda > 500 micron for the origin of the 500 micron excess.

preprint2010arXiv

Dust in the bright supernova remnant N49 in the LMC

We investigate the dust associated with the supernova remnant (SNR) N49 in the Large Magellanic Cloud (LMC) as observed with the Herschel Space Observatory. N49 is unusually bright because of an interaction with a molecular cloud along its eastern edge. We have used PACS and SPIRE to measure the far IR flux densities of the entire SNR and of a bright region on the eastern edge of the SNR where the SNR shock is encountering the molecular cloud. Using these fluxes supplemented with archival data at shorter wavelengths, we estimate the dust mass associated with N49 to be about 10 Msun. The bulk of the dust in our simple two-component model has a temperature of 20-30 K, similar to that of nearby molecular clouds. Unfortunately, as a result of the limited angular resolution of Herschel at the wavelengths sampled with SPIRE, the uncertainties are fairly large. Assuming this estimate of the dust mass associated with the SNR is approximately correct, it is probable that most of the dust in the SNR arises from regions where the shock speed is too low to produce significant X-ray emission. The total amount of warm 50-60 K dust is ~0.1 or 0.4 Msun, depending on whether the dust is modeled in terms of carbonaceous or silicate grains. This provides a firm lower limit to the amount of shock heated dust in N49.

preprint2010arXiv

Herschel Inventory of The Agents of Galaxy Evolution (HERITAGE): the Large Magellanic Cloud dust

The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory&#39;s PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 μm and SPIRE 250, 350 and 500 μm. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 μm emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.

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

Spitzer Spectroscopy of Mass Loss and Dust Production by Evolved Stars in Globular Clusters

We have observed a sample of 35 long-period variables and four Cepheid variables in the vicinity of 23 Galactic globular clusters using the Infrared Spectrograph on the Spitzer Space Telescope. The long-period variables in the sample cover a range of metallicities from near solar to about 1/40th solar. The dust mass-loss rate from the stars increases with pulsation period and bolometric luminosity. Higher mass-loss rates are associated with greater contributions from silicate grains. The dust mass-loss rate also depends on metallicity. The dependence is most clear when segregating the sample by dust composition, less clear when segregating by bolometric magnitude, and absent when segregating by period. The spectra are rich in solid-state and molecular features. Emission from alumina dust is apparent across the range of metallicities. Spectra with a 13-um dust emission feature, as well as an associated feature at 20 um, also appear at most metallicities. Molecular features in the spectra include H_2O bands at 6.4-6.8 um, seen in both emission and absorption, SO_2 absorption at 7.3-7.5 um, and narrow emission bands from CO_2 from 13.5 to 16.8 um. The star Lynga 7 V1 has an infrared spectrum revealing it to be a carbon star, adding to the small number of carbon stars associated with Galactic globular clusters.

preprint2010arXiv

SPITZER-IRS spectral fitting of discs around binary post-AGB stars - Corrigendum

Recently, we have discovered an error in our Monte-Carlo spectral fitting routine, more specifically where the errors on the fluxes were rescaled to get a reduced chi2 of 1. The rescaled errors were too big, resulting in too wide a range of good fits in our 100 step Monte-Carlo routine. This problem affects Figs. 7-9 and Tables A.1, A.2 in Gielen et al. (2008), Table 3 in Gielen et al. (2009a), and Table 4 in Gielen et al. (2009b). We corrected for this error and present the new values and errors in the tables below. The new values and errors nearly all fall within the old error range. Our best chi2 values and overall former scientific results are not affected. With these new errors some possible new trends in the dust parameters might be observed. These will be discussed in an upcoming paper where we extend the sample presented in Gielen et al. (2008) with newly obtained SPITZER-IRS data.

preprint2010arXiv

The Global Gas and Dust budge of the Large Magellanic Cloud --- Importance of Asymptotic Giant Branch stars

It is still an unresolved problem how much AGB stars can contribute to the overall gas and dust enrichment processes in the interstellar medium within galaxies. We start tackling this problem, by using our test case observational data from the Large Magellanic Cloud (LMC), from which we obtain the global gas and dust budget. The photometric data from the LMC is obtained with the Spitzer Space Telescope. We established an infrared colour classification scheme to select AGB stars, which are based on spectroscopically identified AGB stars. We further confirm a correlation between the Spitzer colour and mass-loss rate, which leads to a measurement of the total mass-loss rate from the entire AGB population in the LMC. Indeed, AGB stars are an important gas and dust source.

preprint2010arXiv

The mass-loss return from evolved stars to the Large Magellanic Cloud III. Dust properties for carbon-rich asymptotic giant branch stars

We present a 2Dust model for the dust shell around a LMC long-period variable (LPV) previously studied as part of the OGLE survey. OGLE LMC LPV 28579 (SAGE J051306.40-690946.3) is a carbon-rich asymptotic giant branch (AGB) star for which we have photometry and spectra from the Spitzer SAGE and SAGE-Spec programs along with UBVIJHK_s photometry. By modeling this source, we obtain a baseline set of dust properties to be used in the construction of a grid of models for carbon stars. We reproduce its spectral energy distribution using a mixture of AmC and SiC (15% by mass). The grain sizes are distributed according to the KMH model. The best-fit model has an optical depth of 0.28 for the shell at the peak of the SiC feature, with R_in~1430 R_sun or 4.4 R_star. The temperature at this inner radius is 1310 K. Assuming an expansion velocity of 10 km s^-1, we obtain a dust mass-loss rate of 2.5x10^-9 M_sun yr-1. We calculate a 15% variation in this rate by testing the fit sensitivity against variation in input parameters. We also present a simple model for the molecular gas in the extended atmosphere that could give rise to the 13.7 μm feature seen in the spectrum. We find that a combination of CO and C_2H_2 gas at an excitation temperature of about 1000 K and column densities of 3x10^21 cm^-2 and 10^19 cm^-2 respectively are able to reproduce the observations. Given that the excitation temperature is close to T_dust(R_in), most of the molecular contribution probably arises from the inner shell region. The luminosity corresponding to the first epoch of SAGE observations is 6580 L_sun. For an effective temperature of about 3000 K, this implies a stellar mass of 1.5-2 M_sun and an age of 1-2.5 Gyr. For a gas:dust ratio of 200, we obtain a gas mass-loss rate of 5.0x10^-7 M_sun yr^-1, consistent with the gas mass-loss rates estimated from the period, color and 8 μm flux of the source.

preprint2010arXiv

The SAGE-Spec Spitzer Legacy program: The life-cycle of dust and gas in the Large Magellanic Cloud

The SAGE-Spec Spitzer Legacy program is a spectroscopic follow-up to the SAGE-LMC photometric survey of the Large Magellanic Cloud carried out with the Spitzer Space Telescope. We present an overview of SAGE-Spec and some of its first results. The SAGE-Spec program aims to study the life cycle of gas and dust in the Large Magellanic Cloud, and to provide information essential to the classification of the point sources observed in the earlier SAGE-LMC photometric survey. We acquired 224.6 hours of observations using the InfraRed Spectrograph and the SED mode of the Multiband Imaging Photometer for Spitzer. The SAGE-Spec data, along with archival Spitzer spectroscopy of objects in the Large Magellanic Cloud, are reduced and delivered to the community. We discuss the observing strategy, the specific data reduction pipelines applied and the dissemination of data products to the scientific community. Initial science results include the first detection of an extragalactic &#34;21 um&#34; feature towards an evolved star and elucidation of the nature of disks around RV Tauri stars in the Large Magellanic Cloud. Towards some young stars, ice features are observed in absorption. We also serendipitously observed a background quasar, at a redshift of z~0.14, which appears to be host-less.

preprint2010arXiv

The SAGE-Spec Spitzer Legacy program: The life-cycle of dust and gas in the Large Magellanic Cloud. Point source classification I

We present the classification of 197 point sources observed with the Infrared Spectrograph in the SAGE-Spec Legacy program on the Spitzer Space Telescope. We introduce a decision-tree method of object classification based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership, and variability information, which is used to classify the SAGE-Spec sample of point sources. The decision tree has a broad application to mid-infrared spectroscopic surveys, where supporting photometry and variability information are available. We use these classifications to make deductions about the stellar populations of the Large Magellanic Cloud and the success of photometric classification methods. We find 90 asymptotic giant branch (AGB) stars, 29 young stellar objects, 23 post-AGB objects, 19 red supergiants, eight stellar photospheres, seven background galaxies, seven planetary nebulae, two HII regions and 12 other objects, seven of which remain unclassified.

preprint2009arXiv

Analysis of the infrared spectra of the peculiar post-AGB stars EPLyr and HD52961

Aim: We aim to study in detail the peculiar mineralogy and structure of the circumstellar environment of two binary post-AGB stars, EPLyr and HD52961. Both stars were selected from a larger sample of evolved disc sources observed with Spitzer and show unique solid-state and gas features in their infrared spectra. Moreover, they show a very small infrared excess in comparison with the other sample stars. Methods: The different dust and gas species are identified on the basis of high-resolution Spitzer-IRS spectra. We fit the full spectrum to constrain grain sizes and temperature distributions in the discs. This, combined with our broad-band spectral energy distribution and interferometric measurements, allows us to study the physical structure of the disc, using a self-consistent 2D radiative-transfer disc model. Results: We find that both stars have strong emission features due to CO_2 gas, dominated by ^{12}C^{16}O_2, but with clear ^{13}C^{16}O_2 and even ^{16}O^{12}C^{18}O isotopic signatures. Crystalline silicates are apparent in both sources but proved very hard to model. EP Lyr also shows evidence of mixed chemistry, with emission features of the rare class-C PAHs. Whether these PAHs reside in the oxygen-rich disc or in a carbon-rich outflow is still unclear. With the strongly processed silicates, the mixed chemistry and the low ^{12}C/^{13}C ratio, EP Lyr resembles some silicate J-type stars, although the depleted photosphere makes nucleosynthetic signatures difficult to probe. We find that the disc environment of both sources is, to a first approximation, well modelled with a passive disc, but additional physics such as grain settling, radial dust distributions, and an outflow component must be included to explain the details of the observed spectral energy distributions in both stars.

preprint2009arXiv

Chemical depletion in the Large Magellanic Cloud: RV Tauri stars and the photospheric feedback from their dusty discs

Aims: By studying the photospheric abundances of 4 RV Tauri stars in the LMC, we test whether the depletion pattern of refractory elements, seen in similar Galactic sources, is also common for extragalactic sources. Since this depletion process probably only occurs through interaction with a stable disc, we investigate the circumstellar environment of these sources. Methods: A detailed photospheric abundance study was performed using high-resolution UVES optical spectra. To study the circumstellar environment we use photometric data to construct the spectral energy distributions of the stars, and determine the geometry of the circumstellar environment, whereas low-resolution Spitzer-IRS infrared spectra are used to trace its mineralogy. Results: Our results show that, also in the LMC, the photospheres of RV Tauri stars are commonly affected by the depletion process, although it can differ significantly in strength from source to source. From our detailed disc modelling and mineralogy study, we find that this process, as in the Galaxy, appears closely related to the presence of a stable Keplerian disc. The newly studied extragalactic objects have similar observational characteristics as Galactic post-AGB binaries surrounded by a dusty disc, and are therefore also believed to be part of a binary system. One source shows a very small infrared excess, atypical for a disc source, but still has evidence for depletion. We speculate this could point to the presence of a very evolved disc, similar to debris discs seen around young stellar objects.

preprint2009arXiv

Galactic AGB stars from the IPHAS survey

The INT Photometric H-Alpha Survey (IPHAS) is particular effective in tracing asymptotic giant branch (AGB) stars across the Galactic disk due to its use of broad-band r&#39; and i&#39; filters, and the clear separation of the dwarf and giant branches in the colour-colour diagram. Here we discuss some results from our studies of the Galactic AGB population that has been revealed by IPHAS. This includes a photometric study, our spectroscopic follow-up campaign, and the photometric technique that allows us to separate O-rich, S-type and carbon stars across the entire Galactic disk based on their IPHAS photometry alone. We also discuss the potential impact of this technique on our understanding of the dredge-up process in AGB stars, Galactic chemical evolution, and the structure and properties of the Galactic disk.

preprint2009arXiv

Knots in the Helix Nebula found in H2

We present a deep and wide field-of-view (4&#39;x 7&#39;) image of the planetary nebula (PN) NGC 7293 (the Helix Nebula) in the 2.12 micron H2 v=1-0 S(1) line. The excellent seeing (0.4&#39;&#39;) at the Subaru Telescope, allows the details of cometary knots to be examined. The knots are found at distances of 2.2&#39;-6.4&#39; from the central star (CS). At the inner edge and in the inner ring (up to 4.5&#39; fromthe CS), the knot often show a `tadpole&#39; shape, an elliptical head with a bright crescent inside and a long tail opposite to the CS. In detail, there are variations in the tadpole shapes, such as narrowing tails, widening tails, meandering tails, or multi-peaks within a tail. In the outer ring (4.5&#39;-6.4&#39; from the CS), the shapes are more fractured, and the tails do not collimate into a single direction. The transition in knot morphology from the inner edge to the outer ring is clearly seen. The number density of knots governs the H2 surface brightness in the inner ring: H2 exists only within the knots. Possible mechanisms which contribute to the shaping of the knots are discussed, including photo-ionization and streaming motions. A plausible interpretation of our images is that inner knots are being overrun by a faster wind, but that this has not (yet) reached the outer knots. Based on H2 formation and destruction rates, H2 gas can survive in knots from formation during the late asymptotic giant branch (AGB) phase throughout the PN phase. These observations provide new constraints on the formation and evolution of knots, and on the physics of molecular gas embedded within ionized gas.

preprint2009arXiv

Near-IR spectra of IPHAS extremely red Galactic AGB stars

We present a library of 139 near-IR spectra of cool asymptotic giant branch stars that will be useful for comparison with theoretical model atmosphere calculations and for modeling the integrated emission from intermediate-age stellar populations. The source list was selected from the `extremely red&#39; region of the INT Photometric H Alpha Survey (IPHAS) colour-colour plane that is overwhelmingly dominated by very late-type stars. The spectral library also includes a large fraction of S-type and carbon stars. We present a number of spectral classification sequences highlighting the various molecular features identified and discuss a number of rare features with uncertain identifications in the literature. With its focus on particularly cool photospheres this catalogue serves as a companion to recent spectroscopic atlases of MK standards in the near-IR. Finally the relationship between IPHAS (r&#39;-i&#39;) and (r&#39;-H Alpha) colours and spectroscopically determined properties is discussed and a strong correlation between (r&#39;-H Alpha) colour and the C/O abundance index for S-type and carbon stars is noted. This relation has the potential to separate O-rich, S-type and carbon stars in the Galaxy based on their photometry alone.

preprint2009arXiv

Renormalization of Commensurate Magnetic Peak in Ni-doped La$_{1.85}$Sr$_{0.15}$CuO$_{4}$

We have studied the magnetic excitations in impurity doped La$_{1.85}$Sr$_{0.15}$Cu$_{1-y}$A$_{y}$O$_{4}$ (A=Ni or Zn) by neutron scattering. The dispersion for Zn:$y=0.017$ is similar to that for the impurity free sample: incommensurate peaks with the incommensurability $δ=0.12\pm0.01$ (rlu) do not change their positions up to 21 meV. On the other hand, for Ni:$y=0.029$, two incommensurate peaks observed at low energies suddenly change into a broad commensurate peak at $E_\mathrm{cross}=15$ meV. Compared to the impurity free sample with a similar Sr-concentration $x=0.16$, [B. Vignolle {\it et al.} Nature Physics {\bf 3} (2007) 163], $E_\mathrm{cross}$ for Ni:$y=0.029$ is decreased by nearly the same factor for the reduction in $T_{c}$. This is very similar to the shift of the resonance energy ($E_\mathrm{res}$) in Ni-doped YBa$_{2}$Cu$_{3}$O$_{7}$.[Y. Sidis {\it et al.}: Phys. Rev. Lett. {\bf 84} (2000) 5900]. These common impurity effects on the shift of $E_\mathrm{cross}$ and $E_\mathrm{res}$ suggest the same magnetic origin for the resonance peak in YBa$_{2}$Cu$_{3}$O$_δ$ and that for a crossing point of upward and downward dispersions in the La$_{2-x}$Sr$_{x}$CuO$_{4}$. We propose that the sudden change in the dispersion is better described by a crossover from incommensurate spin fluctuations to a gapped spin wave rather than a hourglass-like dispersion.

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.

preprint2009arXiv

The interface between the stellar wind and interstellar medium around R Cassiopeiae revealed by far-infrared imaging

The circumstellar dust shells of intermediate initial-mass (about 1 to 8 solar masses) evolved stars are generated by copious mass loss during the asymptotic giant branch phase. The density structure of their circumstellar shell is the direct evidence of mass loss processes, from which we can investigate the nature of mass loss. We used the AKARI Infrared Astronomy Satellite and the Spitzer Space Telescope to obtain the surface brightness maps of an evolved star R Cas at far-infrared wavelengths, since the temperature of dust decreases as the distance from the star increases and one needs to probe dust at lower temperatures, i.e., at longer wavelengths. The observed shell structure and the star&#39;s known proper motion suggest that the structure represents the interface regions between the dusty wind and the interstellar medium. The deconvolved structures are fitted with the analytic bow shock structure to determine the inclination angle of the bow shock cone. Our data show that (1) the bow shock cone of 1 - 5 x 10^-5 solar masses (dust mass) is inclined at 68 degrees with respect to the plane of the sky, and (2) the dust temperature in the bow shock cone is raised to more than 20 K by collisional shock interaction in addition to the ambient interstellar radiation field. By comparison between the apex vector of the bow shock and space motion vector of the star we infer that there is a flow of interstellar medium local to R Cas whose flow velocity is at least 55.6 km/s, consistent with an environment conducive to dust heating by shock interactions.

preprint2009arXiv

Unusual Dust Emission from Planetary Nebulae in the Magellanic Clouds

We present a Spitzer Space Telescope spectroscopic study of a sample of 25 planetary nebulae in the Magellanic Clouds. The low-resolution modules are used to analyze the dust features present in the infrared spectra. This study complements a previous work by the same authors where the same sample was analyzed in terms of neon and sulfur abundances. Over half of the objects (14) show emission of polycyclic aromatic hydrocarbons, typical of carbon-rich dust environments. We compare the hydrocarbon emission in our objects to those of Galactic HII regions and planetary nebulae, and LMC/SMC HII regions. Amorphous silicates are seen in just two objects, enforcing the now well-known-fact that oxygen-rich dust is less common at low metallicities. Besides these common features, some planetary nebulae show very unusual dust. Nine objects show a strong silicon carbide feature at 11um and twelve of them show magnesium sulfide emission starting at 25um. The high percentage of spectra with silicon carbide in the Magellanic Clouds is not common. Two objects show a broad band which may be attributed to hydrogenated amorphous carbon and weak low-excitation atomic lines. It is likely that these nebulae are very young. The spectra of the remaining eight nebulae are dominated by the emission of fine-structure lines with a weak continuum due to thermal emission of dust, although in a few cases the S/N in the spectra is low, and weak dust features may not have been detected.

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.

preprint2006arXiv

Detection of a Far-Infrared Bow-Shock Nebula Around R Hya: the First MIRIAD Results

We present the first results of the MIRIAD (MIPS [Multiband Imaging Photometer for Spitzer] Infra-Red Imaging of AGB [asymptotic giant branch] Dustshells) project using the Spitzer Space Telescope. The primary aim of the project is to probe the material distribution in the extended circumstellar envelopes (CSE) of evolved stars and recover the fossil record of their mass loss history. Hence, we must map the whole of the CSEs plus the surrounding sky for background subtraction, while avoiding the central star that is brighter than the detector saturation limit. With our unique mapping strategy, we have achieved better than one MJy/sr sensitivity in three hours of integration and successfully detected a faint (< 5 MJy/sr), extended (~400 arcsec) far-infrared nebula around the AGB star R Hya. Based on the parabolic structure of the nebula, the direction of the space motion of the star with respect to the nebula shape, and the presence of extended H alpha emission co-spatial to the nebula, we suggest that the detected far-IR nebula is due to a bow shock at the interface of the interstellar medium and the AGB wind of this moving star. This is the first detection of the stellar-wind bow-shock interaction for an AGB star and exemplifies the potential of Spitzer as a tool to examine the detailed structure of extended far-IR nebulae around bright central sources. \

preprint2006arXiv

Mid-infrared spectroscopy of carbon stars in the Small Magellanic Cloud

We have observed a sample of 36 objects in the Small Magellanic Cloud (SMC) with the Infrared Spectrometer on the Spitzer Space Telescope. Nineteen of these sources are carbon stars. An examination of the near- and mid-infrared photometry shows that the carbon-rich and oxygen-rich dust sources follow two easily separated sequences. A comparison of the spectra of the 19 carbon stars in the SMC to spectra from the Infrared Space Observatory (ISO) of carbon stars in the Galaxy reveals significant differences. The absorption bands at 7.5 um and 13.7 um due to C2H2 are stronger in the SMC sample, and the SiC dust emission feature at 11.3 um is weaker. Our measurements of the MgS dust emission feature at 26-30 um are less conclusive, but this feature appears to be weaker in the SMC sample as well. All of these results are consistent with the lower metallicity in the SMC. The lower abundance of SiC grains in the SMC may result in less efficient carbon-rich dust production, which could explain the excess C2H2 gas seen in the spectra. The sources in the SMC with the strongest SiC dust emission tend to have redder infrared colors than the other sources in the sample, which implies more amorphous carbon, and they also tend to show stronger MgS dust emission. The weakest SiC emission features tend to be shifted to the blue; these spectra may arise from low-density shells with large SiC grains.

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.

preprint2000arXiv

Infrared Properties of SiO Maser Sources in Late-Type Stars

Results of an SiO maser survey for the late-type stars selected by the IRTS (Infrared Telescope in Space) are presented. We have detected SiO J=1-0, v=1 and/or v=2 lines in 27 stars out of 59 stars. The maser intensity increases with the depth of the H2O absorption in the infrared spectra and redness of the 2.2 and 12 micron color. The column densities of the water vapor in the target stars are estimated from the depth of the water absorption in the IRTS spectra. We found that the SiO maser was detected mostly in the stars with the column density of water vapor higher than 3x10^19-3x10^20 cm^-2. We further estimate the density of hydrogen molecules in the outer atmosphere corresponding to these column densities, obtaining 10^9 - 10^10 cm^-3 as a lower limit. These values are roughly in agreement with the critical hydrogen density predicted by models for the excitation of the SiO masers. It is possible that the SiO masers are excited in clumps with even higher than this density. The present results provide useful information on the understanding of the physical conditions of the outer atmospheres in late-type stars.