Proton-coupled electron transfer at SOFC electrodes.

Understanding the charge transfer processes at solid oxide fuel cell (SOFC) electrodes is critical to designing more efficient and robust materials. Activation losses at SOFC electrodes have been widely attributed to the ambipolar migration of charges at the mixed ionic-electronic conductor-gas interface. Empirical Butler-Volmer kinetics based on the transition state theory is often used to model the current-voltage relationship, where charged particles transfer classically over an energy barrier. However, the hydrogen oxidation/water electrolysis reaction H2(g) + O2- ⇌ H2O(g) + 2e- must be modeled through concerted electron and proton tunneling events, where we unify the theory of the electrostatic surface potential with proton-coupled electron transfer kinetics. We derive a framework for the reaction rate that depends on the electrostatic surface potential, adsorbate dipole moment, the electronic structure of the electron donor/acceptor, and vibronic states of the hydrogen species. This theory was used to study the current-voltage characteristics of the Ni/gadolinium-doped ceria electrode in H2/H2O(g), where we find excellent validation of this novel model. These results yield the first reported quantification of the solvent reorganization energy for an SOFC material and suggest that the three-phase boundary mechanism is the dominant pathway for charge transfer at cermet electrodes.

Inductive Effect Alone Cannot Explain Lewis Adduct Formation and Dissociation at Electrode Interfaces.

Understanding breaking and formation of Lewis bonds at an electrified interface is relevant to a large range of phenomena, including electrocatalysis and electroadsorption. The complexities of interfacial environments and associated reactions often impede a systematic understanding of this type of bond at interfaces. To address this challenge, we report the creation of a main group classic Lewis acid-base adduct on an electrode surface and its behavior under varying electrode potentials. The Lewis base is a self-assembled monolayer of mercaptopyridine and the Lewis acid is BF3, forming a Lewis bond between nitrogen and boron. The bond is stable at positive potentials but cleaves at potentials more negative of approximately -0.3 V vs Ag/AgCl without an associated current. We also show that if the Lewis acid BF3 is supplied from a reservoir of Li+BF4- electrolyte, the cleavage is completely reversible. We propose that the N-B Lewis bond is affected both by the field-induced intramolecular polarization (electroinduction) and by the ionic structures and ionic equilibria near the electrode. Our results indicate that the second effect is responsible for the Lewis bond cleavage at negative potentials. This work is relevant to understanding the fundamentals of electrocatalytic and electroadsorption processes.

Correlation between Electronic Descriptor and Proton-Coupled Electron Transfer Thermodynamics in Doped Graphite-Conjugated Catalysts.

Graphite-conjugated catalysts (GCCs) provide a powerful framework for investigating correlations between electronic structure features and chemical reactivity of single-site heterogeneous catalysts. GCC-phenazine undergoes proton-coupled electron transfer (PCET) involving protonation of phenazine at its two nitrogen atoms with the addition of two electrons. Herein, this PCET reaction is investigated in the presence of defects, such as heteroatom dopants, in the graphitic surface. The proton-coupled redox potentials, EPCET, are computed using a constant potential periodic density functional theory (DFT) strategy. The electronic states directly involved in PCET for GCC-phenazine exhibit the same nitrogen orbital character as those for molecular phenazine. The energy εLUS of this phenazine-related lowest unoccupied electronic state in GCC-phenazine is identified as a descriptor for changes in PCET thermodynamics. Importantly, εLUS is obtained from only a single DFT calculation but can predict EPCET, which requires many such calculations. Similar electronic features may be useful descriptors for thermodynamic properties of other single-site catalysts.

Multilevel Computational Studies Reveal the Importance of Axial Ligand for Oxygen Reduction Reaction on Fe-N-C Materials.

