Chemical kinetic modeling in heterogeneous catalysis is advancing in its ability to provide qualitatively or even quantitatively accurate prediction of real-world behavior because of new advances in the physical and chemical representations of catalytic systems, estimation of relevant kinetic parameters, and capabilities in kinetic modeling. This Perspective describes current trends and future areas of advancement in chemical kinetic modeling, simulation, and parameter estimation: ranging from elementary step calculations to multiscale modeling to the role of advanced statistical methods for incorporating uncertainties in predictions. Multiple new or growing methodologies are covered, examples are provided, and forward-looking topics for advancement are noted.

With the advent of efficient electronic structure methods, effective continuum solvation methods have emerged as a way to, at least partially, include solvent effects into simulations without the need for expensive sampling over solvent degrees of freedom. The multipole moment expansion (MPE) model, while based on ideas initially put forward almost 100 years ago, has recently been updated for the needs of modern electronic structure calculations. Indeed, for an all-electron code relying on localized basis sets and-more importantly-a multipole moment expansion of the electrostatic potential, the MPE method presents a particularly cheap way of solving the macroscopic Poisson equation to determine the electrostatic response of a medium surrounding a solute. In addition to our implementation of the MPE model in the FHI-aims electronic structure theory code [ Blum , V. ; Comput. Phys. Commun. 2009 , 180 , 2175 - 2196 , DOI: 10.1016/j.cpc.2009.06.022 ], we describe novel algorithms for determining equidistributed points on the solvation cavity-defined as a charge density isosurface-and the determination of cavity surface and volume from just this collection of points and their local density gradients. We demonstrate the efficacy of our model on an analytically solvable test case, against high-accuracy finite-element calculations for a set of ≈140000 2D test cases, and finally against experimental solvation free energies of a number of neutral and singly charged molecular test sets [ Andreussi , O. ; J. Chem. Phys. 2012 , 136 , 064102 , DOI: 10.1063/1.3676407 ; Marenich , A. V. ; Minnesota Solvation Database , Version 2012; University of Minnesota : Minneapolis, MN, USA , 2012 . ]. In all test cases we find that our MPE approach compares very well with given references at computational overheads < 20% and sometimes much smaller compared to a plain self-consistency cycle.

We present high-pressure x-ray photoelectron spectroscopy (HP-XPS) and first-principles kinetic Monte Carlo study addressing the nature of the active surface in CO oxidation over Pd(100). Simultaneously measuring the chemical composition at the surface and in the near-surface gas phase, we reveal both O-covered pristine Pd(100) and a surface oxide as stable, highly active phases in the near-ambient regime accessible to HP-XPS. Surprisingly, no adsorbed CO can be detected during high CO(2) production rates, which can be explained by a combination of a remarkably short residence time of the CO molecule on the surface and mass-transfer limitations in the present setup.

The size-modified Poisson-Boltzmann (MPB) equation is an efficient implicit solvation model which also captures electrolytic solvent effects. It combines an account of the dielectric solvent response with a mean-field description of solvated finite-sized ions. We present a general solution scheme for the MPB equation based on a fast function-space-oriented Newton method and a Green's function preconditioned iterative linear solver. In contrast to popular multigrid solvers, this approach allows us to fully exploit specialized integration grids and optimized integration schemes. We describe a corresponding numerically efficient implementation for the full-potential density-functional theory (DFT) code FHI-aims. We show that together with an additional Stern layer correction the DFT+MPB approach can describe the mean activity coefficient of a KCl aqueous solution over a wide range of concentrations. The high sensitivity of the calculated activity coefficient on the employed ionic parameters thereby suggests to use extensively tabulated experimental activity coefficients of salt solutions for a systematic parametrization protocol.

The conceptual idea of degree of rate control (DRC) approaches is to identify the “rate limiting step” in a complex reaction network by evaluating how the overall rate of product formation changes when a small change is made in one of the kinetic parameters. We examine two definitions of this concept by applying it to first-principles kinetic Monte Carlo simulations of the CO oxidation at RuO₂(110). Instead of studying experimental data we examine simulations, because in them we know the surface structure, reaction mechanism, the rate constants, the coverage of the surface and the turn-over frequency at steady state. We can test whether the insights provided by the DRC are in agreement with the results of the simulations thus avoiding the uncertainties inherent in a comparison with experiment. We find that the information provided by using the DRC is non-trivial: It could not have been obtained from the knowledge of the reaction mechanism and of the magnitude of the rate constants alone. For the simulations the DRC provides furthermore guidance as to which aspects of the reaction mechanism should be treated accurately and which can be studied by less accurate and more efficient methods. We therefore conclude that a sensitivity analysis based on the DRC is a useful tool for understanding the propagation of errors from the electronic structure calculations to the statistical simulations in first-principles kinetic Monte Carlo simulations