Computational and Experimental Study of the Mechanism of Hydrogen Generation from Water by a Molecular Molybdenum-Oxo Electrocatalyst

How general is the effect of the bulkiness of organic ligands on the basicity of metal-organic catalysts? H2-evolving Mo oxides/sulphides as case studies

Tailoring the activity of an organometallic catalyst usually requires a targeted ligand design. Tuning the ligand bulkiness and tuning the electronic properties are popular approaches, which are somehow interdependent because substituents of different sizes within ligands can determine inter alia the occurrence of different degrees of inductive effects. Ligand basicity, in particular, turned out to be a key property for the modulation of protonation reactions occurring in vacuo at the metals in complexes bearing organophosphorus ligands; however, when the same reactions take place in a polar organic solvent, their energetics becomes dependent on the trade-off between ligand basicity and bulkiness, with the polarity of the solvent playing a key role in this regard [Bancroft et al., Inorg. Chem., 1986, 25, 3675; Rovaletti et al., J. Phys. Org. Chem., 2018, 31, e3748]. In the present contribution, we carried out molecular dynamics and density functional theory calculations on water-soluble Mo-based catalysts for proton reduction, in order to study the energetics of protonation reactions in complexes where the incipient proton binds a catalytically active ligand (i.e., an oxide or a disulphide). We considered complexes either soaked in water or in a vacuum, and featuring N-based ancillary ligands of different bulkiness (i.e. cages constituted either by pyridine or isoquinoline moieties). Our results show that the energetics of protonation events can be affected by ancillary ligand bulkiness even when the metal center does not play the role of the H⁺ acceptor. In vacuo, protonation at the O or S atom in the α position relative to the metal in complexes featuring the bulky isoquinoline-based ligand is more favored by around 10 kcal mol^-1 when compared to the case of the pyridine-based counterparts, a difference that is almost zero when the same reactions occur in water. Such an outcome is rationalized in light of the different electrostatic properties of complexes bearing ancillary ligands of different sizes. The overall picture from theory indicates that such effects of ligand bulkiness can be relevant for the design of green chemistry catalysts that undergo protonation steps in water solutions.

A combined QTAIM/IRI topological analysis of the effect of axial/equatorial positions of NH2 and CN substituents in the [(PY5Me2)MoO]+ complex

By means of the Interaction Region Indicator (IRI) and Quantum Theory of Atoms in Molecules (QTAIM), the influence exerted by NH₂ (amino) and CN (cyano) as electron donor and electron acceptor substituent groups, respectively, located at para-positions of axial and equatorial pyridine rings of derivatized complexes coming from the [(PY₅Me₂)MoO]⁺ complex during the hydrogen molecular release in the gas phase was analyzed. In any case, a H-H covalent bond is forming at the transition state, with a strengthening of the electron density of 5.5% when the substituent group involved is NH₂ at the para-position of the axial pyridine ring. However, there was no difference between NH₂ and CN when these substituent groups are located at the para-positions of the equatorial pyridine rings. The topological properties of electron densities from the QTAIM are not perturbed by the electron donor and electron acceptor nature of the substituents, even when these substituent groups are located at the axial or equatorial pyridine rings of the Mo-based complex.

Proton reduction by a bimetallic zinc selenolate electrocatalyst

The development of alternative energy sources is the utmost priority of developing society. Unlike many prior homogeneous electrocatalysts that rely on a change in the oxidation state of the metal center and/or electrochemically active ligand, here we report the synthesis and structural characterization of a bimetallic zinc selenolate complex consisting of a redox silent zinc metal ion and a tridentate ligand that catalyzes the reduction of protons into hydrogen gas electrochemically and displays one of the highest reported TOF for a homogeneous TM-metal free ligand centered HER catalyst, 509 s-1. The current-voltage analysis confirms the onset overpotential of 0.86 V vs. Ag/AgCl for the HER process. Constant potential electrolysis (CPE) has been carried out to study the bulk electrolysis of our developed protocol, which reveals that the bimetallic zinc selenolate catalyst is stable under cathodic as well as anodic potentials and generates hydrogen gas with a faradaic efficiency of 75%. Preliminary studies on the heterogeneous catalyst were conducted by depositing the bimetallic zinc selenolate catalyst on the electrode surface.

Aqueous Stability of Zirconium Clusters, Including the Zr(IV) Hexanuclear Hydrolysis Complex [Zr6O4(OH)4(H2O)24]12+, from Density Functional Theory

