Free Access

This article has an erratum: [https://doi.org/10.1051/0004-6361/201732276e]


Issue
A&A
Volume 616, August 2018
Article Number A106
Number of page(s) 34
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201732276
Published online 28 August 2018

© ESO 2018

1 Introduction

The asymptotic giant branch (AGB) phase of stars with an initial mass up to 8 M is characterised by strong mass loss via a stellar outflow or wind. The outflow removes the outer stellar layers and creates an extended circumstellar envelope (CSE). Through their outflows, AGB stars eject nuclear reaction products in the form of gas and dust into the interstellar medium and hence play an important role in the chemical enrichment of the Universe (Habing & Olofsson 2003, and references therein). The AGB outflow can be split up into three dynamically and chemically different regions: the inner region close to the central star in which the materialexperiences shocks caused by stellar pulsations and in which the first chemical processing of photospheric gas takes place; the intermediate region where dust condenses from the gas phase and material is accelerated outward; and the outer region where the material has reached its terminal velocity and the chemistry is dominated by photodissociation caused by the penetration of interstellar UV photons. The molecular composition of the CSE is determined by the elemental carbon-to-oxygen abundance ratio of the AGB star because of the large binding energy of CO. Oxygen-rich stars have C/O < 1, and therefore an oxygen-rich chemistry throughout their outflows. Similarly, carbon-rich stars have C/O > 1 and a carbon-rich chemistry.

The initial molecular composition of the outflow is established by the physical circumstances in its innermost region. Thermodynamic equilibrium (TE) can be assumed at the stellar surface, which leads to a certain molecular composition with a clear O- or C-rich signature (Willacy & Cherchneff 1998). However, several molecular species have been detected that are expected, under the assumption of TE chemistry, either to be absent in the inner wind or to have a smaller abundance than observed in the outer regions. In this paper, we refer to these species as unexpected species for convenience and brevity. In O-rich outflows, there are detections of several unexpected C-rich molecules such as HCN, CS, and CN (Omont et al. 1993; Bujarrabal et al. 1994; Justtanont et al. 1996). Similarly in C-rich outflows, for example, warm H2 O and CH3 CN have been detected close to the central star (Decin et al. 2010a; Neufeld et al. 2011; Agúndez et al. 2015). In both O- and C-rich outflows, NH3 has been detected in greater abundance than predicted by TE (Keady & Ridgway 1993; Menten & Alcolea 1995; Menten et al. 2010; Wong et al. 2018). To reproduce the unexpectedly high observed abundances of these molecules in the outer wind, the chemical models of Willacy & Millar (1997) and Willacy & Cherchneff (1998) relied on the injection of molecules bearing the deficient element. A proposed solution to these discrepancies in abundance is to take into account the non-equilibrium character of the chemistry in the innermost region of the wind. In this region, the chemistry is taken out of TE due to shocks caused by the pulsating AGB star (e.g. Duari et al. 1999, Cherchneff 2006, Gobrecht et al. 2016). These shock-induced non-equilibrium chemistry models are able to reproduce most observed abundances, but fail to account for observed abundances, such as the NH3 abundance in the O-rich star IK Tau (Gobrecht et al. 2016), for example.

Clumpy structures, and non-uniformity in general, have been widely observed in AGB outflows. Both large-scale structures, such as spirals (e.g. Mauron & Huggins 2006, Maercker et al. 2016) and disks (e.g. Kervella et al. 2014), and small-scale inhomogeneities or clumps (e.g. Khouri et al. 2016, Agúndez et al. 2017), have been detected. As an alternative to the shock-induced non-equilibrium models, Agúndez et al. (2010) proposed the penetration of interstellar UV photons into the inner wind due to the clumpiness of the outflow. By allowing UV photons to reach the otherwise shielded inner region of the outflow, the CO and N2 bonds can be broken and the deficient element (C in O-rich outflows, O in C-rich outflows, N in both) is released close to the star. Agúndez et al. (2010) assumed that a certain fraction of UV photons is able to reach the inner wind uninhibited, but did not take the clumpy density distribution of the outflow into account. Their model results in abundances of unexpected molecules well above TE for both O- and C-rich outflows, and produces NH3 in both.

In this paper, we include a clumpy density distribution in our spherically symmetric chemical model by implementing a porosity formalism. The porosity formalism provides us with a mathematical framework that not only describes the effect of a clumpy outflow on its optical depth (and thereby on the penetration of UV photons, i.e. porosity), but also allows us to take the local overdensities (i.e. clumping) into account in the chemical model. We explore the parameter space of the formalism for both O- and C-rich outflows with various mass-loss rates, and study its effects on the abundances of the unexpected species. We consider two scenarios: (i) a simplified outflow where all material is located in mass conserving clumps, and (ii) a perhaps more realistic outflow consisting of both a clump component and an inter-clump component.

In Sect. 2, we describe the chemical model and reaction network used. We elaborate on the analytical description of the clumping and porosity properties of the outflow and their implementation into the model in Sect. 3. The results for both O- and C-rich outflows are presented in Sect. 4. A discussion and conclusions follow in Sects. 5 and 6.

Table 1

Physical parameters of the CSE model.

Table 2

Parent species and their initial abundances relative to H2 for the C- and O-rich CSE.

2 Chemical model of the outflow

The chemical reaction network used is the most recent release of the UMIST Database for Astrochemistry (UDfA), RATE12 (McElroy et al. 2013), a gas-phase only chemical network containing 6173 reactions involving 467 species. Our chemical model is based on the UDfA CSE model. Both the chemical network and CSE model are publicly available1. The one-dimensional model describes a uniformly expanding, spherically symmetric CSE with a constant mass-loss rate. The number density profile n(r) hence falls as 1∕r2, with r the distance from the centre of the star. The velocity of the material is assumed to be constant throughout the outflow. The effect of CO self-shielding is taken into account during the calculation using a single-band approximation (Morris & Jura 1983). H2 is assumed to be completely self-shielded, so that nH2(r)n(r)$n_{\textrm{H}_2}(r) \approx n(r)$. A more detailed description of the model can be found in Millar et al. (2000), Cordiner & Millar (2009), and McElroy et al. (2013).

The initial abundance of all species is set to zero, except for the parent species. These are defined as having formed just above the photosphere and are injected into the CSE at the start of our model at r = 2 R* = 1 × 1014 cm. The adjustments made to the code by introducing the porosity formalism are described in detail in Sect. 3.3. We have also changed the assumed temperature structure throughout the outflow to a power-law, T(r)=T*(R*r)ϵ,\begin{equation*} T(r) = T_* \left( \frac{R_*}{r} \right)^{\epsilon}, \end{equation*}(1)

with T* and R* the stellar temperature and radius, and ϵ the exponent characterising the power-law. The temperature profile has a lower limit of 10 K, so as to prevent unrealistically low temperatures in the outer CSE (Cordiner & Millar 2009).

In order to compare our results to those of Agúndez et al. (2010), we have adopted the same physical parameters (Table 1) and parent species, and their initial abundances for the O- and C-rich CSE (Table 2). The constant outflow velocity ν is the main physical difference with the model of Agúndez et al. (2010), who have assumed an expansion velocity of 5 km s−1 for r ≤ 5 R* and a velocity of 15 km s−1 beyond. The RATE12 network does not include 13CO and its self-shielding, in contrast to the chemical network used by Agúndez et al. (2010).

3 Porosity formalism

By assuming that the non-uniformity of the outflow is caused by a statistical ensemble of small-scale density inhomogeneities or clumps, we are able to use a so-called porosity formalism (Owocki et al. 2004; Owocki & Cohen 2006; Sundqvist et al. 2012, 2014) to describe the effect of a clumpy outflow on its chemistry. See also Feldmeier et al. (2003) and Oskinova et al. (2004, 2006) for an alternative formulation regarding radiative transfer effects in an inhomogeneous medium.

Within this formalism, we assume that the outflow is a stochastic two-component medium consisting of overdense clumps and a rarified inter-clump medium. The clumps are assumed to conserve their mass throughout the outflow. The fraction of the total volume occupied by the clumps is set by fvol, the clump volume filling factor. The mean (laterally averaged) density of the outflow can then be written as ρ=(1fvol)ρic+fvolρcl=M˙4πr2ν,\begin{equation*}\langle \rho \rangle = (1 - f_{\textrm{vol}})\ \rho_{\textrm{ic}} &#x002B; f_{\textrm{vol}}\ \rho_{\textrm{cl}} = \frac{\dot{M}}{4\pi r^2 \nu_{\infty}}, \end{equation*}(2)

where the subscripts “cl” and “ic” refer to the clump and inter-clump density components respectively. The last equality assumes mass-conservation as compared to a corresponding smooth outflow with density ρsm = ⟨ρ⟩. For ease of notation, we drop the radial dependency of all expressions. The density contrast between the inter-clump component and the mean density is set by the parameter (Sundqvist et al. 2014) ficρicρ.\begin{equation*}f_{\textrm{ic}} \equiv \frac{\rho_{\textrm{ic}}}{\langle\rho\rangle}. \end{equation*}(3)

By splitting the outflow into two separate components, their specific over- or under-densities compared to the corresponding smooth outflow can be taken into account in our chemical model. We note that within our statistical formalism, no assumptions on the specific locations of individual clumps need to be made.

Generally, if clumps become optically thick, this leads both to a local self-shielding of opacity within the clumps and to an increased escape of radiation through the porous channels in between the clumps. The overall net effect of a clumpy outflow is therefore a reduction in its mean opacity and hence its optical depth for a given wind mass. A clumpy outflow hence affects the penetration of interstellar UV photons into the circumstellar material, which is important to the chemistry as it induces photodissociation reactions. The extinction of the interstellar UV radiation field in the outflow is given by eτ, with τ the optical depth. For a smooth outflow, the radial optical depth is given by τ=rκ ρsmdr=κrρsm dr,\begin{equation*}\tau = \int_{r^{\prime}}^{\infty} \kappa\ \rho_{\textrm{sm}}\ \textrm{d}r = \kappa \int_{r^{\prime}}^{\infty} \rho_{\textrm{sm}}\ \textrm{d}r, \end{equation*}(4)

where the latter equality assumes that the opacity κ of the outflow is spatially constant. We have here that κ ≈ 1200 cm2 g−1; its calculation can be found in Appendix A. Non-radial rays are accounted for according to Eq. (4) of Morris & Jura (1983). By introducing the radial optical depth down to the stellar surface for a smooth wind as a reference optical depth τ*, given by τ*=κM˙4πνR*,\begin{equation*} \tau_* = \frac{\kappa \dot{M}}{4\pi \nu_{\infty} R_*}, \end{equation*}(5)

we can writethat τ = τ*R*r. The opticaldepth of a smooth outflow is therefore proportional to . For the physical parameters assumed in this paper, we have that τ*≈ 800 for =10−5 M yr−1, τ* ≈ 80 for =10−6 M yr−1 and τ* ≈ 8 for =10−7 M yr−1. The mass-loss rate of a smooth outflow hence strongly influences its optical depth.

The increased interstellar UV radiation field in a clumpy outflow is now taken into account by changing the optical depth τ to an effective optical depth τeff. Equation (4) then changes to τeff=rκeff ρdr,\begin{equation*}\tau_{\textrm{eff}} = \int_{r^{\prime}}^{\infty} \kappa_{\textrm{eff}}\ \langle\rho\rangle \ \textrm{d}r, \end{equation*}(6)

with τeff and κeff the effective optical depth and opacity of the clumpy outflow, respectively.

Within our statistical porosity model, this effective opacity κeff depends on the interaction probability P of a photon with a clump, and thus depends on the optical depth of an individual clump τcl, given by τcl=κρcll,\begin{equation*}\tau_{\textrm{cl}} = \kappa \ \rho_{\textrm{cl}} \ l, \end{equation*}(7)

with l the characteristic length scale of the clumps. The effective extinction κeff ⟨ρ⟩ (per unit length) is then given by κeffρ=P(τcl)nclAcl=P(τcl)h,\begin{equation*}\kappa_{\textrm{eff}}\ \langle \rho \rangle = P\left(\tau_{\textrm{cl}}\right)\ n_{\textrm{cl}}\ A_{\textrm{cl}} = \frac{P(\tau_{\textrm{cl}})}{h}, \end{equation*}(8)

with ncl the clump number density and Acl the clump’s geometric cross section (e.g. Feldmeier et al. 2003), and where the last equality introduces the porosity length h (Owocki et al. 2004; Sundqvist et al. 2012). This paper assumes clumps that are spherical or randomly oriented with an arbitrary shape, which then leads to a statistically isotropic effective opacity (see Sundqvist et al. 2012). The porosity length represents the local mean free path between two clumps. Since the total volume associated with one clump is V tot = 1∕ncl, we can also use V cl = l3, Acl = l2, and fvol = V clV tot to rewrite the porosity length as h=lfvol.\begin{equation*}h = \frac{l}{f_{\textrm{vol}}}. \end{equation*}(9)

The probability that a photon will be absorbed by a given clump is P(τcl)=1eτcl$P\left(\tau_{\textrm{cl}}\right) = 1 - e^{-\tau_{\textrm{cl}}}$. But the relation between the effective and smooth opacity in principal depends also on the properties of the clumpy outflow. Assuming here an exponential distribution of clumps (i.e. a Markovian mixture; see Sundqvist et al. 2012), f(τcl)=eτ/τclτcl,\begin{equation*} f\left(\tau_{\textrm{cl}}\right) = \frac{e^{-\tau/\tau_{\textrm{cl}}}}{\tau_{\textrm{cl}}} ,\end{equation*}(10)

with τ=0τ f(τ)dτ=τcl,\begin{equation*} \langle \tau \rangle = \int_0^{\infty} \tau f\left( \tau \right) \textrm{d}\tau = \tau_{\textrm{cl}}, \end{equation*}(11)

and an outflow where all material is located inside of clumps (a so-called one-component model, see Sect. 3.1), gives upon integration the bridging law for the effective opacity of the outflow (Owocki & Cohen 2006; Sundqvist et al. 2012) κeffκ=11+τcl.\begin{equation*}\frac{\kappa_{\textrm{eff}}}{\kappa} = \frac{1}{1&#x002B;\tau_{\textrm{cl}}}. \end{equation*}(12)

By using this expression in the calculation of the effective optical depth (Eq. (6)), we can now take the effect ofporosity upon the radiation field into account for a one-component outflow. We note that as τcl goes to zero in Eq. (12), the smooth opacity and wind optical depths are recovered. For the other limit of τcl ≫ 1, we obtain instead κeffρ⟩ = 1∕h. That is, for such fully opaque clumps the effective extinction is solely affected by the porosity length, and is independent of κ. For this case, h also represents the actual local mean free path of photons.

This paper further assumes that the clump volume filling factor fvol is spatially constant. Since we also assume mass conserving clumps, this leads to clumps expanding according to l(r)=l*(rR*)2/3,\begin{equation*}l(r) = l_* \left(\frac{r}{R_*}\right)^{2/3}, \end{equation*}(13)

with l* the characteristic length scale near the stellar surface. Following Eq. (9), h then has the same radial dependence, h(r)=h*(rR*)2/3,\begin{equation*}h(r) = h_* \left(\frac{r}{R_*}\right)^{2/3}, \end{equation*}(14)

with h* = l*fvol the porosity length at the stellar surface.

Two cases are considered in this paper: the previously introduced simplified case of a one-component model where all material is located inside clumps and the inter-clump component is effectively void (Sect. 3.1), and a case where material is located inboth the clump and inter-clump components, the two-component model (Sect. 3.2).

3.1 One-component outflow

In a one-component outflow, all material is located inside clumps and the inter-clump component is hence effectively void. Such an extreme outflow gives us an upper limit to the effect of porosity on the chemistry within the outflow, both regarding its effect on the density contrast between a clumpy outflow and a smooth outflow and the increase in UV radiation field throughout the clumpy outflow.

As the inter-clump component is effectively void, Eq. (2) reduces to ρ=fvolρcl.\begin{equation*}\langle \rho \rangle = f_{\textrm{vol}}\ \rho_{\textrm{cl}}. \end{equation*}(15)

Using Eqs. (9), (14), and (15), the clump optical depth (Eq. (7)) can be written as τcl=κρh=τ*R*1/3h*r4/3. \begin{equation*}\begin{split} \tau _{\textrm{cl}} &= \kappa\ \langle\rho\rangle\ h \\ &= \tau_*\ R_*^{1/3}\ h_*\ r^{-4/3}. \end{split} \end{equation*}(16)

Using Eq. (12), the effective optical depth throughout a one-component clumpy outflow is then given by τeff=κrρ 11+τcldr=τ*R*r1r2+τ*R*1/3h*r2/3 dr, \begin{equation*}\begin{split} \tau_{\textrm{eff}} &= \kappa \int_{r^{\prime}}^{\infty} \langle\rho\rangle\ \frac{1}{1&#x002B;\tau_{\textrm{cl}}} \textrm{d}r \\ &= \tau_*\ R_* \int_{r^{\prime}}^{\infty} \frac{1}{r^2 &#x002B;\tau_*\ R_*^{1/3}\ h_*\ r^{2/3}}\ \textrm{d}r, \end{split} \end{equation*}(17)

where Eq. (16) has been used for τcl. If h* = 0 (i.e. τcl = 0), this indeed reduces to the optical depth of a smooth wind τ = τ* R*r. The analytic solution of this integral can be found in Appendix C.1. We note here that the scaling of the effective optical depth of a one-component medium is influenced only by the porosity length h*.

For larger values of h*, both the clump optical depth and the mean free path in between the clumps increase (Fig. 1). As the clump optical depth increases, the ratio of the effective opacity of the clumpy outflow to that of the smooth outflow κeffκ decreases. Alarger fraction of the interstellar radiation field is therefore able to reach the inner and intermediate regions of the outflow: the visual extinction is smaller in these regions compared to that of a smooth outflow. In addition, as the mean free path in between the clumps increases, the onset of the decrease in visual extinction (compared to the visual extinction of a smooth outflow, associated with the position where τcl is no longernegligible) is located further out in the outflow (Fig. 1). Since τcl scales with τ*, the clump optical depth decreases with decreasing mass-loss rate. Following Eq. (12), this leads to a less pronounced effect of porosity on the effective opacity κeff and optical depth τeff for lower mass-loss rates.

thumbnail Fig. 1

Visualextinction AV throughout a one-component outflow with = 10−5 M yr−1 (left panel) and = 10−7 M yr−1 (right panel). Results for a smooth, uniform outflow are shown in black. Results for different values of h* are shown in blue. For reference, 1 R* = 5 × 1013 cm. Figure B.1 shows the corresponding decrease in interstellar UV radiation field.

thumbnail Fig. 2

As Fig. 1, but for a two-component outflow with fic = 0.1. For reference, 1 R* = 5 × 1013 cm.

3.2 Two-component outflow

In a two-component outflow, material is located both inside the clumps and the inter-clump component. Both components will now contribute to the total opacity of the clumpy outflow. The density of a two-component outflow is given by Eq. (2). Using Eqs. (3), (9), and (14), the clump optical depth (Eq. (7)) is then given by τcl=κρh(1(1fvol)fic)=τ*R*1/3h*(1(1fvol)fic)r4/3. \begin{equation*}\begin{split} \tau_{\textrm{cl}} &= \kappa \ \langle\rho\rangle\ h\ \left( 1 - \left( 1 - f_{\textrm{vol}} \right) f_{\textrm{ic}} \right) \\ &= \tau_*\ R_*^{1/3}\ h_*\ \left( 1 - \left( 1 - f_{\textrm{vol}} \right) f_{\textrm{ic}} \right)\ r^{-4/3}. \end{split} \end{equation*}(18)

The bridging law between the effective opacity of the clumpy medium and that of a smooth outflow now needs to be slightly modified. Following Sundqvist et al. (2014), this is given by κeffκ=1+τclfic1+τcl.\begin{equation*}\frac{\kappa_{\textrm{eff}}}{\kappa} = \frac{1 &#x002B;\tau_{\textrm{cl}}\ f_{\textrm{ic}}}{1&#x002B;\tau_{\textrm{cl}}}. \end{equation*}(19)

We note that as τcl goes to zero, the smooth opacity and optical depths are recovered. This is again recovered when fic = 1, corresponding to two equal components. When the inter-clump component is effectively void (fic = 0), both the clump optical depth and the bridging law reduce to their one-component expressions. In the limit where τcl ≫ 1, we find that κeffκ = fic. Unlike a one-component outflow, where the effective extinction κρ⟩ saturates to 1∕h, we now find that the ratio of the effective opacity of the clumpy outflow to that of the smooth outflow κeffκ saturates to fic. The opacity of a two-component outflow hence always depends on κ, implying a larger effective extinction compared to a one-component outflow (Figs. 13).

The effective optical depth throughout a two-component outflow is then given by τeff=κrρ 1+τclfic1+τcldr=τ*R*r1+τ*R*1/3h*(1(1fvol)fic)r4/3r2+τ*R*1/3h*(1(1fvol)fic)r2/3 dr, \begin{equation*}\begin{split} \tau_{\textrm{eff}} &= \kappa \int_{r^{\prime}}^{\infty} \langle\rho\rangle\ \frac{1 &#x002B;\tau_{\textrm{cl}}\ f_{\textrm{ic}}}{1&#x002B;\tau_{\textrm{cl}}} \textrm{d}r\\ &= \tau_*\ R_* \int_{r^{\prime}}^{\infty} \frac{1&#x002B;\tau_*\ R_*^{1/3}\ h_*\ \left( 1 - \left( 1 - f_{\textrm{vol}} \right) f_{\textrm{ic}} \right)\ r^{-4/3}}{r^2 &#x002B;\tau_*\ R_*^{1/3}\ h_*\ \left( 1 - \left( 1 - f_{\textrm{vol}} \right) f_{\textrm{ic}} \right)\ r^{2/3}}\ \textrm{d}r, \end{split} \end{equation*}(20)

where Eq. (18) has been used for τcl. If h* = 0, this again reduces to the optical depth of the smooth wind τ = τ* R*r. The analytic solution of this integral can be found in Appendix C.2. In contrast to the one-component model above, the scaling of the effective optical depth of a two-component outflow is now influenced by two parameters: the porosity length h* and the inter-clump density contrast fic.

Similar to a one-component outflow, larger values of h* cause a larger decrease in visual extinction in the inner and intermediate regions of the outflow and the onset of the decrease in visual extinction (compared to the visual extinction of a smooth outflow) is located further out (Fig. 2). Smaller values of fic further leadto an increase in the clump optical depth and a decrease in the ratio κeffκ (Fig. 3). A larger fraction of the interstellar radiation field is able to reach the inner and intermediate regions of the outflow, leading to a smaller visual extinction in these regions compared to that of a smooth outflow. In analogy with the one-component case, the effect of porosity is less pronounced for lower mass-loss rates.

thumbnail Fig. 3

Visualextinction AV throughout a two-component outflow with h* = 1015 cm and = 10−5 M yr−1 (left panel) and = 10−7 M yr−1 (right panel). Results for a smooth, uniform outflow are shown in black. Results for different values of fic are shown in blue. For reference, 1 R* = 5 × 1013 cm.

3.3 Implementation into the UdfA CSE model

Three extra parameters are thus necessary to fully include the effect of clumping and porosity in the chemical model: fic, fvol, and l*. These are provided as input parameters to the model. A model is calculated for each component separately: one for the clump component, and, if present, one for the inter-clump component.

In all models, the mean density of the corresponding smooth outflow ⟨ρ⟩ is calculated using the mass-loss rate and outflow velocity (Eq. (2)). The mean density ⟨ρ⟩ is then used to calculate the optical depth of the corresponding smooth outflow. This value goes into the calculation of the effective optical depth of the clumpy outflow τeff using Eq. (C.4). In practice, this equation is also used for the one-component outflow, as it corresponds to the special case fic = 0 (see Eq. (C.2)). We note that for a two-component model, our formulation assumes that clumps and inter-clump medium are exposed to the same UV radiation field, so that the same effective optical depth can be used in both components.

The density of the clump component is derived from Eq. (2): ρcl=1(1fvol)ficfvolρρ.\begin{equation*}\rho_{\textrm{cl}} = \frac{1 - \left( 1-f_{\textrm{vol}}\right) f_{\textrm{ic}} }{f_{\textrm{vol}}}\ \langle \rho \rangle\ \geq\ \langle\rho\rangle. \end{equation*}(21)

For a one-component model, this reduces to ρcl = ⟨ρ⟩∕fvol. The density of the inter-clump component is derived from Eq. (3): ρic=ficρρ.\begin{equation*}\rho_{\textrm{ic}} = f_{\textrm{ic}}\ \langle \rho \rangle\ \leq\ \langle\rho\rangle. \end{equation*}(22)

We aldo highlight that the use of τeff (instead of τ) in the calculations affects only the rate of photodissociation reactions induced by UV photons. However, the specific density of the component (ρcl or ρic) is used when calculating the rates of two-body reactions. Since cosmic-ray photons and particles are assumed to have a radially independentincidence, the overall ionisation of the outflow depends on its density distribution.

Finally then, the resulting mean number density of species X ⟨nX ⟩(r) in the clumpy outflow is calculated using fvolncl,X(r)+(1fvol)nic,X(r)$f_{\textrm{vol}}\ n_{\textrm{cl,X}}(r) &#x002B; \left( 1 - f_{\textrm{vol}} \right) n_{\textrm{ic,X}}(r)$, as in Eq. (2). Its fractional abundance relative to H2yX ⟩(r) is therefore given by yX(r)=ycl,X+fic(1fvol)(yic,Xycl,X),\begin{equation*}\langle y_{\textrm{X}} \rangle(r) = y_{\textrm{cl,X}} &#x002B; f_{\textrm{ic}}\ \left( 1-f_{\textrm{vol}} \right)\ \left( y_{\textrm{ic,X}} - y_{\textrm{cl,X}} \right), \end{equation*}(23)

with ycl,X and yic,X its abundance relative to the H2 abundance of the clump and inter-clump component, respectively.

Table 3

Porosity lengths corresponding to the outflows modelled in Sect. 4.

3.4 Relevant parameter ranges

The parameter l* denotes the characteristic length scale of the clumps at the stellar radius. Since it is physically unrealistic to have clumps appearing at the stellar surface with larger volumes than that of the star, we have that l*R*. The parameters fvol and fic are both allowed to range between 0 and 1.

In addition, we here assume that the clumps become optically thin as the outflow merges with the interstellar medium, meaning that porosity effects become negligible. The merger occurs where the density of the outflow becomes comparable to that of the interstellar medium, between 1017 and 1018 cm depending on the mass-loss rate. We find that this assumption imposes a mass-loss rate dependent upper limit on h*. By requiring that the visual extinction of a clumpy outflow is equal to that of its corresponding smooth outflow near the merger, we find that h*≲ 1 × 1015 cm for = 1 × 10−5 M yr−1, h* ≲ 5 × 1015 cm for = 10−6 M yr−1, and h* ≲ 1 × 1016 cm for = 10−7 M yr−1. The assumption of optically thin clumps as the outflow merges with the interstellar medium hence introduces an extra constraint on l* and fvol. By using this maximal value for h*, a lower limit can be found for fvol assuming a certain value of l*, and vice versa.

4 Effect on chemistry

We have studied the effect of clumping and porosity on both an O-rich and C-rich outflow. In both cases, a one-component and a two-component outflow are considered. The only physical variation between the models are the mass-loss rate and the description of the specific clumpiness of the outflow, that is, the values of the parameters fvol, fic, and l*. Despite being widely detected, there are few observational constraints on the specific density structures observed. Imaging of 22 GHz water masers suggests the presence of clouds that are 10–100 times more dense than the outflow, with filling factors of less than one percent (Richards et al. 2012). From images of the molecular shell of IRC+10216, molecular clumps with sizes of 3.9 × 1015 cm up to 15.6× 1015 cm have been detected at a radial distance of about 24 arcsec or ~ 4.6 × 1016 cm, for a distance of 130 pc (Keller et al. 2018). The upper limit of this size range is set by the beam size of these observations, hence smaller clumps may exist as well. Applying our clump expansion assumption (Eq. (13)) on the above size range, we then have 4 × 1013 < l* < 1.5 × 1014 cm. Because the observational constraints on the specific clumpiness and porosity of observed outflows are very limited, we cannot estimate direct values for fvol, fic, and l*. We therefore explore their parameter spaces.

For the one-component outflow, we show the resulting abundance profiles and column densities for a selection of molecules for fvol = 0.05, 0.2, and 0.4, corresponding to an outflow where the clumps take up a small fraction, almost a quarter, and almost half of the total volume of the outflow. The parameter l* has values of 5 × 1012 cm, 1013 cm, and 5 × 1013 cm, corresponding to clumps with a characteristic length scale of 0.1 and 0.2 times the stellar radius, and to the extreme case where l* = R*. The corresponding porosity lengths h* for the different models are listed in Table 3. We note that the outflows with fvol = 0.2, l* = 5 × 1012 cm and fvol = 0.4, l* = 1 × 1013 cm have the same porosity length, but different clump overdensities.

Figures 4 and 6, and Tables 4 and 6 show the results for O- and C-rich one-component outflows, respectively.

For the two-component outflow, we show abundance profiles and column densities for fvol = 0.05, 0.2, and 0.4 and for fic = 0.1, 0.3 and 0.5 (Figs. 5 and 7, and Tables 5 and 7), corresponding to a very rarefied, a moderately rarefied, and an inter-clump component with a density equal to half the density of the corresponding smooth outflow, respectively. The parameter l* is now fixed at 1013 cm to show the effect of the parameter fic. The corresponding porosity lengths h* are listed in Table 3. Figures 5 and 7, and Tables 5 and 7 show the results for respectively O- and C-rich two-component outflows. The visual extinction and decrease in interstellar UV radiation field throughout these one- and two-component outflows are shown in Figs. D.1 and D.2, respectively.

There are several general trends with fvol, fic, and l* visible in the abundance profiles and column densities. Larger values of l*, or similarly h* = l*fvol, increase the production of unexpected species in the inner region, as they lead to a larger UV radiation field in the inner wind, inducing photodissociation of parent species. This results in a larger column density and a larger abundance in the inner region of the unexpected species. A larger value of l* or h* also shifts the photodissociation radius of all species in the outer region more inward. Smaller values of fic and fvol generally result in an increase in column density of the investigated molecules relative to the corresponding smooth outflow. This is due to both clumping and porosity. The increased overdensity of the clumps for smaller values of fic and fvol (Eq. (21)) causes the chemistry to occur at a faster rate and the abundances to increase, since the reaction rates are proportional to the product of the density of the reactants. Smaller values of fic and fvol also correspond to more UV radiation reaching the inner wind (Eq. (19), Figs. 1 and 3). “Highly porous” outflows, with both a large clump overdensity and porosity length, therefore have the largest effect on the circumstellar chemistry.

The effect of clumping and porosity extends of course beyond these general trends, affecting the chemical formation and destruction pathways of all species. This is reflected in the shape of their abundance profiles. In Sects. 4.1 and 4.2, we discuss the effect on the column densities and abundance profiles of some key unexpected molecules in detail for both the O- and C-rich outflows. These species have been chosen as they have been detected with an abundance larger than expected by TE in the inner region (see Sect. 1 for an overview), and as the changes to their chemical formation and destruction pathways induced by clumping and porosity are representative for the other species present in the outflow.

The species discussed in Sects. 4.1 and 4.2 have also been modelled by Agúndez et al. (2010). In order to facilitate a comparison, the column densities and abundance profiles of the other species modelled byAgúndez et al. (2010) can be found in Appendices E.1 and E.2, respectively, together with two other species that showcase the effects of clumping and porosity on the chemistry. The column densities and abundance profiles of so far undetected molecules that experience a large increase in predicted abundance due to a clumpy outflow, making them potentially detectable, are given in Appendix F.

thumbnail Fig. 4

Abundance of NH3 (upper panels), HCN (middle panels), and CS (lower panels) relative to H2 throughout a one-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively. We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

4.1 Oxygen-rich outflows

Under the assumption of TE chemistry, C- and N-bearing molecules are not expected to be abundantly present in the inner region of O-rich outflows (Cherchneff 2006; Agúndez et al. 2010; Gobrecht et al. 2016). In the following, we investigate the effect of clumping and porosity on the chemical pathways involved in the formation and destruction of NH3, HCN, and CS in the inner, intermediate, and outer regions of the stellar wind. The location of these regions depends on the mass-loss rate. The inner region ends roughly around ~ 1015 cm. The intermediate region ends where photodissociation becomes the dominant destruction process, around ~ 3× 1016 cm. These boundaries shift inward with decreasing mass-loss rate.

We compare our results to observed abundances in the O-rich AGBstars IK Tau, TX Cam, and R Dor. IK Tau and TX Cam have high mass-loss rates of ~ 5× 10−6 M yr−1 (Decin et al. 2010b) and ~ 3 × 10−6 M yr−1 (Bujarrabal et al. 1994), respectively. R Dor has a low mass-loss rate of ~ 1× 10−7 M yr−1 (Olofsson et al. 2002). Only a rough comparison to these observations is possible, since the physical parameters and input species of our chemical model are not tailored to these stars. We note that not all molecules have been observed in all three stars.

Table 4

Column density [cm−2] of NH3, HCN, and CS in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

4.1.1 NH3

The NH3 abundance in the inner region can reach up to a few times 10−8 relative to H2 in highly porous one-component outflows, an increase of three orders of magnitude relative to a smooth outflow. In such outflows, the NH3 column density can increase up to four orders of magnitude relative to a smooth outflow. The abundance profiles are shown in the upper panel of Figs. 4 and 5. The column densities are listed in Tables 4 and 5.

For highly porous outflows, the abundance of NH3 reaches its maximum value close to the starting radius (1014 cm) of the model and is constant throughout the inner and intermediate regions. This differs from a smooth outflow, or an outflow where the increase in the UV field in the inner region is not significant enough, where the NH3 abundance increases throughout the inner and intermediate regions.

The observed inner wind abundance of NH3 in IK Tau is 6 × 10−7 relative to H2 (Wong et al. 2018). In the outflows of other oxygen-rich evolved stars, NH3 has a similar abundance of ~ 10−6 (Wong et al. 2018). Our models do not produce such large abundances for any mass-loss rate considered, but including a highly porous clumpy outflow does lead to inner wind abundances about two orders of magnitude larger compared to a smoothoutflow. For high mass-loss rate outflows, a one-component outflow is needed to achieve such an increase.

Inner region

NH3 is predominantly produced via the hydrogenation of N through reactions with H2. In outflows with a well-shielded inner region, N is produced only through N2+He+N++N+e.\begin{equation*} {\textrm{N}_2 &#x002B; \textrm{He}^&#x002B; \rightarrow {\textrm{N}}^&#x002B; &#x002B; \textrm{N} &#x002B; \textrm{e}^-}. \end{equation*}

As He+ is only produced via cosmic-ray ionisation, the formation of N is independent of the UV field. However, for outflows with a larger UV field present in the inner region, either through a lower mass-loss rate or a larger porosity, N is also produced via photodissociation of the parent species N2. This yields a higher abundance of N, and subsequently of NH3. The effect of fvol, fic, and l* on the UV radiation field in the inner region is therefore strongly noticeable in Figs. 4 and 5.

Intermediate region

The sharp decline of NH3 at the end of the intermediate region is caused by photodissociation. In one-component outflows with larger mass-loss rates, the sharp decline shifts inward with increasing h*, equivalent to increasing l* and decreasing fvol. This is due to the strong decrease in extinction in this region (Fig. 1). Low mass-loss rate outflows and two-component outflows do not experience as large a decrease in visual extinction in this region. The onset of photodissociation for these outflows even moves outward with decreasing fvol, as the clumps are self-shielded.

4.1.2 HCN

In highly porous one-component outflows, the HCN abundance can reach up to a few times 10−7 relative to H2 in the inner region, an increase of ten orders of magnitude relative to a smooth outflow. The HCN column density can increase in such outflows up to four orders of magnitude relative to a smooth outflow. The abundance profiles are shown in the middle panel of Figs. 4 and 5. The column densities are listed in Tables 4 and 5.

The overdensity of the clumps increases the overall HCN abundance, but the strength of the UV radiation field in the inner region plays a crucial role in its formation. The importance of the UV radiation field can already be deduced by comparing the HCN abundance in smooth outflows for different mass-loss rates. Its abundance in the inner region is larger for lower mass-loss rates, as these outflows have an overall lower visual extinction.

HCN has been observed in IK Tau with an inner wind abundance of 4.4 × 10−7 by Decin et al. (2010b) and 9.8 × 10−7 by Bujarrabal et al. (1994). In TX Cam, it has been observed with a peak abundance of 2.2 × 10−6 (Bujarrabal et al. 1994). Its inner wind abundance in R Dor is 5.0 × 10−7 (Van de Sande et al. 2018). Our models with = 1 × 10−6 M yr−1 produce an abundance similar to the inner wind abundance of IK Tau and TX Cam only for highly porous one-component outflows. The observed inner wind abundance of R Dor is reproduced within a factor of five with = 1 × 10−7 M yr−1 by both highly porous one- and two-component outflows.

thumbnail Fig. 5

Abundance of NH3 (upper panels), HCN (middle panels), and CS (lower panels) relative to H2 throughout a two-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l*5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

Inner region

In the inner region, HCN is mainly produced via H2+CNHCN+H.\begin{equation*} {\textrm{H}_2 &#x002B; \textrm{CN} \rightarrow {\textrm{HCN}} &#x002B; \textrm{H}}. \end{equation*}

The abundance of CN depends on the UV radiation field in the inner region, as it is mainly produced via C+NOCN+O\begin{equation*} {\textrm{C} &#x002B; \textrm{NO} \rightarrow {\textrm{CN}} &#x002B; \textrm{O}} \end{equation*}

and N+CSCN+S.\begin{equation*} {\textrm{N} &#x002B; \textrm{CS} \rightarrow {\textrm{CN}} &#x002B; \textrm{S}}. \end{equation*}

The abundances of all four reactants increase with increased UV radiation field. C and N are produced in greater abundance due to the contribution of the photodissociation of the parent species CO and N2. NO and CS are produced via reactions between the parent species SO and N and C respectively, and are therefore also more abundant.

Moreover, an increased UV radiation field opens up additional formation pathways that contribute significantly to the formation ofHCN. The abundances of all reactants increase with the UV radiation field. The main additional reactions are N+CH2HCN+H,\begin{equation*} {\textrm{N} &#x002B; \textrm{CH}_2 \rightarrow {\textrm{HCN}} &#x002B; \textrm{H}}, \end{equation*} N+CH3HCN+H2,\begin{equation*} {\textrm{N} &#x002B; \textrm{CH}_3 \rightarrow {\textrm{HCN}} &#x002B; \textrm{H}_2}, \end{equation*}

and H+H2CNHCN+H2,\begin{equation*} {\textrm{H} &#x002B; \textrm{H}_2\textrm{CN} \rightarrow {\textrm{HCN}} &#x002B; \textrm{H}_2}, \end{equation*}

where H2CN is produced via N + CH3. CH2 and CH3 are produced via hydrogenation of C through reactions with H2 and are therefore more abundant as well.

Table 5

Column density [cm−2] of NH3, HCN, and CS in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

Intermediate region

The HCN abundance increases and forms a bump on the abundance profile near the end of the intermediate region. The feature is more prominent for higher mass-loss rates. In this region, HCN is mainly produced via CH+NOHCN+O.\begin{equation*} {\textrm{CH} &#x002B; \textrm{NO} \rightarrow {\textrm{HCN}} &#x002B; {O}}. \end{equation*}

The bump is visible in the abundance profile of CH, and is transferred to the profile of HCN. CH is produced through CH3+$_3^&#x002B;$ + e in this region, where CH3+$_3^&#x002B;$ is the product of hydrogenation of C+ through reactions with H2. In the inner region, C+ is predominantly produced via He+ + CO, and is therefore independent from the UV radiation field. The C+ production reaction shifts to the photoionisation of C near the end of the intermediate region. The shift to photoionisation of C occurs closer to the central star for larger UV radiation fields and increases the C+ abundance. This increase creates the bump feature, which is transferred to the abundance profile of HCN through CH3+$_3^&#x002B;$ and CH.

The end of the feature corresponds to the decrease in NO. In the intermediate reaction, NO is predominantly produced via OH+NNO+H,\begin{equation*} {\textrm{OH} &#x002B; \textrm{N} \rightarrow {\textrm{NO}} &#x002B; \textrm{H}}, \end{equation*}

where OH is produced through photodissociation of the parent species H2 O. The declinein NO abundance corresponds to the photodissociation radius of H2 O, decreasing the abundance of OH and therefore NO. The earlier onset of photodissociation of H2 O in highly porous outflows with a high mass-loss rate and low mass-loss rate outflows shifts the peak in NO abundance inward, and hence makes the bump in the HCN profile less pronounced.

Outer region

The slope of the decline in HCN abundance differs for different values of fvol. At the start of the outer region, reactions with C+ contribute significantly to the destruction of HCN, along with its photodissociation. The C+ abundance in this region is smaller for smaller values of fvol, as the increased self-shielding of the clumps makes the photoionisation of C less efficient. The smaller C+ abundance leads to a slower decline in HCN abundance.

4.1.3 CS

A CS abundance of up to a few times 10−7 relative to H2 is reached in the inner regions of highly porous one-component outflows, an increase of three orders of magnitude relative to a smooth outflow. For such outflows, the CS column density increases up to three orders of magnitude relative to a smooth outflow. The abundance profiles are shown in the lower panel of Figs. 4 and 5. The column densities are listed in Tables 4 and 5.

Both the overdensity of the clumps and the UV radiation field throughout the outflow have a large influence on the abundance of CS throughout the outflow. The overdensity of the clumps increases the overall abundance. For very overdense clumps, the abundance drops significantly in the intermediate region.

The CS abundance in the inner wind of IK Tau lies between 8.1 × 10−8 and 3.0 × 10−7 (Bujarrabal et al. 1994; Decin et al. 2010b; Kim et al. 2010). In TX Cam, its peak abundance lies between 1 × 10−6 and 5.0 × 10−7 (Lindqvist et al. 1988; Olofsson et al. 1991; Bujarrabal et al. 1994). Our models with = 1 × 10−6 M yr−1 have peak abundances similar to these observed values in both one- and two-component highly porous outflows. In models with = 1 × 10−5 M yr−1, such abundances are only produced in highly porous one-component outflows.

thumbnail Fig. 6

Abundance of NH3 (upper panels), H2O (middle panels), and H2S (lower panels) relative to H2 throughout a one-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively.We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

Inner region

The formation of CS in the inner region of the wind occurs predominantly through C+SOCS+O.\begin{equation*} {\textrm{C} &#x002B; \textrm{SO} \rightarrow {\textrm{CS}} &#x002B; \textrm{O}}. \end{equation*}

The abundance of C increases for larger UV fields in the inner region (see Sect. 4.1.2), increasing the CS abundance.

Intermediate region

The large dip in CS abundance is caused by its destruction through OH+CSOCS+H.\begin{equation*} {\textrm{OH} &#x002B; \textrm{CS} \rightarrow {\textrm{OCS}} &#x002B; \textrm{H}}. \end{equation*}

The region where this reaction is the main CS destruction mechanism corresponds to the peak in OH abundance, created by photodissociation of the parent species H2O. The CS abundance increases again as OH is photodissociated, which occurs closer to the central star for a larger porosity. The decrease in abundance is larger for smaller values of fvol, as the reaction can occur more efficiently due to the overdensity of the clumps.

Table 6

Column density [cm−2] of NH3, H2 O, and H2 S in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

Outer region

Before the photodissociation of CS, its abundance first increases. The associated feature is shifted inwards for highly porous outflows.

Near the end of the intermediate region and at the beginning of the outer region, CS is mainly produced via C+SOCS+O.\begin{equation*} {\textrm{C} &#x002B; \textrm{SO} \rightarrow {\textrm{CS}} &#x002B; O}. \end{equation*}

At the location of the trough, the reaction C+HSCS+H\begin{equation*} {\textrm{C} &#x002B; \textrm{HS} \rightarrow {\textrm{CS}} &#x002B; \textrm{H}} \end{equation*}

contributes to the production of CS. This region corresponds to a peak in HS abundance, caused by the reaction H+SiS+HS+Si+,\begin{equation*} {\textrm{H} &#x002B; \textrm{SiS}^&#x002B; \rightarrow {\textrm{HS}} &#x002B; \textrm{Si}^&#x002B;}, \end{equation*}

where SiS+ is mainly produced in this region by C++SiSSiS++C.\begin{equation*} {\textrm{C}^&#x002B; &#x002B; \textrm{SiS} \rightarrow {\textrm{SiS}}^&#x002B; &#x002B; \textrm{C}}. \end{equation*}

The location of the peak in abundance can hence be traced back to the abundance of C+ (see Sect. 4.1.2).

4.2 Carbon-rich outflows

In the inner region of C-rich outflows, O- and N-bearing molecules are not expected to be abundantly present under the assumption of TE chemistry (Cherchneff 2006; Agúndez et al. 2010). In the following, we investigate the effect of clumping and porosity on the chemical pathways involved in the formation and destruction of NH3, H2 O, and H2 S in the the inner, intermediate, and outer regions of the stellar wind. As with the O-rich outflows, the location of these regions depends on the mass-loss rate. The inner region ends roughly around ~ 1015 cm. The end of the intermediate region, around ~ 3 × 1016 cm, corresponds to where photodissociation becomes the dominant destruction process. These boundaries shift inward with decreasing mass-loss rate.

We compare our results to the observed abundances in the well-studied C-rich AGB star IRC+10216 (CW Leo), which has a high mass-loss rate of 2–4 × 10−5 M yr−1 (De Beck et al. 2012; Cernicharo et al. 2015). Only a rough comparison to these observations is possible, since the physical parameters and input species of our chemical model are not tailored to IRC+10216.

4.2.1 NH3

As in O-rich chemistry, the inner region NH3 abundance in C-rich outflows can reach up to a few times 10−8 relative to H2 in highly porous outflows. This corresponds to an increase of three orders of magnitude relative to a smooth outflow. The NH3 column density of such outflows can increase up to four orders of magnitude relative to a smooth outflow. The abundance profiles are shown in the upper panel of Figs. 6 and 7. The column densities are listed in Tables 6 and 7.

For highly porous outflows, the NH3 abundance reaches its maximum value close to the start of the model and is constant throughout the inner and intermediate regions, similar to the O-rich outflow. The destruction of NH3 at the end of the intermediate region is influenced by the overdensity of the clumps and the UV radiation field present.

In IRC+10216, NH3 has an observed inner wind abundance of 1.7 × 10−7 (Keady & Ridgway 1993). Similar to the O-rich outflows, our models do not produce such a large abundance, but including a highly porous clumpy outflow does increase the inner wind abundance significantly. For = 10−5 M yr −1, a one-component outflow is again needed to achieve such an increase.

thumbnail Fig. 7

Abundance of NH3 (upper panels), H2O (middle panels), and H2S (lower panels) relative to H2 throughout a two-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

Inner region

The formation of NH3 in the inner region is similar to that in an O-rich outflow. The extra N-bearing parent species HCN does not influence the formation of NH3 in the inner region, as its photodissociation products are CN and H.

Intermediate region

For higher mass-loss rate outflows with a large porosity, the NH3 abundance declines rapidly near the end of the intermediate region. The decline is due to the onset of the reaction NH3+CNHCN+NH2,\begin{equation*} {\textrm{NH}_3 &#x002B; \textrm{CN} \rightarrow {\textrm{HCN}} &#x002B; \textrm{NH}_2}, \end{equation*}

which contributes significantly to its destruction as the CN abundance peaks in this region. The peak in CN abundance is caused by the photodissociation of HC3N, which is formed mainly through CN+C2H2HC3N+H.\begin{equation*} {\textrm{CN} &#x002B; \textrm{C}_ 2\textrm{H}_2 \rightarrow {\textrm{HC}}_3\textrm{N} &#x002B; \textrm{H}}. \end{equation*}

Since CN is produced through photodissociation of the parent species HCN, the peak in abundance of its daughter species HC3N shifts closer to the central star for larger porosity lengths, as does the photodissociation of HC3N. The decline in NH3 abundance ends when the CN abundance is reduced through the reaction O+CNN+CO,\begin{equation*} { \textrm{O} &#x002B; \textrm{CN} \rightarrow {\textrm{N}} &#x002B; \textrm{CO}}, \end{equation*}

as the O abundance increases due to the photodissociation of the parent species SiO.

Highly porous outflows with a higher mass-loss rate also show an increase in NH3 abundance after its sharp decline near a few times 1016 cm. In this region, NH3 is predominantly produced through reactions with NH4+$_4^&#x002B;$, which originates from the hydrogenation of N+ through reactions with H2. N+ is solely formed through reactions with He+ with N2 and is therefore unaffected by porosity. However, as the UV radiation field increases in this region, the abundances of the intermediary products NH+, NH2+$_2^&#x002B;$, and NH3+$_3^&#x002B;$ increase through photoionisation of their neutral counterparts. This leads to an increased abundance of NH4+$_4^&#x002B;$ and subsequently of NH3.

Table 7

Column density [cm−2] of NH3, H2 O, and H2 S in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

Outer region

At the start of the outer region, reactions with C+ contribute to the destruction of NH3. The different slope of the decline abundance with different values of fvol at the beginning of the outer region has a similar origin as that of HCN in the O-rich outflow (Sect. 4.1.2).

4.2.2 H2O

The H2O abundance can reach up to 10−5 relative to H2 in the inner region in highly porous one-component outflows and in highly porous two-component outflows with a small mass-loss rate, an increase of three orders of magnitude relative to a smooth outflow. In such outflows, the H2 O column density can increase up to three orders of magnitude relative to a smooth outflow. The abundance profiles are shown in the middle panel of Figs. 6 and 7. The column densities are listed in Tables 6 and 7.

The peak in abundance around ~1–3 ×1016 cm, visible in the smooth abundance profile, disappears as fvol decreases. Outflows with very overdense clumps show a steep decline in this region, followed by a plateau-feature.

The inner wind H2O abundance of IRC+10216 lies around 1 × 10−7 (Decin et al. 2010a). Lombaert et al. (2016) find that the inner wind abundance lies between ~ 10−6–10−4 for a sample of C-rich stars with mass-loss rates between ~ 10−7 and 10−5 M yr−1. Only highly porous one-component outflows with = 10−5 M yr−1 produce an inner wind abundance close to these observations. Models with = 10−6 and 10−7 M yr−1 produce inner wind abundances within the range probed by Lombaert et al. (2016) for both one-component and two-component outflows, where the models have a maximum abundance of ~ 1 × 10−5.

Inner region

H2O is formed in the inner region through hydrogenation of O through reactions with H2, similar to the formation ofNH3. O is produced in greater abundance in outflows with a large UV radiation field in the inner region, either through a low mass-loss rate or a large porosity, as the parent species CO and SiO are photodissociated in this region.

Intermediate region

In the intermediate region, the formation of water shifts from a chemistry initiated by the proton transfer of H3+$_3^&#x002B;$ with O atoms to one in which the dissociative recombination of CH2CO+ dominates CH2CO++eH2O+C2.\begin{equation*} {\textrm{CH}_2\textrm{CO}^&#x002B; &#x002B; \textrm{e}^- \rightarrow {\textrm{H}}_2\textrm{O} &#x002B; \textrm{C}_2}. \end{equation*}

However, there are no laboratory measurements of the dissociative recombination of CH2CO+. The channel forming H2O is one of three equal channels suggested by Prasad & Huntress (1980). Given the re-arrangement needed to get H2 O out of the ion, which is structured as H2 CCO+, it may have zero probability of occurring.

Nevertheless, the CH2CO+ abundance peaks near the end of the intermediate region. It is formed by the photoionisation of CH2CO, which also peaks in the same region. The peak in abundance is very sensitive to the temperature of the outflow, as several reactions with a negative temperature dependency play a large role in the chain of reactions forming CH2CO: C2H2+hνC2H2++e,\begin{equation*} {\textrm{C}_2\textrm{H}_2 &#x002B; {\textrm{h}\nu} \rightarrow {\textrm{C}}_2\textrm{H}_2^&#x002B; &#x002B; \textrm{e}^-}, \end{equation*} H2+C2H2+C2H4++hν  (β=2.01),\begin{equation*} {\textrm{H}_2 &#x002B; \textrm{C}_2\textrm{H}_2^&#x002B; \rightarrow {\textrm{C}}_2\textrm{H}_4^&#x002B; &#x002B; {\textrm{h}\nu}}\ \ \ (\beta = -2.01), \end{equation*} C2H4++OCH3++HCO,\begin{equation*} {\textrm{C}_2\textrm{H}_4^&#x002B; &#x002B; \textrm{O} \rightarrow {\textrm{CH}}_3^&#x002B; &#x002B; \textrm{HCO}}, \end{equation*} CH3++COCH3CO++hν   (β=1.2),\begin{equation*} {\textrm{CH}_3^&#x002B; &#x002B; \textrm{CO} \rightarrow {\textrm{CH}}_3\textrm{CO}^&#x002B; &#x002B; {\textrm{h}\nu} }\ \ \ (\beta = -1.2), \end{equation*} CH3CO++eCH2CO+H,\begin{equation*} {\textrm{CH}_3\textrm{CO}^&#x002B; &#x002B; \textrm{e}^- \rightarrow {\textrm{CH}}_2\textrm{CO} &#x002B; \textrm{H}}, \end{equation*}

with β the exponent denoting the temperature dependency of the reaction rate in the modified Arrhenius formula (McElroy et al. 2013). For larger values of the exponent ϵ of the temperature power-law, corresponding to a lower temperatures throughout the outflow, the peak disappears.

For large UV radiation fields, the CH2CO+ peak shifts closer to the central star and becomes comparable in abundance to H3 O+, which removes the peak in H2O abundance.

Outer region

At the start of the outer region, reactions with C+ contribute to the destruction of H2O. The effect of fvol on the slope of the decline in H2O abundance at the beginning of the outer region has a similar origin as that of NH3 and HCN inthe O-rich outflow (Sect. 4.1).

4.2.3 H2S

In highly porous one-component outflows, the inner wind abundance of H2 S can reach up to 10−8 relative to H2, an increase of three orders of magnitude relative to a smooth outflow. The H2 S column density in such outflows can increase up to four orders of magnitude relative to a smooth outflow. The abundance profiles are shown in the lower panel of Figs. 6 and 7. The column densities are listed in Tables 6 and 7.

For one-component outflows with fvol = 0.05, the abundance shows a large dip in the inner and intermediate regions. The H2 S abundance profile also shows a plateau-feature after a sharp decline around 1016 cm.

The inner wind abundance of H2S in IRC+10216 is 4 × 10−9 (Agúndez et al. 2012). In models with = 10−5 M yr−1, such an inner wind abundance is produced in highly porous one-component outflows.

Inner region

In the inner region, H2S is formed through the sequential hydrogenation of S through reactions with H2. In outflows with a well-shielded inner region, S is formed via cosmic-ray induced photodissociation of the parent species CS and SiS. As the UV radiation field increases, the S abundance increases as these parent species are also photodissociated. The increased N abundance through photodissociation of the parent species N2 moreover leads to an increased contribution to S through the reaction N+CSS+CN.\begin{equation*} {\textrm{N} &#x002B; \textrm{CS} \rightarrow {\textrm{S}} &#x002B; \textrm{CN}}. \end{equation*}

For outflows with a large porosity, the reaction N+HSNS+H\begin{equation*} {\textrm{N} &#x002B; \textrm{HS} \rightarrow {\textrm{NS}} &#x002B; \textrm{H}} \end{equation*}

destroys a fraction of the HS after ~ 2 × 1014 cm. The onset of the reaction corresponds to the start of the decrease of HS abundance, and coincides with the increase in N abundance.

Intermediate region

The decrease in H2S abundance before the peak at the end of the intermediate region is also caused by the photodissociation of the parent species CS. It therefore has the same origin as the extended dip in H2 S abundance for outflows with a larger porosity, but is more spatially confined due to the later onset of the CS photodissociation.

The peak in H2S abundance corresponds to the contribution of the reaction H+SiS+HS+Si+\begin{equation*} {\textrm{H} &#x002B; \textrm{SiS}^&#x002B; \rightarrow {\textrm{HS}} &#x002B; \textrm{Si}^&#x002B;} \end{equation*}

to the H2S production, as the abundance of SiS+ reaches its maximal value in this region and forms a peak at the same location as that in the HS abundance, and consequently also the H2 S abundance. The SiS+ peak is caused by reactions between C+ and the parent species SiS, and therefore shifts inward for smaller values of fvol (see Sect. 4.1.2). The SiS+ abundance then decreases through dissociative recombination with electrons, marking the end of the peak in abundance.

The origin of the peak in H2S abundance is similar to that in the H2O abundance (Sect. 4.2.3). The rates of the reactions between SiS and C+ and S+ have a negative temperaturedependency (β = −0.5 for both). Again, the peak disappears for larger values of the exponent ϵ of the temperature power-law, corresponding to lower temperatures.

As the mass-loss rate decreases and/or the porosity of the outflow decreases, the peak in the H2 S abundance disappears. This is due to the earlier onset of its photodissociation in these outflows.

Outer region

The different slope of the decline in H2S abundance with different values of fvol at the beginning of the outer region is again due to the C+ abundance (see e.g. Sect. 4.1.1). In this region, H2 S is formed via reactions between H2 S+ and various anions of C-chains. The different slopes are also present in the decline of the H2 S+ abundance, as it is mainly produced by the reaction H2CO+S+H2S++CO,\begin{equation*} {\textrm{H}_2\textrm{CO}&#x002B; \textrm{S}^&#x002B; \rightarrow {\textrm{H}}_2\textrm{S}^&#x002B; &#x002B; \textrm{CO}}, \end{equation*}

where S+ is produced via the photoionisation of S. Before H2CO is photodissociated, it is destroyed through H2CO+C+CO+CH2+,\begin{equation*} {\textrm{H}_2\textrm{CO}&#x002B; \textrm{C}^&#x002B; \rightarrow {\textrm{CO}} &#x002B; \textrm{CH}_2^&#x002B;}, \end{equation*}

causing different slopes for different values of fvol, which are transferred to the H2S abundance profile.

5 Discussion

In Sect. 5.1, we discuss the general influence of the parameters fvol, fic, and l* on the chemistry throughout the outflow. We compare our results to those of Agúndez et al. (2010), who have also modelled the effect of a clumpy outflow, in Sect. 5.2, and to those of shock-induced non-equilibrium chemistry models in Sect. 5.3. We compare our modelled inner wind abundances to observations in Sect. 5.4. Finally, molecules that experience a large increase in abundance or change in their abundance profiles due to a clumpy outflow are discussed in Sect. 5.5.

5.1 Influence of the parameters fvol, fic, and l*

The porosity of the outflow is determined by the porosity length h* = l*fvol and fic. The porosity length has the largest influence on the decrease in visual extinction, for a certain value of fic. For a certain value of h*, the visual extinction increases with fic. A larger value of fic corresponds to more material in the inter-clump component, hindering radiation travelling through the porous channels between the clumps. An increase in the UV radiation field in the inner region leads to the photodissociation of parent species in this region, which is well-shielded in a smooth outflow. The elements which are normally locked up in stable parent species in smooth outflows (O in C-rich outflows, C in O-rich outflows, and N in both) are therefore released closer to the central star, allowing for an increased formation of the unexpected species in the inner region. These species are, under the assumption of TE chemistry, expected either to be absent in the inner wind or to have a smaller abundance than observed in the outer regions. More generally, an increased UV radiation field influences the chemistry throughout the outflow, as photodissociation of both parent and daughter species occurs closer to the central star, opening up new chemical pathways.

Clumping also influences the overall abundance of the unexpected species. For the same value of h*, smaller values of fvol and fic will generally lead to larger inner region abundances and column densities. The effect of porosity on the increase in abundance is however much larger. Highly porous outflows, characterised by both a large clump overdensity and a large porosity length, therefore have the largest influence on the chemistry throughout the outflow.

Since the inner wind abundances of the unexpected species are mainly determined by the porosity of the outflow, one-component outflows yield a larger increase in abundance than their corresponding two-component outflows. For higher mass-loss rate outflows, the inter-clump component needs to be effectively void in order for porosity to increase the inner wind abundance of the unexpected species.

5.2 Comparison to the chemical model of Agúndez et al. (2010)

The effect of a clumpy outflow on its chemistry has previously been studied by Agúndez et al. (2010). In their model, the outflow is also split into two components: a smooth well-shielded component and a clumpy component, for which a fraction of solid angle of arrival of interstellar photons fΩ is free of matter. The smooth component experiences the same UV radiation field as a uniform outflow, whereas a certain fraction of photons is allowed to reach the inner regions of the clumpy component uninhibited. During the modelling, both components have the same density, namely that of the corresponding smooth outflow. The final abundances are calculated by taking a weighted average of the models of the components using the mass fraction of the total circumstellar mass that is located inside the clump component fm.

The model of Agúndez et al. (2010) hence differs from our model in both its implementation of the density distribution and the alteration of the UV radiation field. Whereas we take the density distribution into account while calculating the chemical network, Agúndez et al. (2010) include it via post-processing. Both components experience the same effective UV radiation field in our model, while the components of Agúndez et al. (2010) experience two distinct UV radiation fields. Finally, the chemical reaction network used by Agúndez et al. (2010) is different to the UDfA RATE12 reaction network.

Agúndez et al. (2010) assumed that fm = 0.1 and fΩ = 0.25, chosen in order to roughly model the clumpy outflow of the C-rich star CW Leo. The smooth component is hence a major contributor to their final abundances, while at the same time a significant fraction of the interstellar radiation is still able to reach the inner wind. This combination is not possible in our model, as a larger density of the inter-clump component leads to a smaller decrease in visual extinction. Translated into our parameters, the fm and fΩ parameters would correspond to fvol ≈ 1 − fΩ = 0.75 and fic ≈ 1 − fm = 0.9. However, in order to have a similar effect on the UV radiation field, we would require fic = 0. It is therefore not possible to make a direct comparison to one of the models shown in Sect. 4.

Nevertheless, we can perform a qualitative comparison and point out several clear differences. Firstly, the shapes of the abundance profiles of Agúndez et al. (2010) are only roughly similar to those shown in this paper, with a large discrepancy for H2CO, C2H2, HCN, and CS in the O-rich outflow, and for OH and H2S in the C-rich outflow. Secondly, the modelled inner wind abundances of Agúndez et al. (2010) can only be obtained with our model by assuming a one-component outflow with a large porosity length. In that case, a large fraction of UV photons is able to reach the inner region of these outflows, which matches the assumption by Agúndez et al. (2010). Moreover, the one-component outflow that reproduces the inner wind abundances of NH3 in the O-rich outflow, yields abundances up to an order of magnitude higher for all other species shown in Agúndez et al. (2010). Similar for the C-rich outflow, where the one-component model that reproduces the inner wind abundances of NH3 and H2 S yields abundances an order of magnitude larger for the other molecules shown.

5.3 Comparison to non-equilibrium chemistry models

In order to explain the presence and abundance of the unexpected species in the inner wind, shock-induced non-equilibrium chemistry in the inner region has been proposed as a possible solution. These models are able to reproduce the observed abundances in the inner regions, except for NH3, which is predicted to have an abundance of ~ 10−19 relative to H2 (Willacy & Cherchneff 1998; Duari et al. 1999; Gobrecht et al. 2016). Although a direct comparison is not possible as we have not specifically modelled IK Tau, we do produce NH3 with a larger abundance than TE. In non-equilibrium models, CH4 is not expected to be present with an abundance larger than predicted by TE chemistry in O-rich outflows (Duari et al. 1999). Our models produce an inner wind abundance up to seven order of magnitude larger than predicted for this species (Figs. E.1 and E.2).

The inner regions of C-rich outflows have not been studied using non-equilibrium chemical models to the same extent as O-rich outflows. However, the recently discovered distribution of CH3CN in the C-rich AGB star CW Leo does not match the predictions of both thermodynamic equilibrium and non-equilibrium chemistry (Agúndez et al. 2015; Millar 2016). Agúndez et al. (2015) suggest a clumpy outflow as a possible origin to this distribution. Similar to their best fitting abundance profile, we find a peak in the CH3CN abundance around 1016 cm (Figs. E.3 and E.4). Moreover, the presence of H2 O in the inner region of CW Leo was not predicted to be present by non-equilibrium chemistry (Cherchneff 2006), and could only be modelled by manually adapting the rate of the reaction between SiO and H (Cherchneff 2012). Although a direct comparison to the observations is again not possible, our models do produce H2 O in larger abundances than TE in the inner region.

5.4 Comparison to observations

For both O-rich and C-rich outflows, we find that the observations are best reproduced by highly porous, one-component outflows. A large UV field is therefore required to be present in the inner region for the deficient elements to be released in order to obtain abundances for the unexpected species similar to those observed. The overdensity of the clumps leads to a further increase in abundance, so that both porosity and clumping are critical in reproducing the observed abundances. For lower mass-loss rates, the requirement of a highly porous, one-component outflow is less strict. The UV field in the inner region is still sufficiently large in these outflows if a fraction of the material resides in the inter-clump component.

The observed NH3 abundances in both O-rich and C-rich AGB stars are not reproduced by our models. Our models also do not produce the upper limit of the observed H2 O abundance range in C-rich AGB stars. Despite this, we do find an increase of up to four orders of magnitude compared to a smooth outflow in the inner wind abundance of NH3, and up to three orders of magnitude in the inner wind abundance of H2 O for highly porous, one-component outflows. However, an effectively void inter-clump medium is likely not realistic for the outflows of AGB stars. The density contrast between the clumped and inter-clump component of IRC+10216 is found to lie between 1.5 and 7.6 (Mauron & Huggins 2000; Dinh-V-Trung & Lim 2008), its inter-clump component can therefore not be considered effectively void.

The requirement of an effectively void inter-clump component might be due to assumption of optically thin clumps at the edge of the outflow, which places an upper limit on the porosity length (Sect. 3.4). If the outer boundary condition were less strict, larger porosity lengths would be possible, potentially allowing for an increased inter-clump density.

5.5 Predictability of models

The clumpy outflow model can also be used to predict the effect of clumping and porosity on the abundances of unobserved molecules, making them potentially observable, and their effect on the abundance profiles of already detected molecules. We have identified a set of molecules whose abundance is strongly influenced by a clumpy outflow by requiring a peak fractional abundance larger than 1 × 10−9 relative to H2 and an increase in column density relative to the smooth outflow of at least two orders of magnitude. Their abundance profiles and column densities are shown in Appendix F. For the O-rich outflow, we find considerably changed abundance profiles and column densities for N2O, C2 N, C3 H, C3 H2 and OCS (Figs. F.1 and F.2, Tables F.1 and F.2). For the C-rich outflow, considerable changes are found for CO2, SO2, HC9N, NO and OCS (Figs. F.3 and F.4, Tables F.3 and F.4).

6 Conclusion

We have determined the effect of a non-uniform, clumpy AGB outflow on its chemistry by implementing a porosity formalism in our one-dimensional chemical model. This theoretical formalism takes into account both the effect of porosity on the UV radiation field throughout the outflow (by adjusting its optical depth), and the effect of clumping in general on the reaction rates for two-body processes (by splitting the outflow into an overdense clumped and a rarified inter-clump component). The specific clumpiness is determined by three parameters: the clump volume filling factor fvol, the density contrast fic, and the characteristic size of the clumps at the stellar surface l*. The first two parameters determine the relative densities of the clumped and inter-clump component. The porosity length h* = l*fvol and fic determine the decrease in visual extinction throughout the outflow, with h* having the largest influence. In this paper, the clumps are assumed to have a spatially constant fvol, which implies that they expand as they move outward. We note that, being a statistical description of a stochastic ensemble of clumps, the porosity formalism does not require any assumptions on the specific location of the clumps.

We have shown results for a one-component outflow, where all material is located inside clumps and the inter-clump component is effectively void, and a two-component outflow, where material is located in both the clumped and inter-clump component. Both the overdensity of the clumps and the increase in interstellar UV radiation in the inner regions of the outflow alter the chemistry throughout the outflow, as they alter the abundances and the chemical formation and destruction pathways. The decrease in visual extinction has the largest influence on the formation of the unexpected species in both O- and C-rich outflows. Porosity, and not clumping, is therefore the main mechanism determining the overall increase in abundance. The larger UV radiation field in the inner region liberates the deficient elements through photodissociation of parent species, enabling the formation of the unexpected species. A highly porous outflow, characterised by a large clump overdensity and a large porosity length, therefore has the largest effect on the chemistry throughout the outflow.

Our chemical model is able to produce all unexpected species in clumpy O- and C-rich outflows with inner region abundances larger than expected by TE chemistry. Their column densities and molecular abundances in the inner regionincrease with clump overdensity and porosity length, and decrease with inter-clump density contrast. Unlike the shock-induced non-equilibrium chemistry models, our models are able to produce NH3 and H2 O in the inner region with a larger abundance than predicted by TE. The differences between our results and those of Agúndez et al. (2010), who have also implemented a clumpy outflow in their model, are most likely due to a different modelling approach. However, both approaches are able to produce all unexpected species in the inner region.

Highly porous, one-component outflows generally reproduce the observationally retrieved values of the unexpected species in higher mass-loss rate outflows. For lower mass-loss rate outflows, the inter-clump component does not need to be effectively void. Nonetheless, such models do not produce the observed inner wind abundance of NH3 in O- and C-rich outflows and the upper limit of the range of H2O in C-rich outflows. Combining both shock-induced non-equilibrium chemistry and a clumpy outflow can potentially solve the discrepancy.

Chemical models that include a non-uniform density structure can lead to an increase in the abundances of the unexpected species in the inner regions of both O- and C-rich outflows, as demonstrated in this paper. This emphasises the importance of including the effect of a clumpy outflow in the chemical modelling of AGB outflows. Even though it cannot account for the observed abundance of all unexpected species investigated in this paper, the non-uniform density structure should be taken into account in order to accurately model the chemistry and to correctly interpret the observed abundance profiles.

Acknowledgements

MVdS, DK, WH, and LD acknowledge support from the ERC consolidator grant 646758 AEROSOL. TJM thanks STFC for support under grant reference ST/P000312/1. FDC acknowledges support from the EPSRC iCASE studentship programme, Intel Corporation and Cray Inc. We thank the anonymous referee for a detailed reading and helping to improve this paper.

Appendix A Opacity of the smooth outflow

The opacity of the smooth outflow is given by κ=1.0862[AUV/AV]1.87×10211μmH1200cm2g1,\begin{equation*} \kappa = {1.086}\ \frac{{2}\ \lbrack A_{UV}/A_V\rbrack}{1.87 \times 10^{21}}\ \frac{1}{\mu\ m_H} \approx {1200}\ \textrm{cm}^2 \ \textrm{g}^{-1} , \end{equation*}(A.1)

with μ = 2.68 the mean molecular mass of the outflow relative to H2, including the fractional abundance of He (Table 2), and mH the atomic massunit. We hence assume that the extinction throughout the outflow is equal to that throughout the ISM of 1.87× 1021 atoms cm−2 mag−1 (Cardelli et al. 1989), which is then scaled according to the ratio of UV and visual extinction, with AUVAV = 4.65 the adopted ratio of UV and visual extinction (Nejad et al. 1984). The factor of 1.086 originates from the conversion of extinction to optical depth.

Appendix B Figures of the decrease in interstellar UV radiation field

Figures B.1B.3 show the decrease in interstellar UV radiation field I(r)∕I0 = eτeff(r) throughout the outflow, with I the UV radiation field throughout the outflow and I0 the interstellar UV radiation field.

thumbnail Fig. B.1

Decrease in interstellar UV radiation field throughout a one-component outflow with = 10−5 M yr−1 (left panel) and = 10−7 M yr−1 (right panel). Results for a smooth, uniform outflow are shown in black. Results for different values of h* are shown in blue. For reference, 1 R* = 5 × 1013 cm.

thumbnail Fig. B.2

As Fig. B.1, but for a two-component outflow with fic = 0.1. For reference, 1 R* = 5 × 1013 cm.

thumbnail Fig. B.3

Decrease in interstellar UV radiation field throughout a two-component outflow with h = 1015 cm and and = 10−5 M yr−1 (left panel) and = 10−7 M yr−1 (right panel). Results for a smooth, uniform outflow are shown in black. Results for different values of fic are shown in blue. For reference, 1 R* = 5 × 1013 cm.

Appendix C Expression of the effective optical depth

C.1 One-component outflow

The integral describing the effective optical depth τeff for the one-component medium (Eq. (17)) can be solved analytically. We have that τeff=τ*R*r1r2+τ*R*1/3h*r2/3 dr=3πτ*1/4R*3/422h*3/4+3τ*1/4R*3/442h*3/4×(2arctan(12r1/3(τ*R*1/3h*)1/4)2arctan(1+2r1/3(τ*R*1/3h*)1/4)+log(τ*R*1/3h*2(τ*R*1/3h*)1/4r1/3+r2/3)log(τ*R*1/3h*+2(τ*R*1/3h*)1/4r1/3+r2/3)). \begin{equation*}\begin{split} \tau_{\textrm{eff}} &= \tau_*\ R_* \int_r^{\infty} \frac{1}{r^2 &#x002B;\tau_*\ R_*^{1/3}\ h_*\ r^{2/3}}\ \textrm{d}r \\ &= \frac{3\pi \ \tau_*^{1/4} R_*^{3/4}}{2\sqrt{2}\ h_*^{3/4}} &#x002B; \frac{3\ \tau_*^{1/4} R_*^{3/4}}{4\sqrt{2}\ h_*^{3/4}} \times \ \Bigg(2\ \textrm{arctan} \left( 1 - \frac{\sqrt{2}\ r^{1/3}}{\left(\tau_*\ R_*^{1/3}h_*\right)^{1/4}} \right) \\ &\quad - 2\ \textrm{arctan} \left( 1 &#x002B; \frac{\sqrt{2}\ r^{1/3}}{\left(\tau_*\ R_*^{1/3}h _*\right)^{1/4}} \right) \\ &\quad&#x002B; \textrm{log}\left( \sqrt{\tau_*\ R_*^{1/3}h_*} - \sqrt{2}\ \left(\tau_*\ R_*^{1/3}h_*\right)^{1/4}\ r^{1/3} &#x002B; r^{2/3} \right) \\ &\quad- \textrm{log}\left( \sqrt{\tau_*\ R_*^{1/3}h_*} &#x002B; \sqrt{2}\ \left(\tau_*\ R_*^{1/3}h_*\right)^{1/4}\ r^{1/3} &#x002B; r^{2/3} \right) \Bigg). \end{split} \end{equation*}(C.1)

Since τ = τ* R*r, this can be rewritten as τeff=3π(τr)1/4R*1/222h*3/4+3(τr)1/4R*1/242h*3/4×(2arctan(12r1/12(τh*R*2/3)1/4)arctan(1+2r1/12(τh*R*2/3)1/4)+log(τh*R*2/3r2(τh*R*2/3)1/4r7/12+r2/3)log(τh*R*2/3r+2(τh*R*2/3)1/4r7/12+r2/3)). \begin{equation*}\begin{split} \tau_{\textrm{eff}} &= \frac{3\pi\ (\tau\ r)^{1/4} R_*^{1/2}}{2\sqrt{2}\ h_*^{3/4}} &#x002B; \frac{3\ (\tau\ r)^{1/4} R_*^{1/2}}{4\sqrt{2}\ h_*^{3/4}}\ \\ &\quad \times\Bigg( 2\ \textrm{arctan} \left( 1 - \frac{\sqrt{2}\ r^{1/12}}{\left(\tau\ h_*\ R_*^{-2/3}\right)^{1/4}} \right) \\ &\quad - \textrm{arctan} \left( 1 &#x002B; \frac{\sqrt{2}\ r^{1/12}}{\left(\tau\ h_*\ R_*^{-2/3}\right)^{1/4}} \right) \\ &\quad &#x002B; \textrm{log}\left( \sqrt{\tau\ h_*\ R_*^{-2/3}\ r} - \sqrt{2}\ \left(\tau\ h_*\ R_*^{-2/3}\right)^{1/4}\ r^{7/12} &#x002B; r^{2/3}\right) \\ &\quad - \textrm{log}\left( \sqrt{\tau\ h_*\ R_*^{-2/3}\ r} &#x002B; \sqrt{2}\ \left(\tau\ h_*\ R_*^{-2/3}\right)^{1/4}\ r^{7/12} &#x002B; r^{2/3}\right) \Bigg). \end{split} \end{equation*}(C.2)

C.2 Two-component outflow

The integral describing the effective optical depth τeff for the two-component medium (Eq. (20)) can be solved analytically. We have that τeff=τ*R*r1+τ*R*1/3h*Br4/3r2+τ*R*1/3h*Br2/3 dr=τ*R*ficr+3πτ*1/4R*3/4(1fic)22(h*B)3/4+3τ*1/4R*3/4(1fic)42(h*B)3/4×(2arctan(12r1/3(τ*R*1/3h*B)1/4)2arctan(1+2r1/3(τ*R*1/3h*B)1/4)+log(τ*R*1/3h*B2(τ*R*1/3h*B)1/4r1/3+r2/3)log(τ*R*1/3h*B+2(τ*R*1/3h*B)1/4r1/3+r2/3)), \begin{equation*}\begin{split} \tau_{\textrm{eff}} &= \tau_*\ R_* \int_r^{\infty} \frac{1&#x002B;\tau_*\ R_*^{1/3}\ h_*\ B\ r^{-4/3}}{r^2 &#x002B;\tau_*\ R_*^{1/3}\ h_*\ B\ r^{2/3}}\ \textrm{d}r \\ &= \frac{\tau_*\ R_*\ f_{\textrm{ic}}}{r} &#x002B; \frac{3\pi \ \tau_*^{1/4}\ R_*^{3/4} \left(1-f_{\textrm{ic}}\right)}{2\sqrt{2}\ (h_*\ B)^{3/4}}\\ &\quad &#x002B; \frac{3\ \tau_*^{1/4}\ R_*^{3/4} \left(1-f_{\textrm{ic}}\right)}{4\sqrt{2}\ (h_*\ B)^{3/4}} \times \ \Bigg(2\ \textrm{arctan} \left( 1 - \frac{\sqrt{2}\ r^{1/3}}{\left(\tau_*\ R_*^{1/3}h_*\ B\right)^{1/4}} \right) \\ &\quad - 2\ \textrm{arctan} \left( 1 &#x002B; \frac{\sqrt{2}\ r^{1/3}}{\left(\tau_*\ R_*^{1/3}h_*\ B\right)^{1/4}} \right) \\ &\quad&#x002B; \textrm{log}\left( \sqrt{\tau_*\ R_*^{1/3}h_*\ B} - \sqrt{2}\ \left(\tau_*\ R_*^{1/3}h_*\ B\right)^{1/4}\ r^{1/3} &#x002B; r^{2/3} \right) \\ &\quad - \textrm{log}\left( \sqrt{\tau_*\ R_*^{1/3}h_*\ B} &#x002B; \sqrt{2}\ \left(\tau_*\ R_*^{1/3}h_*\ B\right)^{1/4}\ r^{1/3} &#x002B; r^{2/3} \right) \Bigg), \end{split} \end{equation*}(C.3)

with B=1(1fvol)fic$B = 1 - \left( 1-f_{\textrm{vol}} \right) f_{\textrm{ic}}$. Since τ = τ* R*r, this can be rewritten as τeff=τfic+3π(τr)1/4(1fic)22(h*R*2/3B)3/4+3(τr)1/4(1fic)42(h*R*2/3B)3/4×(2arctan(12r1/12(τh*R*2/3B)1/4)2arctan(1+2r1/12(τh*R*2/3B)1/4)+log(τh*R*2/3Br2(τh*R*2/3B)1/4r7/12+r2/3)log(τh*R*2/3Br+2(τh*R*2/3B)1/4r7/12+r2/3)). \begin{equation*}\begin{split} \tau_{\textrm{eff}} &= \tau\ f_{\textrm{ic}} &#x002B; \frac{3\pi \ (\tau\ r)^{1/4} \left(1-f_{\textrm{ic}}\right)}{2\sqrt{2}\ (h_*\ R_*^{-2/3} B)^{3/4}}\\ & \quad&#x002B; \frac{3\ (\tau\ r)^{1/4} \left(1-f_{\textrm{ic}}\right)}{4\sqrt{2}\ (h_*\ R_*^{-2/3} B)^{3/4}} \times \ \Bigg( 2\ \textrm{arctan} \left( 1 - \frac{\sqrt{2}\ r^{1/12}}{\left(\tau\ h_*\ R_*^{-2/3} B \right)^{1/4}} \right) \\ &\quad - 2\ \textrm{arctan} \left( 1 &#x002B; \frac{\sqrt{2}\ r^{1/12}}{\left(\tau\ h_*\ R_*^{-2/3} B\right)^{1/4}} \right) \\ &\quad &#x002B; \textrm{log}\left( \sqrt{\tau\ h_*\ R_*^{-2/3} B\ r}\!-\! \sqrt{2}\! \left(\tau\ h_*\ R_*^{-2/3} B\right)^{1/4}\! r^{7/12}\!&#x002B;\!r^{2/3}\right) \\ & \quad- \textrm{log}\left( \sqrt{\tau\ h_*\ R_*^{-2/3} B\ r}\!&#x002B;\!\sqrt{2}\!\left(\tau\ h_*\ R_*^{-2/3} B\right)^{1/4}\! r^{7/12}\!&#x002B;\!r^{2/3}\right)\Bigg). \end{split} \end{equation*}(C.4)

We note that Eqs. (C.3) and (C.4) indeed reduce to the one-component equations for fic = 0.

Appendix D Visual extinction and decrease in interstellar UV radiation field throughout the outflows considered in Sect. sectionlinking 4

Figures D.1 and D.2 show the visual extinction and decrease in interstellar UV radiation field throughout the one-component and two-component outflows, respectively, considered in Sect. 4.

Appendix E Additional molecules

In this Section, abundance profiles and number densities of additional molecules are shown. The effects of clumpiness and porosity on the chemical formation and destruction pathways of these molecules are similar to those presented in Sect. 4. For both the O- and C-rich outflows, we show CH4, H2 CO, and C2 H2 in analogy with Agúndez et al. (2010). We additionally show CN and the parent species SO for the O-rich outflow (Sect. E.1) and the recently detected species CH3CN and C4 H2 for the C-rich outflow (Sect. E.2).

E.1 O-rich outflows

Figures E.1 and E.2 show the abundance profiles of CH4, H2 CO, C2 H2, SO, and CNfor respectively a one-component and two-component O-rich outflow. The corresponding column densities are listed in Tables E.1 and E.2.

Table E.1

Column density [cm−2] of CH4, H2 CO, C2 H2, SO, and CN in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

CH4 and C2 H2 are expected to have very low abundances assuming TE (Agúndez et al. 2010). Similar to the molecules presented in Sect. 4, we find that a clumpy outflow can increase their inner wind abundances significantly. The largest effects are again seen for one-component, highly porous outflows.

H2CO has been detected in IK Tau with an abundance between 10−7–10−8 (Velilla Prieto et al. 2017). Our models with = 10−6 M yr−1 give a peak abundance similar to the observed range for highly porous one-component outflows.

SO is a parent species in the O-rich outflow, and plays an important role in the formation of CS in the inner wind (Sect. 4.1.3). The decrease in SO abundance is linked to the C abundance (Sect. 4.1.2), and shows that parent species can also be affected by the clumpiness of the outflow.

CN plays a large role in the formation of HCN (Sect. 4.1.2). In IK Tau, it is observed in the intermediate wind (at ~ 1016 cm) with an abundance between 2 × 10−10–6 × 10−8, with a peak abundance of 3 × 10−7 around 5 × 1016 cm (Decin et al. 2010b). In TX Cam, CN has a peak abundance of 2 × 10−7 (Olofsson et al. 1991), with an upper limit of 5 × 10−6 (Charnley & Latter 1997). Both one-component and two-component models with = 10−6 M yr−1 have peak abundances within an order of magnitude of those of IK Tau and TX Cam. We note that this does not only hold for highly porous outflows. The location of the peak CN abundance in the outflow of IK Tau is reproduced by one-component and two-component highly porous models with = 10−6 M yr−1.

E.2 C-rich outflows

Figures E.3 and E.4 show the abundance profiles of OH, H2 CO, HC3N, CH3CN, and C4 H2 for a one-component and two-component C-rich outflow, respectively. The corresponding column densities are listed in Tables E.3 and E.4. The effects of clumpiness and porosity on the chemical formation and destruction pathways of these molecules are similar to those presented in Sect. 4.2.

OH, H2CO, HC3N are expected to have a low inner wind abundance assuming TE (Agúndez et al. 2010). Similar to the molecules presented in Sect. 4, their abundance in the inner region of the outflow increases when taking clumpiness into account, with again the largest effects seen in one-component, highly porous outflows. The abundance of OH is linked to that of H2 O (Sect. 4.2.2).

H2CO has been detected in the outflow of IRC+10216 with an abundance of 1.30.8+1.5×108$1.3_{-0.8}^{&#x002B;1.5} \times 10^{-8}$ (Ford et al. 2004). Highly porous models with = 10−5 M yr−1 have a peak value close to the observed abundance, both for a one- and a two-component model.

HC3N is located in a molecular shell around 3 × 1016 cm in the outflow of IRC+10216 (Agúndez et al. 2017). The location of its peak in abundance clearly depends on the specific clumpiness of the outflow.

thumbnail Fig. D.1

Visual extinction and decrease in interstellar UV radiation field throughout the one-component outflows shown in Figs. 4 and 6. We note that the models with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

The observed abundance profile of CH3CN in the outflow of IRC+10216 has a peculiar shape, which cannot be reproduced by a smooth outflow model or clumpy outflow model (using the formalism of Agúndez et al. 2010). Its abundance peaks around 1016 cm at 5.0 × 10−8 (Agúndez et al. 2015). Our one- and two-component highly porous models have peak abundances within an order of magnitude of the observed value, with the peak also occurring around 1016 cm for a highly porous outflow. The one-component outflows have wider distribution compared to the two-component outflows.

C4 H2 has recently been detected in the outflow of IRC+10216 with a column density of (2.4 ± 1.5) × 1016 cm−2 (Fonfría et al. 2018). From the analysis of Fonfría et al. (2018), it is suggested that the photodissociation occurs closer to the central star, due to, for example, a clumpy outer envelope. Our models with = 10−5 M yr−1 yield a column density close to the observed value for highly porous one-component outflows.

Appendix F Predictability of models

We have selected thus far undetected species whose peak abundance and number density change considerably when assuming a clumpy outflow. These species are selected by requiring that a clumpy outflow causes both a peak fractional abundance larger than 1 × 10−9 and an increase in column density relative to the smooth outflow by at least two orders of magnitude. The abundance profiles and column densities of the selected molecules are considerably changed in a clumpy outflow, making them potentially detectable.

As in Sect. 4 and Appendix E, highly porous one-component outflows show the largest increase in both abundance in the inner region and column density.

F.1 O-rich outflows

Figures F.1 and F.2 show the abundance profiles of N2 O, C2 N, C3 H, C3 H2 and OCS for respectively a one- and two-component O-rich outflow. The corresponding column densities are listed in Tables F.1 and F.2.

F.2 C-rich outflows

Figures F.3 and F.4 show the abundance profiles of CO2, SO2, HC9N, NO and OCS for a one- and two-component O-rich outflow, respectively. The corresponding column densities are listed in Tables F.3 and F.4.

thumbnail Fig. D.2

Visual extinction and decrease in interstellar UV radiation field throughout the two-component outflows shown in Figs. 5 and 7. We note that models with fvol = 0.4 have the same porosity length h* = 2.5 × 1013 cm as the one-component models with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

thumbnail Fig. E.1

Abundance of CH4, H2 CO, C2 H2, SO and CN relative to H2 throughout one-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively.We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

thumbnail Fig. E.2

Abundance of CH4, H2 CO, C2 H2, SO, and CN relative to H2 throughout a two-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

thumbnail Fig. E.3

Abundance of OH, H2CO, HC3 N, CH3 CN, and C4H2 relative to H2 throughout a one-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively.We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

thumbnail Fig. E.4

Abundance of OH, H2CO, HC3 N, CH3 CN, and C4H2 relative to H2 throughout a two-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

thumbnail Fig. F.1

Abundance of N2O, C2 N, C3 H, C3 H2 and OCS relative to H2 throughout a one-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively.We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

thumbnail Fig. F.2

Abundance of N2O, C2 N, C3 H, C3 H2 and OCS relative to H2 throughout a two-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

thumbnail Fig. F.3

Abundance of CO2, SO2, HC9 N, NO, and OCS relative to H2 throughout a one-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively.We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

thumbnail Fig. F.4

Abundance of CO2, SO2, HC9 N, NO, and OCS relative to H2 throughout a two-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

Table E.2

Column density [cm−2] of CH4, H2 CO, C2 H2, SO, and CN in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

Table E.3

Column density [cm−2] of OH, H2CO, HC3 N, CH3 CN, and C4 H2 in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

Table E.4

Column density [cm−2] of OH, H2CO, HC3 N, CH3 CN, and C4H2 in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

Table F.1

Column density [cm−2] of N2 O, C2 N, C3 H, C3 H2 and OCS in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

Table F.2

Column density [cm−2] of N2 O, C2 N, C3 H, C3 H2 and OCS in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

Table F.3

Column density [cm−2] of CO2, SO2, HC9 N, NO, and OCS in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

Table F.4

Column density [cm−2] of CO2, SO2, HC9 N, NO, and OCS in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

References

  1. Agundez, M. 2009, PhD Thesis, Departamento de Astrofisica, Centro de Astrobiologia, CSIC-INTA [Google Scholar]
  2. Agúndez, M., Cernicharo, J., & Guélin, M. 2010, ApJ, 724, L133 [NASA ADS] [CrossRef] [Google Scholar]
  3. Agúndez, M., Fonfría, J. P., Cernicharo, J., et al. 2012, A&A, 543, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Agúndez, M., Cernicharo, J., Quintana-Lacaci, G., et al. 2015, ApJ, 814, 143 [NASA ADS] [CrossRef] [Google Scholar]
  5. Agúndez, M., Cernicharo, J., Quintana-Lacaci, G., et al. 2017, A&A, 601, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Bujarrabal, V., Fuente, A., & Omont, A. 1994, A&A, 285, 247 [NASA ADS] [Google Scholar]
  7. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
  8. Cernicharo, J., Marcelino, N., Agúndez, M., & Guélin, M. 2015, A&A, 575, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Charnley, S. B., & Latter, W. B. 1997, MNRAS, 287, 538 [NASA ADS] [CrossRef] [Google Scholar]
  10. Cherchneff, I. 2006, A&A, 456, 1001 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Cherchneff, I. 2012, A&A, 545, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Cordiner, M. A., & Millar, T. J. 2009, ApJ, 697, 68 [NASA ADS] [CrossRef] [Google Scholar]
  13. De Beck, E., Lombaert, R., Agúndez, M., et al. 2012, A&A, 539, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Decin, L., Agúndez, M., Barlow, M. J., et al. 2010a, Nature, 467, 64 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  15. Decin, L., De Beck, E., Brünken, S., et al. 2010b, A&A, 516, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Dinh-V-Trung & Lim, J., 2008, ApJ, 678, 303 [NASA ADS] [CrossRef] [Google Scholar]
  17. Duari, D., Cherchneff, I., & Willacy, K. 1999, A&A, 341, L47 [NASA ADS] [Google Scholar]
  18. Feldmeier, A., Oskinova, L., & Hamann, W.-R. 2003, A&A, 403, 217 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Fonfría, J. P., Agúndez, M., Cernicharo, J., Richter, M. J., & Lacy, J. H. 2018, ApJ, 852, 80 [NASA ADS] [CrossRef] [Google Scholar]
  20. Ford, K. E. S., Neufeld, D. A., Schilke, P., & Melnick, G. J. 2004, ApJ, 614, 990 [NASA ADS] [CrossRef] [Google Scholar]
  21. Gobrecht, D., Cherchneff, I., Sarangi, A., Plane, J. M. C., & Bromley, S. T. 2016, A&A, 585, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Habing, H. J., & Olofsson, H. 2003, Asymptotic Giant Branch Stars (New York, Berlin: Springer) [Google Scholar]
  23. Justtanont, K., de Jong, T., Helmich, F. P., et al. 1996, A&A, 315, L217 [NASA ADS] [Google Scholar]
  24. Keady, J. J., & Ridgway, S. T. 1993, ApJ, 406, 199 [NASA ADS] [CrossRef] [Google Scholar]
  25. Keller, D., de Koter, A., Decin, L., et al. 2018, A&A, submitted [Google Scholar]
  26. Kervella, P., Montargès, M., Ridgway, S. T., et al. 2014, A&A, 564, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Khouri, T., Maercker, M., Waters, L. B. F. M., et al. 2016, A&A, 591, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Kim, H., Wyrowski, F., Menten, K. M., & Decin, L. 2010, A&A, 516, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Lindqvist, M., Nyman, L.-A., Olofsson, H., & Winnberg, A. 1988, A&A, 205, L15 [NASA ADS] [Google Scholar]
  30. Lombaert, R., Decin, L., Royer, P., et al. 2016, A&A, 588, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Maercker, M., Schöier, F. L., Olofsson, H., Bergman, P., & Ramstedt, S. 2008, A&A, 479, 779 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Maercker, M., Danilovich, T., Olofsson, H., et al. 2016, A&A, 591, A44 [Google Scholar]
  33. Mauron, N., & Huggins, P. J. 2000, A&A, 359, 707 [NASA ADS] [Google Scholar]
  34. Mauron, N., & Huggins, P. J. 2006, A&A, 452, 257 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. McElroy, D., Walsh, C., Markwick, A. J., et al. 2013, A&A, 550, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Menten, K. M., & Alcolea, J. 1995, ApJ, 448, 416 [NASA ADS] [CrossRef] [Google Scholar]
  37. Menten, K. M., Wyrowski, F., Alcolea, J., et al. 2010, A&A, 521, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Millar, T. J. 2016, J. Phys. Conf. Ser., 728, 052001 [CrossRef] [Google Scholar]
  39. Millar, T. J., Herbst, E., & Bettens, R. P. A. 2000, MNRAS, 316, 195 [NASA ADS] [CrossRef] [Google Scholar]
  40. Morris, M., & Jura, M. 1983, ApJ, 264, 546 [NASA ADS] [CrossRef] [Google Scholar]
  41. Nejad, L. A. M., Millar, T. J., & Freeman, A. 1984, A&A, 134, 129 [NASA ADS] [Google Scholar]
  42. Neufeld, D. A., González-Alfonso, E., Melnick, G., et al. 2011, ApJ, 727, L29 [NASA ADS] [CrossRef] [Google Scholar]
  43. Olofsson, H., Lindqvist, M., Winnberg, A., Nyman, L.-A., & Nguyen-Q-Rieu. 1991, A&A, 245, 611 [NASA ADS] [Google Scholar]
  44. Olofsson, H., González Delgado, D., Kerschbaum, F., & Schöier, F. L. 2002, A&A, 391, 1053 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Omont, A., Lucas, R., Morris, M., & Guilloteau, S. 1993, A&A, 267, 490 [NASA ADS] [Google Scholar]
  46. Oskinova, L. M., Feldmeier, A., & Hamann, W.-R. 2004, A&A, 422, 675 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Oskinova, L. M., Feldmeier, A., & Hamann, W.-R. 2006, MNRAS, 372, 313 [NASA ADS] [CrossRef] [Google Scholar]
  48. Owocki, S. P.,& Cohen, D. H. 2006, ApJ, 648, 565 [NASA ADS] [CrossRef] [Google Scholar]
  49. Owocki, S. P., Gayley, K. G., & Shaviv, N. J. 2004, ApJ, 616, 525 [NASA ADS] [CrossRef] [Google Scholar]
  50. Prasad, S. S., & Huntress, Jr. W. T. 1980, ApJS, 43, 1 [NASA ADS] [CrossRef] [Google Scholar]
  51. Richards, A. M. S., Etoka, S., Gray, M. D., et al. 2012, A&A, 546, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Schöier, F. L., Olofsson, H., & Lundgren, A. A. 2006, A&A, 454, 247 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Schöier, F. L., Bast, J., Olofsson, H., & Lindqvist, M. 2007, A&A, 473, 871 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Sundqvist, J. O., Owocki, S. P., Cohen, D. H., Leutenegger, M. A., & Townsend, R. H. D. 2012, MNRAS, 420, 1553 [NASA ADS] [CrossRef] [Google Scholar]
  55. Sundqvist, J. O., Puls, J., & Owocki, S. P. 2014, A&A, 568, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Tenenbaum, E. D., Woolf, N. J., & Ziurys, L. M. 2007, ApJ, 666, L29 [NASA ADS] [CrossRef] [Google Scholar]
  57. Teyssier, D., Hernandez, R., Bujarrabal, V., Yoshida, H., & Phillips, T. G. 2006, A&A, 450, 167 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Tsuji, T., Ohnaka, K., Aoki, W., & Yamamura, I. 1997, A&A, 320, L1 [NASA ADS] [Google Scholar]
  59. Van de Sande, M., Decin, L., Lombaert, R., et al. 2018, A&A, 609, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Velilla Prieto, L., Sánchez Contreras, C., Cernicharo, J., et al. 2017, A&A, 597, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Willacy, K., & Cherchneff, I. 1998, A&A, 330, 676 [NASA ADS] [Google Scholar]
  62. Willacy, K., & Millar, T. J. 1997, A&A, 324, 237 [NASA ADS] [Google Scholar]
  63. Wong, K. T., Menten, K. M., Kamiński, T., et al. 2018, A&A, 612, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Ziurys, L. M., Milam, S. N., Apponi, A. J., & Woolf, N. J. 2007, Nature, 447, 1094 [Google Scholar]

All Tables

Table 1

Physical parameters of the CSE model.

Table 2

Parent species and their initial abundances relative to H2 for the C- and O-rich CSE.

Table 3

Porosity lengths corresponding to the outflows modelled in Sect. 4.

Table 4

Column density [cm−2] of NH3, HCN, and CS in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

Table 5

Column density [cm−2] of NH3, HCN, and CS in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

Table 6

Column density [cm−2] of NH3, H2 O, and H2 S in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

Table 7

Column density [cm−2] of NH3, H2 O, and H2 S in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

Table E.1

Column density [cm−2] of CH4, H2 CO, C2 H2, SO, and CN in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

Table E.2

Column density [cm−2] of CH4, H2 CO, C2 H2, SO, and CN in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

Table E.3

Column density [cm−2] of OH, H2CO, HC3 N, CH3 CN, and C4 H2 in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

Table E.4

Column density [cm−2] of OH, H2CO, HC3 N, CH3 CN, and C4H2 in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

Table F.1

Column density [cm−2] of N2 O, C2 N, C3 H, C3 H2 and OCS in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

Table F.2

Column density [cm−2] of N2 O, C2 N, C3 H, C3 H2 and OCS in a smooth O-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

Table F.3

Column density [cm−2] of CO2, SO2, HC9 N, NO, and OCS in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific one-component outflows.

Table F.4

Column density [cm−2] of CO2, SO2, HC9 N, NO, and OCS in a smooth C-rich outflow with different mass-loss rates, together with column density ratios relative to the smooth outflow for specific two-component outflows.

All Figures

thumbnail Fig. 1

Visualextinction AV throughout a one-component outflow with = 10−5 M yr−1 (left panel) and = 10−7 M yr−1 (right panel). Results for a smooth, uniform outflow are shown in black. Results for different values of h* are shown in blue. For reference, 1 R* = 5 × 1013 cm. Figure B.1 shows the corresponding decrease in interstellar UV radiation field.

In the text
thumbnail Fig. 2

As Fig. 1, but for a two-component outflow with fic = 0.1. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. 3

Visualextinction AV throughout a two-component outflow with h* = 1015 cm and = 10−5 M yr−1 (left panel) and = 10−7 M yr−1 (right panel). Results for a smooth, uniform outflow are shown in black. Results for different values of fic are shown in blue. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. 4

Abundance of NH3 (upper panels), HCN (middle panels), and CS (lower panels) relative to H2 throughout a one-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively. We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. 5

Abundance of NH3 (upper panels), HCN (middle panels), and CS (lower panels) relative to H2 throughout a two-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l*5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. 6

Abundance of NH3 (upper panels), H2O (middle panels), and H2S (lower panels) relative to H2 throughout a one-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively.We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. 7

Abundance of NH3 (upper panels), H2O (middle panels), and H2S (lower panels) relative to H2 throughout a two-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. B.1

Decrease in interstellar UV radiation field throughout a one-component outflow with = 10−5 M yr−1 (left panel) and = 10−7 M yr−1 (right panel). Results for a smooth, uniform outflow are shown in black. Results for different values of h* are shown in blue. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. B.2

As Fig. B.1, but for a two-component outflow with fic = 0.1. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. B.3

Decrease in interstellar UV radiation field throughout a two-component outflow with h = 1015 cm and and = 10−5 M yr−1 (left panel) and = 10−7 M yr−1 (right panel). Results for a smooth, uniform outflow are shown in black. Results for different values of fic are shown in blue. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. D.1

Visual extinction and decrease in interstellar UV radiation field throughout the one-component outflows shown in Figs. 4 and 6. We note that the models with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. D.2

Visual extinction and decrease in interstellar UV radiation field throughout the two-component outflows shown in Figs. 5 and 7. We note that models with fvol = 0.4 have the same porosity length h* = 2.5 × 1013 cm as the one-component models with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. E.1

Abundance of CH4, H2 CO, C2 H2, SO and CN relative to H2 throughout one-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively.We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. E.2

Abundance of CH4, H2 CO, C2 H2, SO, and CN relative to H2 throughout a two-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. E.3

Abundance of OH, H2CO, HC3 N, CH3 CN, and C4H2 relative to H2 throughout a one-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively.We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. E.4

Abundance of OH, H2CO, HC3 N, CH3 CN, and C4H2 relative to H2 throughout a two-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. F.1

Abundance of N2O, C2 N, C3 H, C3 H2 and OCS relative to H2 throughout a one-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively.We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. F.2

Abundance of N2O, C2 N, C3 H, C3 H2 and OCS relative to H2 throughout a two-component O-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. F.3

Abundance of CO2, SO2, HC9 N, NO, and OCS relative to H2 throughout a one-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: characteristic clump scale l* = 5 × 1012 cm, porosity length h* = 1 × 1014, 2.5 × 1013, 1.25 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dashed coloured line: l* = 1013 cm, h* = 2 × 1014, 5 × 1013, 2.5 × 1013 cm for fvol = 0.05, 0.2, 0.4, respectively. Dotted coloured line: l* = 5 × 1013 cm, h* = 1 × 1015, 2.5 × 1014, 1.25 × 1014 cm for fvol = 0.05, 0.2, 0.4, respectively.We note that models with fvol = 0.2,  l* = 5 × 1012 cm (green, solid) and fvol = 0.4,  l* = 1 × 1013 cm (red, dashed) have the same porosity length h* = 2.5 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

In the text
thumbnail Fig. F.4

Abundance of CO2, SO2, HC9 N, NO, and OCS relative to H2 throughout a two-component C-rich outflow with different mass-loss rates and clump volume filling factors fvol. The characteristic size of the clumps at the stellar radius is l* = 1013 cm. Blue lines: porosity length h* = 2 × 1014 cm. Green lines: h* = 5 × 1013 cm. Red lines: h* = 2.5 × 1013 cm. Solid black line: calculated abundance for a smooth, uniform outflow. Solid coloured line: density contrast between the inter-clump and smooth outflow fic = 0.1. Dashed coloured line: fic = 0.3. Dotted coloured line: fic = 0.5. We note that the models with fvol = 0.4 (red) have the same porosity length as the one-component outflows with fvol = 0.2,  l* = 5 × 1012 cm and fvol = 0.4,  l* = 1 × 1013 cm. For reference, 1 R* = 5 × 1013 cm.

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.