SEARCH · Engineering Papers

Results for “Modeling Electrochemical Reaction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Electrochemical CO2 Reduction over Metal-/Nitrogen-Doped Graphene Single-Atom Catalysts Modeled Using the Grand-Canonical Density Functional Theory

Renewably driven, electrochemical conversion of carbon dioxide into value-added products is expected to be a critical tool in global decarbonization. However, theoretical studies based on the computational hydrogen electrode largely ignore the nonlinear effects of the applied potential on the calculated results, leading to inaccurate predictions of catalytic behavior or mechanistic pathways. Here, we use grand canonical density functional theory (GC-DFT) to model electrochemical CO2 reduction (CO2R) over metal- and nitrogen-doped graphene catalysts (MNCs) and explicitly include the effects of the applied potential. We used GC-DFT to compute the CO2 to CO reaction intermediate energies at -0.3, -0.7, and -1.2 VSHE catalyzed by MNCs each doped with 1 of the 10 3d block metals coordinated by four pyridinic nitrogen atoms. Our results predict that Sc-, Ti-, Co-, Cu-, and Zn-N4Cs effectively catalyze CO2R at moderate to large reducing potentials (-0.7 to -1.2 VSHE). ZnN4C is a particularly promising electrocatalyst for CO2R to CO both at low and moderate applied potentials based on our thermodynamic analysis. Our findings also explain the observed pH independence of CO production over FeN4C and predict that the rate-determining step of CO2R over FeN4C is not *CO2- formation but rather *CO desorption. Additionally, the GC-DFT-computed density of states analysis illustrates how the electronic states of MNCs and adsorbates change non-uniformly with applied potential, resulting in a significantly increased *CO2- stability relative to other intermediates and demonstrating that the formation of the adsorbed *CO2- anion is critical to CO2R activation. This work demonstrates how GC-DFT paves the way for physically realistic and accurate theoretical simulations of reacting electrochemical systems.

CO2 reduction↗

DOE OSTI2025

Uncovering the True Active Sites in Ni–N–C Catalysts for CO 2 Electroreduction

Understanding and designing active sites in single-atom catalysts (SACs) requires going beyond static models to capture their dynamic evolution under realistic electrochemical conditions. Here, in this work, we develop an integrated theoretical framework that accounts for operational conditions, by combining grand canonical density functional theory (GC-DFT) with machine-learning-accelerated sampling, to uncover structure–activity–stability relationships in Ni–N–C SACs for the CO 2 reduction reaction (CO 2 RR). A library of NiN x C 4–x (x = 0–4) motifs─representing coordination defects likely formed during high-temperature synthesis─was systematically evaluated. Under working conditions, these sites were found to undergo hydrogenation, and NiN 3 C 1_ H 1 was identified as the most probable active site. At reducing potentials, hydrogen adsorbs spontaneously at C–Ni bridge sites rather than Ni top sites, while subsurface hydrogen facilitates bent CO 2 adsorption crucial for activation. High CO 2 RR selectivity toward CO arises from site separation: Ni centers drive CO2RR, while the hydrogen evolution reaction (HER) occurs at the C–Ni bridge or N sites and from thermodynamic suppression of HER at moderate hydrogen coverage. At more negative potentials, a shift in the CO 2 RR rate-determining process (RDP) and Ni out-of-surface displacement induced by coadsorption of H and H 2 O jointly reduce activity and selectivity. Thus, both the high CO2RR selectivity of Ni–N–C catalysts and its reversal with more negative potentials can be rationalized by accounting for hydrogenated surfaces. This highlights the necessity of modeling realistic; in situ conditions. This framework provides generalizable insights into the dynamic behavior of active sites in SACs, offering guidance for the rational design of active and robust catalysts for a wide range of electrochemical reactions.

25 ENERGY STORAGE↗

DOE OSTI2024

Effect of H + Exchange and Surface Impurities on Bulk and Interfacial Electrochemistry of Garnet Solid Electrolytes

