Search

Ammonia is crucial as a fertilizer and in the chemical industry and is considered to be a carbon-free fuel1. Ammonia electrosynthesis from nitrogen under ambient conditions offers an attractive alternative to the Haber–Bosch process2,3, and lithium-mediated nitrogen reduction represents a promising approach to continuous-flow ammonia electrosynthesis, coupling nitrogen reduction with hydrogen oxidation4. However, tetrahydrofuran, which is commonly used as a solvent, impedes long-term ammonia production owing to polymerization and volatility problems. Here we show that a chain-ether-based electrolyte enables long-term continuous ammonia synthesis. We find that a chain-ether-based solvent exhibits non-polymerization properties and a high boiling point (162 °C) and forms a compact solid-electrolyte interphase layer on the gas diffusion electrode, facilitating ammonia release in the gas phase and ensuring electrolyte stability. We demonstrate 300 h of continuous operation in a flow electrolyser with a 25 cm2electrode at 1 bar pressure and room temperature, and achieve a current-to-ammonia efficiency of 64 ± 1% with a gas-phase ammonia content of approximately 98%. Our results highlight the crucial role of the solvent in long-term continuous ammonia synthesis.

Ammonia (NH3) is a key commodity chemical for the agricultural, textile and pharmaceutical industries, but its production via the Haber–Bosch process is carbon-intensive and centralized. Alternatively, an electrochemical method could enable decentralized, ambient NH3production that can be paired with renewable energy. The first verified electrochemical method for NH3synthesis was a process mediated by lithium (Li) in organic electrolytes. So far, however, elements other than Li remain unexplored in this process for potential benefits in efficiency, reaction rates, device design, abundance and stability. In our demonstration of a Li-free system, we found that calcium can mediate the reduction of nitrogen for NH3synthesis. We verified the calcium-mediated process using a rigorous protocol and achieved an NH3Faradaic efficiency of 40 ± 2% using calcium tetrakis(hexafluoroisopropyloxy)borate (Ca[B(hfip)4]2) as the electrolyte. Our results offer the possibility of using abundant materials for the electrochemical production of NH3, a critical chemical precursor and promising energy vector.

Noble metal-based catalysts have numerous industrial uses, and maximum utilization of the precious metals by lowering the metal loading is of significant interest in heterogeneous catalysis research. However, lowering the metal loading could lead to single-atom metal species formation, which may not be active for important reactions like propylene oxidation. We report a way to drastically reduce precious metal loading of catalysts by judiciously choosing an active metal/support pair and using the reversible atomization–nanoparticulate formation of transition metal on a high-surface area support. Here, Pt and MgAl2O4are used as the transition metal and high-surface area support, respectively. Through catalytic testing and characterization using scanning transmission electron microscopy and synchrotron X-ray absorption spectroscopy, a reversible change between atomization and nano-cluster formation under oxidizing and reducing conditions has been found. Via density functional theory, favorable sites for reversible Pt adsorption are identified, including ionic Pt4+sites that can serve to nucleate nanoclusters. Catalytic reaction modeling also rationalizes the catalytic inertness of atomic Pt sites. Finally, a re-activation mechanism for the atomized Pt based on gases present during reaction has been formulated and demonstrated.

Ammonia is a large-scale commodity essential to fertilizer production, but the Haber-Bosch process leads to massive emissions of carbon dioxide. Electrochemical ammonia synthesis is an attractive alternative pathway, but the process is still limited by low ammonia production rate and faradaic efficiency. Herein, guided by our theoretical model, we present a highly efficient lithium-mediated process enabled by using different lithium salts, leading to the formation of a uniform solid-electrolyte interphase (SEI) layer on a porous copper electrode. The uniform lithium-fluoride-enriched SEI layer provides an ammonia production rate of 2.5 ± 0.1 μmol s−1cmgeo−2at a current density of −1 A cmgeo−2with 71% ± 3% faradaic efficiency under 20 bar nitrogen. Experimental X-ray analysis reveals that the lithium tetrafluoroborate electrolyte induces the formation of a compact and uniform SEI layer, which facilitates homogeneous lithium plating, suppresses the undesired hydrogen evolution as well as electrolyte decomposition, and enhances the nitrogen reduction.

