Photodissociation and chemistry of N in the circumstellar envelope of carbon-rich AGB stars
Key Words.:Circumstellar matter – Molecular process – Photodissociation – N – Self-shielding – Nitrogen chemistry – IRC +10216
Context: The envelopes of AGB stars are irradiated externally by ultraviolet photons; hence, the chemistry is sensitive to the photodissociation of \ceN2 and CO, which are major reservoirs of nitrogen and carbon, respectively. The photodissociation of \ceN2 has recently been quantified by laboratory and theoretical studies. Improvements have also been made for CO photodissociation.
Aims: For the first time, we use accurate \ceN2 and CO photodissociation rates and shielding functions in a model of the circumstellar envelope of the carbon-rich AGB star, IRC +10216.
Methods: We use a state-of-the-art chemical model of an AGB envelope, the latest CO and \ceN2 photodissociation data, and a new method for implementing molecular shielding functions in full spherical geometry with isotropic incident radiation. We compare computed column densities and radial distributions of molecules with observations.
Results:The transition of \ceN2 N (also, CO C \ceC+) is shifted towards the outer envelope relative to previous models. This leads to different column densities and radial distributions of N-bearing species, especially those species whose formation/destruction processes largely depend on the availability of atomic or molecular nitrogen, for example, CN (=1, 3, 5), CN (=1, 3, 5), HCN (=1, 3, 5, 7, 9), HCN and CHCN.
Conclusions:The chemistry of many species is directly or indirectly affected by the photodissociation of \ceN2 and CO, especially in the outer shell of AGB stars where photodissociation is important. Thus, it is important to include \ceN2 and CO shielding in astrochemical models of AGB envelopes and other irradiated environments. In general, while differences remain between our model of IRC +10216 and the observed molecular column densities, better agreement is found between the calculated and observed radii of peak abundance.
The asymptotic giant branch (AGB) is the last nuclear-burning phase for low- to intermediate-mass stars (from 0.8 M to 8 M), where the core of the star evolves into an inert and degenerate C-O core (Herwig, 2005). AGB stars can be classified by the elemental C/O ratio, namely, C-rich AGB stars (C/O 1), M-type AGB stars (C/O 1) and S-type AGB stars (C/O 1). Our own Sun will become an M-type AGB star. Generally, the envelopes of C-rich stars contain CO plus C-bearing molecules, such as C, CN, HCN, and CH, while O-rich stars contain CO plus O-bearing molecules such as HO, TiO, and VO (Le Bertre, 1997). AGB stars undergo considerable mass loss and eject dust and molecules into the surrounding regions, creating circumstellar envelopes (CSEs). These gas and dust envelopes eventually merge with the interstellar medium (ISM), enriching molecular clouds in which new stars may be born. Furthermore, CSEs are one of the richest sources for detecting new molecules. The study of AGB stars is of particular interest and importance for our understanding of molecular formation, destruction, and the recycling of material between star birth and star death. The schematic structure of the CSE of an AGB star is shown in Fig. 1.
IRC +10216 (CW Leonis), the brightest object in the sky at mid-infrared wavelengths outside the solar system, is the nearest C-rich AGB star and attracts intensive theoretical and observational study (e.g., Morris, 1975; Bieging & Rieu, 1988; Glassgold, 1996; Mauron & Huggins, 1999; Millar et al., 2000; Cernicharo et al., 2000; Woods et al., 2003; Cordiner & Millar, 2009; De Beck et al., 2012; Agúndez et al., 2012; McElroy et al., 2013). More importantly, it is one of the richest molecular sources in the sky. To date, around 180 molecular species, not counting isotopologues, have been identified in the interstellar medium or CSEs111http://www.astro.uni-koeln.de/cdms/molecules/ while more than 80 of them have been detected in IRC +10216. Among these, more than 30 nitrogen-bearing species have been identified, including \ceHCN, CN, SiN, PN, HNC, MgCN, NaCN, NH, HCN, and CN, amongst others (Wakelam et al., 2010). Even water, which was thought unlikely to form in C-rich AGB envelopes, has been detected here (Melnick et al., 2001; Decin et al., 2010). Some molecular detections in IRC +10216 were their first discoveries in astrophysical environments, for example, the first detection of the cyanide anion CN (Agúndez et al., 2010), and FeCN (Zack et al., 2011).
Most of the detected species in IRC +10216 are found in the cool, expanding outer CSE, where photodissociation processes dominate the destruction of molecules. The investigation of ‘key’ reactions, which might affect the entire chemical network and therefore the comparison of simulations and observations, is an important topic. According to a systematic sensitivity study using a large chemical network, the reaction
is one of the most significant reactions in the outer CSE, which directly affects the abundances of many N-bearing species, such as \ceN, \ceN2, \ceHC2N, \ceC3N, \ceC3N- and \ceC2N (Wakelam et al., 2010). Eq. (1) is the primary destruction route of N in any region where UV photons are present. From an observational point of view, the direct detection of N is very challenging because it has no permanent electric dipole and thus possesses no electric-dipole-allowed pure rotational spectrum. The only reported detection of molecular nitrogen is via its far-UV electronic transitions observed in the interstellar medium (Knauth et al., 2004). One can infer \ceN2 indirectly through the protonated ion, \ceN2H+ (Turner, 1974; Herbst et al., 1977) or its deuterated form, \ceN2D+. However, neither \ceN2H+ nor \ceN2D+ have been identified in IRC +10216. From the simulation point of view, Eq. (1) has usually not been treated properly in models because these have not included \ceN2 self-shielding. The importance of self-shielding of molecules in CSEs was first noticed for CO some 30 years ago (Morris & Jura, 1983). Even with an approximate treatment of CO self-shielding (the ‘one-band approximation’), much better agreement was obtained between simulations and observations.
|2. Mass-loss rate||1.5(-5) M yr|
|3. Envelope expansion velocity||14.5 km s|
|4. Radiation field||Standard interstellar radiation field (Draine, 1978), isotropic incidence|
|5. Dust and gas shells in the outflow||Ignored|
|6. Grain surface reactions||Ignored|
|7. Polycyclic aromatic hydrocarbons (PAHs)||Ignored|
|8. Gas density distribution||Falls as , where is the distance from the central star|
|9. Chemical evolution||Kinetic equations solved as a function of radius as material traverses the CSE|
|10. H||Fully self-shielded, no photodissociation|
|11. Parent species||See Table 2|
|12. , , and gas density||See Fig. 2|
|13. Distance||150 pc (De Beck et al., 2012)|
We here employ the latest reported photodissociation rate and shielding functions for \ceN2 (Li et al., 2013; Heays et al., 2014) to investigate the effects in a chemical model of the CSE of IRC +10216. These are based on a concerted laboratory (e.g., Ajello et al., 1989; Helm et al., 1993; Sprengers et al., 2004; Stark et al., 2008; Lewis et al., 2008a; Heays et al., 2011) and theoretical (e.g., Spelsberg & Meyer, 2001; Lewis et al., 2005, 2008b; Ndome et al., 2008) effort over the last two decades. An update was also made to the photodissociation of CO using the self-shielding functions from Visser et al. (2009) following a similarly large experimental effort over the past decades. While the absolute unshielded rates of both \ceN2 and CO are changed only at the level of 30% with respect to the previous values (van Dishoeck, 1988), it is important to realize that the uncertainty in the photorates is reduced by an order of magnitude. Also, we develop and employ a new fully spherically-symmetric (SS) model that computes the self-shielding of molecules in an isotropic interstellar radiation field, rather than the usual plane-parallel (PP) geometry.
The paper is organised as follows: the CSE model, the improvements in \ceN2 and CO photodissociation, as well as the theory and details for the evaluation of photodissociation rates for the SS model are described in Sect. 2. The results and discussion can be found in Sect. 3, followed by the concluding remarks in Sect. 4. The impact of these improvements on the molecular abundances in the CSE of IRC +10216 are discussed. Special attention is given to those species which have already been detected and those which may be detectable in the near future, e.g., using the Atacama Large Millimeter/submillimeter Array (ALMA). A full description of the SS model used for calculating the photodissociation rates of CO and \ceN2, together with the proposed numerical methods for implementing molecular shielding functions, are included in Appendices A and B.
2.1 CSE Model
The CSE model described in McElroy et al. (2013) was employed and extended in the present work. Specific details can be found in Millar et al. (2000) and Cordiner & Millar (2009). Improvements in the model are summarised in Sect. 2.2. The latest (fifth) release of the UMIST Database for Astrochemistry (UDfA), McElroy et al. (2013), hereafter RATE12, was adopted in all of the calculations. RATE12 contains 6173 gas-phase reactions involving 467 species. McElroy et al. (2013) tested this network in a dark cloud model and in a CSE model of IRC +10216 and compared the results with previous models and observations.
The main assumptions and key parameters adopted in our model of IRC +10216 are summarised in Table . Parent species that are injected at the inner radius of the envelope are listed in Table 2. The simulations start at the inner envelope with the following set of parameters: the radius, molecular hydrogen number density, visual extinction, and kinetic temperature of the gas are = 1.0 10 cm, = 1.3 10cm, = 13.8 mag, and = 575 K, respectively. The outer radius of the envelope is set to = 7 10 cm, where the density has decreased to = 26 cm, the temperature to = 10 K, and the visual extinction to = 0.02 mag. We pay particular attention to the chemistry in the outer CSE where photodissociation is an important process. For all photoreactions, species are destroyed by photons from all directions, i.e., the radiation field is isotropic. For IRC +10216, one has to only take into account photons that come from the interstellar medium. The star, with a temperature of 2330 K (De Beck et al., 2012), is too cool to generate photons that could lead to the photodissociation of molecules.
2.2 What’s new
Our aim is to improve the treatment of N and CO photodissociation. In the work of McElroy et al. (2013), they employed the following parameters and methods:
CO: dust + self-shielding, from Morris & Jura (1983), neither taking into account all lines for CO self-shielding nor shielding from H
N: dust shielding only
Method for implementing self-shielding functions
The ‘one-band approximation’ (Morris & Jura, 1983)
In this work, we employ the most accurate molecular data to date and develop a new method for implementing molecular shielding functions in the SS model:
2.3 Photodissociation rate and shielding function
The definitions of the photodissociation rate, , unshielded rate, , and shielding function, , appropriate for a plane-parallel model (only considering photons from the normal direction) are the same as those described in our recent paper (Li et al., 2013). In this work, we consider photons from all directions in space (SS model) for photoreactions that occur in the outer CSE. Suppose an incident ray, inclined with angle to the outward normal direction, arrives at a specific radius, . The total column density of H integrated from this point along the ray to infinity can be written as (Jura & Morris, 1981; Morris & Jura, 1983)
This expression diverges for close to zero and . In practice, the value of equals 1 if . For the case of , its value is obtained using a linear extrapolation from the preceding data points on our angular grid.
In the CSE of IRC +10216, most hydrogen is locked in H due to self-shielding. Thus, the corresponding visual extinction, , is
with the conversion factor based on the observations of Bohlin et al. (1978) and Rachford et al. (2009). The photodissociation rates of most interstellar species are mainly attenuated by dust, which is characterised by dust shielding functions given by
where is the parameter used to estimate the increased dust extinction at ultraviolet wavelengths. In general, the contribution of a specific ray of interstellar photons to the overall photodissociation rate at may be calculated using
However, the photodissociation rates of N, CO, and H are more complex due to self- and mutual-shielding by H, \ceH2, CO, and other molecules (which are wavelength and column density dependent) and also continuum shielding by dust. A convenient way to take these effects into account in models is by use of molecular shielding functions that are calculated from high accuracy wavelength-dependent absorption cross sections. There is no need to consider H photodissociation in the current study because it is fully shielded everywhere in the CSE. For N and CO, one must consider self-shielding and mutual-shielding by \ceH2. The photodissociation rate of N is calculated by
where is the unshielded photodissociation rate of \ceN2, and is the column density of \ceN2 at radius , integrated from this point along angle to infinity. is the molecular shielding function that is responsible for the multiple shielding effects caused by N (self-) and H. It can be obtained via a bi-linear interpolation of N and H shielding functions that depend only on their column densities.
Correspondingly, the CO photodissociation rate is calculated using
where is the CO unshielded photodissociation rate, and is the column density of CO at radius , integrated along angle to infinity. Atomic hydrogen can also shield \ceN2 and CO, although in practice this is unimportant for the case of CSEs because the majority of hydrogen is locked in \ceH2.
A full description of the SS model and the numerical method for the computation of can be found in Appendix A. The implementation of N and CO shielding functions in the SS model is not straightforward. This is because before starting the evolution of the chemistry, one needs to know in advance the abundances of N and CO at each radius to evaluate their column densities, and to compute the molecular shielding functions. However, these data are generated from the output of the model. One approximate solution is to assume that the N and CO number densities are a constant proportion of H, which falls as (Jura & Morris, 1981), then work out the corresponding column densities. This assumption is inaccurate in the photon-dominated regions where N and CO densities are no longer constant with respect to \ceH2. In this study, we propose a new procedure which gives more accurate and convergent results by iterating the computation of the chemistry. This method is described in Appendix B.
The photodissociation of N and CO becomes considerable at the edge of the CSE, where the visual extinction is less than 1.0 mag and the temperature is around 10 K, as can be seen in Fig. 2. When calculating the N shielding functions, the excitation temperature for both H and \ceN2 are chosen to be 10 K. The Doppler widths of N, H, and H are 0.2, 3, and 5 km s, respectively. The column density of atomic hydrogen is taken to be cm. For CO shielding, the excitation temperature for H and CO are chosen to be 11 and 5 K, respectively. At the low densities in the outer CSE, CO is subthermally excited. The Doppler widths of CO, H, and H are 0.3, 3, and 5 km s. Electronic tables of the \ceN2 and CO shielding functions are available online.444http://home.strw.leidenuniv.nl/ewine/photo/ Note that these shielding functions were calculated by considering (\ceH2, H, and self-) shielding at various temperatures for a few typical column densities. Shielding functions for other temperatures and column densities can be obtained by interpolation. In practice, the uncertainties in the shielding functions over a temperature variation of a few degrees are negligible, see Li et al. (2013). The parameter used to account for the dust extinction at ultraviolet wavelengths for \ceN2 and CO was taken to be 3.9 and 3.5, respectively.
Finally, the radial photodissociation rate in a spherical geometry with isotropic incident radiation is evaluated by
3.1 Improvements in N and CO photodissociation
Figure 3 compares the calculated column densities and photodissociation rates for the spherically-symmetric and plane-parallel models. For the PP case, we consider photons from the radial direction only but the full shielding (dust + self- + H + \ceH2) was included. The envelope was illuminated by a Draine (1978) field, extended to wavelengths longer than 2000 according to van Dishoeck & Black (1988) and with a scaling factor, . In the SS model, photons from all directions in space were considered with and the angular weighting scheme is given by Eq. (8). For the PP model, radiation from one side was taken into account explicitly with and radiation from the far side assumed completely attenuated. Hence, in the absence of the envelope the same unattenuated Draine (1978) radiation field is recovered as for the SS model. Indeed, for , the \ceN2 photodissociation rates in the SS and PP cases are very similar (middle panel Fig. 3), whereas when assuming , the PP rate is up to a factor of two higher at the outer edge. As a result, the PP model with underestimates the \ceN2 column density at this location, but the model results are very close to those of the SS model. The bottom panel of Fig. 3 shows that the PP and SS models start to deviate more in the inner envelope, where the PP model gives a much higher photodissociation rate. However, the absolute photodissociation rates are so small here that they are negligible in the chemistry. Overall, we conclude that for the CSE chemistry of IRC +10216, the PP model with a scaling factor of provides excellent agreement with the full SS model.
In the SS model the abundance of \ceN2 is higher than that in the PP model in the outer envelope. The distribution of CO in the SS and PP models is similar to that for N, also indicating a higher abundance of CO at the edge of the cloud for the SS case.
Including or excluding different shielding effects in the model is another crucial factor in the investigation of N and CO photodissociation. Shielding effects, which act like a ‘smoke screen’, help molecules survive in photon dominated regions. N and CO can be significantly shielded by both dust (wavelength independent) and molecules (wavelength dependent, line-by-line shielding). Fig. 4 presents the various shielding effects on N photodissociation. Self-shielding is the most significant shielding effect, whereas (H + H) shielding contributes the least. Shielding from dust largely depends on the properties of the dust (van Dishoeck et al., 2006). It is necessary to include all sources of shielding for \ceN2 because the combined effect has a strong influence on the abundance and distribution of \ceN2.
Figure 5 shows the fractional abundances of N and N, calculated using various models. The location of the transition zone from N to N is shifted outwards by a factor of 5 when molecular shielding is taken into account in the PP () model. Further changes happen when we use the SS model and PP () model.
The photodissociation of CO has an impact on both carbon and nitrogen chemistry. In particular, reactions with C are important destruction mechanisms for quite a few N-bearing species, for example CN (=1, 3, 5). As shown in Fig. 6, the conversion of CO to C and \ceC+ is affected when CO self-shielding is calculated using the updated shielding functions. The chemistry of other C-containing molecules is also affected by the photodissociation of CO, but are beyond the scope of this work. The SS model with full shielding predicts that CO is relatively abundant up to a radius of cm, which is in better agreement with observations than the PP model (which is described further in Sect. 3.3.2).
3.2 Impacts of N and CO photodissociation on the CSE chemistry
The photodissociation of N and CO significantly affects the chemistry in the outer CSE. This influence has two observable aspects: one leads to changes in the radial distribution of the species’ peak abundances. These distributions will become increasingly observable after exploiting the increased spatial resolution and sensitivity of new facilities such as ALMA. The other effect leads to changes in the peak abundances, which are directly reflected in total column densities. In the following subsections, we discuss fractional abundances of the most interesting (detectable or promisingly detectable) species which are sensitive to the photodissociation of N and CO. Special attention is given to N-bearing species. Unless stated elsewhere, comparisons are made between the SS model including the updated \ceN2 and CO shielding and a recreation of the model of McElroy et al. (2013). Additional comparisons between observations and simulations from this study and McElroy et al. (2013) with regards to the 46 detected species in IRC +10216 are discussed in Sect. 3.3.
The dominant formation and destruction mechanisms of most species vary with physical conditions in the CSE. Here, we mainly focus on the chemistry of the outer CSE where photodissociation processes dominate.
3.2.1 Cn, CN CN
When N photodissociation is treated in an appropriate manner, a higher abundance of N is present in the outer CSE with a corresponding decrease in N. The most direct effects are on those N-bearing species whose formation processes require atomic N. \ceC_nN- (=1, 3, 5) have already been observed in the CSE of IRC +10216 (Agúndez et al., 2010; Thaddeus et al., 2008; Cernicharo et al., 2008). The following discussion of \ceC_nN- chemistry is based on the adopted reaction network from Walsh et al. (2009) which was incorporated into McElroy et al. (2013) and our model. Many of the adopted rate coefficients come from Eichelberger et al. (2007).
The chemistry of the molecular anions, \ceCN- and \ceC3N-, in the CSE of IRC +10216 has been described in Kumar et al. (2013). As shown in Fig 7, CN (=1, 3, 5) are very sensitive to \ceN2 photodissociation and are significantly decreased in the new treatment. This is because their formation processes are strongly affected by the abundance of atomic N in the outer CSE. As an example, at 1.0 10 cm \ceCN- and \ceC3N- are mainly formed through:
The abundances of these anions decrease with the decreasing abundance of N. As shown in Fig. 5, the fractional abundance of N drops from (McElroy et al. (2013) model) to (SS model) at 10 cm, leading to 1 order of magnitude drop in \ceCN- and \ceC3N-. Interestingly, CN has a different formation route. It comes primarily from CN, via
3.2.2 Cn, CN CN
The chemistry of CN (=1, 3, 5) are also sensitive to N photodissociation, but not as sensitive as CN (=1, 3, 5). This is because they are only partly affected by the available abundance of atomic N. Take CN as an example, it is mostly formed by the following reactions,
The photodissociation of HCN dominates CN formation inside a radius of 4 cm. Further out, the reactions involving atomic N play a role, although these reactions do not affect its peak abundance, see Fig. 8. The cases of CN and CN are similar to CN.
There is a shift in the peak abundances of CN, \ceC3N- and CN from the outside of the envelope inwards towards the star. In the case of \ceC5N-, this shift is from 10 to cm. This difference is caused by the change in location of the \ceN2 to N transition zone when full shielding is included.
3.2.3 CN CN
As shown in Fig. 9, the peak abundances of CN and CN are located at radii between 5 10 and 1 10 cm. Formation mechanisms for CN and CN are thought to be very similar to those of CN (=1, 3, 5). They have a formation route via:
These larger CN molecules have not yet been identified in IRC +10216, but they are also sensitive to N photodissociation and may be useful in indirect investigations of the chemistry and abundance of \ceN2. Their predicted peak abundances are around one order of magnitude lower than that of \ceC5N.
3.2.4 HcN ( = 3, 5, 7, 9)
The cyanopolyynes, \ceHC_nN ( = 3, 5, 7, 9), are also influenced by N photodissociation, especially in the outer shell where photodissociation processes are important, as shown in Fig. 10. The fractional abundances of HCN ( = 3, 5, 7, 9) are all affected when the full shielding of \ceN2 is considered. The most significantly affected species is HCN. The changes in their total column densities are rather small because the peak abundances, which contribute most to the column densities, are unchanged. However, the radial distribution of \ceHC5N may be spatially resolvable by sensitive interferometers.
The chemistry of large organic molecules is more complex than for simple species. At the photon-dominated radius 10 cm, the main destruction process for the parent species, HCN, is photodissociation, which gives rise to its daughter molecule CN. Another important process for the destruction of HCN is by C, via
However, there are quite a few reformation reactions at this location. Amongst the most important are two reactions relevant to the present study,
The decrease of HCN in the outer CSE, shown in Fig. 10, is largely due to the enhanced shielding of N, through the decreased abundance of N and, indirectly, \ceCN-.
In our calculations, HCN are synthesised by the reactions
The changes in the radial distributions for HCN (=3, 7, 9) are very similar to those of CN, see Fig. 8. HCN is an exception, and has the largest change among all HCN species in the photodissociation region. This is because the main formation process for HCN is via associative electron detachment,
As discussed previously, the abundance and distribution of \ceC5N- is strongly affected by the inclusion of \ceN2 shielding, thus, the abundance of \ceHC5N is also affected.
3.2.5 HCn ChCn
The HCN and CHCN peak abundances are also significantly decreased when full N shielding is included, as shown in Fig.11. Their main formation processes in the outer CSE are
There has already been a detection of CHCN in IRC +10216 (Agúndez et al., 2008), whereas HCN has only been detected in dark clouds (Ohishi et al., 1994). Since the column density computed here is close to the value observed in dark clouds, this species could also be detectable in IRC +10216. Because the chemistry of the two species are relatively simple and largely depend on the availability of N, we may be able to use them to constrain the fractional abundance of N in space by comparing observations with models.
The chemistry of P-bearing and Si-bearing species in C-rich CSEs are not completely understood (McElroy et al., 2013). The most sensitive P-bearing and Si-bearing species to \ceN2 shielding are HNSi, SiN, and PN, see Fig. 12. Here we concentrate only on the chemistry of PN, which is mainly formed through the following route at a radius of cm:
This process is largely reliant on atomic nitrogen, which is significantly decreased at the PN peak radius when updated shielding is used.
3.3 Comparison with observations
|No.||Species c-CH and l-CH represent cyclic- and linear- CH||ObservedThe compiled observational results were taken from Table 7 of McElroy et al. (2013).||CalculatedThe most complete model (SS model with full shielding) in this work.||Calculated Model similar to that of McElroy et al. (2013) with only a few modifications, e.g.,the photodissociation rates of \ceCN- and \ceC3N- have been updated. These have only a small effect on the calculated column densities of species apart from \ceCH3CCH (old vs. new, 4.9 vs. 1.4 cm) and \ceCH2CHCN (3.6 vs. 5.3 cm).||Percentage changeComparison of computational results between the SS model and the model of McElroy et al. (2013), estimated by .||Radius of peakThe currently predicted radii of peak abundances for the species.||Radius of peakObserved radii of peak emission for a few species, estimated from their interferometric maps. Observational references are given under Fig. 13, with priority given to the newer detections.|
|(This work)||(RATE12 model)||between the||abundance||emission|
|(cm)||(cm)||(cm)||two simulations||(cm & arcsecond, )||(cm & arcsecond, )|
|1||CN||1.1(15)||3.5(15)||4.1(15)||-15||5.0(16) & 22||5.4(16) & 24|
|2||CN||2-4(14)||6.7(14)||5.1(14)||31||4.0(16) & 18|
|3||CN||3-6(12)||3.6(13)||3.5(13)||3||4.5(16) & 20||3.0(16) & 13|
|4||CN||5.0(12)||6.0(10)||7.0(11)||-91||7.9(16) & 35|
|5||CN||2.0(12)||3.2(11)||1.0(12)||-68||4.0(16) & 18|
|6||CN||3.0(12)||3.8(12)||7.7(12)||-51||4.5(16) & 20|
|7||HCN||1-2(15)||5.4(14)||4.8(14)||13||2.5(16) & 11||3.8(16) & 17|
|8||HCN||2-3(14)||1.7(14)||1.4(14)||21||3.2(16) & 14|
|9||HCN||1.0(14)||6.0(13)||4.7(13)||28||3.5(16) & 16|
|10||HCN||3.0(13)||1.8(13)||1.3(13)||38||3.5(16) & 16||4.3(16) & 19|
|11||PN||1.0(13)||2.6(09)||5.4(09)||-52||1.1(16) & 5|
|12||SiN||4.0(13)||3.5(11)||2.4(12)||-85||6.3(16) & 28|
|13||SiNC||2.0(12)||8.7(08)||1.2(09)||-28||4.5(16) & 20|
|14||CHNH||9.0(12)||2.8(11)||2.4(11)||17||2.8(16) & 12|
|15||CHCN||8.4(12)||3.0(11)||4.2(12)||-93||5.6(16) & 25|
|16||CHCHCN||5.0(12)||5.1(11)||5.3(11)||-4||5.6(16) & 25|
|17||CHCN||6-30(12)||4.9(12)||5.6(12)||-13||4.0(16) & 18|
|18||C||1.1(16)||1.9(16)||1.9(16)||0||1.4(17) & 62|
|19||C||7.9(14)||5.2(15)||4.2(15)||24||5.0(16) & 22|
|20||C||1.0(15)||2.0(14)||1.8(14)||11||3.5(16) & 16|
|21||C||1.0(14)||2.0(14)||1.5(14)||33||7.1(16) & 32|
|22||CP||1.0(14)||2.4(12)||2.2(12)||9||4.0(16) & 18|
|23||CP||1.0(12)||1.1(09)||8.3(09)||-87||6.3(16) & 28|
|24||SiC||6.0(13)||1.5(13)||1.0(13)||50||7.1(16) & 32|
|25||SiC||2.0(14)||2.6(15)||2.6(15)||0||4.0(16) & 18||3.6(16) & 16|
|26||SiC||4.0(12)||1.7(12)||1.5(12)||13||4.5(16) & 20|
|27||SiC||7.0(12)||9.6(10)||9.5(10)||1||5.6(16) & 25|
|28||HCS||1.0(13)||4.7(11)||3.8(11)||24||5.6(16) & 25|
|29||CS||9-15(13)||1.7(13)||1.4(13)||21||6.3(16) & 28|
|30||CS||6-11(13)||1.6(13)||1.3(13)||23||4.0(16) & 18|
|31||CO||1.0(12)||6.4(11)||6.0(11)||7||4.5(16) & 20|
|32||HCO||5.0(12)||1.9(11)||1.6(11)||19||5.6(16) & 25|
|33||HCO||3-4(12)||1.4(12)||1.3(12)||8||5.0(16) & 22|
|34||CH||3-5(15)||1.1(16)||9.7(15)||13||3.5(16) & 16||3.4(16) & 15|
|35||CH||3-7(13)||1.5(14)||1.4(14)||7||3.5(16) & 16||3.3(16) & 15|
|36||c-CH||2.0(13)||6.1(13)||5.5(13)||11||3.5(16) & 16|
|37||l-CH||3.0(12)||1.3(13)||1.1(13)||18||2.5(16) & 11|
|38||CHCCH||1.6(13)||1.4(12)||1.3(12)||8||2.5(16) & 11|
|39||CH||2-9(15)||7.2(14)||6.5(14)||11||2.8(16) & 12||3.4(16) & 15|
|40||CH||2-50(13)||4.5(13)||4.1(13)||10||3.5(16) & 16|
|41||CH||7.0(13)||7.3(14)||5.7(14)||28||4.5(16) & 20||3.4(16) & 15|
|42||CH||1-2(12)||8.2(13)||7.3(13)||12||5.6(16) & 25|
|43||CH||5.0(12)||1.6(14)||1.2(14)||33||4.0(16) & 18|
|44||CH||7.0(11)||1.3(13)||1.3(13)||0||2.5(16) & 11|
|45||CH||4.0(12)||9.8(13)||8.2(13)||20||1.3(17) & 58|
|46||CH||2.0(12)||3.2(12)||2.9(12)||10||4.0(16) & 18|
3.3.1 Column density
The radially-integrated column densities of most species are controlled by their peak abundances. After investigation, the abundances of quite a few of the 467 species in the network are affected when we use an appropriate method to calculated the photodissociation of N and CO. A summary of the changes in column densities is an efficient way to extract the effects of these updates on the chemistry in the CSE. Here, we concentrate on those 46 species which have been observed in IRC +10216. The calculated and observed column densities are listed in Table . Of these, the calculated values for 14 species vary by more than 30%, and the values for 6 species vary by more than 80%, when compared with results from McElroy et al. (2013). The most affected N-bearing species with respect to the previous model are \ceCN- (-91%), SiN (-85%), and \ceCH2CN (-93%).
All of the species listed in Table contain nitrogen and/or carbon and their formation and destruction may be influenced by \ceN2 and CO photodissociation. The important reactions for some of these species have been discussed in the previous subsections. The most significantly altered column densities include \ceCN- and \ceC3N-, whose formation is sensitive to the available N, which is reduced by the increased shielding of \ceN2 in our new calculations. A factor-of-10 decrease for SiN and \ceCH2CN column densities can be similarly explained.
Some of these changes reduce the agreement of the CSE chemistry model with observed column densities. There are many possible reasons for the discrepancies. First of all, some of the assumptions in the physical model of the CSE, unrelated to photodissociation (see Table ), may be too simple. Further refinement of these assumptions may improve the agreement. For example, in this study, we ignored the presence of gas and dust shells in the expanding envelope. According to the study of Cordiner & Millar (2009), the inclusion of these shells gives a significant improvement in both modelled column densities and spatial distributions. A study of photodetachment as a destruction mechanism for the N-bearing anions, \ceCN- and \ceC3N-, in the CSE of IRC +10216, concluded that the inclusion of shells with enhanced density, similar to the model of Cordiner & Millar (2009), increases the column densities of the anions by about 20% (Kumar et al., 2013).
One of the challenges in CSE simulations is the uncertainty in the abundance of parent species (see Table 2). In particular the modelled column densities of Si-bearing and P-bearing species in Table are several orders of magnitude lower than that suggested by observations. Possibly, the currently modelled Si-bearing and P-bearing parent species should have a higher abundance or there exist other Si- and P-containing parent species. On the other hand, the uncertainties of some observations are considerable. Take \ceSiC2 as an example. The cited column density for this molecule in Table is 2.0 cm (Thaddeus et al., 1984). However, according to a more recent study by Cernicharo et al. (2010) employing high quality Herschel/HIFI spectra an average column density of 6.4 was found. Our computed value, cm, is in better agreement with the latter observed value.
3.3.2 Radius of peak abundance
The radial location of peak molecular abundance is another observable parameter and is of equal importance to the column densities. Many daughter species have a clear radial peak in their modelled abundances which we may directly compare with spatially resolved observations. A selection of such observations have been reduced to their peak radii (without considering detailed excitation or radiative transfer) and compared with our model calculations, and are summarised in Table and Fig. 13.
The model predicts that most peak radii fall within shells between (4) and cm (89) from the central star. Here, the distance between the Earth and IRC +10216 is assumed to be 150 pc, based on the study by De Beck et al. (2012). Different values for this parameter have been employed in other work which may affect intercomparison, e.g., Huggins et al. (1988) who assumed 200 pc, and Cordiner & Millar (2009) who employed 130 pc.
The outer radius of CO using the SS model is slightly larger, 7 10 cm (), than the previous value, see Fig. 6. Emission from CO in IRC +10216 has been extensively investigated. A good summary can be found in De Beck et al. (2012), who also obtained various physical parameters of the structure of IRC +10216. However, conclusions from these observations are not always the same. For instance, Huggins et al. (1988) detected CO emission out to in the line and 180 in the line. Later, the line was detected out to 190 (Fong et al., 2003, 2006). Recently, Cernicharo et al. (private communication) observed weak emission out to a radius of , in excellent agreement with the present simulations.
Figure 13 compares the simulated and observed peak abundance radii for a few selected species for which interferometric maps are available. The observed peak radii, , were deduced from direct measurement of the peak radius and averaged over those maps for which multiple rotational transitions were available. These observations were done with high angular resolution telescopes, such as the Plateau de Bure Interfrometer (PdBI) (Guélin et al., 1998, 1999; Guélin, 2011) and the Very Large Array (VLA) (Dinh-V-Trung & Lim, 2008). In our models, it is assumed that the peak emission occurs at the peak abundance position, i.e., no excitation or radiative transfer of the model lines is performed.
With these assumptions, it is seen in Fig. 13 that the overall agreement between our simulations and observations is within a factor of two. Interestingly, the simulations predict a peak abundance at 18 for the parent species \ceSiC2 (see Table ). This is very close to the value estimated from observations (Lucas et al., 1995; Guélin, 2011), 16. No other parent species exhibits a similar maximum.
None of the peak radii of observationally-mapped species were found to be significantly altered by the introduction of a more appropriate photodissociation treatment, even if their column densities were substantially changed. However, large differences were deduced for some unmapped carbon chain molecules and anions, specifically \ceC3N-, \ceC5N-, \ceC5N, \ceC7N, and \ceC9N. The peak radii of these species all shift inwards after applying the new photodissociation shielding functions. For the case of \ceC5N-, the peak radius retreats by a factor of 2.5, from (44) to (20) cm.
Our new model predicts that the \ceCN- peaks at a somewhat larger radius than its corresponding neutral CN, whereas the peaks of \ceC3N- and \ceC5N- match those of \ceC3N and \ceC5N, respectively, as shown in Table , Fig. 7, and Fig. 8. This conclusion differs from some previous models, e.g., Kumar et al. (2013) and references therein. Similar differences were found between the modelled peak radii of the neutral and anionic forms of \ceC_nH (n=2, 4, 6) (Table ). Interferometric maps of these species could serve as a test of the model predictions.
3.4 N abundance in IRC +10216
There is some uncertainty in the initial abundance of N (a parent species) in the CSE of IRC +10216. In previous models, its adopted abundance relative to \ceH2 lay between 1 and (e.g., Nejad & Millar, 1987; Millar et al., 2000; MacKay & Charnley, 2001; Agúndez & Cernicharo, 2006). In the present work, we employ the same abundance as that in the model of McElroy et al. (2013), 2 10, which gives an overall good agreement between the modelled and observed N-bearing species’ abundances in Table . However, the observational investigation of Milam et al. (2008) deduced an N abundance of . They used the observation of PN to infer an \ceN2 abundance by assuming approximate equality of the abundance ratios, PN/\ceN2 [P]/[N] and [P]/[N] HCP/HCN 0.001.
To test the sensitivity of daughter species to the initial \ceN2 abundance we recomputed their chemical abundances after lowering the parent \ceN2 abundance to , a factor of 20 lower than in our fiducial models. The results for four species, \ceCN-, \ceC5N-, \ceC5N-, and \ceN2H+ are shown Fig. 14.
N2H+ has not been observed in IRC +10216, but has been detected in other interstellar and circumstellar environments (e.g., Liu et al., 2013; Qi et al., 2013; Codella et al., 2013). This species is thought to be a good tracer for N, from which it forms directly via the reaction,
The utility of \ceN2H+ as an \ceN2 tracer is supported by the calculated abundances plotted in Fig. 14, which are directly proportional to the assumed \ceN2 abundance. The predicted column density of \ceN2H+ assuming initial \ceN2 abundances of and are 9.4 10 and 5.3 10 cm, respectively. The peak abundance of \ceN2H+ is located at a radius of cm (22). Future deeper observations of IRC+10216 with ALMA may determine the column density of \ceN2H+ and allow for a new estimate of the \ceN2 abundance. Note that the abundances of neither \ceN2H+ nor \ceHCO+, formed through a similar reaction between CO and \ceH3+, are significantly affected by the new model treatment.
The abundances of \ceC_nN- anions are also not significantly affected by the altered parent \ceN2 abundance. These molecules are strongly dependent on the availability of atomic nitrogen (Sec. 3.2.1) which in the region of their peak abundance forms primarily from the photodissociation of HCN. Thus, HCN is the important N-bearing parent species affecting their column densities. Note that the main source of atomic N further out in the envelope is via the photodissociation of \ceN2.
4 Concluding remarks
In this work, the effect of self-shielding of \ceN2 and CO, and mutual-shielding by \ceH2, in the outer envelope of the C-rich AGB star, IRC +10216, was studied for the first time. This was performed using the latest available data for \ceN2 and CO photodissociation, and using an extended spherically-symmetric model for molecular self-shielding. The impact of these two improvements on the chemical evolution of the expanding envelope was investigated, with special attention paid to nitrogen chemistry and detectable species.
Key points from this study are:
N and CO are more abundant at the edge of the CSE than predicted by previous models due to molecular shielding from photodissociation.
The photodissociation of \ceN2 and CO affects the chemistry of most species in the outer CSE. Following its improved treatment here, the transition zones of N N and CO \ceC+ C shift towards the outer edge of the envelope. The abundances of some species formed from N are reduced in the outer CSE. This induces large changes in the predicted column densities (factor of 10) of some species (e.g., \ceC_nN and \ceC_nN- carbon chains) and the radii of their peak abundances.
Predictions are made for column densities and peak radii of molecules whose abundances are sensitive to \ceN2 and CO photodissociation and which could conceivably be detected. These may be verified by observations using high resolution and high sensitivity telescopes, especially ALMA and PdBI.
The abundances obtained from models which treat the full angularly-resolved spherically-symmetric radiation field incident on the outer CSE are very similar to those found in plane-parallel models in which the radiation is incident normally but with the intensity reduced by half. The proposed iterative method for implementing molecular shielding functions is most efficient when assuming that the molecule is fully shielded at the beginning of the calculation.
Other astrochemical models simulating the CSEs of AGB stars or other environments may show a similar sensitivity using a more realistic treatment of self-shielding and spherically-symmetric radiation as was found in this model of IRC +10216.
The photodissociation processes for \ceN2 and CO (both the unshielded rate as well as the depth dependence) are now very well understood, so any major discrepancies between observations and models are likely due to other assumptions in the model.
Acknowledgements.The authors would like to thank the anonymous referee for his/her valuable comments and constructive suggestions. X. Li is delighted to thank TJM for his hospitality during a two-week visit to Queen’s University Belfast (QUB) in 2013, and thank Dr. Markus Schmalzl and Prof. Xander Tielens for some stimulating discussions. Astrochemistry in Leiden is supported by the Netherlands Research School for Astronomy (NOVA), by a Spinoza grant and grant 648.000.002 from the Netherlands Organisation for Scientific Research (NWO), and by the European Community’s Seventh Framework Programme FP7/2007-2013 under grant agreements 291141 (CHEMPLAN) and 238258 (LASSIE). Astrophysics at QUB is supported by a grant from the STFC. C. W. acknowledges support from the NWO (program number 639.041.335).
- Agúndez & Cernicharo (2006) Agúndez, M. & Cernicharo, J. 2006, ApJ, 650, 374
- Agúndez et al. (2010) Agúndez, M., Cernicharo, J., Guélin, M., et al. 2010, A&A, 517, L2
- Agúndez et al. (2012) Agúndez, M., Fonfría, J. P., Cernicharo, J., et al. 2012, A&A, 543, A48
- Agúndez et al. (2008) Agúndez, M., Fonfría, J. P., Cernicharo, J., Pardo, J. R., & Guélin, M. 2008, A&A, 479, 493
- Ajello et al. (1989) Ajello, J. M., James, G. K., Franklin, B. O., & Shemansky, D. E. 1989, Phys. Rev. A, 40, 3524
- Bieging & Rieu (1988) Bieging, J. H. & Rieu, N.-Q. 1988, ApJ, 329, L107
- Bohlin et al. (1978) Bohlin, R. C., Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132
- Cernicharo et al. (2008) Cernicharo, J., Guélin, M., Agúndez, M., McCarthy, M. C., & Thaddeus, P. 2008, ApJ, 688, L83
- Cernicharo et al. (2000) Cernicharo, J., Guélin, M., & Kahane, C. 2000, A&AS, 142, 181
- Cernicharo et al. (2010) Cernicharo, J., Waters, L. B. F. M., Decin, L., et al. 2010, A&A, 521, L8
- Codella et al. (2013) Codella, C., Viti, S., Ceccarelli, C., et al. 2013, ApJ, 776, 52
- Cordiner & Millar (2009) Cordiner, M. A. & Millar, T. J. 2009, ApJ, 697, 68
- De Beck et al. (2012) De Beck, E., Lombaert, R., Agúndez, M., et al. 2012, A&A, 539, A108
- Decin et al. (2010) Decin, L., Agúndez, M., Barlow, M. J., et al. 2010, Nature, 467, 64
- Dinh-V-Trung & Lim (2008) Dinh-V-Trung & Lim, J. 2008, ApJ, 678, 303
- Draine (1978) Draine, B. T. 1978, ApJS, 36, 595
- Eichelberger et al. (2007) Eichelberger, B., Snow, T. P., Barckholtz, C., & Bierbaum, V. M. 2007, ApJ, 667, 1283
- Fong et al. (2003) Fong, D., Meixner, M., & Shah, R. Y. 2003, ApJ, 582, L39
- Fong et al. (2006) Fong, D., Meixner, M., Sutton, E. C., Zalucha, A., & Welch, W. J. 2006, ApJ, 652, 1626
- Glassgold (1996) Glassgold, A. E. 1996, ARA&A, 34, 241
- Guélin (2011) Guélin, M. 2011, in EPJ Web of Conferences, 18, 1002
- Guélin et al. (1998) Guélin, M., Neininger, N., & Cernicharo, J. 1998, A&A, 335, L1
- Guélin et al. (1999) Guélin, M., Neininger, N., Lucas, R., & Cernicharo, J. 1999, in The Physics and Chemistry of the Interstellar Medium, ed. V. Ossenkopf, J. Stutzki, & G. Winnewisser, 326
- Heays et al. (2011) Heays, A. N., Dickenson, G. D., Salumbides, E. J., et al. 2011, J. Chem. Phys., 135
- Heays et al. (2014) Heays, A. N., Visser, R., Gredel, R., et al. 2014, A&A, 562, A61
- Helm et al. (1993) Helm, H., Hazell, I., & Bjerre, N. 1993, Phys. Rev. A, 48, 2762
- Herbst et al. (1977) Herbst, E., Green, S., Thaddeus, P., & Klemperer, W. 1977, ApJ, 215, 503
- Herbst & Osamura (2008) Herbst, E. & Osamura, Y. 2008, ApJ, 679, 1670
- Herwig (2005) Herwig, F. 2005, ARA&A, 43, 435
- Huggins et al. (1988) Huggins, P. J., Olofsson, H., & Johansson, L. E. B. 1988, ApJ, 332, 1009
- Jura & Morris (1981) Jura, M. & Morris, M. 1981, ApJ, 251, 181
- Knauth et al. (2004) Knauth, D. C., Andersson, B.-G., McCandliss, S. R., & Warren Moos, H. 2004, Nature, 429, 636
- Kumar et al. (2013) Kumar, S. S., Hauser, D., Jindra, R., et al. 2013, ApJ, 776, 25
- Le Bertre (1997) Le Bertre, T. 1997, Lecture Notes in Physics, Vol. 497, Stellar Atmospheres: Theory and Observations, ed. J. P. Greve, R. Blomme, & H. Hensberge (Berlin, Heidelberg: Springer Berlin Heidelberg), 133
- Lewis et al. (2008a) Lewis, B. R., Baldwin, K. G. H., Sprengers, J. P., et al. 2008a, J. Chem. Phys., 129, 164305
- Lewis et al. (2005) Lewis, B. R., Gibson, S. T., Zhang, W., Lefebvre-Brion, H., & Robbe, J.-M. 2005, J. Chem. Phys., 122, 144302
- Lewis et al. (2008b) Lewis, B. R., Heays, A. N., Gibson, S. T., Lefebvre-Brion, H., & Lefebvre, R. 2008b, J. Chem. Phys., 129, 164306
- Li et al. (2013) Li, X., Heays, A. N., Visser, R., et al. 2013, A&A, 555, A14
- Liu et al. (2013) Liu, X.-L., Wang, J.-J., & Xu, J.-L. 2013, MNRAS, 431, 27
- Lucas et al. (1995) Lucas, R., Guélin, M., Kahane, C., Audinos, P., & Cernicharo, J. 1995, Ap&SS, 224, 293
- MacKay & Charnley (2001) MacKay, D. D. S. & Charnley, S. B. 2001, MNRAS, 325, 545
- Mauron & Huggins (1999) Mauron, N. & Huggins, P. J. 1999, A&A, 349, 203
- McElroy et al. (2013) McElroy, D., Walsh, C., Markwick, A. J., et al. 2013, A&A, 550, A36
- Melnick et al. (2001) Melnick, G. J., Neufeld, D. A., Ford, K. E. S., Hollenbach, D. J., & Ashby, M. L. N. 2001, Nature, 412, 160
- Milam et al. (2008) Milam, S. N., Halfen, D. T., Tenenbaum, E. D., et al. 2008, ApJ, 684, 618
- Millar et al. (2000) Millar, T. J., Herbst, E., & Bettens, R. P. A. 2000, MNRAS, 316, 195
- Morris (1975) Morris, M. 1975, ApJ, 197, 603
- Morris & Jura (1983) Morris, M. & Jura, M. 1983, ApJ, 264, 546
- Ndome et al. (2008) Ndome, H., Hochlaf, M., Lewis, B. R., et al. 2008, J. Chem. Phys., 129, 164307
- Nejad & Millar (1987) Nejad, L. A. M. & Millar, T. J. 1987, A&A, 183, 279
- Ohishi et al. (1994) Ohishi, M., McGonagle, D., Irvine, W. M., Yamamoto, S., & Saito, S. 1994, ApJ, 427, L51
- Qi et al. (2013) Qi, C., Öberg, K. I., & Wilner, D. J. 2013, ApJ, 765, 34
- Rachford et al. (2009) Rachford, B. L., Snow, T. P., Destree, J. D., et al. 2009, ApJS, 180, 125
- Spelsberg & Meyer (2001) Spelsberg, D. & Meyer, W. 2001, J. Chem. Phys., 115, 6438
- Sprengers et al. (2004) Sprengers, J. P., Ubachs, W., Johansson, A., et al. 2004, J. Chem. Phys., 120, 8973
- Stark et al. (2008) Stark, G., Lewis, B. R., Heays, A. N., et al. 2008, J. Chem. Phys., 128, 114302
- Thaddeus et al. (1984) Thaddeus, P., Cummins, S. E., & Linke, R. A. 1984, ApJ, 283, L45
- Thaddeus et al. (2008) Thaddeus, P., Gottlieb, C. A., Gupta, H., et al. 2008, ApJ, 677, 1132
- Turner (1974) Turner, B. E. 1974, Ap&SS, 29, 247
- van Dishoeck (1988) van Dishoeck, E. F. 1988, in Rate Coefficients in Astrochemistry, ed. T. J. Millar, and D. A. Williams, Kluwer Academic Publishers, Dordrecht, 49
- van Dishoeck & Black (1988) van Dishoeck, E. F. & Black, J. H. 1988, ApJ, 334, 771
- van Dishoeck et al. (2006) van Dishoeck, E. F., Jonkheid, B., & van Hemert, M. C. 2006, Faraday Discuss., 303, 231
- Visser et al. (2009) Visser, R., van Dishoeck, E. F., & Black, J. H. 2009, A&A, 503, 323
- Wakelam et al. (2010) Wakelam, V., Smith, I. W. M., Herbst, E., et al. 2010, Space Sci. Rev., 156, 13
- Walsh et al. (2009) Walsh, C., Harada, N., Herbst, E., & Millar, T. J. 2009, ApJ, 700, 752
- Woods et al. (2003) Woods, P. M., Schöier, F. L., Nyman, L.-Å., & Olofsson, H. 2003, A&A, 402, 617
- Zack et al. (2011) Zack, L. N., Halfen, D. T., & Ziurys, L. M. 2011, ApJ, 733, L36
Appendix A Spherically-symmetric model
In this appendix we describe the SS model employed in this work for the computation of the angle-dependent column densities of species at each radius. The SS model is shown in Fig. 15. The star is located at point S and ejects dust and molecules to various shells with radii (not equidistant). Since the star itself is too cool to generate photons that can induce photodissociation of molecules (in our case, IRC +10216, K), one only needs to consider photons from the interstellar radiation field. For a specific position P at radius , suppose is the angle between the direction of the incoming ray and the radius vector from the central star. Generally, the column density, , of a species is defined as
where is the number density in units of cm and is the integral path length in units of cm. In the present study, the total column density, , of a species integrated along a ray from point P to infinity is given by two cases:
, for example, , yielding
where is the number density of the molecule at radius .
, for instance . Since we assume the model is spherically symmetric, one may simplify the calculations by finding a point C that satisfies , then get the column density along with
Both and can be conveniently calculated by case (1). Then, it is necessary to interpolate an abundance between the radially-gridded values in our model for points near the star, such as C.
Figure 16 shows for molecular N. According to our calculations, the major contribution of photons that induce photodissociation of N come from the area within the dashed line ( cm), where cm and . Deeper in the envelope, the fractional abundance of N stays constant because it is fully shielded.
Appendix B Implementing the molecular shielding functions
In this Appendix we introduce the new procedure employed in the current work for the implementation of molecular shielding functions. Here we take N as an example for the demonstration. The method has five steps.
Set up the initial condition. Suppose N is fully shielded everywhere in the CSE, then we can set for each radius, , which means the corresponding shielding functions from dust and molecules are zero, i.e., , and . In all cases, we consider 201 incident rays as a function of between and , where gives the normal direction.
Loop 1: compute the CSE chemistry for the first time, employing prepared from Step 1. The calculated N abundance generates column densities, , and corresponding dust and molecular shielding functions, , and , and new photodissociation rates .
Loop 2: repeat the simulations of the CSE chemistry, but employing produced in Step 2 for each radius. Using the same method, generate another set of photodissociation rates, .
Compare the fractional abundances of N at each radius for Loop 1 and Loop 2, to see if they are the same. If no changes are found, go to Step 5, otherwise go to Step 3, continuing the calculations.
Stop the calculation. The outputs at this step are the final results.
The simulations can also approach convergence by assuming N is unshielded in Step 1, namely at each radius. In this case, the fractional abundance of N will steeply decrease after starting the evolution of the chemical network, then approach the converged results quickly, as illustrated in Fig. 17. The same converged results can be obtained from the two possible assumptions in Step 1. In practice, the fractional abundances of both N and CO (as well as all other species) reach their converged abundances after 15 loops. In practice, starting the calculation assuming full shielding in Loop 1 is most efficient.
The iterative method proposed here may be used in other models, where output variables are also required as inputs.