Quantifying Confidence in DFT Predicted Surface Pourbaix Diagrams and Associated Reaction Pathways for Chlorine Evolution

Abstract

Density functional theory calculations can be used to identify dominant reaction mechanisms. However, the dominant reaction mechanism is sensitive to choice of the exchange correlation functional. Here, we demonstrate using the example case of chlorine evolution reaction on rutile oxides, which can occur through at least three reaction mechanisms each mediated by different surface intermediates and active sites. We utilize Bayesian error estimation capabilities within the BEEF-vdW exchange correlation (XC) functional to quantify the uncertainty associated with predictions of the operative reaction mechanism by systematically propagating the uncertainty originating from DFT-computed adsorption free energies. We construct surface Pourbaix diagrams based on the calculated adsorption free energies for rutile oxides of Ru, Ir, Ti, Pt, V, Sn and Rh. We utilize confidence-value (c-value) to determine the degree of confidence in the predicted surface phase diagrams. Using the scaling relations between the adsorption energies of intermediates we construct a generalized Pourbaix diagram showing the stable surface composition as a function of potential and the oxygen binding energy on the cus site. This is used to consistently determine activity volcano relationships. We incorporate the uncertainty in linear scaling relations to quantify the confidence in generalized Pourbaix diagram and the associated activity. This allows us to compute the expectation limiting potential as a function of $\Delta E_{O^c}$, which provides a more appropriate activity measure incorporating DFT uncertainty. We show that the confidence in the classification problem of identifying the active reaction mechanism is much better than the prediction problem of determining catalytic activity.