اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدA Unified Monte Carlo Treatment of Gas-Grain Chemistry for Large Reaction Networks. II. A Multiphase Gas-Surface-Layered Bulk Model
43
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
The observed gas-phase molecular inventory of hot cores is believed to be significantly impacted by the products of chemistry in interstellar ices. In this study, we report the construction of a full macroscopic Monte Carlo model of both the gas-phase chemistry and the chemistry occurring in the icy mantles of interstellar grains. Our model treats icy grain mantles in a layer-by-layer manner, which incorporates laboratory data on ice desorption correctly. The ice treatment includes a distinction between a reactive ice surface and an inert bulk. The treatment also distinguishes between zeroth and first order desorption, and includes the entrapment of volatile species in more refractory ice mantles. We apply the model to the investigation of the chemistry in hot cores, in which a thick ice mantle built up during the previous cold phase of protostellar evolution undergoes surface reactions and is eventually evaporated. For the first time, the impact of a detailed multilayer approach to grain mantle formation on the warm-up chemistry is explored. The use of a multilayer ice structure has a mixed impact on the abundances of organic species formed during the warm-up phase. For example, the abundance of gaseous HCOOCH3 is lower in the multilayer model than in previous grain models that do not distinguish between layers (so-called two phase models). Other gaseous organic species formed in the warm-up phase are affected slightly. Finally, we find that the entrapment of volatile species in water ice can explain the two-jump behavior of H2CO previously found in observations of protostars.
قيم البحث
اقرأ أيضاً
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 cov
ering 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
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.
Reaction networks are often used to model interacting species in fields such as biochemistry and ecology. When the counts of the species are sufficiently large, the dynamics of their concentrations are typically modeled via a system of differential e
quations. However, when the counts of some species are small, the dynamics of the counts are typically modeled stochastically via a discrete state, continuous time Markov chain. A key quantity of interest for such models is the probability mass function of the process at some fixed time. Since paths of such models are relatively straightforward to simulate, we can estimate the probabilities by constructing an empirical distribution. However, the support of the distribution is often diffuse across a high-dimensional state space, where the dimension is equal to the number of species. Therefore generating an accurate empirical distribution can come with a large computational cost. We present a new Monte Carlo estimator that fundamentally improves on the classical Monte Carlo estimator described above. It also preserves much of classical Monte Carlos simplicity. The idea is basically one of conditional Monte Carlo. Our conditional Monte Carlo estimator has two parameters, and their choice critically affects the performance of the algorithm. Hence, a key contribution of the present work is that we demonstrate how to approximate optimal values for these parameters in an efficient manner. Moreover, we provide a central limit theorem for our estimator, which leads to approximate confidence intervals for its error.
We revisit the relation between the stellar surface density, the gas surface density, and the gas-phase metallicity of typical disk galaxies in the local Universe with the SDSS-IV/MaNGA survey, using the star formation rate surface density as an indi
cator for the gas surface density. We show that these three local parameters form a tight relationship, confirming previous works (e.g., by the PINGS and CALIFA surveys), but with a larger sample. We present a new local leaky-box model, assuming star formation history and chemical evolution is localized except for outflowing materials. We derive closed-form solutions for the evolution of stellar surface density, gas surface density and gas-phase metallicity, and show that these parameters form a tight relation independent of initial gas density and time. We show that, with canonical values of model parameters, this predicted relation match the observed one well. In addition, we briefly describe a pathway to improving the current semi-analytic models of galaxy formation by incorporating the local leaky-box model in the cosmological context, which can potentially explain simultaneously multiple properties of Milky Way-type disk galaxies, such as the size growth and the global stellar mass-gas metallicity relation.
Grain surface chemistry is key to the composition of protoplanetary disks around young stars. The temperature of grains depends on their size. We evaluate the impact of this temperature dependence on the disk chemistry. We model a moderately massive
disk with 16 different grain sizes. We use POLARIS to calculate the dust grain temperatures and the local UV flux. We model the chemistry using the 3-phase astrochemical code NAUTILUS. Photoprocesses are handled using frequency-dependent cross-sections, and a new method to account for self and mutual shielding. The multi-grain model outputs are compared to those of single-grain size models (0.1 $mu$m), with two different assumptions for their equivalent temperature. We find that the Langmuir-Hinshelwood (LH) mechanism at equilibrium temperature is not efficient to form H$_2$ at 3-4 scale heights ($H$), and adopt a parametric fit to a stochastic method to model H$_2$ formation instead. We find the molecular layer composition (1-3 $H$) to depend on the amount of remaining H atoms. Differences in molecular surface densities between single and multi-grain models are mostly due to what occurs above 1.5 $H$. At 100 au, models with colder grains produce H$_2$O and CH$_4$ ices in the midplane, and warmer ones produce more CO$_2$ ices, both allowing efficient depletion of C and O as soon as CO sticks on grain surfaces. Complex organic molecules (COMs) production is enhanced by the presence of warmer grains in the multi-grain models. Using a single grain model mimicking grain growth and dust settling fails to reproduce the complexity of gas-grain chemistry. Chemical models with a single grain size are sensitive to the adopted grain temperature, and cannot account for all expected effects. A spatial spread of the snowlines is expected to result from the ranges in grain temperature. The amplitude of the effects will depend on the dust disk mass.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