The development of inexpensive and abundant catalysts with high activity, selectivity, and stability for the oxygen reduction reaction (ORR) is imperative for the widespread implementation of fuel cell devices. Herein, we present a combined theoretical–experimental approach to discover and design first-row transition metal antimonates as excellent electrocatalytic materials for the ORR. Theoretically, we identify first-row transition metal antimonates─MSb2O6, where M = Mn, Fe, Co, and Ni─as nonprecious metal catalysts with good oxygen binding energetics, conductivity, thermodynamic phase stability, and aqueous stability. Among the considered antimonates, MnSb2O6shows the highest theoretical ORR activity based on the 4e–ORR kinetic volcano. Experimentally, nanoparticulate transition metal antimonate catalysts are found to have a minimum of a 2.5-fold enhancement in intrinsic mass activity (on transition metal mass basis) relative to the corresponding transition metal oxide at 0.7 V vs RHE in 0.1 M KOH. MnSb2O6is the most active catalyst under these conditions, with a 3.5-fold enhancement on a per Mn mass activity basis and 25-fold enhancement on a surface area basis over its antimony-free counterpart. Electrocatalytic and material stability are demonstrated over a 5 h chronopotentiometry experiment in the stability window identified by theoretical Pourbaix analysis. This study further highlights the stable and electrically conductive antimonate structure as a framework to tune the activity and selectivity of nonprecious metal oxide active sites for ORR catalysis.

Understanding hydrogen evolution reaction (HER) behaviors over two-dimensional transition-metal dichalcogenides (2D-TMDs) is critical for the development of nonprecious HER electrocatalysts with better activity. In this work, by combining density functional theory calculations with microkinetic modeling, we thoroughly investigated the HER mechanism on 2D-TMDs. We find an important dependence of simulated cell size on the calculated hydrogen adsorption energy and the activation barrier for MoS2. Distinct from previous “H migration” mechanisms proposed for the Heyrovsky reaction, the rate-determining step for MoS2, we propose that the Mo site only serves as the stabilized transition state rather than H adsorption. In comparison to transition-metal electrocatalysts, we find that the activation barrier of the Heyrovsky reaction on 2D-TMDs scales with the hydrogen adsorption energy exactly as for transition metals except that all activation energies are displaced upward by ca. 0.4 eV. This higher Heyrovsky activation barrier is responsible for the substantially lower activity of 2D-TMDs. We further show that this higher activation barrier stems from the more positively charged adsorbed hydrogen on the chalcogenides interacting repulsively with the incoming proton. Based on these insights, we discuss potential strategies for the design of nonprecious HER catalysts with activity comparable to Pt.

Interfaces between solvents and transition metals play a critical role in many catalytic and electrochemical systems. In this work, we use density functional theory (DFT) to study the interactions between several common nonaqueous solvents and a number of transition metal surfaces. We investigate trends in binding energies among the metals, and the influence of specific surface sites and facets. We examine the electrostatic field dependence of the most stable binding configurations, and shed light on experimentally observed field-dependent or potential-dependent adsorption structures. We demonstrate that the reorientation of DMSO on noble metal surfaces can occur with only a relatively small change in the surface dipole moment; this result is simultaneously consistent with experimental evidence for field-dependent reorientation and measurements of surprisingly low differential capacitance of DMSO on noble metals. These results contribute to understanding competitive adsorption effects in catalysis, and also have implications for fundamental surface properties such as the surface structure, solvated work function, and potential of zero charge.

The dissociation of molecular nitrogen in lithium is of interest for several promising technologies, such as the catalytic synthesis of ammonia in ambient or mild conditions. In this work we simulate nitrogen dissociation in the lithium BCC (110) surface at ambient and elevated temperatures using density functional theory (DFT) metadynamics simulations. The rate constants at temperatures of 300, 400, and 500 K are calculated by statistical analysis of the reaction time distributions from the accelerated simulations. This approach finds and estimates rate constants for transition pathways out of the initial state; the required input is the stable initial state and a reasonable choice of collective variable. A single collective variable is used in this case: the N–N distance. The results are robust to changes in metadynamics parameters, and the reaction time distributions follow the expected exponential distribution. We show that the metadynamics-derived rate constants are in agreement with results from the conventional harmonic approximation approach using a climbing image nudged elastic band (NEB) transition state search. The reaction barriers from metadynamics and the NEB/harmonic approximation agree to within 0.02–0.04 eV at all temperatures studied. This work demonstrates that the harmonic approximation provides an accurate description of the rate constants for nitrogen dissociation in lithium metal, even at temperatures near or above the melting point of lithium, lending credence to previous and future theoretical studies using this approximation. Moreover, this work demonstrates a step toward the automated exploration and discovery of reaction mechanisms and associated rate constants for elementary surface-catalyzed reactions using DFT-based metadynamics.

