Free Access
Issue
A&A
Volume 570, October 2014
Article Number A67
Number of page(s) 18
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201424298
Published online 20 October 2014

© ESO, 2014

1. Introduction

Stars of low- and intermediate-mass (~0.8 to 8 M) populate the asymptotic giant branch (AGB) at the end of their lives. These luminous, extended, and cool objects develop a strong mass loss that controls their evolution from its onset. The driving of the wind is thought to be the result of a combination of pulsations and radiation pressure on dust grains (see e.g. Habing & Olofsson 2003). For carbon-rich AGB stars, this scenario is able to reproduce the observed mass-loss rates as the dust species that form in these environments are opaque enough and can exist close enough to the star to acquire momentum by absorbing infrared photons (Winters et al. 2000). However, for oxygen-rich AGB stars the situation seems more complex. The dust species found in these stars through infrared spectroscopy that would be able to drive the wind by absorption of stellar photons cannot exist close enough to the star to be important for initiating the wind (Woitke 2006). A possible alternative is that the wind is driven through scattering of photons on large (~ 0.3 μm) translucent dust grains (Höfner 2008). Although this alternative scenario seems plausible and many dust species have been identified in oxygen-rich objects, it is still unclear which of these are actually responsible for driving the outflow (e.g. Bladh & Höfner 2012).

Characterizing both the physical structure and molecular and solid state composition of the outflow are crucial to understanding the physics underlying the wind driving. An important step forward in this understanding is the combined analysis of space and ground-based observations of multiple molecular species, covering a range of rotational excitation states of the ground vibrational level as wide as possible. Such data allow the determination of the flow properties from the onset of the wind close to the photosphere, to the outer regions, where molecules are eventually photodissociated by the interstellar radiation field. We have embarked on such an analysis for W Hya, a close-by oxygen-rich AGB star, observed in detail using the Herschel Space Observatory (Pilbratt et al. 2010). In addition to the 12CO lines, the observations reveal a rich HO spectrum with over 150 observed lines, and a broad 28SiO ladder, ranging from 137 to more than 2000 K in upper level excitation energy.

In Khouri et al. (2014, henceforth Paper I), we focussed on the analysis of the 12CO ladder of rotational levels up to Jup = 30. We used these carbon monoxide lines to determine the temperature and velocity structure of the wind, as well as the mass-loss rate. In the study presented here we focus on an analysis of lines from H2O, including its isotopologues, and 28SiO. The modelling of the H2O isotopologues allows us to constrain the 17O/16O and 18O/16O ratios. The analysis of silicon in both the gas phase and solid phase allows us to assess the overall budget of this element and the condensation fraction of this species. This is a crucial issue for understanding the role of silicates for the wind driving mechanism. We also account for the contribution of the dust thermal emission to the radiation field when modelling the molecular emission. The main components of these grains are aluminum-oxides and silicates (e.g. Sharp & Huebner 1990).

Determining the ortho- and para-H2O abundances is important for understanding the chemistry in the outermost layers of the star and the innermost regions of the envelope, where shocks can be important (Cherchneff 2006). For an AGB star with carbon-to-oxygen ratio of 0.75, Cherchneff (2006) finds the abundance ratio between 12CO and HO to be ~12 and ~1.6 for thermal equilibrium and non-equilibrium chemistry respectively. The H2O ortho-to-para ratio is expected to be 3:1, reflecting the ratio of the statistical weights between the species, if the molecules are formed in a high-temperature (T ≫ 30 K) and under local-thermodynamical equilibrium (Decin et al. 2010b). Studies of the HO emission from W Hya have found a low ortho-to-para ratio of around 1 (e.g. Barlow et al. 1996; Justtanont et al. 2005; Zubko & Elitzur 2000). The uncertainties on these obtained ortho-to-para ratios are, however, at least of a factor of two.

The isotopic composition of the outflowing material contains much-needed information on the dredge-up processes that are an important part of giant branch and AGB evolution (e.g. Landre et al. 1990; El Eid 1994; Stoesz & Herwig 2003; Charbonnel & Lagarde 2010; Palmerini et al. 2011). Dredge-ups enrich the surface with isotopic species when convective streams in the star reach down to regions where the composition has been modified because of thermonuclear burning (see e.g. Iben & Renzini 1983; Iben 1975). The characteristics of the isotopic enrichment are expected to vary significantly over the evolution of low- and intermediate-mass stars and are found to be especially sensitive to stellar mass (Boothroyd et al. 1994) – the most important stellar property, which for AGB stars is notoriously difficult to constrain. Unfortunately, the dredge-up process in AGB stars cannot yet be modelled from first principles (Habing & Olofsson 2003), mainly because of the complex (and poorly understood) physics of convective and non-convective mixing (see e.g. Busso et al. 1999; Karakas 2010).

In Sect. 2 we provide general information on W Hya, the 12CO model of Paper I, and the dataset that is used to constrain the HO and 28SiO properties. In Sect. 3, we discuss observed line shapes and provide details on the treatment of HO and 28SiO in our models, and how this treatment relates to that of 12CO. Section 4 is devoted to presenting our model for HO and for the lower abundance isotopologues. The model for 28SiO is covered in Sect. 5. We discuss the results in Sect. 6 and we end with a summary.

2. Basis information, dataset, and model assumptions

2.1. W Hya

An overview of the literature discussing the stellar properties of W Hya, and the uncertainties in these properties, is presented in Paper I. We refer to this paper for details. We adopt a distance to the star of 78 pc (Knapp et al. 2003), implying a luminosity of 5400 L. Assuming that W Hya radiates as a black-body, our 12CO analysis is consistent with an effective temperature of 2500 K, which leads to a stellar radius, R, of 2.93 × 1013 cm or 1.96 AU. Whenever we provide the radial distance in the wind of W Hya in units of R, this is the value we refer to. The 12CO analysis presented in Paper I leads to a mass-loss rate of 1.3 × 10-7M yr-1, consistent with the findings of Justtanont et al. (2005). Interestingly, on a scale that is larger than the 12CO envelope, images of cold dust emission suggest W Hya had a substantially larger mass loss some 103105 years ago (Cox et al. 2012; Hawkins 1990). In this work, we focus on the present-day mass-loss rate by modelling the gas-phase 28SiO and H2O emission, which probe the last 200 years at most. Variations in the mass-loss rate seen on larger timescales will be addressed in a future study.

W Hya features prominent rotational HO emission, first reported by Neufeld et al. (1996) and Barlow et al. (1996), using data obtained with the Infrared Space Observatory (ISO; Kessler et al. 1996). Additional data were obtained by Justtanont et al. (2005), using Odin (Nordh et al. 2003), and by Harwit & Bergin (2002), using SWAS (Melnick et al. 2000). The analysis by the authors mentioned above and by others (Zubko & Elitzur 2000; Maercker et al. 2008, 2009) point to a quite high HO abundance relative to H2, ranging from 10-4 to a few times 10-3. The ortho-to-para H2O ratio reported by most studies is usually in-between 1 and 1.5, which is significantly lower than the expected value of three, for H2O formed at high temperatures (T ≫ 30 K) and in local-thermodynamical equilibrium.