Contact loss and current constriction pose significant challenges at the Li metal interface of solid-state batteries. For garnet-structured Li 7 La 3 Zr 2 O 12 (LLZO), these effects are amplified by Li + /H + exchange and surface contamination reactions, which lead to conductivity losses and poor Li wetting. In this study, we utilize a variety of surface treatment processes across 37 cells to selectively induce proton exchange and contamination reactions in LLZO. The resulting bulk and surface chemistry is systematically characterized and correlated to changes in electrochemical properties. Additionally, we combine impedance analysis and finite element method modeling to deconvolute sources of impedance contributions at the Li metal interface. Specifically, we show that constriction impedance at the Li metal interface arises not solely from voids, but also from ionically-resistive surface contaminants. Further, these findings emphasize the connection between ionic conductivity and constriction, demonstrating that micron-scale ionically-resistive components increase constriction even with identical contact geometries. Finally, we leverage our comprehensive dataset to highlight unstable overpotential growth as a failure mechanism, additionally showing that the phase of a cell’s impedance is a sensitive indicator for the onset of interfacial instability. Overall, this study clarifies the impacts of proton exchange and surface contamination on electrochemical properties at the Li|solid electrolyte interface and elucidates insights that are generalizable to other solid-state battery systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Elucidating the complex interplay between thermodynamics, kinetics, and electrochemistry in battery electrodes through phase-field modeling

Abstract This article highlights applications of phase-field modeling to electrochemical systems, with a focus on battery electrodes. We first provide an overview on the physical processes involved in electrochemical systems and applications of the phase-field approach to understand the thermodynamic and kinetic mechanisms underlying these processes. We employ two examples to highlight how realistic thermodynamics and kinetics can naturally be incorporated into phase-field modeling of electrochemical processes. One is a composite battery cathode with an intercalation compound (Li x FePO 4 ) as the electrochemically active material, and the other is a displacement reaction compound (Li–Cu–TiS 2 ). With the input parameters mostly from atomistic calculations and experimental measurements, phase-field simulations allowed us to untangle the interactions among transport, reaction, electricity, chemistry, and thermodynamics that lead to highly complex evolution of the materials within battery electrodes. The implications of these observations for battery performance and degradation are discussed. Graphical abstract

Andrews, W. Beck (ORCID:0000000287824621)↗

DOE OSTI2026

Multi-Scale Modeling and Prototype Development for Electrochemical CO2 Reduction (CRADA Final Report)

In this CRADA project, Lawrence Livermore National Laboratory, Stanford University, SLAC National Laboratory, and TotalEnergies collaboratively executed a multidisciplinary investigation of electrochemical reduction of CO2 to produce sustainable fuels and chemicals. Overall, the project led to an increased understanding of the fundamental processes involved in CO2 electrolysis, from the atomistic scale to the full electrolyzer device scale, ultimately leading to design guidelines for CO2 electrolyzers that will help in their future commercialization. As the model systems, Ag- and Cu-based catalysts were investigated in various forms depending on the electrochemical platform that was utilized to study the activity, selectivity, and durability towards electrochemical CO2 reduction. By employing experimental, theoretical, and computational techniques, the project team experimentally validated multi-physics models, evaluated the experimental levers that lead to increased electrolyzer reaction selectivity and energy efficiency, and used computational optimization to design higher performance electrodes. The learnings of this project were extensively documented in publicly available peer-reviewed journal publications and conference presentations, which serve as a foundation for further work to build from.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Exploring non-adiabaticity to CO reduction reaction through ab initio molecular dynamics simulation

Non-adiabatic chemical reaction refers to the electronic excitation during reactions. This effect cannot be modeled by the ground-state Born-Oppenheimer molecular dynamics (BO-MD), where the electronic structure is at the ground state for every step of ions' movement. Although the non-adiabatic effect has been explored extensively in gas phase reactions, its role in electrochemical reactions, such as water splitting and CO 2 reduction, in electrolyte has been rarely explored. On the other hand, electrochemical reactions usually involve electron transport; thus, a non-adiabatic process can naturally play a significant role. In this work, using one-step CO 2 reduction as an example, we investigated the role of the non-adiabatic effect in the reaction. The reaction barriers were computed by adiabatic BO-MD and non-adiabatic real-time time dependent density functional theory (rt-TDDFT). We found that by including the non-adiabatic effect, rt-TDDFT could increase the reaction barrier up to 6% compared to the BO-MD calculated barrier when the solvent model is used to represent water. Simulations were carried out using explicit water molecules around the reaction site under different overpotentials, and similar non-adiabatic effects were found.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Intrinsic chemical reactivity of solid-electrolyte interphase components in silicon–lithium alloy anode batteries probed by FTIR spectroscopy

