Grain surface chemistry and its treatment in gas-grain chemical models is an area of large uncertainty. Whilst laboratory experiments are making progress, there is still much that is unknown about grain surface chemistry. Further, the results and parameters produced by experiment are often not easily translated to the rate equation approach most commonly used in astrochemical modelling. It is possible that statistical methods can reduce the uncertainty in grain surface chemical networks. In this work, a simple model of grain surface chemistry in a molecular cloud is developed and a Bayesian inference of the reactions rates is performed through MCMC sampling. Using observational data of the solid state abundances of major chemical species in molecular clouds, the posterior distributions for the rates of seven reactions producing CO, CO$_2$, CH$_3$OH and H$_2$O are calculated, in a form that is suitable for rate equation models. This represents a vital first step in the development of a method to infer reaction rates from observations of chemical abundances in astrophysical environments.

We present a publicly available, open source version of the time dependent gas-grain chemical code UCLCHEM. UCLCHEM propagates the abundances of chemical species through a large network of chemical reactions in a variety of physical conditions. The model is described in detail along with its applications. As an example of possible uses, UCLCHEM is used to explore the effect of protostellar collapse on commonly observed molecules and to study the behaviour of molecules in C-type shocks. We find the collapse of a simple Bonnor-Ebert sphere successfully reproduces most of the behaviour of CO,CS and NH$_3$ from cores observed by \citet{Tafalla2004} but cannot predict the behaviour of N$_2$H$^+$. In the C-shock application, we find that molecules can be categorized so that they can be useful observational tracers of shocks and their physical properties. Whilst many molecules are enhanced in shocked gas, we identify two groups of molecules in particular. A small number of molecules are enhanced by the sputtering of the ices as the shock propagates and then remain high in abundance throughout the shock. A second, larger set are also enhanced by sputtering but then are destroyed as the gas temperature rises. Through these applications the general applicability of UCLCHEM is demonstrated.

Understanding the physical conditions of dark molecular clouds and star forming regions is an inverse problem subject to complicated chemistry that varies non-linearly with time and the physical environment. In this paper we apply a Bayesian approach based on a Markov Chain Monte Carlo (MCMC) method for solving the non-linear inverse problems encountered in astrochemical modelling. We use observations for ice and gas species in dark molecular clouds and a time dependent, gas grain chemical model to infer the values of the physical and chemical parameters that characterize quiescent regions of molecular clouds. We show evidence that in high dimensional problems, MCMC algorithms provide a more efficient and complete solution than more classical strategies. The results of our MCMC method enable us to derive statistical estimates and uncertainties for the physical parameters of interest as a result of the Bayesian treatment.