28SiO lines of the vibrational excited v = 1 and v = 2 state show maser emission and have been intensively studied (see e.g. Imai et al. 2010; Vlemmings et al. 2011). Ground vibrational lines do not suffer from strong amplification and are better probes of the silicon abundance. Studies of this molecule (e.g. González Delgado et al. 2003; Bieging et al. 2000; Bujarrabal et al. 1986) suggest 28SiO is depleted from the gas, probably due to the formation of silicate grains. Spatially resolved data of the v = 0,J = 2 − 1 line have been presented by Lucas et al. (1992). The authors determined the half-intensity radius of this transition to be 0.9 arcsec. From model calculations, González Delgado et al. (2003) conclude that the half-intensity radius determined by Lucas et al. (1992) is approximately three times smaller than the radius where the 28SiO abundance has decreased to 37% of its initial value.

The dust envelope of W Hya has been imaged by Norris et al. (2012) using aperture-masking polarimetric interferometry. They discovered a close-in shell of large (~0.3 μm) and translucent grains which might be responsible for driving the outflow through scattering (Höfner 2008). Unfortunately, the composition of the observed grains could, not be unambiguously identified. W Hya has also been observed in the near-infrared using MIDI/VLTI (Zhao-Geisler et al. 2011). This revealed that the silicate dust emission must come from an envelope with an inner radius of 50 AU (or 28 R). Zhao-Geisler et al. (2011) also argue that Al2O3 grains and H2O molecules close to the star are responsible for the observed increase in diameter at wavelengths longer than 10 μm.

2.2. Dataset

W Hya was observed by all three instruments onboard Herschel in the context of the guaranteed-time key programmes HIFISTARS (Menten et al. 2010) and MESS (Groenewegen et al. 2011). These are the Heterodyne Instrument for the Far Infrared (HIFI; de Graauw et al. 2010), the Spectral and Photometric Imaging Receiver Fourier-Transform Spectrometer (SPIRE; Griffin et al. 2010), and the Photodetector Array Camera and Spectrometer (PACS; Poglitsch et al. 2010). The data reduction procedure of the PACS and SPIRE data is outlined in Paper I ; that of HIFI is presented by Justtanont et al. (2012).

We applied two methods to identify the HO and 28SiO lines. The 28SiO lines were identified in very much the same way as the 12CO lines presented in Paper I: we inspected the spectra at the wavelengths of the 28SiO transitions, identified the lines and measured their fluxes. HO, however, has a plethora of transitions that, moreover, are not regularly spaced in frequency, as are the ones of 12CO and 28SiO. Our approach in this case was to calculate spectra with different values for the HO-envelope parameters and compare these to the observations. Transitions that were predicted to stand out from the noise were identified and extracted.

In order to extract the integrated fluxes from the PACS spectrum, we fitted Gaussians to the identified transitions using version 11.0.1 of Herschel interactive processing environment (HIPE1). For extracting transitions measured by SPIRE, we used the script Spectrometer Line Fitting available in HIPE. The script simultaneously fits a power law to the continuum and a cardinal sine function to the lines in the unapodized SPIRE spectrum. The spectral resolution of both instruments is smaller than twice the terminal velocity of W Hya (υ = 7.5 km s-1; see Table 1). Therefore, what may appear to be a single observed line might be a blend of two or more transitions. We removed such blended lines from our analysis by flagging them as blends whenever two or more transitions were predicted to be formed closer together than the native full width at half-maximum (FWHM) of the instrument, or when the FWHM of the fitted Gaussian was 20% or more larger than the expected FWHM for single lines. For PACS, a total of 50 ortho-HO, 24 para-HO, and only 1 28SiO unblended transitions were extracted. The 28SiO transitions that lie in the spectral region observed by PACS are high excitation lines (Jup> 35). These are weak compared to the much stronger HO lines and are below the detection limit of the PACS spectrum for Jup> 38. For SPIRE, 15 ortho-HO, 20 para-HO and 24 28SiO lines were extracted and not flagged as blends. Properties of all HO transitions measured by PACS and those of 28SiO measured by all instruments are given in Appendix B.

The transitions observed with HIFI are spectrally resolved. This gives valuable information on the velocity structure of the flow. Because of the high spectral resolution the lines are easily identified and their total fluxes can be measured accurately. In Table 2, we list the integrated main beam brightness temperatures and the excitation energy of the upper level for the observed transitions. HIFI detected ten ortho-HO transitions, two of which are clearly masering, six para-HO transitions, one of which is clearly masering, five 28SiO transitions, three of which are from vibrationally excited states and appear to be masering, two transitions of ortho-HO and 29SiO, and one transition each of ortho-HO, para-HO, para-HO and 30SiO. We have not included the rarer isotopologues of SiO, 29SiO and 30SiO, in our model calculations. The vibrationally excited lines from 28SiO probe high temperature gas (T ~ 2000 K) which is very close to the star and were not included in our analysis. The upper level excitation energies of the 28SiO ground-vibrational lines range from 137 to 1462 K and these transitions probe the gas temperatures in which silicates are expected to condense, T ~ 1000 K (e.g. Gail & Sedlmayr 1999).

Table 1

Model parameters of W Hya as derived in Paper I.

2.2.1. HO masers

The already challenging task of modelling HO emission is further complicated by the fact that this molecule has the predisposition to produce maser emission. Matsuura et al. (2013) found that for the oxygen-rich supergiant VY CMa 70% of the HO lines in the PACS and SPIRE spectra are affected by population inversions. In our line emission code we artificially suppress strong stimulated emission and we have therefore excluded from the present analysis lines that show masering in at least one of our models. In the PACS data, four transitions were excluded. Those are the ones at 133.549, 159.051 and 174.626 μm of ortho-HO and 138.528 μm of para-HO. The intensity of the ortho-HO transition at 73.415 was under-predicted by more than a factor of three by all models. Although this line is not masering in any of our models, the discrepancies are so large that we suspect missing physics to be responsible. We have also excluded this line from the analysis. In four of the remaining PACS lines, we do find modest population inversions for some of the models. These, however, do not appear to have a strong impact on the strength of the lines and the transitions were not excluded. For the SPIRE lines, the case is worse. About 70 % of the lines – in agreement with the results reported for VY CMA by Matsuura et al. (2013) – present population inversions. The intensity of most of these lines was strongly under-predicted by our models. This is a much higher fraction than the less than 10% seen for PACS. To be on the safe side, we decided not to include SPIRE data at all in our calculations for HO. As a result of this mismatch a comprehensive identification of the HO lines in the SPIRE spectrum was not possible and we do not present extracted line fluxes for HO lines obtained by this instrument.

Table 2

H2O and 28SiO transitions observed by HIFI for W Hya.

2.2.2. Comparison between the total line fluxes measured by HIFI, SPIRE and PACS

Some of the lines were observed by more than one instrument, either by HIFI and SPIRE or by HIFI and PACS. In Table 3 we present a comparison of the observed total line fluxes. The antennae temperatures observed by HIFI were converted to flux units assuming that the emission is not spatially resolved by the telescope. We find that the fluxes of HO lines as detected by HIFI are systematically lower than those by both SPIRE and PACS.

SPIRE lines are typically 25% stronger than those of HIFI. This difference may be due to uncertainties in the flux calibration of both instruments and in the baseline fitting of the spectra. Small blends and pointing errors can also account for some of the observed differences. A flux loss of less than 10% is expected from the pointing errors in the HIFI observations. The different shapes of the response function of the instruments may also have a significant effect.

