Towards accurate prediction of catalytic activity in IrO 2 nanoclusters via first principles-based variable charge force field

• Published in: Journal of materials chemistry. A, Materials for energy and sustainability Vol. 3; no. 37; pp. 18970 - 18982
• Main Authors: Sen, F. G., Kinaci, A., Narayanan, B., Gray, S. K., Davis, M. J., Sankaranarayanan, S. K. R. S., Chan, M. K. Y.
• Format: Journal Article
• Language: English
• Published: United States Royal Society of Chemistry 2015
• ISSN: 2050-7488; 2050-7496
• DOI: 10.1039/C5TA04678E

IrO 2 is one of the most efficient electrocatalysts for the oxygen evolution reaction (OER), and also has other applications such as in pH sensors. Atomistic modeling of IrO 2 is critical for understanding the structure, chemistry, and nanoscale dynamics of IrO 2 in these applications. Such modeling has remained elusive due to the lack of an empirical force field (EFF) for IrO 2 . We introduce a first-principles-based EFF that couples the Morse (MS) potential with a variable charge equilibration method, QEq. The EFF parameters are optimized using a genetic algorithm (GA) on a density functional theory (DFT)-based training set. The resultant Morse plus QEq EFF, “MS-Q” in short, successfully reproduces the lattice parameters, elastic constants, binding energies, and internal coordinates of various polymorphs of IrO 2 from DFT calculations. More importantly, MS-Q accurately captures key metrics for evaluating structural and chemical properties of catalysts such as surface energetics, equilibrium shape, electrostatic charges, oxygen vacancy formation energies, relative stability of low index rutile IrO 2 surfaces, and pressure-induced phase transformations. The MS-Q EFF is used to predict the oxygen binding energy ( E ad ), a well-known descriptor for OER activity, on various sites of a nanocatalyst. We find E ad to be more favorable at low coordination sites, i.e. edges and corners, compared to planar facets; E ad is also correlated with charge transfer between the adsorbed O and nanocrystal, highlighting the importance of variable charge electrostatics in modeling catalysis on metal oxide surfaces. Our variable charge force field offers encouraging prospects for carrying out large-scale reactive simulations to evaluate catalytic performance of IrO 2 surfaces and nanostructures.