The systematic improvement of Fe-N-C materials for fuel cell applications has proven challenging, due in part to an incomplete atomistic understanding of the oxygen reduction reaction (ORR) under electrochemical conditions. Herein, a multilevel computational approach, which combines ab initio molecular dynamics simulations and constant potential density functional theory calculations, is used to assess proton-coupled electron transfer (PCET) processes and adsorption thermodynamics of key ORR intermediates. These calculations indicate that the potential-limiting step for ORR on Fe-N-C materials is the formation of the FeIII-OOH intermediate. They also show that an active site model with a water molecule axially ligated to the iron center throughout the catalytic cycle produces results that are consistent with the experimental measurements. In particular, reliable prediction of the ORR onset potential and the Fe(III/II) redox potential associated with the conversion of FeIII-OH to FeII and desorbed H2O requires an axial H2O co-adsorbed to the iron center. The observation of a five-coordinate rather than four-coordinate active site has significant implications for the thermodynamics and mechanism of ORR. These findings highlight the importance of solvent-substrate interactions and surface charge effects for understanding the PCET reaction mechanisms and transition-metal redox couples under realistic electrochemical conditions.

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer.

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. A vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential-dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

Understanding the Solid-State Electrode-Electrolyte Interface of a Model System Using First-Principles Statistical Mechanics and Thin-Film X-ray Characterization.

Intermixing of atomic species at the electrode-electrolyte boundaries can impact the properties of the interfaces in solid-state batteries. Herein, this work uses first-principles statistical mechanics along with experimental characterization to understand intermixing at the electrode-electrolyte interface. For the model presented in this work, lithium manganese oxide (LiMn2O4, LMO) and lithium lanthanum titanate (Li3xLa2/3-xTiO3, LLTO) are employed as the cathode and electrolyte, respectively. The results of the computational work show that Ti-Mn intermixing at the interface is significant at synthesis temperatures. The experimental results in this work find that, at some critical temperatures between 600 and 700 °C for material preparation, the interface of LLTO-LMO becomes blurred. Calculations predict that the interface is unstable with regard to Ti-Mn intermixing starting at 0 K, suggesting that the critical temperature found in the experiment is related to kinetics. The work overall suggests that, in designing a solid-state battery, the fundamental reactions such as intermixing need to be considered.

Proton-Coupled Defects Impact O-H Bond Dissociation Free Energies on Metal Oxide Surfaces.

Proton-coupled electron transfer (PCET) reactions on metal oxides require coupling between proton transfer at the solid-liquid interface and electron transfer involving defects at or near the band edge. Herein, hybrid functional periodic density functional theory is used to elucidate the impact of proton-coupled defects on the bond dissociation free energies (BDFEs) of O-H bonds on anatase TiO2 surfaces. These O-H BDFEs are directly related to interfacial PCET thermochemistry. Comparison between geometrically similar O-H bonds associated with different defect types, namely conduction d-band electrons or valence p-band holes, reveals that the BDFEs differ by ∼81 kcal/mol (3.50 eV), comparable to the wide TiO2 band gap. These differences are shown to be determined primarily by differences in electron transfer driving forces, which are analyzed by using band energies and inner-sphere reorganization energies within a Marcus theory framework. These fundamental insights about the impact of proton-coupled defects on PCET thermochemistry at semiconductor surfaces have broad implications for electrocatalysis.

Mechanistic Insights about Electrochemical Proton-Coupled Electron Transfer Derived from a Vibrational Probe.

Proton-coupled electron transfer (PCET) is a fundamental step in a wide range of electrochemical processes, including those of interest in energy conversion and storage. Despite its importance, several mechanistic details of such reactions remain unclear. Here, we have combined a proton donor (tertiary ammonium) with a vibrational Stark-shift probe (benzonitrile), to track the process from the entry of the reactants into the electrical double layer (EDL), to the PCET reaction associated with proton donation to the electrode, and the formation of products. We have used operando vibrational spectroscopy and periodic density functional theory under electrochemical bias to assign the reactant and product peaks and their Stark shifts. We have identified three main stages for the progress of the PCET reaction as a function of applied potential. First, we have determined the potential necessary for desolvation of the reactants and their entry into the polarizing environment of the EDL. Second, we have observed the appearance of product peaks prior to the onset of steady state electrochemical current, indicating formation of a stationary population of products that does not turn over. Finally, more negative of the onset potential, the electrode attracts additional reactants, displacing the stationary products and enabling steady state current. This work shows that the integration of a vibrational Stark-shift probe with a proton donor provides critical insight into the interplay between interfacial electrostatics and heterogeneous chemical reactions. Such insights cannot be obtained from electrochemical measurements alone.