Herein we report the solid reaction products from the chemical reaction of aprotic battery electrolyte and three purported components of the Si-based anode SEI : SiO 2 nanoparticles (NPs), lithium silicate (Li x SiO y ) powders, and Si NPs. We use FTIR and classical molecular dynamics/density functional perturbation theory to assess the solid products remaining with these model materials after exposure to electrolyte. The absence of electrochemical bias provides a view of the chemical speciation resulting from early-stage chemical reactivity during battery assembly as well as under open circuit storage conditions. We believe these species represent the initial stages of SEI growth and predict they likely drive subsequent chemical and electrochemical reactions by controlling molecular interactons at the Si active material interface. We find that nominally equivalent materials react differently even before any electrochemistry is performed (e.g., acidic SiO 2 dissolves whereas alkaline SiO 2 is relatively robust), and derive new understanding of the chemical species that could and could not form stable SEI components in Si-based anodes. These results can be used to inform how to passivate Si anode surfaces and potentially generate an artificially engineered SEI that would be stable and enable next-generation battery anodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

BEAST DB: Grand-Canonical Database of Electrocatalyst Properties

We present BEAST DB, an open-source database comprised of ab initio electrochemical data computed using grand-canonical density functional theory in implicit solvent at consistent calculation parameters. The database contains over 20,000 surface calculations and covers a broad set of heterogeneous catalyst materials and electrochemical reactions. Calculations were performed at self-consistent fixed potential as well as constant charge to facilitate comparisons to the computational hydrogen electrode. This article presents common use cases of the database to rationalize trends in catalyst activity, screen catalyst material spaces, understand elementary mechanistic steps, analyze the electronic structure, and train machine learning models to predict higher fidelity properties. Users can interact graphically with the database by querying for individual calculations to gain a granular understanding of reaction steps or by querying for an entire reaction pathway on a given material using an interactive reaction pathway tool. BEAST DB will be periodically updated, with planned future updates to include advanced electronic structure data, surface speciation studies, and greater reaction coverage.

database↗

DOE OSTI2020

Thermally Driven Interfacial Degradation between Li 7 La 3 Zr 2 O 12 Electrolyte and LiNi 0.6 Mn 0.2 Co 0.2 O 2 Cathode

Solid-state batteries offer higher energy density and enhanced safety compared to the present lithium-ion batteries using liquid electrolytes. A challenge to implement them is the high resistances, especially at the solid electrolyte interface with the cathode. Sintering at elevated temperature is needed in order to get good contact between the ceramic solid electrolyte and oxide cathodes and thus to reduce contact resistances. Many solid electrolyte and cathode materials react to form secondary phases. It is necessary to find out which phases arise as a result of interface sintering and evaluate their effect on electrochemical properties. Herein, we assessed the interfacial reactions between LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622) and Li 7 La 3 Zr 2 O 12 (LLZO) as a function of temperature in air. We prepared model systems by depositing thin-film NMC622 cathode layers on LLZO pellets. The thin-film cathode approach enabled us to use interface-sensitive techniques such as X-ray absorption spectroscopy in the near-edge as well as the extended regimes and identify the onset of detrimental reactions. We found that the Ni and Co chemical environments change already at moderate temperatures, on-setting from 500 °C and becoming especially prominent at 700 °C. By analyzing spectroscopy results along with X-ray diffraction, we identified Li 2 CO 3 , La 2 Zr 2 O 7 , and La(Ni,Co)O 3 as the secondary phases that formed at 700 °C. The interfacial resistance for Li transfer, measured by electrochemical impedance spectroscopy, increases significantly upon the onset and evolution of the detected interface chemistry. Our findings suggest that limiting the bonding temperature and avoiding CO 2 in the sintering environment can help to remedy the interfacial degradation.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Plating and Stripping of Lithium Metal Stabilized by a Block Copolymer Electrolyte: Local Current Density Measurement and Modeling

