How the Bioinspired Fe2Mo6S8 Chevrel Breaks Electrocatalytic Nitrogen Reduction Scaling Relations.

The nitrogen reduction reaction (NRR) is a renewable alternative to the energy- and CO2-intensive Haber-Bosch NH3 synthesis process but is severely limited by the low activity and selectivity of studied electrocatalysts. The Chevrel phase Fe2Mo6S8 has a surface Fe-S-Mo coordination environment that mimics the nitrogenase FeMo-cofactor and was recently shown to provide state-of-the-art activity and selectivity for NRR. Here, we elucidate the previously unknown NRR mechanism on Fe2Mo6S8 via grand-canonical density functional theory (GC-DFT) that realistically models solvated and biased surfaces. Fe sites of Fe2Mo6S8 selectively stabilize the key *NNH intermediate via a narrow band of free-atom-like surface d-states that selectively hybridize with p-states of *NNH, which results in Fe sites breaking NRR scaling relationships. These sharp d-states arise from an Fe-S bond dissociation during N2 adsorption that mimics the mechanism of the nitrogenase FeMo-cofactor. Furthermore, we developed a new GC-DFT-based approach for calculating transition states as a function of bias (GC-NEB) and applied it to produce a microkinetic model for NRR at Fe2Mo6S8 that predicts high activity and selectivity, in close agreement with experiments. Our results suggest new design principles that may identify effective NRR electrocatalysts that minimize the barriers for *N2 protonation and *NH3 desorption and that may be broadly applied to the rational discovery of stable, multinary electrocatalysts for other reactions where narrow bands of surface d-states can be tuned to selectively stabilize key reaction intermediates and guide selectivity toward a target product. Furthermore, our results highlight the importance of using GC-DFT and GC-NEB to accurately model electrocatalytic reactions.

X-ray absorption spectroscopy insights on the structure anisotropy and charge transfer in Chevrel Phase chalcogenides.

The electronic structure and local coordination of binary (Mo6T8) and ternary Chevrel Phases (MxMo6T8) are investigated for a range of metal intercalant and chalcogen compositions. We evaluate differences in the Mo L3-edge and K-edge X-ray absorption near edge structure across the suite of chalcogenides MxMo6T8 (M = Cu, Ni, x = 1-2, T = S, Se, Te), quantifying the effect of compositional and structural modification on electronic structure. Furthermore, we highlight the expansion, contraction, and anisotropy of Mo6 clusters within these Chevrel Phase frameworks through extended X-ray absorption fine structure analysis. Our results show that metal-to-cluster charge transfer upon intercalation is dominated by the chalcogen acceptors, evidenced by significant changes in their respective X-ray absorption spectra in comparison to relatively unaffected Mo cations. These results explain the effects of metal intercalation on the electronic and local structure of Chevrel Phases across various chalcogen compositions, and aid in rationalizing electron distribution within the structure.

Machine Learning Guided Synthesis of Multinary Chevrel Phase Chalcogenides.

The Chevrel phase (CP) is a class of molybdenum chalcogenides that exhibit compelling properties for next-generation battery materials, electrocatalysts, and other energy applications. Despite their promise, CPs are underexplored, with only ∼100 compounds synthesized to date due to the challenge of identifying synthesizable phases. We present an interpretable machine-learned descriptor (Hδ) that rapidly and accurately estimates decomposition enthalpy (ΔHd) to assess CP stability. To develop Hδ, we first used density functional theory to compute ΔHd for 438 CP compositions. We then generated >560 000 descriptors with the new machine learning method SIFT, which provides an easy-to-use approach for developing accurate and interpretable chemical models. From a set of >200 000 compositions, we identified 48 501 CPs that Hδ predicts are synthesizable based on the criterion that ΔHd < 65 meV/atom, which was obtained as a statistical boundary from 67 experimentally synthesized CPs. The set of candidate CPs includes 2307 CP tellurides, an underexplored CP subset with a predicted preference for channel site occupation by cation intercalants that is rare among CPs. We successfully synthesized five of five novel CP tellurides attempted from this set and confirmed their preference for channel site occupation. Our joint computational and experimental approach for developing and validating screening tools that enable the rapid identification of synthesizable materials within a sparse class is likely transferable to other materials families to accelerate their discovery.

Stabilizing Hydrogen Adsorption through Theory-Guided Chalcogen Substitution in Chevrel-Phase Mo6X8 (X=S, Se, Te) Electrocatalysts.

