Showing papers by "Pierre Royer published in 2014"

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

Abstract: 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.

Herschel view of the large-scale structure in the Chamaeleon dark clouds

Abstract: Context. The Chamaeleon molecular cloud complex is one of the nearest star-forming sites and encompasses three molecular clouds (Cha I, II, and III) that have a different star-formation history, from quiescent (Cha III) to actively forming stars (Cha II), and one that reaches the end of star-formation (Cha I).Aims. We aim at characterising the large-scale structure of the three sub-regions of the Chamaeleon molecular cloud complex by analysing new far-infrared images taken with the Herschel Space Observatory.Methods. We derived column density and temperature maps using PACS and SPIRE observations from the Herschel Gould Belt Survey and applied several tools, such as filament tracing, power-spectra, Δ-variance, and probability distribution functions (PDFs) of the column density, to derive the physical properties.Results. The column density maps reveal a different morphological appearance for each of the three clouds, with a ridge-like structure for Cha I, a clump-dominated regime for Cha II, and an intricate filamentary network for Cha III. The filament width is measured to be about 0.12 ± 0.04 pc in the three clouds, and the filaments are found to be gravitationally unstable in Cha I and II, but mostly subcritical in Cha III. Faint filaments (striations ) are prominent in Cha I and are mostly aligned with the large-scale magnetic field. The PDFs of all regions show a lognormal distribution at low column densities. For higher densities, the PDF of Cha I shows a turnover indicative of an extended higher density component and culminates in a power-law tail. Cha II shows a power-law tail with a slope characteristic of gravity. The PDF of Cha III can be best fit by a single lognormal.Conclusions. The turbulence properties of the three regions are found to be similar, pointing towards a scenario where the clouds are impacted by large-scale processes. The magnetic field might possibly play an important role for the star formation efficiency in the Chamaeleon clouds if it can be proven that it can effectively channel material onto Cha I, and possibly Cha II, but this is probably less efficient on the quiescent Cha III cloud.

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

Abstract: Context. 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. Aims. We characterize 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 on board the Herschel Space Observatory and ground-based telescopes. 12CO and 13CO lines are used to constrain the intrinsic 12C/13C ratio from resolved HIFI lines. Methods. We combined a state-of-the-art molecular line emission code and a dust continuum radiative transfer code to model the CO lines and the thermal dust continuum. Results. The acceleration of the outflow up to about 5.5 km s-1 is quite slow and can be represented by a β-type velocity law with index β = 5. Beyond this point, acceleration up the terminal velocity of 7 km s-1 is faster. Using the J = 10–9, 9–8, and 6–5 transitions, we find an intrinsic 12C/13C ratio of 18 ± 10 for W Hya, where the error bar is mostly due to uncertainties in the 12CO abundance and the stellar flux around 4.6 μm. To match the low-excitation CO lines, these molecules need to be photo-dissociated at ~500 stellar radii. The radial dust emission intensity profile of our stellar wind model matches PACS images at 70 μm out to 20′′ (or 800 stellar radii). For larger radii the observed emission is substantially stronger than our model predicts, indicating that at these locations there is extra material present. Conclusions. The initial slow acceleration of the wind may imply inefficient dust formation or dust 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⊙ or more. However, the uncertainty in the isotopologic ratio is large, which makes it difficult to set reliable limits on W Hya’s main-sequence mass.

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

Abstract: 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 $\mu$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 $\mu$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.

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

