UCLCHEMCMC: A MCMC Inference tool for Physical Parameters of Molecular Clouds

Abstract

We present the publicly available, open source code UCLCHEMCMC, designed to estimate physical parameters of an observed cloud of gas by combining Monte Carlo Markov Chain (MCMC) sampling with chemical and radiative transfer modeling. When given the observed values of different emission lines, UCLCHEMCMC runs a Bayesian parameter inference, using a MCMC algorithm to sample the likelihood and produce an estimate of the posterior probability distribution of the parameters. UCLCHEMCMC takes a full forward modeling approach, generating model observables from the physical parameters via chemical and radiative transfer modeling. While running UCLCHEMCMC, the created chemical models and radiative transfer code results are stored in an SQL database, preventing redundant model calculations in future inferences. This means that the more UCLCHEMCMC is used, the more efficient it becomes. Using UCLCHEM and RADEX, the increase of efficiency is nearly two orders of magnitude, going from 5185.33 \pm 1041.96 s for ten walkers to take one thousand steps when the database is empty, to 68.89 \pm 45.39 s when nearly all models requested are in the database. In order to demonstrate its usefulness we provide an example inference of UCLCHEMCMC to estimate the physical parameters of mock data, and perform two inferences on the well studied prestellar core, L1544, one of which show that it is important to consider the substructures of an object when determining which emission lines to use.

arXiv comment14 pages, 8 figures, 4 tables accepted by ApJ

References in corpus (16)

• emcee: The MCMC Hammer
• A computer program for fast non-LTE analysis of interstellar line spectra
• An atomic and molecular database for analysis of submillimetre line observations
• The UMIST database for astrochemistry 2012
• A KInetic Database for Astrochemistry (KIDA)
• The origin of complex organic molecules in prestellar cores
• Molecular ions in L1544. I. Kinematics
• Multilayer modeling of porous grain surface chemistry I. The GRAINOBLE model
• UCLCHEM: A Gas-Grain Chemical Code
• Sulphur chemistry in the L1544 pre-stellar core
• RUN DMC: An efficient, parallel code for analyzing Radial Velocity Observations using N-body Integrations and Differential Evolution Markov chain Monte Carlo
• Seeds of Life in Space (SOLIS). III. Zooming into the methanol peak of the pre-stellar core L1544
• Observations of CH${}_3$OH and CH${}_3$CHO in a Sample of Protostellar Outflow Sources
• Molecular transitions as probes of the physical conditions of extragalactic environments
• Incorporating Astrochemistry into Molecular Line Modelling via Emulation
• Molecular-Cloud-Scale Chemical Composition III: Constraints of Average Physical Properties through Chemical Models

Cited by in corpus (4)

• Machine learning-accelerated chemistry modeling of protoplanetary disks
• Neural network-based emulation of interstellar medium models
• Beetroots: spatially-regularized Bayesian inference of physical parameter maps -- Application to Orion
• The Si+SO${}_2$ collision and an extended network of neutral-neutral reactions between silicon and sulphur bearing species