Electrolyte properties in general and electrolyte ions in particular have been shown to have significant effects on the kinetics of electrochemical CO2and CO reduction at metal electrodes, but these effects have not yet been fully understood. We investigate the effects of cations at the aqueous electrolyte–Cu(211) interface on adsorbate binding energies and the electrolyte structure using density functional theory (DFT). Charging the interface via explicit Na+has systematic effects on adsorbate–electrolyte interactions and conformations. We describe specific local adsorbate–ion interactions, including direct alkali ion–adsorbate coordination and hydrogen bonding via ion-coordinated water molecules. The relative importance of these specific interactions and purely electrostatic field–adsorbate interactions is investigated by comparing the DFT-calculated ion effects to those predicted by purely electrostatic models of the interface and the adsorbates. We find that the trend in ion effects among different adsorbates at constant surface charge density is well explained by a purely electrostatic interaction model. The binding energy of OCCHO is found to depend strongly on the surface charge density as well as the spatial distribution of ions at constant surface charge density. These effects are also explained by a purely electrostatic local electrostatic field–adsorbate model. This indicates that alkali ion effects can be mainly attributed to purely electrostatic field interactions and that the local field at ion-stabilized active sites can depend significantly on both the overall charge density and on the spatial distribution of ions at constant charge density. This work provides new insight on alkali ion effects on a variety of adsorbates relevant to the CO2and CO reduction reactions by describing specific local ion–adsorbate interactions, systematic changes in adsorbate–electrolyte interactions induced by explicit surface charge, and comparisons of the calculated effects on binding energies to simple electrostatic field–adsorbate models to explain trends among various adsorbates and at varying surface charge densities and supercell configurations.

Catalytic systems whose properties can be systematically tuned via changes in synthesis conditions are highly desirable for the next-generation catalyst design and optimization. Herein, we present a two-dimensional (2D) conductive metal–organic framework consisting of M-N4units (M = Ni, Cu) and a hexaaminobenzene (HAB) linker as a catalyst for the oxygen reduction reaction. By varying synthetic conditions, we prepared two Ni-HAB catalysts with different crystallinities, resulting in catalytic systems with different electric conductivities, electrochemical activity, and stability. We show that crystallinity has a positive impact on conductivity and demonstrate that this improved crystallinity/conductivity improves the catalytic performance of our model system. Additionally, density functional theory simulations were performed to probe the origin of M-HAB’s catalytic activity, and they suggest that M-HAB’s organic linker acts as the active site with the role of the metal being to modulate the linker sites’ binding strength.

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 IrO2and IrO3. 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 AB2and AB3prototypes in existing materials databases (more than 38000). Next, using an active learning approach, we find 196 IrO2polymorphs within the thermodynamic amorphous synthesizability limit and reaffirm the global stability of the rutile structure. We find 75 synthesizable IrO3polymorphs and report a previously unknown FeF3-type structure as the most stable, termed α-IrO3. 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 IrO2and IrO3with 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–H2O analysis shows that α-IrO3is the globally stable solid phase under acidic OER conditions and supersedes the stability of rutile IrO2. Calculation of theoretical OER surface activities reveal ideal weaker binding of the OER intermediates on α-IrO3than 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.

In this work, we develop a new microkinetic model for the oxygen reduction reaction (ORR) which incorporates field effects into the established computational hydrogen electrode model. We find that the field can significantly alter the binding energy of ORR adsorbates, particularly *OOH, *O2, and *H2O2. By considering these effects, we find that the model gives polarization curves and rotating ring disk electrode currents consistent with those found experimentally on the Pt(111), Au(111), and Au(100) electrodes, allowing us to predict onset potentials for both overall ORR activity and activity toward hydrogen peroxide as well as reproduce experimental Tafel slopes. In particular, we resolve the peculiar behavior observed on the Au(100) surface, where increasing pH leads to greatly enhanced 4e-ORR activity, without strongly affecting 2e-ORR activity. We then use these predictions to better understand activity across all ORR catalysts, showing that weak-binding catalysts are more influenced by pH because of the strong field effects on the *OOH adsorbate. Finally, we argue that considering field effects can expand the search for improved ORR catalysts by allowing us to deconvolute the activity dependence of catalysts on the standard hydrogen electrode and the reversible hydrogen electrode.