Framework materials constitute a broad family of solids that range from zeolites and metal-organic frameworks (MOFs) to coordination polymers. The synthesis of such network structures typically rely on precursor molecular building blocks. As an example, the UiO-66 MOF series is constructed of hexanuclear [Zr₆O₄(OH)₄(CO₂)₁₂] cluster nodes and linear carboxylate linkers. Unfortunately, these Zr MOF cluster nodes cannot currently be manufactured in a sustainable way, motivating a search for "green" alternative synthesis methods. Stabilizing the hexanuclear Zr(IV) cluster (i.e., the hexamer, {Zr₆¹²⁺}) without the use of organic ligation would enable the use of environmentally friendly solvents such as water. The Zr(IV) tetranuclear cluster (i.e., the tetramer, {Zr₄⁸⁺}) can be stabilized in solution with or without organic ligands, yet the hexamer has yet to be synthesized without supporting ligands. The reasons why certain zirconium clusters are favored in aqueous solution over others are not well understood. This study reports the relative thermodynamic instability of the hypothetical hexamer {Zr₆¹²⁺} compared to the ubiquitous {Zr₄⁸⁺} tetramer. Density functional theory calculations were performed to obtain the hydrolysis Gibbs free energy of these species and used to construct Zr Pourbaix diagrams that illustrate the effects of electrochemical potential, pH, and Zr(IV) concentration. It was found that the aqueous {Zr₆¹²⁺} hexamer is ∼17.8 kcal/mol less stable than the aqueous {Zr₄⁸⁺} tetramer at pH = 0, V = 0, and [Zr(IV)] = 1 M, which is an energy difference on the order of counterion interactions. Electronic structure analyses were used to explore trends in the highest occupied molecular orbital-lowest unoccupied molecular orbital gap, frontier molecular orbitals, and electrostatic potential distribution of these clusters. The evidence suggests that the aqueous {Zr₆¹²⁺} hexamer may be promoted with more strategic syntheses incorporating minimal ligands and counterions.

The density polarization reveals directions of electron displacements due to the substituent effect: Analysis performed on a metal-organic Mo-Oxo catalyst

Some Mo-oxo complexes bearing pyridine rings have the capability for dihydrogen production from water. However, energy barrier and overall energy vary depending on the effect exerted by several substituent groups located at different positions around one or more pyridine rings which are ligands of these compounds. Based on the Karunadasa and coworkers investigation where the para-position was experimentally tested in compounds derivatised from the 2,6-bis[1,1-bis(2-pyridil)ethyl]-pyridine oxo-molybdenum complex synthesized (Karunadasa et al., Nature, 2010, 464, 1329), we tested the combined effect of electron-withdrawing and electron-donating groups simulated as perturbations represented by point-charges. Then, we used the density polarization concept, δρ(r), a local reactivity descriptor corresponding to the partially integrated linear response function, χ(r, r') (a non-local reactivity descriptor), which is able to reveal different displacements of π-electrons on molecular structures. We perturbed the para-positions in the pentadentate ligand 2,6-bis[1,1-bis(2-pyridil)ethyl]-pyridine in the Mo-based complex by means of point-charges. They were located in three different configurations of the organic ligand (trans, geminal, and cis) which could help to explain energy barriers and overall energy of reactions catalyzed by this type of Mo-complexes. Our results indicate that the trans configuration of point-charges induces the most amount of fraction of electron shifted on the complex. A Mo-based complex bearing the same trans configuration for electron-withdrawing and electron-donating substituent groups (cyano and amino, respectively), leads to a kinetically more favorable H₂ release than the cis or geminal configuration of the substituent groups aforementioned.

Rationalizing Photo-Triggered Hydrogen Evolution Using Polypyridine Cobalt Complexes: Substituent Effects on Hexadentate Chelating Ligands

Four novel polypyridine cobalt(II) complexes were developed based on a hexadentate ligand scaffold bearing either electron-withdrawing (-CF₃ ) or electron-donating (-OCH₃ ) groups in different positions of the ligand. Experiments and theoretical calculations were combined to perform a systematic investigation of the effect of the ligand modification on the hydrogen evolution reaction. The results indicated that the position, rather than the type of substituent, was the dominating factor in promoting catalysis. The best performances were observed upon introduction of substituents on the pyridine moiety of the hexadentate ligand, which promoted the formation of the Co(II)H intermediate via intramolecular proton transfer reactions with low activation energy. Quantum yields of 11.3 and 10.1 %, maximum turnover frequencies of 86.1 and 76.6 min^-1 , and maximum turnover numbers of 5520 and 4043 were obtained, respectively, with a -OCH₃ and a -CF₃ substituent.

Electronic and geometric structure effects on one-electron oxidation of first-row transition metals in the same ligand framework

Developing new transition metal catalysts requires understanding of how both metal and ligand properties determine reactivity. Since metal complexes bearing ligands of the Py5 family (2,6-bis-[(2-pyridyl)methyl]pyridine) have been employed in many fields in the past 20 years, we set out here to understand their redox properties by studying a series of base metal ions (M = Mn, Fe, Co, and Ni) within the Py5OH (pyridine-2,6-diylbis[di-(pyridin-2-yl)methanol]) variant. Both reduced (M^II) and the one-electron oxidized (M^III) species were carefully characterized using a combination of X-ray crystallography, X-ray absorption spectroscopy, cyclic voltammetry, and density-functional theory calculations. The observed metal-ligand interactions and electrochemical properties do not always follow consistent trends along the periodic table. We demonstrate that this observation cannot be explained by only considering orbital and geometric relaxation, and that spin multiplicity changes needed to be included into the DFT calculations to reproduce and understand these trends. In addition, exchange reactions of the sixth ligand coordinated to the metal, were analysed. Finally, by including published data of the extensively characterised Py5OMe (pyridine-2,6-diylbis[di-(pyridin-2-yl)methoxymethane])complexes, the special characteristics of the less common Py5OH ligand were extracted. This comparison highlights the non-innocent effect of the distal OH functionalization on the geometry, and consequently on the electronic structure of the metal complexes. Together, this gives a complete analysis of metal and ligand degrees of freedom for these base metal complexes, while also providing general insights into how to control electrochemical processes of transition metal complexes.

Enthalpy-Controlled Insertion of a "Nonspectator" Tricoordinate Phosphorus Ligand into Group 10 Transition Metal-Carbon Bonds

Insertion of a tricoordinate phosphorus ligand into late metal-carbon bonds is reported. Metalation of a P^P-chelating ligand (L1), composed of a nontrigonal phosphorous (i.e., P(III)) triamide moiety, P(N(o-N(Ar)C6H4)2, tethered by a phenylene linker to a -PiPr2 anchor, with group 10 complexes L2M(Me)Cl (M = Ni, Pd) results in insertion of the nontrigonal phosphorus site into the metal-methyl bond. The stable methylmetallophosphorane compounds thus formed are characterized spectroscopically and crystallographically. Metalation of L1 with (cod)PtII(Me)(Cl) does not lead to a metallophosphorane but rather to the standard bisphosphine chelate (κ2-L1)Pt(Me)(Cl). These divergent reactivities within group 10 are rationalized by reference to periodic variation in M-C bond enthalpies.

KID Procedure Applied on the [(PY5Me2)MoO]+ Complex

The KID (Koopmans in DFT) protocol usually applies in organic molecules of the closed-shell type. We used the KID procedure on an open-shell Mo-based system for the first time to choose the most suitable density functional to compute global and local reactivity descriptors obtained from the conceptual density-functional theory (DFT). From a set of 18 density functionals, spread from the second until the fourth rung of Jacob's ladder: BLYP, BP86, B97-D, MN12-L, MN15-L, M06-L, M11-L, CAM-B3LYP, PBE0, B3LYP, N12-SX, M06-2X, MN15, MN12-SX, ωB97X-D, M11, LC-ωHPBE, and APFD, we concluded that CAM-B3LYP provides the best outcome, and in the second place, M06-2X. Because the vertical first ionization potential and vertical first electron affinity in the ground state (gs) are defined as follows I = E _gs(N - 1) - E _gs(N) and A = E _gs(N) - E _gs(N + 1), where E _gs(N - 1), E _gs(N), and E _gs(N + 1) correspond to energies of the system bearing N, N + 1, and N - 1 electrons, along with Koopmans' theorem (KT) given by I ≈ -ε_HOMO (ε_HOMO, highest occupied molecular orbital energy) and A ≈ -ε_LUMO (ε_LUMO, lowest unoccupied molecular orbital energy), the deviation from the KT was performed by the use of the index, such that J _I = |E _gs(N - 1) - E _gs(N) + ε_HOMO| and J _A = |E _gs(N) - E _gs(N + 1) + ε_LUMO|, which are absolute deviations from the perspective of I and A, respectively. Furthermore, the ε_SOMO (SOMO: singly-occupied molecular orbital energy) leads us to another index given by |ΔSL| = |ε_SOMO - ε_LUMO|. Therefore, J _HL and |ΔSL| are indexes defined to evaluate the quality of the KT when employed within the context of quantum chemical calculations based on DFT and not the Hartree-Fock theory. We propose the index that could be more suitable to choose the most proper density functional because the J _HL and |ΔSL| are independent indexes.

Artificial photosynthesis: opportunities and challenges of molecular catalysts

Molecular catalysis plays an essential role in both natural and artificial photosynthesis (AP). However, the field of molecular catalysis for AP has gradually declined in recent years because of doubt about the long-term stability of molecular-catalyst-based devices. This review summarizes the development history of molecular-catalyst-based AP, including the fundamentals of AP, molecular catalysts for water oxidation, proton reduction and CO2 reduction, and molecular-catalyst-based AP devices, and it provides an analysis of the advantages, challenges, and stability of molecular catalysts. With this review, we aim to highlight the following points: (i) an investigation on molecular catalysis is one of the most promising ways to obtain atom-efficient catalysts with outstanding intrinsic activities; (ii) effective heterogenization of molecular catalysts is currently the primary challenge for the application of molecular catalysis in AP devices; (iii) development of molecular catalysts is a promising way to solve the problems of catalysis involved in practical solar fuel production. In molecular-catalysis-based AP, much has been attained, but more challenges remain with regard to long-term stability and heterogenization techniques.

Kinetics and mechanisms of catalytic water oxidation

The kinetics and mechanisms of thermal and photochemical oxidation of water with homogeneous and heterogeneous catalysts, including conversion from homogeneous to heterogeneous catalysts in the course of water oxidation, are discussed in this review article. Molecular and homogeneous catalysts have the advantage to clarify the catalytic mechanisms by detecting active intermediates in catalytic water oxidation. On the other hand, heterogeneous nanoparticle catalysts have advantages for practical applications due to high catalytic activity, robustness and easier separation of catalysts by filtration as compared with molecular homogeneous precursors. Ligand oxidation of homogeneous catalysts sometimes results in the dissociation of ligands to form nanoparticles, which act as much more efficient catalysts for water oxidation. Since it is quite difficult to identify active intermediates on the heterogeneous catalyst surface, the mechanism of water oxidation has hardly been clarified under heterogeneous catalytic conditions. This review focuses on the kinetics and mechanisms of catalytic water oxidation with homogeneous catalysts, which may be converted to heterogeneous nanoparticle catalysts depending on various reaction conditions.

Design and photophysical studies of iridium(iii)–cobalt(iii) dyads and their application for dihydrogen photo-evolution

We report several new dyads constituted of cationic iridium(iii) photosensitizers and cobalt(iii) catalyst connected via free pendant pyridine on the photosensitizers.

Design and reactivity of pentapyridyl metal complexes for ammonia oxidation

Computational and experimental work shows that Mo pentapyridal complexes can oxidize ammonia in the presence of a chemical mediator and evolve N₂.

Scrutinizing the substituent effect on Mo-based electrocatalysts for molecular hydrogen release through axial–equatorial decomposition: a DFT study

The effect of electron-donating and -withdrawing groups is studied in terms of axial and equatorial locations at the p-pyridine rings in the [(PY₅Me₂)MoO]²⁺ for the H₂ release in aqueous media. It opens a way to modify its kinetics and thermodynamics.

Formation of Mono Oxo Molybdenum(IV) PNP Pincer Complexes: Interplay between Water and Molecular Oxygen

The synthesis of cationic mono oxo MoIV PNP pincer complexes of the type [Mo(PNPMe-iPr)(O)X]+ (X = I, Br) from [Mo(PNPMe-iPr)(CO)X2] is described. These compounds are coordinatively unsaturated and feature a strong Mo≡O triple bond. The formation of these complexes proceeds via cationic 14e intermediates [Mo(PNPMe-iPr)(CO)X]+ and requires both molecular oxygen and water. ESI MS measurements with 18O labeled water (H2 18O) and molecular oxygen (18O2) indicates that water plays a crucial role in the formation of the Mo≡O bond. A plausible mechanism based on DFT calculations is provided. The X-ray structure of [Mo(PNPMe-iPr)(O)I]SbF6 is presented.

Reaction mechanism of hydrogen evolution catalysed by Co and Fe complexes containing a tetra-dentate phosphine ligand - a DFT study

The reaction mechanism of the electro-catalytic proton reduction in neutral phosphate buffer enabled by mononuclear cobalt and iron complexes containing a tetra-dentate phosphine ligand (MP₄N₂, M = Fe, Co) has been elucidated by density functional calculations. The phosphate from the buffer was found to play a crucial role by coordinating to the metal and delivering a proton to the metal hydride in the H-H bond formation. For the more efficient cobalt catalyst, the starting species is a Co^II complex with a hydrogen phosphate and a water molecule ligated at the two vacant coordination sites. Two sequential proton-coupled electron transfer reductions lead to the formation of a Co^II-H intermediate with a dihydrogen phosphate ligand, and the reduction potentials for these two steps were calculated to be -0.58 V and -0.72 V, respectively. Subsequently, the H-H bond formation takes place via coupling of the Co^II-H and the proton from the dihydrogen phosphate ligand. The total barrier was calculated to be 18.2 kcal mol^-1 with an applied potential of -0.5 V, which can further decrease to only 11.2 kcal mol^-1 with an applied potential of -0.8 V. When the phosphate is displaced by a water molecule, the total barrier for the dihydrogen formation increases by 7.3 kcal mol^-1. For the iron catalyst, the overall mechanism is essentially the same; however, the first reduction (Fe^II/Fe^I, potential of -1.13 V) is likely the rate-limiting step. The calculated results are in good agreement with the experimental data, which showed an onset potential of -0.50 V for the cobalt complex and -1.03 V for the iron complex.

Electrochemical behaviors and relativistic DFT calculations to understand the terminal ligand influence on the [Re6(μ3-Q)8X6]4− clusters

A Born–Haber thermodynamic cycle was used to determine the redox potential in a series of rhenium(iii) clusters theoretical analysis at DFT level was considered to estimate the free energy of the reversible process ReIII6/ReIII5Re^IV.

A 3D visualization of the substituent effect : A brief analysis of two components of the operational formula of dual descriptor for open-shell systems

Six organometallic compounds coming from a basic Mo-based complex were analyzed from the perspective of the dual descriptor in order to detect subtle influences that a substituent group could exert on the reactive core at a long range. Since the aforementioned complexes are open-shell systems, the used operational formula for the dual descriptor is that one defined for those aforementioned systems, which was then compared with spin density. In addition, dual descriptor was decomposed into two terms, each of which was also applied on every molecular system. The obtained results indicated that components of dual descriptor could become more useful than the operational formula of dual descriptor because differences exerted by the substituents at the para position were better detected by components of dual descriptor rather than the dual descriptor by itself.

Synthesis, characterization and solution behavior of a systematic series of pentapyridyl-supported RuII complexes: comparison to bimetallic analogs

A series of RuII complexes stabilized with the pentapyridyl ligand Py5Me2 (Py5Me2 = 2,6-bis(1,1-bis(2-pyridyl)ethyl)pyridine) and with an axial X ligand (X = Cl-, H2O, N3-, MeCN) were prepared and characterized in the solid state and in non-aqueous solution. The cyclic voltammograms of these complexes in MeCN reflect a reversible substitution of the axial X ligand with MeCN. Irreversible ligand substitution of [(Py5Me2)RuN3]+ is also observed in propylene carbonate, but only at oxidizing potentials that decompose the azide ligand. The monometallic chloride and azide species are compared with analogous Ru2 metal-metal bonded complexes, which have been reported to undergo irreversible chloride dissociation upon reduction.

Proton reduction reaction catalyzed by homoleptic nickel bis-1,2-dithiolate complexes: Experimental and theoretical mechanistic investigations

A series of homoleptic monoanionic nickel dithiolene complexes [Ni(bdt)2](NBu4), [Ni(tdt)2](NBu4), and [Ni(mnt)2](NBu4) containing the ligands benzene-1,2-dithiolate (bdt2-), toluene-3,4-dithiolate (tdt2-) and maleonitriledithiolate (mnt2-), respectively, have been employed as electrocatalysts in the hydrogen evolution reaction with trifluoroacetic acid as proton source in acetonitrile. All complexes were active catalysts with TONs reaching 113, 158 and 6 for [Ni(bdt)2](NBu4), [Ni(tdt)2](NBu4), and [Ni(mnt)2](NBu4), respectively. Faradaic yield for hydrogen evolution reaction reaches 88 % for 2- , which also displays the minimal overpotential requirement value (467 mV) within the series. Two pathways for H2 evolution can be hypothesized that differ on on the sequence of protonation and reduction steps. DFT calculations are in agreement with experimental data and indicate that protonation at sulfur follows reduction to the dianion. Hydrogen evolves from the direduced-diprotonated form via a highly distorted nickel hydride intermediate. The effects of acid strength and concentration in the hydrogen-evolving mechanism are also discussed.

μ-Pyridine-bridged copper complex with robust proton-reducing ability

The interconversion of the binuclear copper complex [Cu(DQPD)]2 to mononuclear [Cu(DQPD)]⁺ has been studied and their catalytic behaviour towards proton reduction has been reported.

Hydrogenase Enzymes and Their Synthetic Models: The Role of Metal Hydrides

Hydrogenase enzymes efficiently process H2 and protons at organometallic FeFe, NiFe, or Fe active sites. Synthetic modeling of the many H2ase states has provided insight into H2ase structure and mechanism, as well as afforded catalysts for the H2 energy vector. Particularly important are hydride-bearing states, with synthetic hydride analogues now known for each hydrogenase class. These hydrides are typically prepared by protonation of low-valent cores. Examples of FeFe and NiFe hydrides derived from H2 have also been prepared. Such chemistry is more developed than mimicry of the redox-inactive monoFe enzyme, although functional models of the latter are now emerging. Advances in physical and theoretical characterization of H2ase enzymes and synthetic models have proven key to the study of hydrides in particular, and will guide modeling efforts toward more robust and active species optimized for practical applications.

The secondary coordination sphere and axial ligand effects on oxygen reduction reaction by iron porphyrins: a DFT computational study

Oxygen reduction reaction (ORR) catalyzed by a bio-inspired iron porphyrin bearing a hanging carboxylic acid group over the porphyrin ring, and a tethered axial imidazole ligand was studied by DFT calculations. BP86 free energy calculations of the redox potentials and pK a's of reaction components involved in the proton coupled electron transfer (PCET) reactions of the ferric-hydroxo and -superoxo complexes were performed based on Born-Haber thermodynamic cycle in conjunction with a continuum solvation model. The comparison was made with iron porphyrins that lack either in the hanging acid group or axial ligand, suggesting that H-bond interaction between the carboxylic acid and iron-bound hydroxo, aquo, superoxo, and peroxo ligands (de)stabilizes the Fe-O bonding, resulting in the increase in the reduction potential of the ferric complexes. The axial ligand interaction with the imidazole raises the affinity of the iron-bound superoxo and peroxo ligands for proton. In addition, a low-spin end-on ferric-hydroperoxo intermediate, a key precursor for O-O cleavage, can be stabilized in the presence of axial ligation. Thus, selective and efficient ORR of iron porphyrin can be achieved with the aid of the secondary coordination sphere and axial ligand interactions.

Uranium-mediated electrocatalytic dihydrogen production from water

Depleted uranium is a mildly radioactive waste product that is stockpiled worldwide. The chemical reactivity of uranium complexes is well documented, including the stoichiometric activation of small molecules of biological and industrial interest such as H2O, CO2, CO, or N2 (refs 1 - 11), but catalytic transformations with actinides remain underexplored in comparison to transition-metal catalysis. For reduction of water to H2, complexes of low-valent uranium show the highest potential, but are known to react violently and uncontrollably forming stable bridging oxo or uranyl species. As a result, only a few oxidations of uranium with water have been reported so far; all stoichiometric. Catalytic H2 production, however, requires the reductive recovery of the catalyst via a challenging cleavage of the uranium-bound oxygen-containing ligand. Here we report the electrocatalytic water reduction observed with a trisaryloxide U(III) complex [(((Ad,Me)ArO)3mes)U] (refs 18 and 19)--the first homogeneous uranium catalyst for H2 production from H2O. The catalytic cycle involves rare terminal U(IV)-OH and U(V)=O complexes, which have been isolated, characterized, and proven to be integral parts of the catalytic mechanism. The recognition of uranium compounds as potentially useful catalysts suggests new applications for such light actinides. The development of uranium-based catalysts provides new perspectives on nuclear waste management strategies, by suggesting that mildly radioactive depleted uranium--an abundant waste product of the nuclear power industry--could be a valuable resource.

Hydrogen production from water using a bis(imino)pyridine molybdenum electrocatalyst

Reduction of [(^Ph2PPrPDI)MoO][PF₆]₂ affords an unusual Mo(ii) oxo compound that mediates the electrocatalytic reduction of water.

Self-assembly and photocatalytic H2 evolution activity of two unprecedented polytantalotungstates based on the largest {Ta18} and {Ta18Yb2} clusters

Two unprecedented polytantalotungstates, 1 and 2 based on the largest {Ta₁₈} and {Ta₁₈Yb₂} clusters, respectively, were synthesized. 1 and 2 exhibit significant UV photocatalytic water splitting activity.

Proton-Coupled Electron Transfer: Moving Together and Charging Forward

Proton-coupled electron transfer (PCET) is ubiquitous throughout chemistry and biology. This Perspective discusses recent advances and current challenges in the field of PCET, with an emphasis on the role of theory and computation. The fundamental theoretical concepts are summarized, and expressions for rate constants and kinetic isotope effects are provided. Computational methods for calculating reduction potentials and pKa's for molecular electrocatalysts, as well as insights into linear correlations and non-innocent ligands, are also described. In addition, computational methods for simulating the nonadiabatic dynamics of photoexcited PCET are discussed. Representative applications to PCET in solution, proteins, electrochemistry, and photoinduced processes are presented, highlighting the interplay between theoretical and experimental studies. The current challenges and suggested future directions are outlined for each type of application, concluding with an overall view to the future.

Metal-polypyridyl catalysts for electro- and photochemical reduction of water to hydrogen

Climate change, rising global energy demand, and energy security concerns motivate research into alternative, sustainable energy sources. In principle, solar energy can meet the world's energy needs, but the intermittent nature of solar illumination means that it is temporally and spatially separated from its consumption. Developing systems that promote solar-to-fuel conversion, such as via reduction of protons to hydrogen, could bridge this production-consumption gap, but this effort requires invention of catalysts that are cheap, robust, and efficient and that use earth-abundant elements. In this context, catalysts that utilize water as both an earth-abundant, environmentally benign substrate and a solvent for proton reduction are highly desirable. This Account summarizes our studies of molecular metal-polypyridyl catalysts for electrochemical and photochemical reduction of protons to hydrogen. Inspired by concept transfer from biological and materials catalysts, these scaffolds are remarkably resistant to decomposition in water, with fast and selective electrocatalytic and photocatalytic conversions that are sustainable for several days. Their modular nature offers a broad range of opportunities for tuning reactivity by molecular design, including altering ancillary ligand electronics, denticity, and/or incorporating redox-active elements. Our first-generation complex, [(PY4)Co(CH3CN)2](2+), catalyzes the reduction of protons from a strong organic acid to hydrogen in 50% water. Subsequent investigations with the pentapyridyl ligand PY5Me2 furnished molybdenum and cobalt complexes capable of catalyzing the reduction of water in fully aqueous electrolyte with 100% Faradaic efficiency. Of particular note, the complex [(PY5Me2)MoO](2+) possesses extremely high activity and durability in neutral water, with turnover frequencies at least 8500 mol of H2 per mole of catalyst per hour and turnover numbers over 600 000 mol of H2 per mole of catalyst over 3 days at an overpotential of 1.0 V, without apparent loss in activity. Replacing the oxo moiety with a disulfide affords [(PY5Me2)MoS2](2+), which bears a molecular MoS2 triangle that structurally and functionally mimics bulk molybdenum disulfide, improving the catalytic activity for water reduction. In water buffered to pH 3, catalysis by [(PY5Me2)MoS2](2+) onsets at 400 mV of overpotential, whereas [(PY5Me2)MoO](2+) requires an additional 300 mV of driving force to operate at the same current density. Metalation of the PY5Me2 ligand with an appropriate Co(ii) source also furnishes electrocatalysts that are active in water. Importantly, the onset of catalysis by the [(PY5Me2)Co(H2O)](2+) series is anodically shifted by introducing electron-withdrawing functional groups on the ligand. With the [(bpy2PYMe)Co(CF3SO3)](1+) system, we showed that introducing a redox-active moiety can facilitate the electro- and photochemical reduction of protons from weak acids such as acetic acid or water. Using a high-throughput photochemical reactor, we examined the structure-reactivity relationship of a series of cobalt(ii) complexes. Taken together, these findings set the stage for the broader application of polypyridyl systems to catalysis under environmentally benign aqueous conditions.

Hydrogen evolution catalyzed by cobalt diimine-dioxime complexes

Mimicking photosynthesis and producing solar fuels is an appealing way to store the huge amount of renewable energy from the sun in a durable and sustainable way. Hydrogen production through water splitting has been set as a first-ranking target for artificial photosynthesis. Pursuing that goal requires the development of efficient and stable catalytic systems, only based on earth abundant elements, for the reduction of protons from water to molecular hydrogen. Cobalt complexes based on glyoxime ligands, called cobaloximes, emerged 10 years ago as a first generation of such catalysts. They are now widely utilized for the construction of photocatalytic systems for hydrogen evolution. In this Account, we describe our contribution to the development of a second generation of catalysts, cobalt diimine-dioxime complexes. While displaying similar catalytic activities as cobaloximes, these catalysts prove more stable against hydrolysis under strongly acidic conditions thanks to the tetradentate nature of the diimine-dioxime ligand. Importantly, H2 evolution proceeds via proton-coupled electron transfer steps involving the oxime bridge as a protonation site, reproducing the mechanism at play in the active sites of hydrogenase enzymes. This feature allows H2 to be evolved at modest overpotentials, that is, close to the thermodynamic equilibrium over a wide range of acid-base conditions in nonaqueous solutions. Derivatization of the diimine-dioxime ligand at the hydrocarbon chain linking the two imine functions enables the covalent grafting of the complex onto electrode surfaces in a more convenient manner than for the parent bis-bidentate cobaloximes. Accordingly, we attached diimine-dioxime cobalt catalysts onto carbon nanotubes and demonstrated the catalytic activity of the resulting molecular-based electrode for hydrogen evolution from aqueous acetate buffer. The stability of immobilized catalysts was found to be orders of magnitude higher than that of catalysts in the bulk. It led us to evidence that these cobalt complexes, as cobaloximes and other cobalt salts do, decompose under turnover conditions where they are free in solution. Of note, this process generates in aqueous phosphate buffer a nanoparticulate film consisting of metallic cobalt coated with a cobalt-oxo/hydroxo-phosphate layer in contact with the electrolyte. This novel material, H2-CoCat, mediates H2 evolution from neutral aqueous buffer at low overpotentials. Finally, the potential of diimine-dioxime cobalt complexes for light-driven H2 generation has been attested both in water/acetonitrile mixtures and in fully aqueous solutions. All together, these studies hold promise for the construction of molecular-based photoelectrodes for H2 evolution and further integration in dye-sensitized photoelectrochemical cells (DS-PECs) able to achieve overall water splitting.

Bioinspired Tungsten Dithiolene Catalysts for Hydrogen Evolution: A Combined Electrochemical, Photochemical, and Computational Study

Bis(dithiolene)tungsten complexes, W(VI)O2 (L = dithiolene)2 and W(IV)O (L = dithiolene)2, which mimic the active site of formate dehydrogenases, have been characterized by cyclic voltammetry and controlled potential electrolysis in acetonitrile. They are shown to be able to catalyze the electroreduction of protons into hydrogen in acidic organic media, with good Faradaic yields (75-95%) and good activity (rate constants of 100 s(-1)), with relatively high overpotentials (700 mV). They also catalyze proton reduction into hydrogen upon visible light irradiation, in combination with [Ru(bipyridine)3](2+) as a photosensitizer and ascorbic acid as a sacrificial electron donor. On the basis of detailed DFT calculations, a reaction mechanism is proposed in which the starting W(VI)O2 (L = dithiolene)2 complex acts as a precatalyst and hydrogen is further formed from a key reduced W-hydroxo-hydride intermediate.

Earth-abundant metal complexes as catalysts for water oxidation; is it homogeneous or heterogeneous?

This minireview focuses on the aspects that determine whether particular catalysts for the oxidation of water are homogeneous or heterogeneous.

A well-defined terminal vanadium(III) oxo complex

The ubiquity of vanadium oxo complexes in the V+ and IV+ oxidation states has contributed to a comprehensive understanding of their electronic structure and reactivity. However, despite being predicted to be stable by ligand-field theory, the isolation and characterization of a well-defined terminal mononuclear vanadium(III) oxo complex has remained elusive. We present the synthesis and characterization of a unique terminal mononuclear vanadium(III) oxo species supported by the pentadentate polypyridyl ligand 2,6-bis[1,1-bis(2-pyridyl)ethyl]pyridine (PY5Me2). Exposure of [V(II)(NCCH3)(PY5Me2)](2+) (1) to either dioxygen or selected O-atom-transfer reagents yields [V(IV)(O)(PY5Me2)](2+) (2). The metal-centered one-electron reduction of this vanadium(IV) oxo complex furnishes a stable, diamagnetic [V(III)(O)(PY5Me2)](+) (3) species. The vanadium(III) oxo species is unreactive toward H- and O-atom transfer but readily reacts with protons to form a putative vanadium hydroxo complex. Computational results predict that further one-electron reduction of the vanadium(III) oxo species will result in ligand-based reduction, even though pyridine is generally considered to be a poor π-accepting ligand. These results have implications for future efforts toward low-valent vanadyl chemistry, particularly with regard to the isolation and study of formal vanadium(II) oxo species.

Computational electrochemistry: prediction of liquid-phase reduction potentials

This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car-Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.

Proton-coupled electron transfer in molecular electrocatalysis: theoretical methods and design principles

Molecular electrocatalysts play an essential role in a wide range of energy conversion processes. The objective of electrocatalyst design is to maximize the turnover frequency and minimize the overpotential for the overall catalytic cycle. Typically, the catalytic cycle is dominated by key proton-coupled electron transfer (PCET) processes comprised of sequential or concerted electron and proton transfer steps. Theoretical methods have been developed to investigate the mechanisms, thermodynamics, and kinetics of PCET processes in electrocatalytic cycles. Electronic structure methods can be used to calculate the reduction potentials and pKa's and to generate thermodynamic schemes, free energy reaction pathways, and Pourbaix diagrams, which indicate the most stable species under certain conditions. These types of calculations have assisted in identifying the thermodynamically favorable mechanisms under specified experimental conditions, such as acid strength and overpotential. Such calculations have also revealed linear correlations among the thermodynamic properties, which can be used to predict the impact of modifying the ligands, substituents, or metal centers. The thermodynamic properties can be tuned with electron-withdrawing or electron-donating substituents. Ligand modification can exploit the role of noninnocent ligands. For example, ligand protonation typically decreases the overpotential. Calculations of rate constants for electron and proton transfer, as well as concerted PCET, have assisted in identifying the kinetically favorable mechanisms under specified conditions. The concerted PCET mechanism is thought to lower the overpotential required for catalysis by avoiding high-energy intermediates. Rate constant calculations have revealed that the concerted mechanism involving intramolecular proton transfer will be favored by designing more flexible ligands that facilitate the proton donor-acceptor motion while also maintaining a sufficiently short equilibrium proton donor-acceptor distance. Overall, theoretical methods have assisted in the interpretation of experimental data and the design of more effective molecular electrocatalysts.

Production of hydrogen by electrocatalysis: making the H–H bond by combining protons and hydrides

Electrocatalytic production of hydrogen by nickel complexes is reviewed, with an emphasis on heterocoupling of protons and hydrides.

An unsupported metal hydroxide for the design of molecular μ-oxo bridged heterobimetallic complexes

A terminal and unsupported chromium(iii) hydroxide is reported.