Interfacial Field-Driven Proton-Coupled Electron Transfer at Graphite-Conjugated Organic Acids.

Interfacial proton-coupled electron transfer (PCET) reactions are central to the operation of a wide array of energy conversion technologies, but molecular-level insights into interfacial PCET are limited. At carbon surfaces, designer sites for interfacial PCET can be incorporated by conjugating organic acid functional groups to graphite edges though aromatic phenazine linkages. At these graphite-conjugated catalysts (GCCs) bearing organic acid moieties, PCET is driven by complex interfacial electrostatic and field gradients that are difficult to probe experimentally. Herein, the spatially inhomogeneous interfacial electrostatic potentials and electric fields of GCC organic acids are computed as functions of applied potential. The calculated proton-coupled redox potentials for the PCET reactions at the GCC phenazine bridges and organic acid sites are in agreement with cyclic voltammetry measurements for a series of GCC acids. The trends in these redox potentials are explained in terms of the acidity of the molecular analogues and continuous conjugation between the acid and the graphite surface. The calculations illustrate that this conjugation is interrupted in a GCC acetic acid system, providing an explanation for the absence of a cyclic voltammetry peak corresponding to PCET at this acid site. This combined theoretical and experimental study demonstrates the critical role of continuous conjugation and strong electronic coupling between the GCC acid site and the graphite to enable interfacial field-driven PCET at the acid site. Understanding the connection between the atomic structure of the surface and the interfacial electrostatic potentials and fields that govern PCET thermochemistry may guide heterogeneous catalyst design.

Oriented LiMn2O4 Particle Fracture from Delithiation-Driven Surface Stress.

The insertion and removal of Li+ ions into Li-ion battery electrodes can lead to severe mechanical fatigue because of the repeated expansion and compression of the host lattice during electrochemical cycling. In particular, the lithium manganese oxide spinel (LiMn2O4, LMO) experiences a significant surface stress contribution to electrode chemomechanics upon delithiation that is asynchronous with the potentials where bulk phase transitions occur. In this work, we probe the stress evolution and resulting mechanical fracture from LMO delithation using an integrated approach consisting of cyclic voltammetry, electron microscopy, and density functional theory (DFT) calculations. High-rate electrochemical cycling is used to exacerbate the mechanical deficiencies of the LMO electrode and demonstrates that mechanical degradation leads to slowing of delithiation and lithiation kinetics. These observations are further supported through the identification of significant fracturing in LMO using scanning electron microscopy. DFT calculations are used to model the mechanical response of LMO surfaces to electrochemical delithiation and suggest that particle fracture is unlikely in the [001] direction because of tensile stresses from delithiation near the (001) surface. Transmission electron microscopy and electron backscatter diffraction of the as-cycled LMO particles further support the computational analyses, indicating that particle fracture instead tends to preferentially occur along the {111} planes. This joint computational and experimental analysis provides molecular-level details of the chemomechanical response of the LMO electrode to electrochemical delithiation and how surface stresses may lead to particle fracture in Li-ion battery electrodes.

A Long-Cycle-Life Lithium-CO2 Battery with Carbon Neutrality.