Here, plating and stripping of lithium protrusions in lithium metal symmetric cells containing a solid block copolymer electrolyte was studied as a function of time in 3D using time-resolved X-ray tomography. These measurements enabled determination of the spatial variation in current densities at the plating and stripping electrodes. The initial interelectrode distance was 27 μm. Correlation functions were calculated to reveal the relationships between current densities at the two electrodes and local electrolyte thickness. Current densities at opposing electrode locations during protrusion growth is uncorrelated until the local interelectrode distance decreases to less than 6 μm, just before the cell shorts. Mass balance was used to determine the area from which lithium ions that form a protrusion were stripped. Computational modeling of the plating and stripping process reveals the interplay between electrochemical and mechanical driving forces and their effect on nonuniform current distribution. Model predictions were compared with experiments without resorting to any adjustable parameters. The computed correlation functions were in qualitative agreement with experiments. Finally, the model was used to calculate contour plots of electrochemical potential within the electrolyte, shedding light on how geometry, salt concentration, interelectrode distance, and mechanical stress influence local rates of electrochemical reaction.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Minimum conditions for accurate modeling of urea production via co-electrolysis

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Ni coarsening under humid atmosphere in the electrode of solid oxide cells: A combined study of density-functional theory and phase-field modeling

Ni coarsening in solid oxide cell (SOC) electrodes is known to be significantly faster under a humid atmosphere. The underlying mechanisms, though, have not been fully understood. In this work, we examine the surface diffusion of Ni(OH) x by a combination of density-functional theory and phase-field modeling. Here, density-functional theory is used to evaluate the adsorption and diffusion of the Ni(OH) x species on Ni (111) surface, and the results are used to obtain the effective surface diffusivity of Ni versus steam and hydrogen gas partial pressures. This diffusivity is used as an input for a phase-field model to investigate the Ni coarsening in SOC electrodes. It is found that Ni(OH) x formed through chemical reaction cannot accelerate coarsening unless an extremely high steam to hydrogen ratio is reached. However, if Ni(OH) x is formed through electrochemical reactions near triple phase boundaries (TPBs), surface diffusion of Ni(OH) x may cause faster Ni coarsening in fuel cell mode under a large overpotential. Specifically, surface diffusion of Ni(OH) may be comparable to or faster than that of Ni under an overpotential that is large but still possible under fuel cell operating conditions.

25 ENERGY STORAGE↗

DOE OSTI2022

A Critical Review of Radiolysis Issues in Water-Cooled Fission and Fusion Reactors: Part I, Assessment of Radiolysis Models

A critical review is presented on modeling of the radiolysis of the coolant water in nuclear power reactors with emphasis on ITER. The review is presented in two parts: In Part I, we assess previous work in terms of compliance with important chemical principles and conclude that no model proposed to date is completely satisfactory, in this regard. Thus, some reactions that have been proposed in various radiolysis models are not elementary in nature and can be decomposed into two or more elementary reactions, some of which are already included in the models. These reactions must be removed in formulating a viable model. Furthermore, elementary reactions between species of like charge are also commonly included, but they can be discounted upon the basis of Coulombic repulsion under the prevailing conditions (T < 350 °C) and must also be removed. Likewise, it is concluded that the current state of knowledge with respect to radiolytic yields (i.e., G-values) is also unsatisfactory. More work is required to ensure that the yields used in radiolysis models are truly “primary” yields corresponding to a time scale of nanoseconds or less. This is necessary to ensure that the impact of the reactions that occur outside of the spurs (ionizing particle tracks in the medium) are not counted twice. In Part II, the authors review the use of the radiolysis models coupled with electrochemical models to predict the water chemistry, corrosion potential, crack growth rate in Type 304 SS, and accumulated damage in the coolant circuits of boiling water reactors, pressurized water reactors, and the test fusion reactor, ITER. Based on experience with fission reactors, the emphasis should be placed on the control of the electrochemical corrosion potential because it is the parameter that best describes the state of corrosion in coolant circuits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Process model for multilayer slide coating of polymer electrolyte membrane fuel cells

