Volume 521, October 2010
|Number of page(s)||11|
|Section||Interstellar and circumstellar matter|
|Published online||18 October 2010|
Subsurface chemistry of mantles of interstellar dust grains in dark molecular cores
J. Kalvans1 - I. Shmeld2
1 - Institute of Astronomy, University of Latvia, Raina 19, Riga, Latvia
2 - Ventspils International Radioastronomy Centre of Ventspils University College, In enieru iela 101, Ventspils, Latvia
Received 3 February 2010 / Accepted 13 May 2010
Context. The abundances of many observed compounds in interstellar molecular clouds still lack an explanation, despite extensive research that includes both gas and solid (dust-grain surface) phase reactions.
Aims. We aim to qualitatively prove the idea that a hydrogen-poor subsurface chemistry on interstellar grains is responsible for at least some of these chemical ``anomalies''. This chemistry develops in the icy mantles when photodissociation reactions in the mantle release free hydrogen, which escapes the mantle via diffusion. This results in serious alterations of the chemical composition of the mantle because pores in the mantle provide surfaces for reactions in the new, hydrogen-poor environment.
Methods. We present a simple kinetic model, using existing astrochemical reaction databases. Gas phase, surface and subsurface pore reactions are included, as are physical transformations of molecules.
Results. Our model produces significantly higher abundances for various oxidized species than most other models. We also obtain quite good results for some individual species that have adequate reaction network. Thus, we consider that the hydrogen-poor mantle chemistry may indeed play a role in the chemical evolution of molecular clouds.
Conclusions. The significance of outward hydrogen diffusion has to be proved by further research. A huge number of solid phase reactions between many oxidized species is essential to obtain good, quantitative modeling results for a comparison with observations. We speculate that a variety of unobservable hydrogen-poor sulfur oxoacid derivatives may be responsible for the ``disappearance'' of sulfur in dark cloud cores.
Key words: ISM: abundances - ISM: clouds - astrochemistry - molecular processes
It is widely accepted that molecular hydrogen and many other interstellar molecules form on interstellar dust grains. There has been wide research in the field of gas-grain chemistry occurring in the dark, dense cores of interstellar molecular clouds, including observations and calculations. Various desorption mechanisms and chemical reactions on the surfaces of the grains are investigated.
Several molecules in the interstellar medium at least are known whose abundances are not easily explained by gas-phase and grain surface-phase chemistry. These include OCS, HCN, SO, cyanopolyynes, and others. In the warming star-formation regions there appear highly oxidized organic compounds that indicate that a different chemistry occurs in these regions or, as we believe, the ejection of heavily processed interstellar grain mantles into the gas phase.
The grain surface reactions in interstellar clouds by definition are subjected to heavy hydrogenation. The proper production of heavier and hydrogen-poor species is difficult to reproduce by calculations at least in some cases (Hasegawa & Herbst 1993a; van Weeren et al. 2009; Hatchell et al. 1998). We present an alternative explanation of this problem by considering the possibility that within the grain mantles below the surface layer there is a hydrogen-poor, chemically active material. The aim of this article is to qualitatively answer this question through existing knowledge and means of astrochemical problem solving.
There are only few such models that take into account more than one layer of the accreted species. There are some serious researches done which insist that the chemistry of the species frozen onto grains does not end with the formation of the next accreted layer. Shalabiea & Greenberg (1994) examine photon-induced processes inside the mantle. They conclude that photoprocessing of grain mantles is the start of the synthesis of many species. Hasegawa & Herbst (1993b) use a 3-phase model to examine the formation and composition of the inner mantle. They found that radicals like OH, etc. are absorbed into mantles in large numbers. We argue that the chemistry below the outer surface of the icy mantle has to be common and is an important path in interstellar molecule synthesis. Besides, they also note that grain surface reactions tend to overproduce hydrogenated species, which is the main problem we attempt to tackle in this work. Schutte & Greenberg (1991) examine the possibility of molecule desorption from the grains by chemical explosions within the grain mantle. Naturally, these reactions can be expected to alter the chemical composition of the mantle itself. Freund & Freund (2006) present a dust grain model based on the principle of solid solutions, producing results that explain important features in the molecular cloud composition.
We present a model of the processes on the surface and inside the ``frozen'' grain mantles with a basic concept that chemical reactions occur on the surface of pores (or cracks) inside the mantles. The point that makes the difference between the surface and mantle reactions is that the mantle is not directly exposed to the ocean of hydrogen in a nebula, and thus a significantly different chemistry may develop. This gives an opportunity to present an explanation for some of the astrochemical mysteries. We do not use advanced calculation techniques or new reactions; we evaluate the importance of hydrogen diffusion through the grain mantle on the chemical composition of the mantles in dark molecular cloud cores. Thus, the model includes various physiochemical processes but otherwise keeps a rather conservative approach.
2 The model
2.1 Model considerations
Current models, e.g. Das et al. (2008), Goldsmith et al. (2008), van Weeren et al. (2009) insist that the collapse of an interstellar cloud to densities around 105 cm-3 occurs at a time around 106 years and that molecules form simultaneously. Deposition onto grains proceeds at timescales comparable to cloud cooling and collapse, and most of the accreted matter will accumulate at late stages of the cloud evolution. Taking this into consideration we developed a steady-phase chemical model for timescales much longer than the cloud formation.
In order to investigate the molecular abundances in dark cloud cores under chemical equilibrium we used a chemical kinetics model with 352 molecules. The interstellar UV radiation (molecule photodissociation and photodesorption) has been neglected, the integration time was taken to be 1016 s, comparable to an interstellar cloud entire lifetime, and the cloud density cm-3. The gas temperature is taken to be 15 K, the dust grain temperature - 10 K. We used the UDFA06 dipole (udfa06) astrochemistry database (Woodall et al. 2006) is used to provide the gas-phase chemical reaction set. We adopted elemental abundances provided by Jenkins (2009). The elements permitted are H, He, C, N, O, Na, Mg, Si, S, Fe.
The grain chemistry is described in terms of surface reactions between species that are located on the outer grain surface or on a surface of pores inside the grain mantle itself. It might be possible to even more adequately describe the grain mantle chemistry with the concept of a solid solution or nanoporous matter. For the grain processes the surface reaction set by Hasegawa et al. 1992 and Hasegawa & Herbst (1993a) is used (neutral molecules only). Species without gas phase udfa06 reactions were excluded from calculations, because they also lack the important photodissociation sink reactions (Sect. 2.6). These include , HOC, NaOH, etc. In order to provide a more adequate reaction set for each molecule, organic molecules are included up to only in the grain surface and mantle model. This is because larger molecules tend to have less surface reactions included, especially when one takes into account the number of atoms they contain. Notably, there is an almost complete lack of oxidation reactions for the more complex carbon species. These reactions are extremely important because our model shows that a highly oxidative environment is possible in the grain mantles.
Because of the use of a huge gas-phase reaction database and a limited solid-phase reaction database, there are many gas phase molecules that do not take part in the surface and mantle processes. Since our aim is to investigate exclusively the solid phase of interstellar molecules, we decided for the sake of completeness to leave the udfa06 as intact as possible in our model, with some 260 gas-only species. Most of the molecules with a gas phase only are not relevant to the solid phase (ions), and the remainder are expected to have negligible abundances anyway (complicated organic species). When one compares similar molecules with and without a solid phase, the difference in abundances is about one order of magnitude. The full list of results for gas-phase species is given in Table 4.
We used a chemical model integrated in a 3-phase system similar that of Hasegawa & Herbst (1993b). The three phases are gas, grain surface, and grain mantle. The grain surface consists of reactive species that accrete from gas, may be easily desorbed by several mechanisms, and are subjected to photodissociation. In this context the surface consists of the first few layers of a mantle. The mantle itself is formed by buried surface molecules. H atoms and especially molecules, which are created in the mantle by photolysis can migrate away from their parent molecule, and the mantle becomes enriched with free hydrogen. This excess hydrogen is can escape to the outer surface by diffusion. Thus, with photolysis a mostly hydrogen-poor chemical environment forms below the surface.
Molecules adsorbed onto grains are divided into two layers - the surface and the mantle. We employed a model of an equilibrium state, where surface and gas phase molecules are in a dynamic equilibrium, while mantle molecules are almost permanently locked away in a frozen state. Thus, ``surface'' represents a few (rugged) top layers that are subjected to desorption and surface reactions and may be brought to the very top layer by these processes. ``Mantle'', by definition, is never exposed to the surface and is only slowly returned to the gas phase by direct ejection caused by cosmic ray hits. In order to properly describe the physical and chemical processes we avoided the division of the mantle in layers. Surface and mantle molecules are treated as solid species uniformly dispersed in volume with abundances expressed in cm-3.
2.2 Molecule accretion onto grains
Molecule accretion onto grains happens in an equilibrium with the
various desorption processes. In our model only neutral molecules
accrete. The rate coefficient (s-1) is calculated according to Willacy & Williams (1993) and Nejad & Wagenblast (1999) by the formula
where Si is the sticking coefficient, Mi is molecular mass of species i in atomic mass units, and is total number of H atoms per cm3, is gas temperature (K).
The rate of accretion, cm-3 s-1, is thus
where ni is the gas phase abundence of species i.
The sticking coefficient used by previous authors (i.e. Willacy & Williams 1993; Nejad & Wagenblast 1999; Roberts et al. 2007; Aikawa et al. 1997; Turner 1998b; Willacy & Millar 1998; Brown & Charnley 1990) for heavy species is 0.1 to 1, usually around 0.33, and that for light species is 0 to 1. We took an approximate average path with Si = 0.33 for all heavy species and 0.1 for hydrogen atoms and molecules.
2.3 Thermal evaporation
The thermal evaporation rate is calculated from the evaporation times given in Hasegawa & Herbst (1993a).
2.4 Direct cosmic-ray heating desorption
Desorption by heating, caused by Fe nuclei of cosmic rays, is calculated with the rate coefficients given by Hasegawa & Herbst (1993a).
2.5 Cosmic ray induced photodesorption
Rate coefficient (s-1) for desorption by cosmic ray
induced photons for surface species is calculated by the formula
adapted from Willacy & Williams (1993)
where the photon hit rate for a single grain
Yield Y is taken to be 0.1, is photon flux, taken 4875 cm-2 s-1 from Roberts et al. (2007). is the average cross section of a grain ( cm2). Like for other desorption mechanisms, the desorption rate is
is the average abundance (cm-3) of species i residing on the outer surface of a grain.
2.6 Dissociation by cosmic ray induced photons in solid stateBesides the surface binary reactions of Hasegawa et al. (1992) and Hasegawa & Herbst (1993a), we also include a limited set of cosmic ray photon induced photoreactions on grain surfaces.
We pretend that the individual molecules on and in the amorphous mantles with mixed composition essentially keep their UV absorption properties. The gas phase reaction coefficients we obtain from udfa06. H and formed in photodissociation on the outer surface of grains are always allowed to escape into gas phase, while all other species remain on the grain surfaces (their desorption by CR induced photons is described above in Sect. 2.5). We consider it a reasonable approximation given that the surface species in the model of Hasegawa & Herbst (1993b) are probably too intensively hydrogenated, and Andersson & van Dischoeck (2008) concludes that photodissociation of amorphous water ice mostly results in an escaping hydrogen atom.
Several sources (Shalabiea & Greenberg 1994; Andersson & van Dischoeck 2008; Öberg et al. 2009) indicate that photoreactions can and do occur in the layers below the surface, although they usually do not lead to desorption.
Thus, the formula for photodissociation rate coefficient is
Generally, the quantum yield is taken to be 0.1 for any given species on grain surface . We assume that within the mantle only molecules with access to a surface can be effectively dissociated, so that the products do not recombine again. We expect 1/100 of the mantle to be exposed to inner pore surfaces. The dissociation yield for mantle molecules is
The dissociation products of these inert molecules are then assumed to be reactive until they become frozen again (see Sect. 2.8).
2.7 Desorption resulting from H2 formation
We calculated the rate coefficient (s-1) for desorption of surface molecules by heat released by
molecule formation on grains according to Roberts et al. (2007)
where the abundance (cm-3) of dust grains is
is the formation rate of molecules and is number of molecules desorbed per act of molecule formation. The species desorbed are those with binding energies below a treshold energy. According to the discussion by Willacy & Millar (1998) and Roberts et al. (2007), the efficience of desorption by heat released by formation is a subject of discussion and suggestions. Following Roberts et al. (2007) we take . Like in other models, this selective desorption mechanism has a strong effect on the abundances of major species in both, gas and solid phases. We choose to diminish its importance by taking K (binding energies from Hasegawa & Herbst 1993a).
2.8 Binary reactions on grains
For chemical processes on grain surfaces, i.e. reactions between
accreted species, we use the reaction set provided by Hasegawa
et al. (1992) and Hasegawa & Herbst (1993a). The rate coefficient is
where the diffusion rate is
Because the roughness and inhomogeneity of the mantle are key features of this model, we used the higher quantity for number surface adsorption sites instead of 106, which characterizes a perfect spherical grain. is the probability for the reaction to happen upon an encounter, is the hopping time (s) for mobile species.
We assume that physical transformations of a fully formed grain mantle
in steady conditions inside dark cores are driven by cosmic rays
passing through the grain. When a heavy cosmic ray particle hits the
grain, the mantle is heated, shaken, and at least locally rearranged.
New surfaces appear and molecules and some radicals emerge, migrate,
and react, while other species become locked in ice again. Thus the
mantle (and, perhaps, the grain) is not a frozen chunk of molecules but
undergoes changes in timescales comparable to cloud lifetimes.
Generally, reactions occurring within the icy mantles should be of
great importance because most of the frozen molecules are accumulated
in the volume, not the surface. Cosmic-ray-induced alternations are
generally regarded as first-order reactions. We keep this approach
here. Molecules in a mantle are activated with a rate
where is the average time between two successive strikes by a Fe cosmic ray nucleus (calculated s by Hasegawa & Herbst 1993a) and u is the number of strikes needed to fully reorganize the grain mantle, i.e. to expose to an inner surface an amount of molecules equivalent to the total number of molecules per mantle. u is dependent on the porousness and thickness of an average mantle. The thickness is estimated to be 100 monolayers by Hasegawa & Herbst (1993b) and 25 layers by Turner (1998a). The former value is generally more consistent with our model. We assumed u = 200.
The reactions of the activated (exposed to a pore surface) reactants proceed with the rate coefficient (Hasegawa et al. 1992)
There is no knowledge about pores of interstellar dust; their size and numbers should be affected by the exact conditions in the dark core. However, they are limited by the number of mantle layers and size of the grain (see above). Undoubtedly, most of the pores are too small to provide a functioning reaction surface, while there should be only a few larger pores and with a higher possibility to be connected to the outside, thus indeed becoming at least partially a gully of the outer surface. For the calculation of the diffusion rate the number of adsorption sites per average pore in Eq. (10) is now taken to be 103. We take the same number for the pores per grain, . We consider this choice a possibly neutral estimate.
The factor that changes the reaction rate taken from Hasegawa et al. (1992) and Hasegawa & Herbst (1993a)
is inversely the number of isolated surfaces in a grain multiplied by
the number of adsorption sites on the available surface. For outer
surface (as given in Hasegawa et al. 1992, Eqs. (4) and (9)) this factor is
For pores we take
That is, we assume that there are thousand pores in an average grain with a surface of thousand adsorption sites each.
The molecules are assumed to be in the activated state from a Fe cosmic
ray hit to about the time of the next cosmic ray hit. The interval of
hits is assumed to be
(s). The rate coefficient (s-1) is thus
The chosen rate coefficients ensure that at any given time instant roughly 0.5% of mantle species are ``activated'' and are taking part in reactions on the pore surfaces. It is half the number of molecules assumed to be exposed to pore surface. The remaining half is assumed to be inactive for some reason (molecules which reside in pores too small for reactions and temporarily blocked sites). The molecular abundances (gas, surface, and mantle) produced by the model are only very slightly dependent on this percentage. It is because the chemical reactions on surface are anyway much faster than the radical production on surfaces by CR induced photons.
2.9 Hydrogen diffusion through mantle
Hydrogen diffusion from the surface to the mantle pores and from the
mantle to the surface is included in our model. We calculated the
diffusion rate coefficient assuming that H atoms and
molecules mostly reside on outer or inner surfaces. The hydrogen within
the mantle lattice is only a relatively rare intermediate state. Thus
where D is the diffusion coefficient, is the probability for diffusing species to find an appropriate surface. L is the length of the diffusion path (we take it cm), approximately half the thickness of the mantle. We adopt D for H atoms cm2 s-1 estimated by Awad et al. (2005). For the molecules we use the coefficient given by Strauss & Chen (1994) for diffusion in hexagonal ice extrapolated to 10 K ( cm2 s-1). This may sound rather tentative as the diffusion in amorphous ice is expected to be significantly slower. We could not find any estimates for in amorphous ice, however, in our model the real hydrogen content is rather independent of the exact value of as long as and the real diffusion rate is much shorter than the photodissociation rate - both of which are expected to be true.
We use the value for outward diffusion (assuming that on average 50% of all diffusion directions from pores lead outwards) and for inward diffusion because in our model we expect only 1/100 of the mantle volume to be occupied by the larger pores, able to host enough reactants on their surfaces for reactions to occur. This means that 99% of surface hydrogen diffusing inward returns to the surface and only 1% reach an inner pore to reside in. These probabilities result in an eqilibrium where the outward flux of hydrogen dominates. A release of hydrogen and the formation of less hydrogenated species is observed in numerous experiments involving photolysis of hydrogenated species in vacuum, e.g. Gerakines et al. (1996), Andersson & van Dischoeck (2008).
2.10 Direct ejection by cosmic rays of mantle molecules
Certainly there must be a process besides the transfer of gas molecules
to the surface and then to the mantle that works in the opposite
direction or a complete freeze-out will result, which is not observed.
We propose a simple cosmic-ray-driven mechanism which consistently
returns a proportion of the inner mantle into the gas phase. These are
the molecules believed to be directly, unselectively ejected by a hit
of a Fe cosmic ray nucleus. The rate of ejection (cm-3 s-1) is calculated by
is the proportion of species i in the grain mantle. is the number of mantle molecules ejected by this process in each hit. It should be a negligible number compared to the number of molecules desorbed by the cosmic-ray heating mechanism, estimated to be at least 106 by Hasegawa & Herbst (1993a). We take % of this value, that is 100 molecules per hit of cosmic-ray Fe nuclei. Gas and solid phase abundances are unaffected by the value of over several orders of magnitude. The returning of mantle molecules into the gas means that gas, surface, and mantle phases are in an interconnected equilibrium in our model. Molecule abundances in both surface and mantle directly affect the abundances in gas, where the species are easily observable.
2.11 Surface-to-mantle transition
We calculated the transformation of the surface species into mantle
species in the local thermodynamic equilibrium with a constant rate
is the total circulation rate of a surface molecule at an equilibrium state. The molecule accretes from the gas phase, undergoes surface diffusion, reactions, desorption acts, and is eventually buried by other molecules, and finally ends up in the mantle. We treat this coefficient as an unknown variable and use the model with the given, often assumed, input conditions (Sects. 2.2 to 2.5 and 2.7, and 2.10) to choose . We took s-1 to produce a ratio 30:1 between the total abundances of all atoms in the mantle and on the surface (except hydrogen). Gas, surface, and mantle abundances are very strongly affected by the value of , consequently it has to be chosen with a precision to decimals. H and are not incorporated into the mantle by this mechanism.
It is not a physically adequate approach to calculate the rate of surface-to-mantle transition. However, with so many assumptions and poorly known values used (e.g. Si, , u etc.) we consider it satisfactory for the aims of this paper. This approach also implies that the program might produce somewhat biased results in terms of absolute abundances (cm-3).
2.12 Desorption by chemical explosions
We do not include desorption by chemical explosions (Schutte Greenberg 1991) of grain mantles. In our current model the radicals released by photoprocess do not react rapidly and violently as is required by the explosion theory. They are entirely consumed by the reactions on the pore surfaces. Less than 1/100 of the mantle species are radicals in calculation results of our model, partially due to the incomplete reaction network. However, the photon flux in dark cores is significantly lower than in experiments producing explosion effects. Fully researching them requires a more sophisticated model, because the explosive reactions also directly affect the chemical composition of the icy mantle.
2.13 Model credibility
To focus our research on chemical processes at the local thermodynamic equilibrium we investigated a standard case of a very long-existing interstellar dark cloud core, completely isolated from the interstellar UV radiation field. The parameters regarding the various physical and chemical transformations of solid phase species (Sects. 2.1-2.11) were chosen in accordance with our understanding of the grain structure and processes. These parameters are able to influence the model output (calculated fractional abundances) in a wide range.
Observational results found in the literature were used to evaluate the rate of ejection of the mantle molecules by cosmic rays (Sect. 2.10) and surface-to-mantle transition (Sect. 2.11). These two rates essentially determine the abundances of gas phase species. We opted for long-lived (presumably LTE) interstellar cloud cores - dark, cold, and with metals heavily depleted on grains (Hollenbach et al. 2009; Turner 2000; Tafalla et al. 2004).
There are several aspects or assumptions in the model that require a credibility assessment. Most important is, how close the pore surface model represents the real conditions on interstellar grain mantles. Under certain (given) input conditions the model can produce feasible results. We assume that small pores are insignificant for the mantle chemistry, but they may affect the hydrogen content within the mantle to an unknown extent. Also, we do not take into account the ability of hydrogen to react while diffusing through the solid amorphous ice. Generally these effects should lead to an increased abundance of hydrogenated species at the expense of radicals (perhaps reducing the efficiency of desorption by explosions).
A second aspect regarding the validity of the model is the yields of various photoprocesses, especially the yield of the photodissociation inside the mantle. The model works best with the assumption that below the few surface layers dissociation essentially occurs on pore surfaces only. This is certainly not true, as dissociation basically is not a surface-related process. However, the assumption is justified by the higher possibility of recombination of dissociation products unable to diffuse away in a (amorphous) lattice. However, some the hydrogen released within the lattice may be able to diffuse, again increasing the radical content and supplying atomic hydrogen for reactions and diffusion to the surface. Si, Na, Mg and Fe have a very limited or nonexistent set of solid phase reactions and species and we do not include these elements in our discussion about results. Silicon is kept in the surface-mantle model in order to increase the total number of available reactions and thus improve the general accuracy of results for other species.
Table 1: Calculated fractional abundances of species in gas, surface, and mantle phases, and the normalized mantle-to-surface ratio.
We recall here that the simplified approach in the model is a quick solution to promote the discussion about the importance of hydrogen diffusion effects on grain mantle chemical composition. We emphasize that an important feature is the limited number of reactions available for each species in the mantle. A typical solid heavy molecule, for example, or SO, has one or two dissociation reactions by cosmic ray photons. Also, there are usually two to three reactions where the species is a product. There are usually no more than six reactions (production and destruction) involving a complicated molecule. It is a poor set of chemical transformations, especially in the mantles where hydrogenation reactions play a role of relatively little importance. Thus the calculated abundances for several important species are highly biased, as noted below.
The calculation results are summarised in Tables 1-3 and in Fig. 1.
The overall results of the model are shown in Table 1. Results (Table 2)
show that the mantle-to-surface (M/S) ratios for heavy elements are
different from each other. This means that, for example, M/S of
molecules is not directly comparable, because the M/S ratio of carbon itself is higher than that of oxygen (cf. Table 2)
and there is a larger abundance of carbon species overall in the
mantle. To be able to compare the mantle-to-surface ratios of different
molecules, we calculate a weighed average elemental M/S ratio (WR) for each individual specie according to
a (the index in a molecular formula) is the number of atoms of elements i to j in the molecule m. RE is the average M/S ratio for each element in the molecule m (C, N, O, Na, Mg, Si, S or Fe) from i through j. Because hydrogen wanders rather freely between the surface and the mantle, it is not counted here. We name the relation
Table 2: Total elemental fractional abundances in surface and mantle phases, calculated from Table 1 data.
Table 3: Fractional abundances for hydrogen species in surface and mantle phases.
The normalized mantle-to-surface (nM/S) ratio for selected species, obtained from Table 1 data.
|Open with DEXTER|
the normalized mantle-to-surface ratio (nM/S). For species whose fractional abundance is equal in the mantle and on the surface, nM/S is unity. If the concentration in mantle is higher this value is higher than unity and vice versa. The nM/S for selected species are graphically shown in Fig. 1.
The cosmic-ray-induced photodissociation yield in lower mantle layers should be rather low, because with yields larger than 10-3 in our simulations we obtain a mixture of compounds with water not the dominant molecule, which is inconsistent with the present knowledge about interstellar grain mantle composition. In this case the remaining oxygen is lost mostly to highly oxidized carbon-bearing species.
Hydrogen in its atomic, molecular, and chemically bound states is a special case because it is the only element capable of crossing the boundary between mantle and surface (Sect. 2.9). The data regarding ``solid'' phases of H are shown in Table 3. As a standard comparison data for this paper (Table 1) we use the results given by Hasegawa & Herbst (1993b), Table 1, t = 1.0E+8 yr. To properly compare the results we also calculated the molecular nM/S ratios for these data. Both works aim to investigate the solid phase chemistry in molecular clouds, and both have a 3-phase system (gas, surface, mantle). There are several important differences between the two models. Our model includes a reduced number of solid state species (see Sect. 2.1). The model of Hasegawa & Herbst lacks cosmic-ray-induced photodissociation of solid state species (Sect. 2.6), mantle-phase reactions (Sect. 2.8) and hydrogen diffusion (Sect. 2.9).
The main difference in the chemical modeling results is that we obtain more oxidized compounds, like in the mantle. They are mainly , HCOOH and less for carbon, for sulfur, more HCN, and less for nitrogen. Generally one can say the improvements of our solid phase model diminishes the dominance of hydrogenated species and leads to a much greater chemical diversity with high abundances of fully or partially oxidized species. Another substantial difference in our model is the high concentration of atomic and other radical species, both on the surface and in the mantle. It can be attributed to the inclusion of photodissociation reactions for all solid species. Our model produces a much better overall fit to observations for gas phase species. It is because molecules dominant in the mantle phase are slowly returned into gas and partially reprocessed, and a total freeze-out for any species cannot occur. The main results of our investigation are:
- Important molecules on the surface are , , . Those in the mantle are , , HCOOH, (Table 1). The most abundant molecule among both the surface and mantle species is water. However, the molecule containing most of the oxygen in the mantle is . In our opinion this represents the general shift from highly to poorly hydrogenated species by the photochemical processing of the mantle.
- The nM/S ratio for highly hydrogenated species is usually rather low, in the range 0.1-0.01 (Fig. 1). The most important species here are and other saturated carbon molecules. Methane has a high concentration on the surface, which significantly decreases in the mantle, releasing huge amounts of carbon now available for various oxidized and chain-like compounds.
- The highest nM/S is seen for cyanopolyyne related molecules, showing that mantle processing efficiently transforms saturated hydrocarbons to these species. The related HNC is consumed in the mantle, and the HCN/HNC gas abundance ratio of 1.9 is adequate for dark cores. Remarkable are the high nM/S values for oxidized sulfur species (OCS, SO, ), although these molecules have a very limited reaction set. The abundance of S oxides in the mantle exceeds that of . We expect that hydrogen depletion in the mantle should essentially explain the high abundance of OCS and sulfur oxides in hot core regions noted by many observers (e.g. Mookerjea et al. 2007; Hatchell et al. 1998). We predict that 3-phase models with a more complete reaction network will produce OCS and other sulfur species in an amount more consistent with the observations.
- According to calculations, the hydrogen depletion in the mantle has different effects on the abundances of carbon-oxygen compounds. Molecules HCOOH and HCO have large nM/S, for , , it is mediocre, and for , nM/S is less than unity. With a limited degree of certainty one can conclude that weakly hydrogenated C-O species are those favored by the conditions inside grain mantles, which is confirmed by observations of hot cores (e.g. Mookerjea et al. 2007). Molecules , and have a very deficient chemical reaction set and are not considered here.
- An increased radical content in the mantle. Species like C, O, CH, S and others (notably excluding nitrogen species) show nM/S ratios moderately higher than 1. The model-based explanation is that only a limited reaction network is available and that not all species produced by the dissociation of CR induced photons are able to readily react and generate stable molecules. However, the real radical content depends on photodissociation yields and on the exact conditions in the mantle (see Sect. 2.13).
Table 4: Calculated gas phase fractional abundances of molecules.
There is some experimental evidence available in the literature that backs up the results. Ferrante et al. (2008) and Garozzo et al. (2010) show experimentally that OCS is readily formed in laboratory-simulated interstellar conditions, while our model shows OCS to be a molecule with a very high nM/S ratio. However, we note again that OCS has a very limited reaction set. Several sources (e.g. Weaver et al. 2005, and references therein) note high abundances of organic molecules with H, C and O in regions of grain mantle disruption or evaporation, indicating that these species are formed on grains. This agrees with our model results, where most of these species have nM/S higher (much higher for some molecules) than unity (Fig. 1).
Our model shows that the hydrogen diffusion outwards from the mantle is much more significant than the inward diffusion. The species initially included in the mantle are hydrogen-rich due to the surface reactions. Photoprocessing releases the chemically bound hydrogen, it diffuses, escapes to the surface and, eventually, the gas phase. According to the model, this process is not overcome by the hydrogen diffusing into the mantle from outside.
Although the hydrogen diffusion coefficient in amorphous, dirty ice is most probably smaller than the one we used (in hexagonal ice), the model is insensitive to the exact value of D over several orders of magnitude. This is the case as long as the time of the diffusion is much shorter than the heavy molecule residing time in the mantle. Also, the model is rather insensitive to the exact value of inward diffusion probability , if it is kept lower (for example, by a factor of 2/3) than the outward . The solid phase model results are essentially determined by the reasonable estimate (see Sect. 2.9) that . With a degree of certainty one can say that the total direction of hydrogen diffusion at 10 K (in, immobile, or out) is determined by the structure of the mantle, not by the exact diffusion speed.
The hydrogen atoms released in mantle photodissociation are very slow to diffuse and are very reactive. However, among other molecules, there always form molecules that are generally much less reactive and able to easily diffuse and leave the mantle. While hydrogen is depleted, methane and other highly saturated hydrocarbons transform to cyanopolyynes and other hydrogen-poor carbon chain compounds (Fig. 1). Water slowly loses hydrogen and the remaining oxygen forms compounds with other metals (in our model it is mostly ).
The most serious counter-argument produced by our model is the overabundance of the molecule in calculation results. However, like many others, this important species is badly affected by a poor reaction set, its only sink reaction is the slow dissociation by cosmic-ray-induced photons. This is also the exact reason why chains with six and more carbon atoms were excluded from our calculations. The concentration of overproduced molecules would be reduced by interaction with radicals (including atomic hydrogen) in the mantle, meanwhile reducing the amount of radicals themselves. In the reaction set by Hasegawa et al. (1992) and Hasegawa & Herbst (1993a) most important radicals that lack solid phase reactions have the general formula where M is C, N, O, S or Si. They are generated by photodissociation in our model. The reactions between any complicated species and a variety of radicals are desirable for future mantle chemical models.
The differences in elemental abundance within mantle and surface (Table 2) can be clearly accounted for selective desorption effects by direct cosmic-ray-heating and by formation. Carbon is accumulated because it tends to form chain-like heavy molecules, while sulfur molecules generally have high binding energies and also do not easily desorb.
The hydrogen-poor conditions within the grain mantles may promote a much more colorful chemistry than is usually thought about interstellar matter. In the model it can be seen best with the C-O compounds. Reasoning from the calculation results we expect that the sub-surface mantle conditions would favor an assortment of (1) rather complex; and (2) oxidized species. If we attribute these properties to the heavier elements, we can conclude that a significant proportion of Si, P, and S should be highly oxidized and in molecules forming many bonds. Strong inorganic or semi-organic acids, their salts and esters with the organic alcohols, peroxy acids, etc., can be some of the compounds whose formation can be permitted by the lack of hydrogen in the mantle. To estimate the abundance of these compounds, one needs a greatly expanded solid-phase reaction set.
Last but not least, we offer our explanation to the observable ``lack'' of sulfur molecules in dark clouds (see e.g. Goicoechea et al. 2007; Wakelam et al. 2004; Palumbo et al. 1997). The calculated large abundance and n/MS of sulfur in highest oxidation state, namely , in our model output allow us to suppose that the conditions within the mantles are suitable for many oxidized sulfur compounds that can essentially be described as derivatives of the sulfuric and sulfurous acids. These would include esters, sulfinic and sulfonic acids, sulfoxides, sulfones, peroxysulfuric acids, sulfamides, sulfimides, sulfamic acids, their derivatives and acid salts. Plenty of these molecules are known to be stable at 273 K and should be durable in the mantles at 10 K if they are able to form. Organic species include unsaturated thiols, thioethers and thioketones. Other species would be carbon disulfide and those containing an S-S bond, already mentioned by Garozzo et al. (2010). Our idea is that a significant amount (say, some 40%) of sulfur is dissipated over a great variety of minor species (including many simpler radicals) and thus is hardly observable. An important assumption to explain why the derivatives are not observable in hot cores, is that when these species evaporate, they are reduced by hydrogen and dissociated in fragments by radiation. That is, these molecules lack gas production pathways. We suspect that a higher abundance of oxidized sulfur (mostly SO, , CS) bearing species should be observable in the middle evolution stages of a hot core. Indeed, these observation results are presented, for example, by Mookerjea et al. (2007), Wakelam et al. (2003), Wakelam et al. (2004), Chandler et al. (2005), Jiménez-Serra et al. (2005) with abundances of and, especially, SO typically larger than those of and comparable to OCS. A less well pronounced but similar effect can also be expected for the phosphorus family of compounds, as this element also is known to have a rich chemistry. One may argue that the ``disappeared'' sulfur in dark cores resides on grains in the form of oxides. This is however not confirmed by observations, e.g. Palumbo et al. (1997). It is well known that sulfur oxides tend to combine with water, which is thought to be the dominant molecule on interstellar grains. In our opinion mantle sulfur chemistry can be expected to lead to the various acidic derivatives, besides other more conventional compounds. Some investigations of the sulfur chemistry (Palumbo et al. 1997; Ferrante et al. 2008) insist that the most important sulfur molecule known in dark cores OCS forms in water-deficient ices, so sulfur acid derivatives perhaps would be what one can expect from water-rich ices.
- According to our model, which includes several basic assumptions (see Sect. 2), it is highly possible that chemical processes below the outer surface of interstellar grain mantles play an important role in the chemistry of dark molecular cloud cores.
- An important factor that should be further investigated is the hydrogen diffusion through the grain mantles. This, combined with the continued dissociation of molecules by cosmic ray induced photons, leads to an overall outward flux of hydrogen from the mantle. According to our model results, the more dense the grain mantle, the more efficient is the outward diffusion of hydrogen. A more diverse, H-poor chemistry is encouraged below the surface, explaining the abundance of at least some species observed in dark molecular clouds and hot molecular cores.
- A model based on the concept of pore surface reactions can at least partially describe the transformations occurring within the icy mantle.
- Fe nuclei of cosmic-rays can be a cause of physical alteration of the mantle structure, but other possibilities, like light cosmic-ray particles, CR induced UV photons, and slow thermal diffusion may provide alteration with generally similar chemical consequences.
- A combination of outer and inner mantle surface chemistry is able to produce a wide set of mantle species. It may ultimately lead to more accurate calculations of the composition of molecular clouds.
- Further research is required to clarify many factors that are only approximately estimated. These include photodissociation and desorption yields, H and diffusion rate in amorphous dirty ices, the thickness of the mantle, the properties and number of the pores, reaction mechanism inside the grain mantles, effectiveness of selective desorption mechanisms, etc.
- A chemically active and hydrogen poor environment in the mantle may explain the difficulties of observing sulfur in dark molecular cores. In grain mantles the rich chemistry of sulfur permits the formation of many various S molecules (mostly oxoacid derivatives) low on hydrogen, most of them with abundances too low to be observed. The large abundance of sulfur oxides in hot star-forming cores may be a direct consequence of these compounds being ejected into the gas phase.
- Aikawa, Y., Umebayashi, T., Nakano, T., & Miyama, S. M. 1997 ApJ, 486, L51 [Google Scholar]
- Andersson, S., & van Dishoeck, E. F. 2008, A&A, 491, 907 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Awad, Z., Chigai, T., Kimura, Y., Shalabiea, O. M., & Yamamoto, T. 2005, ApJ, 626, 262 [NASA ADS] [CrossRef] [Google Scholar]
- Baragiola, R. A., Bahr, D. A., & Vidal, R. A. 2003 in Astrophysics of Dust, ed. A. N. Witt [Google Scholar]
- Brown, P. D., Charnley, S. B. 1990, MNRAS, 244, 432 [NASA ADS] [Google Scholar]
- Chandler, C. J., Brogan, C. L., Shirley, Y. L., & Loinard, L. 2005, ApJ, 632, 371 [NASA ADS] [CrossRef] [Google Scholar]
- Chang, Q., Cuppen, H. M., & Herbst. E. 2007, A&A, 469, 973 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dalgarno, A. 2006, PNAS, 103, 12269 [NASA ADS] [CrossRef] [Google Scholar]
- Das, A., Chakrabarti, S. K., Acharyya, K., & Chakrabarti, S. 2008, New Astron., 13, 457 [NASA ADS] [Google Scholar]
- Ferrante, R. F., Moore, M. H., Spiliotis, M. M., & Hudson, R. L. 2008, ApJ, 684, 1210 [NASA ADS] [CrossRef] [Google Scholar]
- Freund, M. M., & Freund, F. T. 2006, ApJ, 639, 210 [NASA ADS] [CrossRef] [Google Scholar]
- Garozzo, M., Fulvio, D., Kanuchova, Z., Palumbo, M. E., & Strazzulla, G. 2010, A&A, 509, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gerakines, P. A., Schutte, W. A., & Ehrenfreund, P. 1996, A&A, 312, 289 [NASA ADS] [Google Scholar]
- Gerin, M., Falgarone, E., Joulain, K., et al. 1997, A&A, 318, 579 [NASA ADS] [Google Scholar]
- Goicoechea, J. R., Pety, J., Gerin, M., et al. 2007 [arXiv:astro-ph/0703393] [Google Scholar]
- Goldsmith, P. F., Li, D., & Krco, M. 2007, ApJ, 654, 273 [NASA ADS] [CrossRef] [Google Scholar]
- Hartquist, T. W., & Williams, D. A. 1990 Mon. Not. R. Astron. Soc., 247, 343 [Google Scholar]
- Hasegawa, T. I., & Herbst, E. 1993a, MNRAS, 261, 83 [NASA ADS] [CrossRef] [Google Scholar]
- Hasegawa, T. I., & Herbst, E. 1993b, MNRAS, 263, 589 [NASA ADS] [CrossRef] [Google Scholar]
- Hasegawa, T. I., Herbst, E., & Leung, C. E. 1992, ApJS, 82, 167 [NASA ADS] [CrossRef] [Google Scholar]
- Hatchell, J., Thompson, M. A., Millar, T. J., & MacDonald, G. H. 1998, A&A, 338, 713 [NASA ADS] [Google Scholar]
- Hollenbach, D., Kaufman, M. J., Bergin, E. A., & & Melnick, G. J. 2009, ApJ, 690, 1497 [CrossRef] [Google Scholar]
- Jenkins, E. B. 2009, ApJ, 700, 1299 [NASA ADS] [CrossRef] [Google Scholar]
- Jiménez-Serra, I., Martín-Pintado, J., Rodríguez-Franco, A., & Martín, S. 2005, ApJ, 627, L121 [NASA ADS] [CrossRef] [Google Scholar]
- Mookerjea, B., Casper, E., Mundy, L. G., & Looney, L. W. 2007, ApJ, 659, 447 [NASA ADS] [CrossRef] [Google Scholar]
- Nejad, L. A. M., & Wagenblast R. 1999, A&A, 350, 204 [NASA ADS] [Google Scholar]
- Öberg, K. I., van Dishoeck, E. F., & Linnartz, H. 2009, A&A, 496, 281 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Padovani, M., Galli, D., & Glassgold, A. E. 2009, A&A, 501, 619 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Palumbo, M. E., Geballe, T. R., & Tielens, A. G. G. M. 1997, ApJ, 479, 839 [NASA ADS] [CrossRef] [Google Scholar]
- Roberts, J. F., Rawlings, J. M. C., Viti, S., & Williams, D. A. 2007, MNRAS, 382, 733 [NASA ADS] [CrossRef] [Google Scholar]
- Sánchez-Salcedo, F. J., & Vázquez-Semadeni, E., & Gazol, A. 2002, ApJ, 577, 768 [NASA ADS] [CrossRef] [Google Scholar]
- Schutte, W. A., & Greenberg, J. M. 1991, A&A, 244, 190 [NASA ADS] [Google Scholar]
- Shalabiea, O. M., & Greenberg, J. M. 1994, A&A, 290, 266 [NASA ADS] [Google Scholar]
- Strauss, H. L., Chen, Z., & Loong., C.-K. 1994, J. Chem. Phys., 101, 7177 [NASA ADS] [CrossRef] [Google Scholar]
- Tafalla, M., Myers, P. C., Caselli, P., & Walmsley, C. M. 2004, A&A, 416, 191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Turner, B. E. 1998a, ApJ, 495, 804 [NASA ADS] [CrossRef] [Google Scholar]
- Turner, B. E. 1998b, ApJ, 501, 731 [NASA ADS] [CrossRef] [Google Scholar]
- Turner, B. E. 2000, ApJ, 542, 837 [NASA ADS] [CrossRef] [Google Scholar]
- van Weeren, R. J., Brinch, C., & Hogerheijde, M. R. 2009, A&A, 497, 773 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wakelam, V., Castets, C., Ceccarelli, C., et al. 2003, in Proc. Conf. Chemistry as a Diagnostic of Star Formation, ed. C. L. Curry, & M. Fich, 445 [Google Scholar]
- Wakelam, V., Castets, A., Ceccarelli, C., et al. 2004, A&A, 413, 609 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Widicus Weaver, S. L., Kelley, M. J., & Blake, G. A. 2005, Proc. 231st Symp. IAU, 223 [Google Scholar]
- Willacy, K., & Millar, T. J. 1998, MNRAS, 298, 562 [NASA ADS] [CrossRef] [Google Scholar]
- Willacy, K., & Williams, D. A. 1993, MNRAS, 260, 635 [NASA ADS] [Google Scholar]
- Woodall, J., Agúndez, M., Markwick-Kemper, A. J., & Millar, T. J. 2006, A&A, 466, 1197 [Google Scholar]
Table 1: Calculated fractional abundances of species in gas, surface, and mantle phases, and the normalized mantle-to-surface ratio.
Table 2: Total elemental fractional abundances in surface and mantle phases, calculated from Table 1 data.
Table 3: Fractional abundances for hydrogen species in surface and mantle phases.
Table 4: Calculated gas phase fractional abundances of molecules.
The normalized mantle-to-surface (nM/S) ratio for selected species, obtained from Table 1 data.
|Open with DEXTER|
|In the text|
Copyright ESO 2010
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.