Computational catalyst screening has the potential to significantly accelerate heterogeneous catalyst discovery. Typically this involves developing microkinetic reactor models that are based on parameters obtained from density functional theory (DFT) and transition state theory (TST). To reduce the large computational cost involved in computing various adsorption and transition state energies of all possible surface states on a large number of catalyst models, linear scaling relations for surface intermediates and transition states have been developed that only depend on a few, typically one or two descriptors such as the carbon atom adsorption energy. As a result, only the descriptor values have to be computed for various active site models to generate volcano curves in activity or selectivity. Unfortunately, for more complex chemistries the predictability of linear scaling relations is unknown. Also, the selection of descriptors is essentially a trial and error process. Here, using a database of adsorption energies of the surface species involved in the decarboxylation and decarbonylation of propionic acid over eight mono-metalic transition metal catalyst surfaces (Ni, Pt, Pd, Ru, Rh, Re, Cu, Ag), we tested if non-linear machine learning (ML) models can outperform the linear scaling relations in prediction accuracy when predicting the adsorption energy for various species on a metal surface based on data from the rest of the metal surfaces. We found linear scaling relations to hold well for predictions across metals with a mean absolute error (MAE) of 0.12 eV and ML methods being unable to outperform linear scaling relations. Only when predicting random species energies on a random metal surface (based on random training data), a currently unlikely scenario for catalyst screening, do kernel based ML models significantly outperform linear scaling relations. Finally, we propose an approach for automatic discovery of appropriate metal descriptors using principal component analysis (PCA).

This paper presents a systematic approach to quantify uncertainties of various quantities of interest (QoIs) in catalysis determined by microkinetic models developed from first principles. One of the main sources of uncertainty in any microkinetic simulation is attributed to the exchange-correlation approximations in density functional theory (DFT) used to calculate the rate constants for all elementary reaction steps within transition state theory. These DFT approximations are at the core of significant discrepancies between computational simulations and experimental measurements. Therefore, any model calculation should be accompanied by a measure of uncertainty. This work uses probability to represent uncertainties and latent variable models to develop probabilistic models that account for errors and correlations in DFT energies. These probabilistic models are further constrained to known reaction thermodynamics, and then propagated to QoIs such as turnover frequency (TOF), apparent activation barrier, and reaction orders. The proposed uncertainty quantification (UQ) framework is applied on the water-gas shift reaction (WGS: CO+H_2 O⇌CO_2+H_2). Specifically, this WGS study models a Pt/TiO2 catalyst as a Pt8 cluster supported on a rutile TiO2 (110) surface, where DFT energies are obtained using four separate functionals PBE, RPBE, HSE and M06L that each have their own justification for being appropriate for this study. In this way, information from three different classes of functionals, GGA (generalized gradient approximation), meta-GGA, and hybrid functionals, are used to generate a free energy probabilistic model. Although the uncertainty in model results spans orders of magnitude, a new approach is introduced to identify the dominant catalytic cycle under uncertainty. Overall, we find that our model captures various experimental kinetic data; however, the probability densities for TOF, apparent activation barrier and reaction orders are relatively wide due to different flavors of DFT predicting a wide variation of transition state and oxygen vacancy formation energies. Nevertheless, we can conclude with high certainty that a CO-promoted redox cycle is the dominant mechanism over the temperature range of 473 K to 600 K and that formate and carboxyl pathways are not playing any role for the investigated active site model.