Slide coating is a precision method suitable for depositing multiple liquid-film layers simultaneously. Originally developed in the photographic film industry, it has been deployed for manufacturing of other products that benefit from multilayer coatings. One emerging application is the manufacture of polymer electrolyte membrane fuel cells (PEMFCs), which are used to produce electricity through electrochemical reactions of hydrogen and oxygen gas. The membrane-electrode assembly (MEA), in which key electrochemical reactions occur, consists of three layers that are typically deposited separately in serial fashion and then laminated together to form the three-layer MEA, i.e., three sequential steps of coat and dry. Adapting the process to simultaneous, multilayer slide coating of all three layers will save equipment cost and space while minimizing possible exposure to contaminants during transition between the steps. We are developing a three-layer slide coating model to aid the manufacturing process design of PEMFC. The model accounts for rheology of each layer, which typically exhibit shear thinning behavior. Model predictions are used to investigate simultaneous coatability of catalyst inks and to determine the best layer-by-layer ink selection.

25 ENERGY STORAGE↗

DOE OSTI2024

From Fundamental Interfacial Reaction Kinetics to Macroscopic Current–Voltage Characteristics: Case Study of Solid Acid Fuel Cell Limitations and Possibilities

The unique properties of solid acid electrolytes, in particular CsH 2 PO 4 , are in many ways ideal for fuel cell operation. However, the technology is constrained by high cathode overpotentials. Here a simplified cathode geometry is employed to obtain the fundamental electrochemical parameters (exchange current density and charge transfer coefficient) describing the oxygen reduction reaction (ORR) at the CsH 2 PO 4 -Pt-gas interface. The parameters are incorporated into a 1D model of the voltage–current characteristics of realistic SAFC cathodes, which reproduced the measured polarization behavior of such cathodes without recourse to fitting adjustable parameters. Following this validation, the model is utilized to evaluate the impact of changes to cathode properties, microstructure, and operating conditions. Of these, the charge transfer coefficient, measured to have a value of ≈0.6 for ORR on Pt in the SAFC cathode environment, is found to have the greatest impact on power output. Nevertheless, even without material modifications, a combination of microstructural and operational modifications are identified with projected performance metrics meeting Department of Energy targets (0.8 V at 300 mA cm –2 , and peak power density of 1 W cm –2 ), albeit at high Pt loadings. However, the analysis indicates that truly meaningful advances will likely necessitate the discovery of alternative ORR catalysts.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Predictive Models and Novel Accelerated Tests for the Reliability of Cell Metallization in Photovoltaic Modules (Final Report)

Studies of metallization corrosion m photovoltaics have mainly been limited to comparisons of modules placed in accelerated chambers to fielded modules [1]. Damp Heat accelerated tests and phenomenological equations [2] are used to assess metallization corrosion without understanding the effect of UV light and temperature and humidity cycles on encapsulant adhesion degradation. The roles of encapsulant in-and out-diffusions of moisture and encapsulant impurities are important. Furthermore, few photovoltaic metallization corrosion studies included the role of bias and leakage currents, which are crucial in the electrochemical reaction of metallization. Leakage currents can highly accelerate the corrosion mechanism and are important to include in the studies of corrosion. In our research plan, we will address the following gaps in the PV community's understanding of metallization corrosion: (1) metallization corrosion with bias, humidity, and impurities in the encapsulant or metallization; (2) humidity diffusion through fresh and degraded encapsulants; and (3) comparison of model predictions with outdoor field modules and SunPower's extensive data for its back-contact and front-contact fleets [2] along with NREL's store of >20 year old modules. Our goal is to build models and accelerated tests to predict long term degradation of metallization corrosion of photovoltaic modules in the field. Our studies will include metals used in c-Si solar cells (Cu, Ag, and Al) and commonly used encapsulants (EV A (ethylene vinyl acetate), TPO (thermoplastic olefin), and silicone).