In this work, we implement a facile microwave-assisted synthesis method to yield three binary Chevrel-Phase chalcogenides (Mo6X8; X = S, Se, Te) and investigate the effect of increasing chalcogen electronegativity on hydrogen evolution catalytic activity. Density functional theory predictions indicate that increasing chalcogen electronegativity in these materials will yield a favorable electronic structure for proton reduction. This is confirmed experimentally via X-ray absorption spectroscopy as well as traditional electrochemical analysis. We have identified that increasing the electronegativity of X in Mo6X8 increases the hydrogen adsorption strength owing to a favorable shift in the p-band position as well as an increase in the Lewis basicity of the chalcogen, thereby improving hydrogen evolution reaction energetics. We find that Mo6S8 exhibits the highest hydrogen evolution activity of the Mo6X8 series of catalysts, requiring an overpotential of 321 mV to achieve a current density of 10 mA cm-2ECSA, a Tafel slope of 74 mV per decade, and an exchange current density of 6.01 × 10-4 mA cm-2ECSA. Agreement between theory and experiment in this work indicates that the compositionally tunable Chevrel-Phase chalcogenide family is a promising framework for which electronic structure can be predictably modified to improve catalytic small-molecule reduction reactivity.

Computational and Experimental Evaluation of Peroxide Oxidants for Amine-Peroxide Redox Polymerization.

Amine-peroxide redox polymerization (APRP) is the prevalent method for producing radical-based polymers in the many industrial and medical applications where light or heat activation is impractical. We recently developed a detailed description of the APRP initiation process through a combined computational and experimental effort to show that APRP proceeds through SN2 attack by the amine on the peroxide, followed by the rate-determining homolysis of the resulting intermediate. Using this new mechanistic understanding, a variety of peroxides were computationally predicted to initiate APRP with fast kinetics. In particular, the rate of APRP initiation can be improved by radical and anion stabilization through increased π-electron conjugation or by increasing the electrophilicity of the peroxy bond through the addition of electron-withdrawing groups. On the other hand, the addition of electron-donating groups lowered the initiation rate. These design principles enabled the computational prediction of several new peroxides that exhibited improved initiation rates over the commonly used benzoyl peroxide. For example, the addition of nitro groups (NO2) to the para positions of benzoyl peroxide resulted in a theoretical radical generation rate of 1.9 × 10-9 s-1, which is ~150 times faster than the 1.3 × 10-11 s-1 radical generation rate observed with unsubstituted benzoyl peroxide. These accelerated kinetics enabled the development of a redox-based direct-writing process that exploited the extremely rapid reactivity of an optimized redox pair with a custom inkjet printer, capable of printing custom shapes from polymerizing resins without heat or light. Furthermore, the application of more rapid APRP kinetics could enable the acceleration of existing industrial processes, make new industrial manufacturing methods possible, and improve APRP compatibility with biomedical applications through reduced initiator concentrations that still produce rapid polymerization rates.

Rational Design of Efficient Amine Reductant Initiators for Amine-Peroxide Redox Polymerization.

Amine-peroxide redox polymerization (APRP) has been highly prevalent in industrial and medical applications since the 1950s, yet the initiation mechanism of this radical polymerization process is poorly understood so that innovations in the field are largely empirically driven and incremental. Through a combination of computational prediction and experimental analysis, we elucidate the mechanism of this important redox reaction between amines and benzoyl peroxide for the ambient production of initiating radicals. Our calculations show that APRP proceeds through SN2 attack by the amine on the peroxide but that homolysis of the resulting intermediate is the rate-determining step. We demonstrate a correlation between the computationally predicted initiating rate and the experimentally measured polymerization rate with an R2 = 0.80. The new mechanistic understanding was then applied to computationally predict amine reductant initiators with faster initiating kinetics. This led to our discovery of N-(4-methoxyphenyl)pyrrolidine (MPP) as amine reductant, which we confirmed significantly outperforms current state-of-the-art tertiary aromatic amines by ∼20-fold, making it the most efficient amine-peroxide redox initiator to date. The application of amines with superior kinetics such as MPP in APRP could greatly accelerate existing industrial processes, facilitate new industrial manufacturing methods, and improve biocompatibility in biomedical applications conducted with reduced initiator concentrations yet higher overall efficiency.