The three lines measured jointly by PACS and HIFI are more than 2–3 times stronger in the PACS traces than in the HIFI spectra. The cross-calibration problem of PACS and HIFI is known, but is not yet resolved. It is often reported as a flux mismatch between the 12CO lines observed by both instruments (Puga et al., in prep.). For W Hya we find that the mismatch between predicted and PACS line fluxes correlates with wavelength, becoming larger for longer wavelengths. At the long wavelength end of the PACS spectral range (which is where PACS and HIFI spectra overlap) the discrepancy strongly increases. Since the HIFI observations are much less susceptible to line blending and the HIFI line fluxes in this overlap region agree much better with the trends predicted by our models than the line fluxes measured by PACS, we decided to exclude these PACS lines from our analysis, and to use the HIFI lines instead.

Table 3

Observed total line fluxes of HO and 28SiO transitions observed by two different Herschel instruments, HIFI and SPIRE or PACS.

2.3. Modelling strategy and 12CO model

Here, we focus on modelling the emission of ortho-HO, para-HO, ortho-HO, para-HO, ortho-HO, para-HO and 28SiO. The molecular data used in our calculations are described in Appendix A. Our model is based on the envelope structure and dust model obtained in Paper I, which was obtained from the analysis of 12CO lines, observed with Herschel, APEX, SMT and SEST, and of the dust emission. The 12CO transitions modelled in Paper I probe a large range in excitation temperature and, therefore, the full extent of the outflowing 12CO envelope, from the regions close to the star where the wind is accelerated, to the outer regions where 12CO is photo-dissociated. The dust properties were constrained by modelling the thermal emission spectrum observed by ISO with the continuum radiative-transfer code MCMax (Min et al. 2009). The composition and radial distribution of the dust are used as inputs to the modelling of the molecular species, for which we employ GASTRoNOoM (Decin et al. 2006, 2010a). The coupling between dust and gas is treated as described by Lombaert et al. (2013). The main parameters of the 12CO model are listed in Table 1.

In Paper I, we found that the 12CO envelope has to be smaller than predicted for our model to better fit the high- (Jup ≥ 6) and low-excitation (Jup< 6) 12CO transitions simultaneously. The size of the 12CO envelope is set in our model by the parameter R1/2 (Mamon et al. 1988), which represents the radius at which the 12CO abundance has decreased by half. We note that this discrepancy between model and observations in the outer wind is not expected to affect the modelled HO or 28SiO lines, since these molecules occupy a much smaller part of the envelope than 12CO, see Fig. 1.

Furthermore, to properly represent the PACS and HIFI transitions excited in the inner part of the wind, our best model requires a value of 5.0 for the exponent of the β-type velocity law,(1)This corresponds to a slow acceleration of the flow in this part of the envelope. However, this high value of β underpredicts the width of the 12CO J = 6 − 5 transition, while predicting well the width of lower excitation transitions (Jup< 6). This indicates a more rapid acceleration between the formation region of transitions J = 10 − 9 and J = 6 − 5 than considered in our β = 5.0 model. This could be due to the addition of extra opacity in the wind at distances beyond where the J = 10 − 9 transition forms, which is roughly inside 50 R. As the HO and 28SiO envelopes are smaller than the excitation region of the 12CO Jup = 6 level, emission from these two molecules comes from a wind region that is well described by a β = 5.0 velocity law (see also Fig. 1).

thumbnail Fig. 1

Normalized abundance profile of HO (long-dashed green), 28SiO (double-short-dashed yellow) and 12CO (solid red) are compared to the normalized populations of levels J = 16 (dashed-dotted brown), J = 10 (dotted purple) and J = 6 (short-dashed blue) of 12CO.

Open with DEXTER

Our dust model is motivated by the work of Justtanont et al. (2004, 2005) and adopts a dust mass-loss rate of 2.8 × 10-10M yr-1. Astronomical silicates, amorphous aluminum oxide (Al2O3), and magnesium-iron oxide (MgFeO) account for 58, 34, and 8% of the dust mass, respectively. The optical constants for astronomical silicates are from Justtanont & Tielens (1992), those for amorphous aluminum oxide and magnesium-iron oxide were retrieved from the University of Jena database and are from the works of Begemann et al. (1997) and Henning et al. (1995).

3. Comparisons between 12CO, HO and 28SiO

3.1. Photodissociation radii

12CO, 28SiO, and HO are formed in the atmosphere of the AGB star and the fate of these molecules is to become dissociated in the outer envelope by interstellar UV photons, causing their abundances to decrease sharply. The radius where dissociation sets in is different for each molecule and depends on molecular and circumstellar parameters and on the spectral shape of the interstellar radiation field.

Dust condensation might play a role in shaping the abundance profile of a molecule if the conditions for condensation are met before dissociation sets in. For W Hya we expect 28SiO to condense and form silicate grains. 12CO and HO abundances should not be affected by such depletion and we assume their abundances to be constant up to the point where dissociation starts. For 28SiO we consider different abundance profiles, mimicking depletion due to silicate formation (see Sect. 5).

Dissociation radii are usually poorly constrained. In our 12CO model presented in Paper I, we find that for W Hya the 12CO envelope is likely smaller than expected for a standard interstellar radiation field. Our results point to an envelope with 12CO being fully dissociated at roughly 800 R. We have assumed dissociation profiles for HO and 28SiO of the type (2)where the standard value of the exponent p is 2 and re is the e-folding radius, the radius at which the abundance has decreased by a factor e from its initial value. HO is expected to dissociate closer to the star than 12CO (e.g. Netzer & Knapp 1987). Groenewegen (1994) argues that HO emission should come from within roughly 100 R, from considering studies of OH density profiles in AGB stars. This limit corresponds to an abundance profile with e-folding radius of 65 R, or 1.8 × 1015 cm. For 28SiO, there are no theoretical estimates that give its dissociation radius in terms of envelope parameters. That said, González Delgado et al. (2003) modelled low-excitation 28SiO transitions and determined the e-folding radius to be re = 2.4 × 1015 cm for W Hya, which corresponds to roughly 85 R in the context of our model. At about 200 R all of the 28SiO will then have disappeared. We adopt these values. They show that the 28SiO envelope is comparable in size to the H2O envelope and that both are considerably smaller than the CO envelope. Because the 28SiO transitions modelled by us are formed mostly deep inside the 28SiO envelope, we expect the problem of constraining the 28SiO dissociation radius to have little impact on our results.

In Fig. 1, we plot the normalized abundance profiles of 12CO, 28SiO and HO compared to the excitation region of the 12CO transitions observed by HIFI. This shows that the HO and 28SiO emissions probe a relatively small region of the envelope compared to 12CO emission.

thumbnail Fig. 2

Observed line shapes of the transitions J = 16 − 15 (dotted purple) and J = 10 − 9 (solid red) of 12CO, J = 14 − 13 of 28SiO (long-dashed green), and 11,0-00,1 of ortho-HO (short-dashed blue). The vertical line marks the adopted υLSR = 40.4 km s-1.

Open with DEXTER

3.2. Observed line shapes: blue wing absorption

thumbnail Fig. 3

Effects of varying the HO abundance, the turbulence velocity, the temperature power law exponent and the velocity power law exponent on the normalized line strength are shown for o-HO (11,0 − 10,1). The standard model is represented by the red line in each panel, it has its values given in Table 1 and fH2O = 4 × 10-4. The short-dashed purple and the long-dashed green represent models with respectively: fH2O = 1 × 10-4 and 1.6 × 10-3 in the first panel; υturb = 0.8 and 2.0 km s-1 in the second panel; ϵ = 0.2 and 1.1 in the third panel; and β = 1.5 and 10 in the fourth panel. All other parameters have the same value as the standard model. We note that varying υturb has the effect of a line shift.

