RESUMO
The demand for green hydrogen has raised concerns over the availability of iridium used in oxygen evolution reaction catalysts. We identify catalysts with the aid of a machine learning-aided computational pipeline trained on more than 36,000 mixed metal oxides. The pipeline accurately predicts Pourbaix decomposition energy (Gpbx) from unrelaxed structures with a mean absolute error of 77 meV per atom, enabling us to screen 2070 new metallic oxides with respect to their prospective stability under acidic conditions. The search identifies Ru0.6Cr0.2Ti0.2O2 as a candidate having the promise of increased durability: experimentally, we find that it provides an overpotential of 267 mV at 100 mA cm-2 and that it operates at this current density for over 200 h and exhibits a rate of overpotential increase of 25 µV h-1. Surface density functional theory calculations reveal that Ti increases metal-oxygen covalency, a potential route to increased stability, while Cr lowers the energy barrier of the HOO* formation rate-determining step, increasing activity compared to RuO2 and reducing overpotential by 40 mV at 100 mA cm-2 while maintaining stability. In situ X-ray absorption spectroscopy and ex situ ptychography-scanning transmission X-ray microscopy show the evolution of a metastable structure during the reaction, slowing Ru mass dissolution by 20× and suppressing lattice oxygen participation by >60% compared to RuO2.
RESUMO
This paper introduces WhereWulff, a semiautonomous workflow for modeling the reactivity of catalyst surfaces. The workflow begins with a bulk optimization task that takes an initial bulk structure and returns the optimized bulk geometry and magnetic state, including stability under reaction conditions. The stable bulk structure is the input to a surface chemistry task that enumerates surfaces up to a user-specified maximum Miller index, computes relaxed surface energies for those surfaces, and then prioritizes those for subsequent adsorption energy calculations based on their contribution to the Wulff construction shape. The workflow handles computational resource constraints such as limited wall-time as well as automated job submission and analysis. We illustrate the workflow for oxygen evolution reaction (OER) intermediates on two double perovskites. WhereWulff nearly halved the number of Density Functional Theory (DFT) calculations from â¼240 to â¼132 by prioritizing terminations, up to a maximum Miller index of 1, based on surface stability. Additionally, it automatically handled the 180 additional resubmission jobs required to successfully converge 120+ atoms systems under a 48-h wall-time cluster constraint. There are four main use cases that we envision for WhereWulff: (1) as a first-principles source of truth to validate and update a closed-loop self-sustaining materials discovery pipeline, (2) as a data generation tool, (3) as an educational tool, allowing users (e.g., experimentalists) unfamiliar with OER modeling to probe materials they might be interested in before doing further in-domain analyses, (4) and finally, as a starting point for users to extend with reactions other than the OER, as part of a collaborative software community.
Assuntos
Oxigênio , Software , Fluxo de Trabalho , Adsorção , Fatores de TempoRESUMO
Machine learning (ML) methods have shown promise for discovering novel catalysts but are often restricted to specific chemical domains. Generalizable ML models require large and diverse training data sets, which exist for heterogeneous catalysis but not for homogeneous catalysis. The tmQM data set, which contains properties of 86,665 transition metal complexes calculated at the TPSSh/def2-SVP level of density functional theory (DFT), provided a promising training data set for homogeneous catalyst systems. However, we find that ML models trained on tmQM consistently underpredict the energies of a chemically distinct subset of the data. To address this, we present the tmQM_wB97MV data set, which filters out several structures in tmQM found to be missing hydrogens and recomputes the energies of all other structures at the ωB97M-V/def2-SVPD level of DFT. ML models trained on tmQM_wB97MV show no pattern of consistently incorrect predictions and much lower errors than those trained on tmQM. The ML models tested on tmQM_wB97MV were, from best to worst, GemNet-T > PaiNN ≈ SpinConv > SchNet. Performance consistently improves when using only neutral structures instead of the entire data set. However, while models saturate with only neutral structures, more data continue to improve the models when including charged species, indicating the importance of accurately capturing a range of oxidation states in future data generation and model development. Furthermore, a fine-tuning approach in which weights were initialized from models trained on OC20 led to drastic improvements in model performance, indicating transferability between ML strategies of heterogeneous and homogeneous systems.
Assuntos
Complexos de Coordenação , Redes Neurais de Computação , Aprendizado de Máquina , Hidrogênio , TermodinâmicaRESUMO
Ultra-high-density single-atom catalysts (UHD-SACs) present unique opportunities for harnessing cooperative effects between neighboring metal centers. However, the lack of tools to establish correlations between the density, types, and arrangements of isolated metal atoms and the support surface properties hinders efforts to engineer advanced material architectures. Here, this work precisely describes the metal center organization in various mono- and multimetallic UHD-SACs based on nitrogen-doped carbon (NC) supports by coupling transmission electron microscopy with tailored machine-learning methods (released as a user-friendly web app) and density functional theory simulations. This approach quantifies the non-negligible presence of multimers with increasing atom density, characterizes the size and shape of these low-nuclearity clusters, and identifies surface atom density criteria to ensure isolation. Further, it provides previously inaccessible experimental insights into coordination site arrangements in the NC host, uncovering a repulsive interaction that influences the disordered distribution of metal centers in UHD-SACs. This observation holds in multimetallic systems, where chemically-specific analysis quantifies the degree of intermixing. These fundamental insights into the materials chemistry of single-atom catalysts are crucial for designing catalytic systems with superior reactivity.
RESUMO
Heterogeneous catalysts in the form of atomically dispersed metals on a support provide the most efficient utilization of the active component, which is especially important for scarce and expensive late transition metals. These catalysts also enable unique opportunities to understand reaction pathways through detailed spectroscopic and computational studies. Here, we demonstrate that atomically dispersed iridium sites on indium tin oxide prepared via surface organometallic chemistry display exemplary catalytic activity in one of the most challenging electrochemical processes, the oxygen evolution reaction (OER). In situ X-ray absorption studies revealed the formation of IrV=O intermediate under OER conditions with an Ir-O distance of 1.83 Å. Modeling of the reaction mechanism indicates that IrV=O is likely a catalyst resting state, which is subsequently oxidized to IrVI enabling fast water nucleophilic attack and oxygen evolution. We anticipate that the applied strategy can be instrumental in preparing and studying a broad range of atomically dispersed transition metal catalysts on conductive oxides for (photo)electrochemical applications.
RESUMO
Understanding metal oxide MO2 (M = Ti, Ru, and Ir)-water interfaces is essential to assess the catalytic behavior of these materials. The present study analyzes the H2O-MO2 interactions at the most abundant (110) and (011) surfaces, at two different water coverages: isolated water molecules and full monolayer, by means of Perdew-Burke-Ernzerhof-D2 static calculations and ab initio molecular dynamics (AIMD) simulations. Results indicate that adsorption preferably occurs in its molecular form on (110)-TiO2 and in its dissociative form on (110)-RuO2 and (110)-IrO2. The opposite trend is observed at the (011) facet. This different behavior is related to the kind of octahedral distortion observed in the bulk of these materials (tetragonal elongation for TiO2 and tetragonal compression for RuO2 and IrO2) and to the different nature of the vacant sites created, axial on (110) and equatorial on (011). For the monolayer, additional effects such as cooperative H-bond interactions and cooperative adsorption come into play in determining the degree of deprotonation. For TiO2, AIMD indicates that the water monolayer is fully undissociated at both (110) and (011) surfaces, whereas for RuO2, water monolayer exhibits a 50% dissociation, the formation of H3O2 - motifs being essential. Finally, on (110)-IrO2, the main monolayer configuration is the fully dissociated one, whereas on (011)-IrO2, it exhibits a degree of dissociation that ranges between 50 and 75%. Overall, the present study shows that the degree of water dissociation results from a delicate balance between the H2O-MO2 intrinsic interaction and cooperative hydrogen bonding and adsorption effects.