Abstract: Context. 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. Aims. 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 including an extensive set of H 2 Oa nd 28 SiO lines. It is the first time such a comprehensive study is performed for this source. 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. Methods. We model the molecular emission observed by the three instruments on board Herschel Space Observatory using a stateof-the-art molecular excitation and radiative transfer code. We also account for the dust radiation field in our calculations. Results. We fi nd an H 2O 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� . Although the ortho- and para-H2O lines observed by HIFI appear to trace gas of slightly different physical properties, we find that a turbulence velocity of 0.7 ± 0. 1k m s −1 fits the HIFI lines of both spin isomers and those of 28 SiO well. Conclusions. The modelling of H2 Oa nd 28 SiO confirms the properties of the envelope model of W Hya, as derived from 12 CO lines, and allows us to constrain the turbulence velocity. The ortho- and para-H 16 Oa nd 28 SiO abundances relative to H2 are (6 +3

Detailed modelling of the circumstellar molecular line emission of the S-type AGB star W Aquilae

Abstract: Context. S-type AGB stars have a C/O ratio which suggests that they are transition objects between oxygen-rich M-type stars and carbon-rich C-type stars. As such, their circumstellar compositions of gas and dust are thought to be sensitive to their precise C/O ratio, and it is therefore of particular interest to examine their circumstellar properties. Aims. We present new Herschel HIFI and PACS sub-millimetre and far-infrared line observations of several molecular species towards the S-type AGB star W Aql. We use these observations, which probe a wide range of gas temperatures, to constrain the circumstellar properties of W Aql, including mass-loss rate and molecular abundances. Methods. We used radiative transfer codes to model the circumstellar dust and molecular line emission to determine circumstellar properties and molecular abundances. We assumed a spherically symmetric envelope formed by a constant mass-loss rate driven by an accelerating wind. Our model includes fully integrated H2O line cooling as part of the solution of the energy balance. Results. We detect circumstellar molecular lines from CO, H2O, SiO, HCN, and, for the first time in an S-type AGB star, NH3. The radiative transfer calculations result in an estimated mass-loss rate for W Aql of 4.0 x 10(-6) M-circle dot yr(-1) based on the (CO)-C-12 lines. The estimated (CO)-C-12/(CO)-C-13 ratio is 29, which is in line with ratios previously derived for S-type AGB stars. We find an H2O abundance of 1.5 x 10(-5), which is intermediate to the abundances expected for M and C stars, and an ortho/para ratio for H2O that is consistent with formation at warm temperatures. We find an HCN abundance of 3 x 10(-6), and, although no CN lines are detected using HIFI, we are able to put some constraints on the abundance, 6 x 10(-6), and distribution of CN in W Aql's circumstellar envelope using ground-based data. We find an SiO abundance of 3 x 10(-6), and an NH3 abundance of 1.7 x 10(-5), confined to a small envelope. If we include uncertainties in the adopted circumstellar model - in the adopted abundance distributions, etc. - the errors in the abundances are of the order of factors of a few. The data also suggest that, in terms of HCN, S-type and M-type AGB stars are similar, and in terms of H2O, S-type AGB stars are more like C-type than M-type AGB stars. We detect excess blue-shifted emission in several molecular lines, possibly due to an asymmetric outflow. Conclusions. The estimated abundances of circumstellar HCN, SiO and H2O place W Aql in between M-and C-type AGB stars, i.e., the abundances are consistent with an S-type classification.

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

Abstract: 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.

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

Abstract: Context. We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67 72 m. This wavelength range covers the 69 m band of crystalline olivine (Mg 2 2x Fe(2x)SiO4). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, ˙ M 10 5 M /yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. Aims. The goal of this study is to exploit the spectral properties of the 69 m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. Methods. We fit the 69 m band and use its width and position to probe the composition and temperature of the crystalline olivine. Results. For 27 sources in the sample, we detected the 69 m band of crystalline olivine (Mg(2 2x)Fe(2x)SiO4). The 69 m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. The disk sources in our sample with crystalline olivine are very diverse. They show either no 69 m band, a moderately strong band, or a very strong band, together with a temperature for the crystalline olivine in their disk that is either very warm ( 600 K), moderately warm ( 200 K), or cold ( 120 K), respectively.

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

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

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

Abstract: Context. The Helix nebula (NGC 7293) is our closest planetary nebulae. Therefore, it is an ideal template for photochemical studies at small spatial scales in planetary nebulae. Aims. 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. Methods. We used five SPIRE FTS pointing observations to make atomic and molecular spectral maps. We analyzed the molecular gas by modeling the CO rotational lines using a non-local thermodynamic equilibrium (non-LTE) radiative transfer model. Results. 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), [N ii] at 1461 GHz from ionized gas, and [C i] (3P2–3P1), which together with the OH+ lines trace extended CO photodissociation regions along the rims. The estimated OH+ column density is ~ 1012 − 1013 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) × 105 cm-3) and the gas temperature (Tk ~ 20−40 K). Conclusions. The SPIRE spectral-maps suggest that CO arises from dense and shielded clumps in the western rims of the Helix nebula, whereas OH+ and [C i] lines trace the diffuse gas and the UV and X-ray illuminated clump surfaces where molecules reform after CO photodissociation. The [N ii] line traces a more diffuse ionized gas component in the interclump medium.

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

Abstract: 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.

Detailed modelling of the circumstellar molecular line emission of the S-type AGB star W Aquilae

Abstract: S-type AGB stars have a C/O ratio which suggests that they are transition objects between oxygen-rich M-type stars and carbon-rich C-type stars. As such, their circumstellar compositions of gas and dust are thought to be sensitive to their precise C/O ratio, and it is therefore of particular interest to examine their circumstellar properties. We present new Herschel HIFI and PACS sub-millimetre and far-infrared line observations of several molecular species towards the S-type AGB star W Aql. We use these observations, which probe a wide range of gas temperatures, to constrain the circumstellar properties of W Aql, including mass-loss rate and molecular abundances. We used radiative transfer codes to model the circumstellar dust and molecular line emission to determine circumstellar properties and molecular abundances. We assumed a spherically symmetric envelope formed by a constant mass-loss rate driven by an accelerating wind. Our model includes fully integrated H2O line cooling as part of the solution of the energy balance. We detect circumstellar molecular lines from CO, H2O, SiO, HCN, and, for the first time in an S-type AGB star, NH3. The radiative transfer calculations result in an estimated mass-loss rate for W Aql of 4.0e-6 Msol yr-1 based on the 12CO lines. The estimated 12CO/13CO ratio is 29, which is in line with ratios previously derived for S-type AGB stars. We find an H2O abundance of 1.5e-5, which is intermediate to the abundances expected for M and C stars, and an ortho/para ratio for H2O that is consistent with formation at warm temperatures. We find an HCN abundance of 3e-6, and, although no CN lines are detected using HIFI, we are able to put some constraints on the abundance, 6e-6, and distribution of CN in W Aql's circumstellar envelope using ground-based data. We find an SiO abundance of 3e-6, and an NH3 abundance of 1.7e-5, confined to a small envelope.