No Arabic abstract
Advanced telescopes, such as ALMA and JWST, are likely to show that the chemical universe may be even more complex than currently observed, requiring astrochemical modelers to improve their models to account for the impact of new data. However, essential input information for gas-grain models, such as binding energies of molecules to the surface, have been derived experimentally only for a handful of species, leaving hundreds of species with highly uncertain estimates. We present in this paper a systematic study of the effect of uncertainties in the binding energies on an astrochemical two-phase model of a dark molecular cloud, using the rate equations approach. A list of recommended binding energy values based on a literature search of published data is presented. Thousands of simulations of dark cloud models were run, and in each simulation a value for the binding energy of hundreds of species was randomly chosen from a normal distribution. Our results show that the binding energy of H$_{2}$ is critical for the surface chemistry. For high binding energy, H$_{2}$ freezes out on the grain forming an H$_{2}$ ice. This is not physically realistic and we suggest a change in the rate equations. The abundance ranges found are in reasonable agreement with astronomical ice observations. Pearson correlation coefficients revealed that the binding energy of HCO, HNO, CH$_{2}$, and C correlate most strongly with the abundance of dominant ice species. Finally, the formation route of complex organic molecules was found to be sensitive to the branching ratios of H$_{2}$CO hydrogenation.
We study the effects of grain surface reactions on the chemistry of protoplanetary disks where gas, ice surface layers and icy mantles of dust grains are considered as three distinct phases. Gas phase and grain surface chemistry is found to be mainly driven by photo-reactions and dust temperature gradients. The icy disk interior has three distinct chemical regions: (i) the inner midplane with low FUV fluxes and warm dust ($gtrsim 15$K) that lead to the formation of complex organic molecules, (ii) the outer midplane with higher FUV from the ISM and cold dust where hydrogenation reactions dominate and, (iii) a molecular layer above the midplane but below the water condensation front where photodissociation of ices affects gas phase compositions. Some common radicals, e.g., CN and C$_2$H, exhibit a two-layered vertical structure and are abundant near the CO photodissociation front and near the water condensation front. The 3-phase approximation in general leads to lower vertical column densities than 2-phase models for many gas-phase molecules due to reduced desorption, e.g., H$_2$O, CO$_2$, HCN and HCOOH decrease by $sim$ two orders of magnitude. Finally, we find that many observed gas phase species originate near the water condensation front; photo-processes determine their column densities which do not vary significantly with key disk properties such as mass and dust/gas ratio.
Theoretical studies have revealed that dust grains are usually moving fast through the turbulent interstellar gas, which could have significant effects upon interstellar chemistry by modifying grain accretion. This effect is investigated in this work on the basis of numerical gas-grain chemical modeling. Major features of the grain motion effect in the typical environment of dark clouds (DC) can be summarised as follows: 1) decrease of gas-phase (both neutral and ionic) abundances and increase of surface abundances by up to 2-3 orders of magnitude; 2) shifts of the existing chemical jumps to earlier evolution ages for gas-phase species and to later ages for surface species by factors of about ten; 3) a few exceptional cases in which some species turn out to be insensitive to this effect and some other species can show opposite behaviors too. These effects usually begin to emerge from a typical DC model age of about 10^5 yr. The grain motion in a typical cold neutral medium (CNM) can help overcome the Coulomb repulsive barrier to enable effective accretion of cations onto positively charged grains. As a result, the grain motion greatly enhances the abundances of some gas-phase and surface species by factors up to 2-6 or more orders of magnitude in the CNM model. The grain motion effect in a typical molecular cloud (MC) is intermediate between that of the DC and CNM models, but with weaker strength. The grain motion is found to be important to consider in chemical simulations of typical interstellar medium.
In the denser and colder ($leq$20 K) regions of the interstellar medium (ISM), near-infrared observations have revealed the presence of sub-micron sized dust grains covered by several layers of Htextsubscript{2}O-dominated ices and dirtied by the presence of other volatile species. Whether a molecule is in the gas or solid-phase depends on its binding energy (BE) on ice surfaces. Thus, BEs are crucial parameters for the astrochemical models that aim to reproduce the observed evolution of the ISM chemistry. In general, BEs can be inferred either from experimental techniques or by theoretical computations. In this work, we present a reliable computational methodology to evaluate the BEs of a large set (21) of astrochemical relevant species. We considered different periodic surface models of both crystalline and amorphous nature to mimic the interstellar water ice mantles. Both models ensure that hydrogen bond cooperativity is fully taken into account at variance with the small ice cluster models. Density functional theory adopting both B3LYP-D3 and M06-2X functionals was used to predict the species/ice structure and their BE. As expected from the complexity of the ice surfaces, we found that each molecule can experience multiple BE values, which depend on its structure and position at the ice surface. A comparison of our computed data with literature data shows agreement in some cases and (large) differences in others. We discuss some astrophysical implications that show the importance of calculating BEs using more realistic interstellar ice surfaces to have reliable values for inclusion in the astrochemical models.
Abridged: We detail and benchmark two sophisticated chemical models developed by the Heidelberg and Bordeaux astrochemistry groups. The main goal of this study is to elaborate on a few well-described tests for state-of-the-art astrochemical codes covering a range of physical conditions and chemical processes, in particular those aimed at constraining current and future interferometric observations of protoplanetary disks. We consider three physical models: a cold molecular cloud core, a hot core, and an outer region of a T Tauri disk. Our chemical network (for both models) is based on the original gas-phase osu_03_2008 ratefile and includes gas-grain interactions and a set of surface reactions for the H-, O-, C-, S-, and N-bearing molecules. The benchmarking is performed with the increasing complexity of the considered processes: (1) the pure gas-phase chemistry, (2) the gas-phase chemistry with accretion and desorption, and (3) the full gas-grain model with surface reactions. Using atomic initial abundances with heavily depleted metals and hydrogen in its molecular form, the chemical evolution is modeled within 10^9 years. The time-dependent abundances calculated with the two chemical models are essentially the same for all considered physical cases and for all species, including the most complex polyatomic ions and organic molecules. This result however required a lot of efforts to make all necessary details consistent through the model runs, e.g. definition of the gas particle density, density of grain surface sites, the strength and shape of the UV radiation field, etc. The reference models and the benchmark setup, along with the two chemical codes and resulting time-dependent abundances are made publicly available in the Internet: http://www.mpia.de/homes/semenov/Chemistry_benchmark/home.html
The abundance of interstellar ice constituents is usually expressed with respect to the water ice because, in denser regions, a significant portion of the interstellar grain surface would be covered by water ice. The binding energy (BE), or adsorption energy of the interstellar species regulates the chemical complexity of the interstellar grain mantle. Due to the high abundance of water ice, the BE of surface species with the water is usually provided and widely used in astrochemical modeling. However, the hydrogen molecules would cover some part of the grain mantle in the denser and colder part of the interstellar medium. Even at around ~ 10K, few atoms and simple molecules with lower adsorption energies can migrate through the surface. The BE of the surface species with H2 substrate would be very different from that of a water substrate. However, adequate information regarding these differences is lacking. Here, we employ the quantum chemical calculation to provide the BE of 95 interstellar species with H2 substrate. These are representative of the BEs of species to a H2 overlayer on a grain surface. On average, we notice that the BE with the H2 monomer substrate is almost ten times lower than the BE of these species reported earlier with the H2 O c-tetramer configuration. The encounter desorption of H and H2 was introduced (with ED (H, H2 ) =45 K and ED (H2 , H2 ) =23 K) to have a realistic estimation of the abundances of the surface species in the colder and denser region. Our quantum chemical calculations yield higher adsorption energy of H2 than that of H (ED (H, H2 ) = 23 - 25 K and ED (H2, H2 ) =67 - 79 K). We further implement an astrochemical model to study the effect of encounter desorption with the resent realistic estimation. The encounter desorption of the N atom (calculations yield ED (N, H2 ) =83 K) is introduced to study the differences with its inclusion.