Lithium-CO2 batteries are attractive energy-storage systems for fulfilling the demand of future large-scale applications such as electric vehicles due to their high specific energy density. However, a major challenge with Li-CO2 batteries is to attain reversible formation and decomposition of the Li2 CO3 and carbon discharge products. A fully reversible Li-CO2 battery is developed with overall carbon neutrality using MoS2 nanoflakes as a cathode catalyst combined with an ionic liquid/dimethyl sulfoxide electrolyte. This combination of materials produces a multicomponent composite (Li2 CO3 /C) product. The battery shows a superior long cycle life of 500 for a fixed 500 mAh g-1 capacity per cycle, far exceeding the best cycling stability reported in Li-CO2 batteries. The long cycle life demonstrates that chemical transformations, making and breaking covalent CO bonds can be used in energy-storage systems. Theoretical calculations are used to deduce a mechanism for the reversible discharge/charge processes and explain how the carbon interface with Li2 CO3 provides the electronic conduction needed for the oxidation of Li2 CO3 and carbon to generate the CO2 on charge. This achievement paves the way for the use of CO2 in advanced energy-storage systems.

Understanding the Role of Overpotentials in Lithium Ion Conversion Reactions: Visualizing the Interface.

Oxide conversion reactions are known to have substantially higher specific capacities than intercalation materials used in Li-ion batteries, but universally suffer from large overpotentials associated with the formation of interfaces between the resulting nanoscale metal and Li2O products. Here we use the interfacial sensitivity of operando X-ray reflectivity to visualize the structural evolution of ultrathin NiO electrodes and their interfaces during conversion. We observe two additional reactions prior to the well-known bulk, three-dimensional conversion occurring at 0.6 V: an accumulation of lithium at the buried metal/oxide interface (at 2.2 V) followed by interfacial lithiation of the buried NiO/Ni interface at the theoretical potential for conversion (at 1.9 V). To understand the mechanisms for bulk and interfacial lithiation, we calculate interfacial energies using density functional theory to build a potential-dependent nucleation model for conversion. These calculations show that the additional space charge layer of lithium is a crucial component for reducing energy barriers for conversion in NiO.

New Class of Electrocatalysts Based on 2D Transition Metal Dichalcogenides in Ionic Liquid.

The optimization of traditional electrocatalysts has reached a point where progress is impeded by fundamental physical factors including inherent scaling relations among thermokinetic characteristics of different elementary reaction steps, non-Nernstian behavior, and electronic structure of the catalyst. This indicates that the currently utilized classes of electrocatalysts may not be adequate for future needs. This study reports on synthesis and characterization of a new class of materials based on 2D transition metal dichalcogenides including sulfides, selenides, and tellurides of group V and VI transition metals that exhibit excellent catalytic performance for both oxygen reduction and evolution reactions in an aprotic medium with Li salts. The reaction rates are much higher for these materials than previously reported catalysts for these reactions. The reasons for the high activity are found to be the metal edges with adiabatic electron transfer capability and a cocatalyst effect involving an ionic-liquid electrolyte. These new materials are expected to have high activity for other core electrocatalytic reactions and open the way for advances in energy storage and catalysis.

Thermodynamic Stability of Low- and High-Index Spinel LiMn2O4 Surface Terminations.

Density functional theory calculations are performed within the generalized gradient approximation (GGA+U) to determine stable terminations of both low- and high-index spinel LiMn2O4 (LMO) surfaces. A grand canonical thermodynamic approach is employed, permitting a direct comparison of off-stoichiometric surfaces with previously reported stoichiometric surface terminations at various environmental conditions. Within this formalism, we have identified trends in the structure of the low-index surfaces as a function of the Li and O chemical potentials. The results suggest that, under a range of chemical potentials for which bulk LMO is stable, Li/O and Li-rich (111) surface terminations are favored, neither of which adopts an inverse spinel structure in the subsurface region. This thermodynamic analysis is extended to identify stable structures for certain high-index surfaces, including (311), (331), (511), and (531), which constitute simple models for steps or defects that may be present on real LMO particles. The low- and high-index results are combined to determine the relative stability of each surface facet under a range of environmental conditions. The relative surface energies are further employed to predict LMO particle shapes through a Wulff construction approach, which suggests that LMO particles will adopt either an octahedron or a truncated octahedron shape at conditions in which LMO is thermodynamically stable. These results are in agreement with the experimental observations of LMO particle shapes.