Open with DEXTER

Figure 2 shows the spectrally resolved profiles of the lowest excitation transitions of HO and 28SiO, as observed by HIFI, and of the J = 16 − 15 and 10 − 9 transitions of 12CO. The adopted local-standard-of-rest velocity of 40.4 km s-1 was derived by Khouri et al. (2014) by modelling the 12CO and 13CO transitions. The plotted HO and 28SiO lines are expected to form in about the same region of the envelope as the pure rotational J = 10 − 9 line of 12CO, as they are low-excitation transitions expected to be formed in the outer parts of the 28SiO and HO envelopes. The 12CO J = 10 − 9 transition shows emission up to the same expansion velocity in the red- and blue-shifted wings. The HO and 28SiO transitions, however, show very asymmetrical line profiles with emission extending to larger velocities in the red-shifted wing than on the blue-shifted wing. To show the effect of differences in line forming region (which cannot be the cause of this behaviour, as the lines are selected to form in about the same part of the flow), we also show the 12CO J = 16−15 line profile. This line forms more to the base of the outflow. The profile is narrower on both sides of line centre.

Huggins & Healy (1986) studied the line shapes of transitions J = 2−1 and J = 1−0 of 12CO in IRC +10216 in comparison to the line shape of the optically thinner 13CO J = 2−1 line. The authors observed asymmetries between the 12CO and 13CO lines and concluded that these were due to the higher optical depths in the 12CO lines. They found the effect to be strongly dependent on the turbulence velocity in the line formation region. In Fig. 2, the differences between H2O and 28SiO on the one hand and 12CO J = 10 − 9 on the other hand also arise as a result of differences in line optical depth, and the presence of a turbulence velocity field. The 12CO radial optical depth at line centre is about unity for the transitions discussed. The Einstein coefficient for spontaneous emission of the HO and 28SiO transitions modelled here are typically at least an order of magnitude larger than that of 12CO. The 28SiO J = 16 − 15 and 14 − 13 lines have even two orders of magnitude higher values. As a result, higher optical depths are built-up for 28SiO and HO transitions than for 12CO transitions, even when the different abundances are taken into account. The effect is important even in the wind of a low mass-loss rate AGB star such as W Hya, and it reveals itself by the shift in the observed line centre velocity with respect to the υLSR of the source.

In Fig. 3, we show that υturb indeed affects strongly the observed behaviour of HO and 28SiO. The o-HO (11,0 − 10,1) transition is shown for models based on the best 12CO model obtained in Paper I. We varied either the HO abundance, the turbulence velocity in the envelope, υturb, the exponent of the temperature power law, ϵ, or the exponent of the velocity power law, β. Although all these parameters have an impact on the HO line shapes or central peak position, it is the turbulence velocity that causes a systematic shift of the line centre position. We stress that varying these parameters does have a significant effect on the line strength. To bring out the effects on profile shape and profile shift, we have normalized the lines. As pointed out by Huggins & Healy (1986), the magnitude of the shift is mainly set by and can be used to determine the value of the turbulence velocity in the line formation region.

thumbnail Fig. 4

Normalized profiles of ortho- and para-HO transitions observed by HIFI compared to models with different values of the turbulence velocity and all other parameters kept fixed. The solid red line represents the data, and models with turbulence velocities of 0.5, 0.8 and 1.4 km s-1 are represented by dotted purple, long-dashed green and short-dashed blue lines, respectively. The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines.

Open with DEXTER

4. HO model

Table 4

Turbulence velocities and abundances relative to H2 considered in the models studying the effects of HO.

We calculated models using the parameters found in Paper I for different values of the ortho- and para-HO abundance and the turbulence velocity in the wind. The turbulence velocity was included as a free parameter since it was clear from the first calculated HO profiles that the originally assumed value (υturb = 1.4 km s-1; see Table 1) was too large to match the line-centre shifts seen in the line profiles obtained with HIFI (see Sect. 3.2). The HO abundance at a radius r is given by Eq. (2). The value adopted by us for the H2O dissociation radius is that given by Groenewegen (1994). Our model calculations, however, show that the derived HO abundance is not strongly affected by the assumed dissociation radius.

The turbulence velocity and the abundance of each spin isomer constitute a degeneracy. We can, however, determine with a good accuracy the value of the turbulence velocity from matching the line-centre shifts seen by HIFI (see Sect. 3.2). The parameters used in the calculation of the HO models are listed in Table 4. The line profiles computed for different values of the turbulence velocity are compared to the line shapes observed by HIFI in Fig. 4. The values of the turbulence velocity that best match the shift seen for ortho-HO and para-H2O are 0.8 and 0.6 km s-1, respectively. Although the difference in the derived turbulence velocity is small (~0.2 km s-1), the high-quality HIFI data seems to suggest that the main line-formation region for the o-H2O lines has a slightly different turbulence velocity than for p-H2O. However, as an uncertainty of 0.1 km s-1 in the turbulence velocity, propagates to an uncertainty of only 5% at most for the predicted line profiles, we adopt a value of 0.7 km s-1 for both spin isomers.

After fixing the turbulence velocity, we calculated the reduced-χ2 of the fits to the PACS and HIFI line fluxes for models with different values of the ortho-HO and para-HO abundances. For ortho-HO, we get the best result for an abundance of 6 × 10-4 for both HIFI and PACS data. For para-HO, a value of 3 × 10-4 fits better the PACS data, while a value between 4 × 10-4 and 3 × 10-4 is the best match for the HIFI observations. Taking the whole dataset into account, the best fits for ortho- and para-HO are achieved with abundances of and , respectively. These models are compared to the lines observed by HIFI in Figs. 5 and 6 for ortho- and para-HO, respectively.

thumbnail Fig. 5

Best fit model (dashed-blue line) to the ortho-HO lines compared to the lines observed by HIFI (solid-red line). The model has an ortho-HO abundance of 6 × 10-4 and a turbulence velocity of 0.7 km s-1. The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines.

Open with DEXTER

thumbnail Fig. 6

Best fit model (dashed-blue line) to the para-HO lines compared to the lines observed by HIFI (solid-red line). The model has a para-HO abundance of 3 × 10-4 and a turbulence velocity of 0.7 km s-1. The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines.

Open with DEXTER

4.1. H2O isotopologues

We calculated models for each of the four isotopologues of H2O, ortho-HO and HO and para-HO and HO, considering both a turbulence velocity of 0.7 km s-1 and different values for the abundance of each isotopologue. The optical depths of the transitions are considerably smaller than the corresponding ones of the main isotopologues, causing the derived abundances to be roughly independent of υturb. The uncertainty on υturb propagates to an uncertainty of a few percent in the line fluxes of the rarer isotopologues.

thumbnail Fig. 7

The solid red line represents the HIFI data, the dashed green, short-dashed blue, dotted purple, dashed-dotted light-blue and short-dashed-dotted brown lines represent, respectively, models with abundances relatively to H2 of: (4, 2, 1.2, 0.4, 0.2) × 10-6 for for ortho-HO; (1.2, 0.8, 0.4, 0.2, 0.12) × 10-4 for ortho-HO; (8, 4, 2, 1.2, 0.4) × 10-7 for para-HO; and (8, 4, 2, 1.2, 0.4) × 106 for para-HO. We adopted a value of 0.7 km s-1 for the turbulence velocity in these calculations. The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines.

