Using Statistical Emulation and Knowledge of Grain-Surface Diffusion for Bayesian Inference of Reaction Rate Parameters: An Application to a Glycine Network

Abstract

There exists much uncertainty surrounding interstellar grain-surface chemistry. One of the major reaction mechanisms is grain-surface diffusion for which the the binding energy parameter for each species needs to be known. However, these values vary significantly across the literature which can lead to debate as to whether or not a particular reaction takes place via diffusion. In this work we employ Bayesian inference to use available ice abundances to estimate the reaction rates of the reactions in a chemical network that produces glycine. Using this we estimate the binding energy of a variety of important species in the network, by assuming that the reactions take place via diffusion. We use our understanding of the diffusion mechanism to reduce the dimensionality of the inference problem from 49 to 14, by demonstrating that reactions can be separated into classes. This dimensionality reduction makes the problem computationally feasible. A neural network statistical emulator is used to also help accelerate the Bayesian inference process substantially. The binding energies of most of the diffusive species of interest are found to match some of the disparate literature values, with the exceptions of atomic and diatomic hydrogen. The discrepancies with these two species are related to limitations of the physical and chemical model. However, the use of a dummy reaction of the form H + X -> HX is found to somewhat reduce the discrepancy with the binding energy of atomic hydrogen. Using the inferred binding energies in the full gas-grain version of UCLCHEM results in almost all the molecular abundances being recovered.

arXiv comment21 pages, 7 figures, accepted for publication in ApJ

References in corpus (33)

• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
• X-ray spectral modelling of the AGN obscuring region in the CDFS: Bayesian model selection and catalogue
• Observations of the Icy Universe
• The UMIST database for astrochemistry 2012
• Importance Nested Sampling and the MultiNest Algorithm
• Our astrochemical heritage
• A three-phase chemical model of hot cores: the formation of glycine
• On the Formation of CO2 and Other Interstellar Ices
• Binding energies: new values and impact on the efficiency of chemical desorption
• Dust as interstellar catalyst I. Quantifying the chemical desorption process
• Atom addition reactions in interstellar ice analogues
• Modeling Complex Organic Molecules in dense regions: Eley-Rideal and complex induced reaction
• Formation of Complex Organic Molecules in Cold Interstellar Environments through non-diffusive grain-surface and ice-mantle chemistry
• UCLCHEM: A Gas-Grain Chemical Code
• Binding energies of interstellar molecules on crystalline and amorphous models of water ice by ab-initio calculations
• A non-energetic mechanism for glycine formation in the interstellar medium
• Sensitivity analysis of grain surface chemistry to binding energies of ice species
• Gas-grain chemistry in cold interstellar cloud cores with a microscopic Monte Carlo approach to surface chemistry
• Chemical modelling of complex organic molecules with peptide-like bonds in star-forming regions
• Emulation of reionization simulations for Bayesian inference of astrophysics parameters using neural networks
• Fast cosmological parameter estimation using neural networks
• Direct measurement of desorption and diffusion energies of O and N atoms physisorbed on amorphous surfaces
• Binding Energy of Molecules on Water Ice: Laboratory Measurements and Modeling
• Bayesian emulator optimisation for cosmology: application to the Lyman-alpha forest
• Cosmological parameter estimation via iterative emulation of likelihoods
• ECoPANN: A Framework for Estimating Cosmological Parameters using Artificial Neural Networks
• A novel framework to study the impact of binding energy distributions on the chemistry of dust grains
• Chemulator: Fast, accurate thermochemistry for dynamical models through emulation
• Incorporating Astrochemistry into Molecular Line Modelling via Emulation
• Bayesian Inference of Reaction Rates in Icy Mantles
• Understanding the Formation and Evolution of Interstellar Ices: A Bayesian Approach
• Exploiting Network Topology for Accelerated Bayesian Inference of Grain Surface Reaction Networks