Molybdenum nitride (Mo–N) catalysts have shown promising activity and stability for the oxygen reduction reaction (ORR) in acid. However, the effect of oxygen (O) incorporation (from synthesis, catalysis, or exposure to air) on their activity remains elusive. Here, we use reactive sputtering to synthesize three compositions of thin-film catalysts and use extensive materials characterization to investigate the depth-dependent structure and incorporated O. We show that the as-deposited Mo–N films are highly oxidized both at the surface (>30% O) and in the bulk (3–21% O) and that the ORR performance is strongly correlated with the bulk structure and composition. Activity for 4e–ORR is highest for compositions with the highest N/O and N/Mo ratio. Furthermore, H2O2production for the films with moderate O content is comparable to or higher than the most H2O2-selective nonprecious metal catalysts in acidic electrolyte, on a moles per mass or surface area of catalyst basis. Density functional theory provides insight into the energetics of O incorporation and vacancy formation, and we hypothesize that activity trends with O/N ratios can be traced to the varying crystallite phases and their interactions with ORR adsorbates. This work demonstrates the prevalence and significance of O in metal nitride electrocatalysts and motivates further investigation into the role of O in other nonprecious metal materials.

We estimate the rate of electron transfer to CO2at the Au (211)|water interface during adsorption in an electrochemical environment under reducing potentials. On the basis of density functional theory calculations at the generalized gradient approximation and hybrid levels of theory, we find electron transfer to the adsorbed *CO2to be very facile. This high rate of transfer is estimated by the energy distribution of the adsorbate-induced density of states as well as from the interaction between diabatic states representing neutral and negatively charged CO2. Up to 0.62 electrons is transferred to CO2, and this charge adiabatically increases with the bending angle to a lower limit of 137°. We conclude that this rate of electron transfer is extremely fast compared to the time scale of the nuclear degrees of freedom, that is, the adsorption process.

The dissociative adsorption of carboxylic acids on oxide surfaces is important for understanding adsorbed carboxylates, which are important as intermediates in catalytic reactions, for the organo-functionalization of oxide surfaces, and in many other aspects of oxide surface chemistry. We present here the first direct experimental measurement of the heat of dissociative adsorption of any carboxylic acid on any single-crystal oxide surface. The enthalpy of the dissociative adsorption of formic acid, the simplest carboxylic acid, to produce adsorbed formate and hydrogen (as a surface hydroxyl) on a (2 × 2)-NiO(111) surface is measured by single crystal adsorption calorimetry. The differential heat of adsorption decreases with formic acid coverage from 202 to 99 kJ/mol at saturation (0.25 ML). The structure of the adsorbed products is clarified by density functional theory (DFT) calculations, which provide energies in reasonable agreement with the calorimetry. These calculations show that formic acid readily dissociates on both the oxygen and Ni terminations of the octapolar NiO(111) surfaces, donating its acid H to a surface lattice oxygen, while HCOO adsorbs preferentially with bridging-type geometry near the M-O3/O-M3sites. The calculated energetics at low coverages agrees well with experimental data, while larger differences are observed at high coverage (0.25 ML). The large decrease in experimental heat of adsorption with coverage can be brought into agreement with the DFT energies if we assume that both types of octapolar surface terminations (O- and Ni-) are present on the starting surface.

The oxygen reduction reaction (ORR) plays a crucial role in various energy devices such as proton-exchange membrane fuel cells (PEMFCs) and metal–air batteries. Owing to the scarcity of the current state-of-the-art Pt-based catalysts, cost-effective Pt-free materials such as transition metal nitrides and their derivatives have gained overwhelming interest as alternatives. In particular, cobalt nitride (CoN) has demonstrated a reasonably high ORR activity. However, the nature of its active phase still remains elusive. Here, we employ density functional theory calculations to study the surface reactivity of rocksalt (RS) and zincblend (ZB) cobalt nitride. The performances of the catalysts terminated by the facets of (100), (110), and (111) are studied for the ORR. We demonstrate that the cobalt nitride surface is highly susceptible to oxidation under ORR conditions. The as-formed oxide overlayer on the facets of CoNRS(100) and CoNZB(110) presents a significant promotional effect in reducing the ORR overpotential, thereby increasing the activity in comparison with those of the pure CoNs. The results of this work rationalize a number of experimental reports in the literature and disclose the nature of the active phase of cobalt nitrides for the ORR. Moreover, they offer guidelines for understanding the activity of other transition metal nitrides and designing efficient catalysts for future generation of PEMFCs.