Open with DEXTER

The models are compared to the observations in Fig. 7 and the best abundance values of each isotopologue relative to H2 are given in Table 5. The best models were selected by comparing the integrated line fluxes with the observations. Then, we computed the different isotopologic ratios given in Table 6, which allow us to derive the 16O/17O, the 16O/18O isotopic ratios and the ortho-to-para ratio. We note, however, that the line shapes of the two ortho-H2O isotopologues are not well reproduced, particularly that of o-HO. The single observed line of o-HO is shifted with respect to the modelled transitions and very high values of the o-H218O abundance are needed to predict the observed line strength. Therefore, the values derived for the 16O/18O ratio and for the H2O ortho-to-para ratio based on this line do not agree with those derived using any other pair of observed lines. Because of this mismatch, we have not included the values derived from the o-H218O line in the ortho-to-para and isotopic ratios calculations. Modifying the dissociation radius and turbulence velocity does not improve the fitting of this transition. Discarding the o-H218O line, we obtain: ortho-to-para = 2.5 , 16O/17O = 1250  and 16O/18O = 190  (see also Table 6).

Table 5

Best fit results for the abundances of the isotopologues of H2O.

Table 6

Derived isotopologic ratios.

4.2. Consequences of a lower υturb for the 12CO lines

The turbulence velocity that resulted from the 12CO analysis in Paper I was 1.4 km s-1. However, this value was not strongly constrained. Although this higher value is better at reproducing the 12CO line wings, the differences between models for 0.60.8 and 1.4 km s-1 are small. In Fig. 8 we show a comparison between the 12CO model from Paper I and a model with υturb = 0.7 km s-1. The change in the value of the turbulence velocity only has an impact on the total line flux of the transitions having Jup< 6, most notably that of J = 2 − 1. When a lower value of the turbulence velocity is used, the model predictions for the J = 1 − 0 and 2 − 1 transitions get somewhat stronger, while those for the 3 − 2 to 6 − 5 transitions get somewhat weaker. The emission from line J = 2 − 1 could be decreased by considering an even smaller 12CO dissociation radius than the one obtained in Paper I, but this will not resolve the discrepancy in the J = 1 − 0 line. The poor fit to the second line is the main shortcoming of our 12CO model. We conclude that a lower value for the turbulence velocity does not affect significantly the quality of our fit, except for the very low-excitation lines J = 2 − 1 and J = 1 − 0, the second of which was also poorly reproduced by our original 12CO model. Furthermore, we can expect the turbulence velocity to be different in the formation regions of HO and 28SiO and of the low-excitation 12CO lines, since the 12CO envelope is significantly larger than that of HO and 28SiO. However, since it is not possible to determine a precise value for the turbulence velocity from the 12CO lines, we cannot draw any conclusion on changes in this parameter between the HO and 28SiO envelopes and the 12CO outer envelope.

thumbnail Fig. 8

A model for the 12CO lines with υturb = 0.7 km s-1 (short-dashed blue line) is compared to the best model found in Paper I (υturb = 1.4 km s-1; dotted-dashed brown line). The observed 12CO lines are shown in red. The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines. The lines J = 16 − 15, 10 − 9 and 6 − 5 were observed with HIFI, J = 4 − 3 and 3 − 2 with APEX, J = 2 − 1 with SMT and J = 1 − 0 with SEST.

Open with DEXTER

thumbnail Fig. 9

PACS spectrum (represented by the grey-filled histogram) compared to our best model, with the parameters given in Table 1 but with υturb of 0.7 km s-1 and with ortho- and para-HO abundances of 6 × 10-4 and 3 × 10-4, respectively. The 12CO model is shown by the full-blue line, the ortho-HO, by the full-red line and the para-HO, by the full-yellow line.

Open with DEXTER

4.3. Reproducing the PACS spectrum

In Fig. 9, we compare our 12CO and ortho- and para-HO models to the PACS spectrum of W Hya. The 28SiO lines modelled do not contribute significantly to the PACS spectrum and were not included in the plot. The vast majority of the prominent lines seen in the spectrum can be accounted for by our HO model. A small fraction of the strong lines, however, is not predicted. These lines have peaks around 61.52, 72.84, 78.47, 79.12, 86.52, 89.78, 154.88 and 163.12 μm. The lines observed at 79.12 and 163.12 μm can be associated with OH transitions (Sylvester et al. 1997). Those at 78.47, 86.52 and 89.78 μm might be explained by 28SiO maser emission (Decin et al., in prep.). We were not able to identify the lines observed at 61.52, 72.84 and 154.88 μm.

5. Model for gas-phase 28SiO emission

In Paper I, we assumed a standard value for the 12CO abundance relative to H2 of 4 × 10-4. Were all carbon and silicon used to form 12CO and 28SiO, respectively, and solar composition assumed, then the 12CO-to-28SiO abundance ratio would be roughly 8.3, corresponding to an 28SiO photospheric abundance of 4.8 × 10-5.

In order to establish the location in the outflow where silicate dust particles condense, and the fraction of the 28SiO that is converted from the molecular phase to the solid phase in the inner wind, we modelled the 28SiO line emission detected by Herschel. The broad spectral coverage of PACS and SPIRE provides a series of lines with upper-level energies ranging from 137 to 1462 K above the ground state, covering the region where silicates are expected to condense (Gail & Sedlmayr 1999) and, consequently, where a decrease in the 28SiO abundance should be seen.

We compare these data to a grid of models in which the photospheric 28SiO abundance relative to H2, , the condensation fraction of 28SiO in solid material, fcond, and radius at which this happens, Rcond, are varied. For we apply (10, 8, 6, 5, 4, 2) × 10-5. For fcond we adopt 0, 0.35, 0.65, and 0.90. The first value implies no condensation of solids. For each of the abundance profiles in which condensation was considered, we take Rcond to be either 5, 10, or 20 R. This adds up to a total of sixty models. We note that models without condensation are equivalent to models with a higher 28SiO abundance in which condensation occurs deep in the envelope, at about 2 R, as our models are not sensitive to the abundance at radii smaller than about 5 R. To give an example, a model with an initial abundance of 4 × 10-5 is analogous to a model with and fcond = 0.33 in which condensation occurs close to the star.

As shown in Table 7, the models were ranked based on the calculated reduced-χ2 fit to the line fluxes obtained with SPIRE, HIFI and PACS, listed in Table B.1. The ten best models have 28SiO abundances between 2.5 × 10-5 and 4.0 × 10-5 relative to H2 in the region between 10 and 100 R, i.e. in the region where the 28SiO emission originates (see Fig. 1). The condensation radius of 28SiO is not strongly constrained. However, condensation at 10 R or less seems preferred over condensation at 20 R. In the top-fifteen-ranked models, those with a condensation fraction of 0.65 all have higher photospheric 28SiO abundance than expected on the basis of a solar carbon-to-silicon ratio. If the photospheric 28SiO abundance is assumed to be solar ( between 4.0 and 6.0 × 10-5), a condensation fraction of 0.35 or less is preferred. Furthermore, the slope seen in Fig. 10 for the observed line fluxes distribution in terms of J is much better reproduced by models with no or very little condensation.

The fit to the J = 14 − 13 and 16 − 1528SiO lines observed by HIFI, shown in Fig. 11, is consistent with υturb = 0.7 ± 0.1 km s-1.

