Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA hybrid moment equation approach to gas-grain chemical modeling
368
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
[Context] The stochasticity of grain chemistry requires special care in modeling. Previously methods based on the modified rate equation, the master equation, the moment equation, and Monte Carlo simulations have been used. [Aims] We attempt to develop a systematic and efficient way to model the gas-grain chemistry with a large reaction network as accurately as possible. [Methods] We present a hybrid moment equation approach which is a general and automatic method where the generating function is used to generate the moment equations. For large reaction networks, the moment equation is cut off at the second order, and a switch scheme is used when the average population of certain species reaches 1. For small networks, the third order moments can also be utilized to achieve a higher accuracy. [Results] For physical conditions in which the surface reactions are important, our method provides a major improvement over the rate equation approach, when benchmarked against the rigorous Monte Carlo results. For either very low or very high temperatures, or large grain radii, results from the rate equation are similar to those from our new approach. Our method is faster than the Monte Carlo approach, but slower than the rate equation approach. [Conclusions] The hybrid moment equation approach with a cutoff and switch scheme is applicable to large gas-grain networks, and is accurate enough to be used for astrochemistry studies. The layered structure of the grain mantle could also be incorporated into this approach, although a full implementation of the grain micro-physics appears to be difficult.
rate research
Read More
We present a method to compute optical spectra and exciton binding energies of molecules and solids based on the solution of the Bethe-Salpeter equation (BSE) and the calculation of the screened Coulomb interaction in finite field. The method does not require the explicit evaluation of dielectric matrices nor of virtual electronic states, and can be easily applied without resorting to the random phase approximation. In addition it utilizes localized orbitals obtained from Bloch states using bisection techniques, thus greatly reducing the complexity of the calculation and enabling the efficient use of hybrid functionals to obtain single particle wavefunctions. We report exciton binding energies of several molecules and absorption spectra of condensed systems of unprecedented size, including water and ice samples with hundreds of atoms.
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
In this work we investigate the effects of ion accretion and size-dependent dust temperatures on the abundances of both gas-phase and grain-surface species. While past work has assumed a constant areal density for icy species, we show that this assumption is invalid and the chemical differentiation over grain sizes are significant. We use a gas-grain chemical code to numerically demonstrate this in two typical interstellar conditions: dark cloud (DC) and cold neutral medium (CNM). It is shown that, although the grain size distribution variation (but with the total grain surface area unchanged) has little effect on the gas-phase abundances, it can alter the abundances of some surface species by factors up to $sim2-4$ orders of magnitude. The areal densities of ice species are larger on smaller grains in the DC model as the consequence of ion accretion. However, the surface areal density evolution tracks are more complex in the CNM model due to the combined effects of ion accretion and dust temperature variation. The surface areal density differences between the smallest ($sim 0.01mu$m) and the biggest ($sim 0.2mu$m) grains can reach $sim$1 and $sim$5 orders of magnitude in the DC and CNM models, respectively.
Interstellar complex organic molecules (iCOMs) are assumed to be mainly formed on dust-grain surfaces. However, neutral gas-phase reactions in the interstellar medium (ISM) can play an important role. In this paper, by investigating the reaction between aldehydes and the cyano radical, we show that both formaldehyde (CH$_2$O) and acetaldehyde (CH$_3$CHO) can lead to the formation of formyl cyanide (HCOCN). Owing to accurate quantum-chemical computations followed by rate constant evaluations, we have been able to suggest and validate an effective mechanism for the formation of HCOCN, one of the molecules observed in the ISM. Quite interestingly, the mechanism starting from CH$_2$O is very effective at low temperature, while that involving CH$_3$CHO becomes more efficient at temperatures above 200 K.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
