The adaptive kinetic Monte Carlo method was used to calculate the decomposition dynamics of a methanol molecule on Cu(100) at room temperature over a time scale of minutes. Mechanisms of reaction were found using minimum mode following saddle point searches based on forces and energies from density functional theory. Rates of reaction were calculated with harmonic transition state theory. The dynamics followed a pathway from , , , HCO, and finally to CO. Our calculations confirm that methanol decomposition starts with breaking the O–H bond followed by breaking C–H bonds in the dehydrogenated intermediates until CO is produced. The bridge site on the Cu(100) surface is the active site for scissoring chemical bonds. Reaction intermediates are mobile on the surface which allows them to find this active reaction site. This study illustrates how the adaptive kinetic Monte Carlo method can model the dynamics of surface chemistry from first principles.
Skip Nav Destination
Research Article| December 31 2009
Adaptive kinetic Monte Carlo simulation of methanol decomposition on Cu(100)
Lijun Xu, Donghai Mei, Graeme Henkelman; Adaptive kinetic Monte Carlo simulation of methanol decomposition on Cu(100). J. Chem. Phys. 28 December 2009; 131 (24): 244520. https://doi.org/10.1063/1.3281688
Download citation file: