Metal-free hydrides are of increasing research interest due to their roles in recent scientific advances in catalysis, such as hydrogen activation with frustrated Lewis pairs and electrocatalytic CO2 reduction with pyridinium and other aromatic cations. The structural design of hydrides for specific applications necessitates the correct description of their thermodynamic and kinetic prowess using reliable parameters - thermodynamic hydricity (ΔGH-) and nucleophilicity (N). This review summarizes reported experimental and calculated hydricity values for more than 200 metal-free hydride donors, including carbon-, boron-, nitrogen- and silicon-based hydrides. We describe different experimental and computational methods used to obtain these thermodynamic and kinetic parameters. Furthermore, tabulated data on metal-free hydrides are discussed in terms of structure-property relationships, relevance to catalysis and contemporary limitations for replacing transition-metal hydride catalysts. Finally, several selected applications of metal-free hydrides in catalysis are described, including photosynthetic CO2 reduction and hydrogen activation with frustrated Lewis pairs.

We report a novel metal-free chemical reduction of CO2 by a recyclable benzimidazole-based organo-hydride, whose choice was guided by quantum chemical calculations. Notably, benzimidazole-based hydride donors rival the hydride-donating abilities of noble-metal-based hydrides such as [Ru­(tpy)­(bpy)­H]+ and [Pt­(depe)2H]+. Chemical CO2 reduction to the formate anion (HCOO–) was carried out in the absence of biological enzymes, a sacrificial Lewis acid, or a base to activate the substrate or reductant. 13CO2 experiments confirmed the formation of H13COO– by CO2 reduction with the formate product characterized by 1H NMR and 13C NMR spectroscopy and ESI-MS. The highest formate yield of 66% was obtained in the presence of potassium tetrafluoroborate under mild conditions. The likely role of exogenous salt additives in this reaction is to stabilize and shift the equilibrium toward the ionic products. After CO2 reduction, the benzimidazole-based hydride donor was quantitatively oxidized to its aromatic benzimidazolium cation, establishing its recyclability. In addition, we electrochemically reduced the benzimidazolium cation to its organo-hydride form in quantitative yield, demonstrating its potential for electrocatalytic CO2 reduction. These results serve as a proof of concept for the electrocatalytic reduction of CO2 by sustainable, recyclable, and metal-free organo-hydrides.

This communication reports a combined experimental and computational study of mechanisms by which biomimetic NADH analogs can be electrochemically regenerated.

Development of efficient electrocatalysts for the CO2 reduction reaction (CO2RR) to multicarbon products has been constrained by high overpotentials and poor selectivity. Here, we introduce iron phosphide (Fe2P) as an earth-abundant catalyst for the CO2RR to mainly C2–C4 products with a total CO2RR Faradaic efficiency of 53% at 0 V vs RHE. Carbon product selectivity is tuned in favor of ethylene glycol formation with increasing negative bias at the expense of C3–C4 products. Both Grand Canonical-DFT (GC-DFT) calculations and experiments reveal that *formate, not *CO, is the initial intermediate formed from surface phosphino-hydrides and that the latter form ionic hydrides at both surface phosphorus atoms (H@Ps) and P-reconstructed Fe3 hollow sites (H@P*). Binding of these surface hydrides weakens with negative bias (reactivity increases), which accounts for both the shift to C2 products over higher C–C coupling products and the increase in the H2 evolution reaction (HER) rate. GC-DFT predicts that phosphino-hydrides convert *formate to *formaldehyde, the key intermediate for C–C coupling, whereas hydrogen atoms on Fe generate tightly bound *CO via sequential PCET reactions to H2O. GC-DFT predicts the peak in CO2RR current density near −0.1 V is due to a local maximum in the binding affinity of *formate and *formaldehyde at this bias, which together with the more labile C2 product affinity, accounts for the shift to ethylene glycol and away from C3–C4 products. Consistent with these predictions, addition of exogenous CO is shown to block all carbon product formation and lower the HER rate. These results demonstrate that the formation of ionic hydrides and their binding affinity, as modulated by the applied potential, controls the carbon product distribution. This knowledge provides new insight into the influence of hydride speciation and applied bias on the chemical reaction mechanism of CO2RR that is relevant to all transition metal phosphides.

We report the results of modeling CO2 reduction (CO2R) to CO over Ag(110) and Cu(211) surfaces at different applied potentials using grand-canonical density functional theory (GC-DFT), a method specifically designed to accurately model electrochemical systems. In addition to demonstrating GC-DFT’s ability to accurately model electrochemical processes, we also compare it with the computational hydrogen electrode (CHE) approach. GC-DFT predicts that the geometries of these reacting systems strongly depend on the applied potential, and the Helmholtz free energies vary nonlinearly with the applied potential, which contradicts a central assumption of the CHE approach. The CHE approach neglects the change in the number of electrons on the electrode surface at different applied potentials, which reduces its accuracy as the potential changes from the potential of zero charge. Our results further demonstrate that the grand free energies of the reaction intermediates not only depend on the value of the applied potential but also on the metal surface type, adsorption site, and adsorbate. GC-DFT’s ability to predict the effect of the applied potential on adsorbate geometry enables it to evaluate different possible reaction mechanisms at different applied potentials. For instance, GC-DFT predicts that the first step of CO2R likely switches from proton-coupled electron transfer to sequential electron transfer and then proton transfer at more reducing potentials, a result that cannot be determined by the CHE because it assumes that all electron transfers are coupled to proton transfers and neglects the effect of the applied potential on the adsorbate geometry.

This work studies the electrocatalytic disproportionation of CO2 into CO and CO3 2– mediated by a [Mn-2,6-bis­(1-(alkyl)­imidazol-2-ylidene)­pyridine]Br pincer complex. We identify three mechanistic scenarios involving one or two catalytically active metal centers in the activation of carbon dioxide and use density functional theory to map out the energy landscape for each of the mechanistic steps. Experimentally determined second-order kinetics in CO2 consumption, the formation of CO and CO3 2– together with an observed order in catalyst of 0.5 determined by cyclic voltammetry using the Burés normalized timescale method, suggest a turnover-limiting CO3 2– loss through the scission of a bimetallic species. Faradaic efficiencies for the reduction of CO2 to CO of 86 ± 4% are observed. Given the reaction mass balance, we estimate a Faradaic efficiency from the conversion of CO2 to CO3 2– of 93 ± 4%. Our combined experimental−theoretical approach suggests that two sequential CO2 insertions are followed by a rearrangement to produce carbon monoxide and carbonate from two molecules of CO2 at singly reduced catalyst molecules. The disproportionation reaction reported herein combines the power of catalytic CO2 conversion with sequestration and provides a new chemical avenue for the conversion of carbon dioxide.

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.

Selective reduction of CO2 to formate represents an ongoing challenge in photoelectrocatalysis. To provide mechanistic insights, we investigate the kinetics of hydride transfer (HT) from a series of metal-free hydride donors to CO2. The observed dependence of experimental and calculated HT barriers on the thermodynamic driving force was modeled by using the Marcus hydride transfer formalism to obtain the insights into the effect of reorganization energies on the reaction kinetics. Our results indicate that even if the most ideal hydride donor were discovered, the HT to CO2 would exhibit sluggish kinetics (<100 turnovers per second at −0.1 eV driving force), indicating that the conventional HT may not be an appropriate mechanism for solar conversion of CO2 to formate. We propose that the conventional HT mechanism should not be considered for CO2 reduction catalysis and argue that the orthogonal HT mechanism, previously proposed to address thermodynamic limitations of this reaction, may also lead to lower kinetic barriers for CO2 reduction to formate.