Table 7

Calculated reduced-χ2 of the fit to the line fluxes of the set of 28SiO lines listed in Table B.1.

6. Discussion

6.1. The turbulence velocity

The turbulence velocity probed by a given line is that of the region where the line is excited and from which photons can escape. The lines observed by HIFI for the ortho- and para-HO transitions have similar excitation energies, hence they may be expected to form in a similar part of the outflow if the spin isomers themselves occupy the same region.

We found a small difference in the turbulent velocity value that predicts best the observed line shapes of the ortho- and para-H2O, 0.8 km s-1 and 0.6 km s-1 respectively. The υturb derived by Khouri et al. (2014) based on the 12CO line shapes was 1.4 ± 1.0 km s-1. Since the low-excitation CO lines probe the outer parts of the wind, this could be another indication of a turbulent velocity gradient. However, the uncertainty in the υturb values from CO are sizable. Moreover, the observed line shapes of both H2O spin isomers are not perfectly fitted for any value of υturb. Although these diagnostics suggest a gradient, GASTRoNOoM can only calculate models with a constant turbulent velocity and therefore we cannot test this possibility at the moment. We adopt υturb = 0.7 ± 0.1 km s-1.

Maercker et al. (2009) found the dissociation radius and the shape of the abundance profile to also have an impact on the line shapes. We discuss our assumptions for these properties in Sect. 3.1. We have tested the impact of decreasing the dissociation radius on the line shapes. The lines become narrower, the red-wing being more strongly affected. The dissociation radii and the abundance profile have an impact on the line shapes and their peak position and add to the uncertainty on the determined turbulence velocity. However, if the abundance profile of ortho- and para-HO are similar, there should be no relative difference between the line shapes of the two spin isomers.

6.2. Isotopic ratios and evolutionary status of W Hya

During the evolution of low- and intermediate-mass stars leading up to the AGB and on the AGB, the abundances of the two minor oxygen isotopes are expected to vary considerably, especially due to the first- and third-dredge-up events. The effect of the first-dredge-up on the 16O/17O surface isotopic ratio is found to depend quite sensitively on the initial stellar mass, the values of this ratio after this first-dredge-up event is a steeply decreasing function of stellar mass for stars with main sequence mass between 1 and 3 M (e.g. Boothroyd et al. 1994; Lattanzio & Boothroyd 1997; Palmerini et al. 2011; Charbonnel & Lagarde 2010). During the third-dredge-up, the surface oxygen isotopic ratios are expected to change only if hot bottom burning is active (e.g. Busso et al. 1999; Charbonnel & Lagarde 2010; Karakas 2011). In this case, the abundance of 18O drops strongly while that of 17O remains almost unchanged. Justtanont et al. (2013) recently lend support to this prediction by presenting observations that clearly show that the abundance of H217O is considerably higher than that of H218O for a set of OH/IR stars observed with HIFI.

The isotopic ratios determined for W Hya, 16O/17O = 1250  and 16O/18O = 190 , are lower than the solar values of (16O/17O) = 2600 and (16O/18O) = 500. The observed value of 1250 implies that W Hya had an initial mass of about 1.5 M. Evolutionary models show the 16O/17O ratio after the first dredge-up to be independent of metallicity. Since the 16O/17O surface ratio is such a steep function of initial mass, the initial mass of W Hya would be constrained to be within 1.3 and 1.7 M. Such a star would reach the AGB phase in about 3 gigayears. If W Hya has less metals than the Sun, the determined value for the 16O/17O surface ratio would also be consistent with it having an initial mass of more than 4 M. It is unlikely, however, that W Hya is either metal-poorer than the Sun or so massive.

The value for the 16O/18O ratio does not agree with what is found with available evolutionary models. The observed value of 190 is lower than that of the Sun. All models predict this ratio to be a weak function of mass and to increase during evolution, therefore the observed value cannot be reconciled with predictions. Interestingly, Decin et al. (2010b) determine the 16O/18O ratio of the also oxygen-rich AGB star IK Tau to be 200, a value that is very close to the value determined by us for W Hya.

thumbnail Fig. 10

Best model for the 28SiO line emission (red line and crosses), with = 4 × 10-5 and fcond = 0, is compared to the line fluxes observed by SPIRE (green), PACS (blue) and HIFI (purple).

Open with DEXTER

thumbnail Fig. 11

Best 28SiO model (green-dashed line) is compared to the line shapes observed by HIFI (solid red line). The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines.

Open with DEXTER

One solution to this problem may be that both W Hya and IK Tau are richer in metals than the Sun, since the 16O/18O ratio is expected to be inversely proportional to metal content (Timmes et al. 1995). The observed isotopic ratio then requires that W Hya and IK Tau are about twice as rich in metals than is the Sun. If that is not the case, and the low 16O/18O ratio is confirmed for these two objects, our findings would imply that the evolution of the 18O surface abundance up to the AGB stage is not yet well understood. However, the uncertainty associated with the 16O/18O measurement does not allow one to draw a firm conclusion on this matter at the moment.

6.3. 28SiO condensation

If we consider a photospheric 28SiO abundance expected for solar composition (), a dust mass-loss rate of 2.8 × 10-10M yr-1 and the silicon-bearing grains to consist of olivine silicates (MgFeSiO4; see Paper I), our model requires about one-third of the silicon atoms in the wind of W Hya to be in dust grains. We thus expect to see a decrease in the 28SiO abundance from 4.8 × 10-5 to 3.2 × 10-5 over the region where silicates condense. Furthermore, our dust model predicts that silicates are formed in W Hya’s wind as close as ~5 stellar radii (or 10 AU) to the surface. Observations carried out by Zhao-Geisler et al. (2011) with MIDI/VLTI set a lower limit on the silicate formation radius at 28 photospheric radii (equivalent to 50 AU or 0.5 arcsec). That corresponds to 22 R (equivalent to ~40 AU) in our model, when we correct for the different distance adopted by them. Furthermore, aperture-masked polarimetric interferometry observations carried out by Norris et al. (2012) reveal a close-in halo of large transparent grains in W Hya. The composition of the grains could not be determined, but if this material also contains silicon it has to be considered in the silicon budget.

Our molecular-line-emission calculations indicate that 28SiO does not suffer from strong depletion in W Hya’s wind. Furthermore, an 28SiO abundance of (3.3 ± 0.8) × 10-5 between 10 and 100 R is required in order to reproduce the observed 28SiO lines. Adding to that the silicon that is in our dust model, which corresponds to an abundance of 1.6 × 10-5, we reach a total silicon abundance of 4.9 × 10-5. This is very close to the abundance expected based on a solar silicon-to-carbon ratio, i.e. 4.8 × 10-5 for 28SiO. Models with condensation occurring at radii equal to, or smaller than, 10 stellar radii are preferred to those with condensation at 20 stellar radii but we are not able to determine the condensation radius based in our data. One could expect condensation to happen over a few or even tens of stellar radii, a scenario not explored in our calculations. Despite the good agreement regarding the silicon budget, we note that our dust model was based on the one obtained by Justtanont et al. (2005) and that we have not studied the dust envelope in detail, as considering different dust species and/or distribution. Furthermore, the present dust model does not agree with the observations carried out by Zhao-Geisler et al. (2011), which show that silicates do not condense closer than about 40 AU. We will analyse in depth the dust envelope of W Hya under the light of the gas-phase wind model presented here in an upcoming study.

