A combined molecular dynamics (MD) and Monte Carlo approach was used to bridge time scales, enabling calculations of surface recombination rates for hydrogen on silica. MD was used for temperatures between 10 and 600 K at a high pressure of 10 atm, yielding recombination coefficients between 0.1 and 1. For the lower pressures more common in applications, low recombination rates make the corresponding calculations intractably expensive. A Monte Carlo technique, informed by the MD simulations, was designed to bridge the essential time scales. Distinct weak and strong surface binding sites for atomic hydrogen with densities of approximately 10 nm-2 were found using grand canonical Monte Carlo (GCMC) simulations, which, in turn, were used to obtain Eley-Rideal rate constants based on semiequilibrium theory. Monte Carlo variational transition state theory (MCVTST) was used to calculate Langmuir-Hinshelwood and thermal desorption rate constants for hydrogen atoms in strong and weak adsorption sites. Calculated reaction rates were used in a Langmuir kinetics model to estimate the recombination coefficient for T = 10-2000 K at gas-phase radical densities between 1012 and 1016 cm-3, yielding values of = 10-4-0.9.

Original languageEnglish (US)
Pages (from-to)24137-24147
Number of pages11
JournalJournal of Physical Chemistry C
Issue number42
StatePublished - Oct 27 2016

ASJC Scopus subject areas

  • Electronic, Optical and Magnetic Materials
  • General Energy
  • Physical and Theoretical Chemistry
  • Surfaces, Coatings and Films


Dive into the research topics of 'Hydrogen Recombination Rates on Silica from Atomic-Scale Calculations'. Together they form a unique fingerprint.

Cite this