A graph-theoretical kinetic Monte Carlo framework for on-lattice chemical kinetics.

There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

Abstract

Existing kinetic Monte Carlo (KMC) frameworks for the simulation of adsorption, desorption, diffusion, and reaction on a lattice often assume that each participating species occupies a single site and represent elementary events involving a maximum of two sites. However, these assumptions may be inadequate, especially in the case of complex chemistries, involving multidentate species or complex coverage and neighboring patterns between several lattice sites. We have developed a novel approach that employs graph-theoretical ideas to overcome these challenges and treat easily complex chemistries. As a benchmark, the Ziff-Gulari-Barshad system is simulated and comparisons of the computational times of the graph-theoretical KMC and a simpler KMC approach are made. Further, to demonstrate the capabilities of our framework, the water-gas shift chemistry on Pt(111) is simulated.

Related collections

Most cited references 33

Record: found
Abstract: found
Article: found

Is Open Access

The SIESTA method for ab initio order-N materials simulation

Alberto García, Daniel Sánchez-Portal, Pablo Ordejon … (2002)

We have developed and implemented a self-consistent density functional method using standard norm-conserving pseudopotentials and a flexible, numerical LCAO basis set, which includes multiple-zeta and polarization orbitals. Exchange and correlation are treated with the local spin density or generalized gradient approximations. The basis functions and the electron density are projected on a real-space grid, in order to calculate the Hartree and exchange-correlation potentials and matrix elements, with a number of operations that scales linearly with the size of the system. We use a modified energy functional, whose minimization produces orthogonal wavefunctions and the same energy and density as the Kohn-Sham energy functional, without the need of an explicit orthogonalization. Additionally, using localized Wannier-like electron wavefunctions allows the computation time and memory, required to minimize the energy, to also scale linearly with the size of the system. Forces and stresses are also calculated efficiently and accurately, thus allowing structural relaxation and molecular dynamics simulations.