Regarding the outer 28SiO envelope, González Delgado et al. (2003) modelled the 28SiO pure rotational emission of the ground-vibrational state and obtained an 28SiO abundance of 1.5 × 10-5. The authors, however, studied the 28SiO abundance relative to 12CO mainly in a statistical way. They compared their models to low-excitation transitions, J = 2 − 1, 3 − 2, 5 − 4 and 6 − 5, which trace mostly the outer parts of the 28SiO envelope. The value González Delgado et al. (2003) derive for the 28SiO abundance is a factor of two lower than the 3.3 × 10-5 abundance found by us. However, if the lower mass-loss rate, 8 × 10-8 and the smaller distance, 65 parsecs, considered by them are taken into account the derived abundance should be even smaller, in the context of our model. Our calculations overpredict the emission seen in these low-excitation transitions by a factor of four, consistent with the difference in the abundances that are obtained. Lucas et al. (1992) determined the half-intensity angular radius for the 28SiO transition J = 2 − 1 to be 0.9 ± 0.1 arcsec. This value is substantially smaller than the value derived by González Delgado et al. (2003) for the e-folding radius of the 28SiO abundance, 2.4 × 1015 cm or 2 arcsec in the context of our model. The inconsistency is only apparent, as González Delgado et al. point out that these two radii are indeed expected to differ. The authors find that the e-folding radii determined by them are about three times larger than the half-intensity radius of the J = 2 − 128SiO transition for model envelopes. Schöier et al. (2004) observed the 28SiO J = 2 − 1 transition of R Dor and L2 Pup with the Australia Telescope Compact Array. By modelling the interferometric data, the authors found that the 28SiO abundance is better described by a two-component profile, a high abundance (fSiO ≈ 4 × 10-5) inner component and a lower abundance (fSiO ≈ 2 − 3 × 10-6) extended component. The radius where the abundance drops is found to be between 1 and 3 × 1015 cm.

The low-excitation transitions probe mostly the outer envelope, where dissociation occurs. The population of level J = 6 of 28SiO reaches its maximum at 60 R (1.5 arcsec or 1.8 × 1015 cm in the context of our model). Therefore, the abundances derived based on transitions from J = 6 and lower levels will depend on the assumed dissociation profile. We have calculated models with smaller dissociation radii and shallower dissociation profiles. We have done so by decreasing, respectively, the value of the e-folding radius re and of the exponent p, initially kept at p = 2, in the expression for the 28SiO abundance profile (see Eq. (2)). For reasonable values of these two parameters, the models are still unable to fit the low-excitation transitions. An alternative possibility to explain simultaneously the high- and low-excitation lines may be that 28SiO suffers from further depletion from the gas-phase in-between the region where the lines observed by Herschel and those observed by SEST are excited. At such large distances from the star, however, condensation and dissociation are indistinguishable on the basis of 28SiO line emission modelling. The Jup> 10 lines modelled by us, are all produced closer to the star and, therefore, trace the depletion of 28SiO independently of dissociation. Thanks to the apparent complex nature of the 28SiO dissociation region and since the choice of dissociation radius does not have a significant impact on the derived value for the 28SiO abundance, we do not attempt to fit the low-excitation 28SiO lines in detail. The value for the 28SiO depletion obtained by us is representative of the inner wind, for r< 1.5 × 1015 cm.

7. Summary

We present an analysis of the ortho-HO, para-HO, ortho-HO, para-HO, ortho-HO, para-HO and 28Si16O emission from the wind of the nearby oxygen-rich AGB star W Hya, as measured by the three instruments on board Herschel. The work builds on the structure model of Khouri et al. (2014), derived on the basis of 12CO lines, and is the first combined 12CO, HO and 28SiO analysis of this source.

The original structure model poorly constrained the turbulence component of the velocity field in the outflow. HO and 28SiO lines put much firmer constraints on the value of the turbulence velocity, essentially because they are much optically thicker than the 12CO lines. The presence of turbulence motions causes the H2O and 28SiO lines to shift to longer wavelengths, which, when compared to 12CO profiles that form in roughly the same part of the wind appear to imply a blue-wing absorption. We find slightly different values for υturb for ortho-HO and para-HO, 0.8 and 0.6 km s-1 respectively, but as our code is not able to calculate models with a gradient in υturb, we have not explored this further.

The abundance of ortho-HO and para-HO relative to H2 are and . We also place constraints on the abundances of ortho-HO and para-HO, and find an ortho-to-para ratio of 2.5  – in agreement with the value of three expected for AGB stars. The 16O/17O ratio is found to be 1250  and suggests that W Hya has an initial mass of about 1.5 M. We find an 16O/18O ratio of 190 , which cannot be explained by the current generation of evolutionary models. It might be reconciled with predictions if W Hya is richer in metals than the Sun, but no firm conclusions can be drawn on this matter given the large uncertainties on the abundance determination.

We find an 28SiO abundance between 10 and 100 R of 3.3 ± 0.8 × 10-5 relative to H2. Adding to this gas-phase abundance the abundance needed by our dust model, equivalent to 1.6 × 10-5, we can account for all silicon in the wind of W Hya if a solar silicon-to-carbon ratio is assumed.


1

HIPE is available for download at http://herschel.esac.esa.int/hipe/

Acknowledgments

HIFI has been designed and built by a consortium of institutes and university departments from across Europe, Canada, and the United States under the leadership of SRON Netherlands Institute for Space Research, Groningen, The Netherlands and with major contributions from Germany, France, and the US. Consortium members are Canada: CSA, U. Waterloo; France: CESR, LAB, LERMA, IRAM; Germany: KOSMA, MPIfR, MPS; Ireland, NUI Maynooth; Italy: ASI, IFSI-INAF, Osservatorio Astrofisico di Arcetri-INAF; Netherlands: SRON, TUD; Poland: CAMK, CBK; Spain: Observatorio Astronómico Nacional (IGN), Centro de Astrobiología (CSIC-INTA). Sweden: Chalmers University of Technology Đ MC2, RSS & GARD; Onsala Space Observatory; Swedish National Space Board, Stockholm University Đ SStockholm Observatory; Switzerland: ETH Zurich, FHNW; USA: Caltech, JPL, NHSC. PACS has been developed by a consortium of institutes led by MPE (Germany) and including UVIE (Austria); KUL, CSL, IMEC (Belgium); CEA, OAMP (France); MPIA (Germany); IFSI, OAP/AOT, OAA/CAISMI, LENS, SISSA (Italy); IAC (Spain). This development has been supported by the funding agencies BMVIT (Austria), ESA-PRODEX (Belgium), CEA/CNES (France), DLR (Germany), ASI (Italy), and CICYT/MCYT (Spain). SPIRE has been developed by a consortium of institutes led by Cardiff Univ. (UK) and including Univ. Lethbridge (Canada); NAOC (China); CEA, LAM (France); IFSI, Univ. Padua (Italy); IAC (Spain); Stockholm Observatory (Sweden); Imperial College London, RAL, UCL-MSSL, UKATC, Univ. Sussex (UK); Caltech, JPL, NHSC, Univ. Colorado (USA). This development has been supported by national funding agencies: CSA (Canada); NAOC (China); CEA, CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB (Sweden); STFC (UK); and NASA (USA). HIFISTARS: The physical and chemical properties of circumstellar environments around evolved stars, (P.I. Bujarrabal, is a Herschel/HIFI guaranteed time key program (KPGT_vbujarra_1) devoted to the study of the warm gas and water vapour contents of the molecular envelopes around evolved stars: AGB stars, red super- and hyper-giants; and their descendants: pre-planetary nebulae, planetary nebulae, and yellow hyper-giants. HIFISTARS comprises 366 observations, totalling 11 186 min of Herschel/HIFI telescope time. See http://hifistars.oan.es; and Key_Programmes.shtml and UserProvidedDataProducts.shtml in the Herschel web portal (http://herschel.esac.esa.int/) for additional details. T.Kh. gratefully acknowledges the support from NWO grant 614.000.903. F.K. is supported by the FWF project P23586 and the ffg ASAP project HIL. M.Sc. and R.Sz. acknowledge support by the National Science Center under grant (N 203 581040).

References

Online material

Appendix A: Molecular models

When modelling the HO transitions for all isotopologues, we include the 45 lowest levels of the ground and first vibrational states (i.e. the bending mode ν2 = 1 at 6.3 μm). For the two spin isomers of the main isotopologue, we have also included excitation to the first excited vibrational state of the asymmetric stretching mode (ν3 = 1). The difference on the model line fluxes due to the inclusion of the ν3 = 1 level is found to be 20% at maximum (Decin et al. 2010a). The frequencies, level energies and Einstein A coefficients were retrieved from the HITRAN HO line list (Rothman et al. 2009). The collisional rates between HO and H2 were extracted from Faure et al. (2007).

Following Decin et al. (2010a), we consider the 40 lowest rotational levels of the ground and first vibrationally excited states when modelling the 28SiO transitions. The Einstein A coefficients, energy levels and frequencies were taken from Langhoff & Bauschlicher (1993). The collisional rates between 28SiO and H2 were retrieved from the LAMBDA-database (Schöier et al. 2005).

Appendix B: Observed 28SiO and HO line fluxes

Table B.1

Integrated line fluxes and uncertainties for the 28SiO transitions observed with Herschel for W Hya.

Table B.2

Extracted ortho-HO lines from the PACS observations.

Table B.3

Extracted p-HO lines from the PACS observations.

All Tables

Table 1

Model parameters of W Hya as derived in Paper I.

Table 2

H2O and 28SiO transitions observed by HIFI for W Hya.

Table 3

Observed total line fluxes of HO and 28SiO transitions observed by two different Herschel instruments, HIFI and SPIRE or PACS.

Table 4

Turbulence velocities and abundances relative to H2 considered in the models studying the effects of HO.

Table 5

Best fit results for the abundances of the isotopologues of H2O.

Table 6

Derived isotopologic ratios.

Table 7

Calculated reduced-χ2 of the fit to the line fluxes of the set of 28SiO lines listed in Table B.1.

Table B.1

Integrated line fluxes and uncertainties for the 28SiO transitions observed with Herschel for W Hya.

Table B.2

Extracted ortho-HO lines from the PACS observations.

Table B.3

Extracted p-HO lines from the PACS observations.

All Figures

thumbnail Fig. 1

Normalized abundance profile of HO (long-dashed green), 28SiO (double-short-dashed yellow) and 12CO (solid red) are compared to the normalized populations of levels J = 16 (dashed-dotted brown), J = 10 (dotted purple) and J = 6 (short-dashed blue) of 12CO.

Open with DEXTER
In the text
thumbnail Fig. 2

Observed line shapes of the transitions J = 16 − 15 (dotted purple) and J = 10 − 9 (solid red) of 12CO, J = 14 − 13 of 28SiO (long-dashed green), and 11,0-00,1 of ortho-HO (short-dashed blue). The vertical line marks the adopted υLSR = 40.4 km s-1.

Open with DEXTER
In the text
thumbnail Fig. 3

Effects of varying the HO abundance, the turbulence velocity, the temperature power law exponent and the velocity power law exponent on the normalized line strength are shown for o-HO (11,0 − 10,1). The standard model is represented by the red line in each panel, it has its values given in Table 1 and fH2O = 4 × 10-4. The short-dashed purple and the long-dashed green represent models with respectively: fH2O = 1 × 10-4 and 1.6 × 10-3 in the first panel; υturb = 0.8 and 2.0 km s-1 in the second panel; ϵ = 0.2 and 1.1 in the third panel; and β = 1.5 and 10 in the fourth panel. All other parameters have the same value as the standard model. We note that varying υturb has the effect of a line shift.

Open with DEXTER
In the text
thumbnail Fig. 4

Normalized profiles of ortho- and para-HO transitions observed by HIFI compared to models with different values of the turbulence velocity and all other parameters kept fixed. The solid red line represents the data, and models with turbulence velocities of 0.5, 0.8 and 1.4 km s-1 are represented by dotted purple, long-dashed green and short-dashed blue lines, respectively. The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines.

Open with DEXTER
In the text
thumbnail Fig. 5

Best fit model (dashed-blue line) to the ortho-HO lines compared to the lines observed by HIFI (solid-red line). The model has an ortho-HO abundance of 6 × 10-4 and a turbulence velocity of 0.7 km s-1. The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines.

Open with DEXTER
In the text
thumbnail Fig. 6

Best fit model (dashed-blue line) to the para-HO lines compared to the lines observed by HIFI (solid-red line). The model has a para-HO abundance of 3 × 10-4 and a turbulence velocity of 0.7 km s-1. The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines.

Open with DEXTER
In the text
thumbnail Fig. 7

The solid red line represents the HIFI data, the dashed green, short-dashed blue, dotted purple, dashed-dotted light-blue and short-dashed-dotted brown lines represent, respectively, models with abundances relatively to H2 of: (4, 2, 1.2, 0.4, 0.2) × 10-6 for for ortho-HO; (1.2, 0.8, 0.4, 0.2, 0.12) × 10-4 for ortho-HO; (8, 4, 2, 1.2, 0.4) × 10-7 for para-HO; and (8, 4, 2, 1.2, 0.4) × 106 for para-HO. We adopted a value of 0.7 km s-1 for the turbulence velocity in these calculations. The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines.

Open with DEXTER
In the text
thumbnail Fig. 8

A model for the 12CO lines with υturb = 0.7 km s-1 (short-dashed blue line) is compared to the best model found in Paper I (υturb = 1.4 km s-1; dotted-dashed brown line). The observed 12CO lines are shown in red. The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines. The lines J = 16 − 15, 10 − 9 and 6 − 5 were observed with HIFI, J = 4 − 3 and 3 − 2 with APEX, J = 2 − 1 with SMT and J = 1 − 0 with SEST.

Open with DEXTER
In the text
thumbnail Fig. 9

PACS spectrum (represented by the grey-filled histogram) compared to our best model, with the parameters given in Table 1 but with υturb of 0.7 km s-1 and with ortho- and para-HO abundances of 6 × 10-4 and 3 × 10-4, respectively. The 12CO model is shown by the full-blue line, the ortho-HO, by the full-red line and the para-HO, by the full-yellow line.

Open with DEXTER
In the text
thumbnail Fig. 10

Best model for the 28SiO line emission (red line and crosses), with = 4 × 10-5 and fcond = 0, is compared to the line fluxes observed by SPIRE (green), PACS (blue) and HIFI (purple).

Open with DEXTER
In the text
thumbnail Fig. 11

Best 28SiO model (green-dashed line) is compared to the line shapes observed by HIFI (solid red line). The adopted value for the υLSR of 40.4 km s-1 was subtracted from the observed lines.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.