14 SOLAR ENERGY↗

DOE OSTI2023

Chemical Kinetic Method for Active-Site Quantification in Fe-N-C Catalysts and Correlation with Molecular Probe and Spectroscopic Site-Counting Methods

Mononuclear Fe ions ligated by nitrogen (FeN x ) dispersed on nitrogen-doped carbon (Fe-N-C) serve as active centers for electrocatalytic O 2 reduction and thermocatalytic aerobic oxidations. Despite their promise as replacements for precious metals in a variety of practical applications, such as fuel cells, the discovery of new Fe-N-C catalysts has relied primarily on empirical approaches. In this context, the development of quantitative structure–reactivity relationships and benchmarking of catalysts prepared by different synthetic routes and by different laboratories would be facilitated by the broader adoption of methods to quantify atomically dispersed FeN x active centers. In the present study, we develop a kinetic probe reaction method that uses the aerobic oxidation of a model hydroquinone substrate to quantify the density of FeN x centers in Fe- N-C catalysts. The kinetic method is compared with low-temperature Mössbauer spectroscopy, CO pulse chemisorption, and electrochemical reductive stripping of NO derived from NO 2 – on a suite of Fe-N-C catalysts prepared by diverse routes and featuring either the exclusive presence of Fe as FeN x sites or the coexistence of aggregated Fe species in addition to FeN x . The FeN x site densities derived from the kinetic method correlate well with those obtained from CO pulse chemisorption and Mössbauer spectroscopy. The broad survey of Fe-N-C materials also reveals the presence of outliers and challenges associated with each site quantification approach. Furthermore, the kinetic method developed here does not require pretreatments that may alter active-site distributions nor specialized equipment beyond reaction vessels and analytical instrumentation (e.g., NMR).

36 MATERIALS SCIENCE↗

DOE OSTI2020

Active Learning Accelerated Discovery of Stable Iridium Oxide Polymorphs for the Oxygen Evolution Reaction

The discovery of high-performing and stable materials for sustainable energy applications is a pressing goal in catalysis and materials science. Understanding the relationship between a material’s structure and functionality is an important step in the process, such that viable polymorphs for a given chemical composition need to be identified. Machine-learning-based surrogate models have the potential to accelerate the search for polymorphs that target specific applications. Herein, we report a readily generalizable active-learning (AL) accelerated algorithm for identification of electrochemically stable iridium oxide polymorphs of IrO 2 and IrO 3 . The search is coupled to a subsequent analysis of the electrochemical stability of the discovered structures for the acidic oxygen evolution reaction (OER). Structural candidates are generated by identifying all 956 structurally unique AB2 and AB3 prototypes in existing materials databases (more than 38000). Next, using an active learning approach, we find 196 IrO 2 polymorphs within the thermodynamic amorphous synthesizability limit and reaffirm the global stability of the rutile structure. We find 75 synthesizable IrO 3 polymorphs and report a previously unknown FeF 3 -type structure as the most stable, termed α-IrO 3 . To test the algorithms performance, we compare to a random search of the candidate space and report at least a 2-fold increase in the rate of discovery. Additionally, the AL approach can acquire the most stable polymorphs of IrO 2 and IrO 3 with fewer than 30 density functional theory optimizations. Analysis of the structural properties of the discovered polymorphs reveals that octahedral local coordination environments are preferred for nearly all low-energy structures. Subsequent Pourbaix Ir–H 2 O analysis shows that α-IrO3 is the globally stable solid phase under acidic OER conditions and supersedes the stability of rutile IrO 2 . Calculation of theoretical OER surface activities reveal ideal weaker binding of the OER intermediates on α-IrO 3 than on any other considered iridium oxide. We emphasize that the proposed AL algorithm can be easily generalized to search for any binary metal oxide structure with a defined stoichiometry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