Computational Characterization of Redox Non-Innocence in Cobalt-Bis(Diaryldithiolene)-Catalyzed Proton Reduction

Precise Fingerprint Determination of Vibrational Infrared Spectra in a Series of Co(II) Clathrochelates through Experimental and Theoretical Analyses

The energetic demands of modern society for clean energy vectors, such as H2, have caused a surge in research associated with homogeneous and immobilized electrocatalysts that may replace Pt. In particular, clathrochelates have shown excellent electrocatalytic properties for the hydrogen evolution reaction (HER). However, the actual mechanism for the HER catalyzed by these d-metal complexes remains an open debate, which may be addressed via Operando spectroelectrochemistry. The prediction of electrochemical properties via density functional theory (DFT) needs access to thermodynamic functions, which are only available after Hessian calculations. Unfortunately, there is a notable lack in the current literature regarding the precise evaluation of vibrational spectra of such complexes, given their structural complexity and the associated tangled IR spectra. In this work, we have performed a detailed theoretical and experimental analysis in a family of Co(II) clathrochelates, in order to establish univocally their IR pattern, and also the calculation methodology that is adequate for such predictions. In summary, we have observed the presence of multiple common bands shared by this clathrochelate family, using the B3LYP functional, the LANL2DZ basis, and effective core potentials (ECP) for heavy atoms. The most important issue addressed in this article was therefore related to the detailed assignment of the fingerprint associated with cobalt(II) clathrochelates, which is a challenging endeavor due to the crowded nature of their spectra.

Oxidation State Localized Orbitals: A Method for Assigning Oxidation States Using Optimally Fragment-Localized Orbitals and a Fragment Orbital Localization Index

Oxidation states represent the ionic distribution of charge in a molecule and are significant in tracking redox reactions and understanding chemical bonding. While effective algorithms already exist based on formal Lewis structures as well as using localized orbitals, they exhibit differences in challenging cases where effects such as redox noninnocence are at play. Given a density functional theory (DFT) calculation with chosen total charge and spin multiplicity, this work reports a new approach to obtaining fragment-localized orbitals that is termed oxidation state localized orbitals (OSLO), together with an algorithm for assigning the oxidation state using the OSLOs and an associated fragment orbital localization index (FOLI). Evaluating the FOLI requires fragment populations, and for this purpose a new version of the intrinsic atomic orbital (IAO) scheme is introduced in which the IAOs are evaluated using a reference minimal basis formed from on-the-fly superposition of atomic density (IAO-AutoSAD) calculations in the target basis set and at the target level of theory. The OSLO algorithm is applied to a range of challenging cases including high valent metal oxide complexes, redox noninnocent NO and dithiolate transition metal complexes, a range of carbene-containing TM complexes, and other examples including the potentially inverted ligand field in [Cu(CF₃)₄]^-. Across this range of cases, OSLO produces generally satisfactory results. Furthermore, in borderline cases, the OSLOs and associated FOLI values provide direct evidence of the emergence of covalent interactions between fragments that nicely complements existing approaches.

H2 Evolution from Electrocatalysts with Redox-Active Ligands: Mechanistic Insights from Theory and Experiment vis-à-vis Co-Mabiq

Electrocatalytic hydrogen production via transition metal complexes offers a promising approach for chemical energy storage. Optimal platforms to effectively control the proton and electron transfer steps en route to H₂ evolution still need to be established, and redox-active ligands could play an important role in this context. In this study, we explore the role of the redox-active Mabiq (Mabiq = 2-4:6-8-bis(3,3,4,4-tetramethlyldihydropyrrolo)-10-15-(2,2-biquinazolino)-[15]-1,3,5,8,10,14-hexaene1,3,7,9,11,14-N₆) ligand in the hydrogen evolution reaction (HER). Using spectro-electrochemical studies in conjunction with quantum chemical calculations, we identified two precatalytic intermediates formed upon the addition of two electrons and one proton to [Co^II(Mabiq)(THF)](PF₆) (Co_Mbq). We further examined the acid strength effect on the generation of the intermediates. The generation of the first intermediate, Co_Mbq-H¹, involves proton addition to the bridging imine-nitrogen atom of the ligand and requires strong proton activity. The second intermediate, Co_Mbq-H², acquires a proton at the diketiminate carbon for which a weaker proton activity is sufficient. We propose two decoupled H₂ evolution pathways based on these two intermediates, which operate at different overpotentials. Our results show how the various protonation sites of the redox-active Mabiq ligand affect the energies and activities of HER intermediates.

Priority of Mixed Diamine Ligands in Cobalt Dithiolene Complex-Catalyzed H2 Evolution: A Theoretical Study

Redox non-innocent metal dithiolene or diamine complexes are potential alternative catalysts in hydrogen evolution reaction and have been incorporated into 2D metal-organic frameworks to obtain unexpected electrocatalytic activity. According to an experimental study, Co-bis(dithiolene), Co-bis(diamine), and Co-dithiolene-diamine portions are considered as active sites where the generation of H₂ occurs and a diamine ligand is necessary for high catalytic efficiency. We are interested in the difference between these catalytic active sites, and mechanistic studies on extracted Co-bis(dithiolene), Co-bis(diamine), and Co-dithiolene-diamine complex-catalyzed hydrogen evolution reactions are carried out by using density functional methods. Our calculated results indicate that the priority of ligand mixed complexes resulted from the readily occurring protonation of diamine ligands and large electron affinity of dithiolene ligands as well as the lowest overall barrier for H₂ evolution.

Improving and Understanding the Hydrogen Evolving Activity of a Cobalt Dithiolene Metal-Organic Framework

Despite the promising previous reports on the development of electrocatalytic dithiolene-based metal-organic frameworks (MOFs) for the hydrogen evolution reaction (HER), these materials often display poor reproducibility of the HER performance because of their poor bulk properties upon integration with electrode materials. We demonstrate here an in-depth investigation of the electrocatalytic HER activity of a cobalt 2,3,6,7,10,11-triphenylenehexathiolate (CoTHT) MOF. To enhance the durability and charge transport properties of the constructed CoTHT/electrode architecture, CoTHT is deposited as an ink composite (1) composed of Nafion and carbon black. We leverage here the well-established use of catalyst inks in the literature to increase adhesion of the catalyst to the electrode surface and to improve the overall electrical conductivity of the integrated catalyst/electrode. The utilization of the composite 1 leads to a significant improvement in the overpotential (η) to reach a current density of 10 mA/cm² (η = 143 mV) compared to prior reports, resulting in the most active MOF-based electrocatalyst for the HER that contains only earth-abundant elements. Extensive density functional theory (DFT) calculations were applied to understand the structure of CoTHT and the mechanistic pathways of the HER. The computational results suggest that an AB stacking geometry is energetically favorable, where one layer is slipped by 1.6 Å relative to the neighboring one along the a and b vectors. Additionally, the DFT calculations indicate that the catalytic cycle likely involves a Volmer discharge step to generate a cobalt hydride, followed by a Heyrovsky step to form a cobalt-H₂ intermediate, and finally the dissociation of H₂.

Computational Modeling Predicts the Stability of Both Pd+ and Pd2+ Ion-Exchanged into H-CHA

Passive NOx adsorbers (PNA) using Pd/zeolites have emerged as a promising solution for the reduction of cold-start emissions from vehicle exhaust. However, the nature of the active sites and the mechanisms underlying NOx adsorption in Pd/zeolites remain a subject of ongoing investigation. In this study, we employ quantum chemical simulations to investigate the structure of Pd species in cation-exchange sites at isolated Al and Al pairs in the 6-ring and 8-ring of the CHA framework, before the introduction of NOx. Our calculations show that the speciation of Pd in these exchange sites strongly depends on the precise Al arrangement within the framework, as well as the operating conditions. Ionically dispersed Pd is found to be the most favorable species over a wide range of oxidizing and reducing conditions. Small oligomers of PdO and metallic Pd do not appear to be competitive at either isolated Al or Al pairs. Notably, our calculations show that ion exchange sites other than next-next-nearest neighbor Al pairs in the 6-ring will be preferentially occupied by Pd+ instead of Pd2+. The stability of Pd+ in the zeolite environment is an interesting contrast with its rareness in molecular Pd compounds. Nonetheless, a detailed analysis of the electronic structure shows that predicted Pd oxidation states are consistent with chemical intuition for all complexes investigated in this study. We also discuss the potential ambiguity in Pd characterization provided by typical experimental techniques such as XANES, EXAFS and UV-VIS, and highlight the need for additional EPR spectroscopy studies to further elucidate the initial Pd speciation in zeolites for PNA applications.

Ni(ii) dithiolate anion composites with two-dimensional materials for electrochemical oxygen evolution reactions (OERs)

A redox active anionic nickel dithiolate complex was synthesized and its composites with GO, rGO and GN were prepared and used as electrocatalyts in the OER.

Electrocatalytic Proton Reduction by a Cobalt Complex Containing a Proton-Responsive Bis(alkylimdazole)methane Ligand: Involvement of a C-H Bond in H2 Formation

Homogeneous electrocatalytic proton reduction is reported using cobalt complex [1](BF₄)₂. This complex comprises two bis(1‐methyl‐4,5‐diphenyl‐1H‐imidazol‐2‐yl)methane (HBMIMPh2 ) ligands that contain an acidic methylene moiety in their backbone. Upon reduction of [1](BF₄)₂ by either electrochemical or chemical means, one of its HBMIMPh2 ligands undergoes deprotonation under the formation of dihydrogen. Addition of a mild proton source (acetic acid) to deprotonated complex [2](BF₄) regenerates protonated complex [1](BF₄)₂. In presence of acetic acid in acetonitrile solvent [1](BF₄)₂ shows electrocatalytic proton reduction with a k_obs of ≈200 s⁻¹ at an overpotential of 590 mV. Mechanistic investigations supported by DFT (BP86) suggest that dihydrogen formation takes place in an intramolecular fashion through the participation of a methylene C−H bond of the HBMIMPh2 ligand and a Co^II−H bond through formal heterolytic splitting of the latter. These findings are of interest to the development of responsive ligands for molecular (base)metal (electro)catalysis.

The rupture mechanism of rubredoxin is more complex than previously thought

The surprisingly low rupture force and remarkable mechanical anisotropy of rubredoxin have been known for several years. Exploiting the first combination of steered molecular dynamics and the quantum chemical Judgement of Energy DIstribution (JEDI) analysis, the common belief that hydrogen bonds between neighboring amino acid backbones and the sulfur atoms of the central FeS₄ unit in rubredoxin determine the low mechanical resistance of the protein is invalidated. The distribution of strain energy in the central part of rubredoxin is elucidated in real-time with unprecedented detail, giving important insights into the mechanical unfolding pathway of rubredoxin. While structural anisotropy as well as the contribution of angle bendings in the FeS₄ unit have a significant influence on the mechanical properties of rubredoxin, these factors are insufficient to explain the experimentally observed low rupture force. Instead, the rupture mechanism of rubredoxin is far more complex than previously thought and requires more than just a hydrogen bond network.

Effect of Stacking Interactions on the Translation of Structurally Related Bis(thiosemicarbazonato)nickel(II) HER Catalysts to Modified Electrode Surfaces

A series of crystalline nickel(II) complexes (1-3) based on inexpensive bis(thiosemicarbazone) ligands diacetylbis(4-methyl-3-thiosemicarbazone) (H₂ATSM), diacetylbis(4,4-dimethyl-3-thiosemicarbazone) (H₂ATSDM), and diacetylbis[4-(2,2,2-trifluoroethyl)-3-thiosemicarbazone] (H₂ATSM-F₆) were synthesized and characterized by single-crystal X-ray diffraction and NMR, UV-visible, and Fourier transform infrared spectroscopies. Modified electrodes GC-1-GC-3 were prepared with films of 1-3 deposited on glassy carbon and evaluated as potential hydrogen evolution reaction (HER) catalysts. HER studies in 0.5 M aqueous H₂SO₄ (10 mA cm^-2) revealed dramatic shifts in the overpotential from 0.740 to 0.450 V after extended cycling for 1 and 2. The charge-transfer resistances for GC-1-GC-3 were determined to be 270, 160, and 630 Ω, respectively. Characterization of the modified surfaces for GC-1 and GC-2 by scanning electron microscopy and Raman spectroscopy revealed similar crystalline coatings before HER that changed to surface-modified crystallites after conditioning. The surface of GC-3 had an initial glasslike appearance before HER that delaminated after HER. The differences in the surface morphology and the effect of conditioning are correlated with crystal-packing effects. Complexes 1 and 2 pack as columns of interacting complexes in the crystallographic a direction with short interplanar spacings between 3.37 and 3.54 Å. Complex 3 packs as columns of isolated molecules in the crystallographic b direction with long-range interplanar spacings of 9.40 Å.

Thermodynamic Properties of Hydrogen-Producing Cobaloxime Catalysts: A Density Functional Theory Analysis

Density functional theory calculations were carried out to study the electrochemical properties including reduction potentials, pK _a values, and thermodynamic hydricities of three prototypical cobaloxime complexes, Co(dmgBF₂)₂ (dmgBF₂ = difluoroboryl-dimethylglyoxime), Co(dmgH)₂ (dmgH = dimethylglyoxime), and Co(dmgH)₂(py)(Cl) (py = pyridine) in the acetonitrile (AN)-water solvent mixture. The electrochemical properties of Co(dmgBF₂)₂ in pure AN and pure water were also considered for comparison to reveal the key roles of the solvent on the catalytic reaction. In agreement with previous studies, hydrogen production pathways starting from reduction of the resting state of Co^II and involving formation of the Co^IIIH and Co^IIH intermediates are the favorable ones for both bimetallic and monometallic pathways. However, we found that in pure AN, both the Co^IIIH and Co^IIH intermediates can react with a proton to produce H₂. In the presence of water in the solvent, the reduction of Co^IIIH to Co^IIH is necessary for the reaction with a proton to occur to form H₂. This suggests that it is possible to design catalytic systems by suitably tuning the composition of the AN-water mixture. We also identified the key role of axial coordination of the solvent molecules in affecting the catalytic reaction, which allows further catalyst design strategy. The highest hydride donor ability of Co(dmgH)₂(py)(Cl) indicates that this complex displays the best catalytic hydrogen-producing performance among the three cobaloximes studied in this work.

Approach to the Mechanism of Hydrogen Evolution Electrocatalyzed by a Model Co Clathrochelate: A Theoretical Study by Density Functional Theory

The hydrogen evolution reaction (HER) has attracted much attention within the scientific community because of increasing demands of modern society for clean and renewable energy sources. Molecular complexes of 3d-transition metals, such as cobalt, hold potential to replace platinum for the HER in acidic media. Among these, cage complexes such as tris-glyoximate metal clathrochelates, have demonstrated promising catalytic properties towards the HER. However, it is not clear whether the catalytic activity of this molecule stems from metal-centered activation of H⁺ , due to a low oxidation state of the metal stabilized by the surrounding organic cage, or if it is the organic cage playing a further cooperative role in bringing protons together. Herein, we report on a density functional theory study of two possible mechanisms for the HER catalyzed by a model Co clathrochelate. To assess the putative ligand involvement in the mechanism, several combinations of single and double protonation sites were investigated. The structural and energetic analysis of relevant intermediates suggests that the electrocatalytic mechanism is not based on the cooperation between the ligand and the metal. Instead, it is mainly due to the activation of H⁺ by the Co metallocenter. Our calculations further suggest that the last step in the mechanism is a proton coupled electron transfer step.

Reactivity of Two-Electron-Reduced Boron Formazanate Compounds with Electrophiles: Facile N-H/N-C Bond Homolysis Due to the Formation of Stable Ligand Radicals

The reactivity of a boron complex with a redox-active formazanate ligand, LBPh₂ [L = PhNNC(p-tol)NNPh], was studied. Two-electron reduction of this main-group complex generates the stable, nucleophilic dianion [LBPh₂]^2–, which reacts with the electrophiles BnBr and H₂O to form products that derive from ligand benzylation and protonation, respectively. The resulting complexes are anionic boron analogues of leucoverdazyls. N–C and N–H bond homolysis of these compounds was studied by exchange NMR spectroscopy and kinetic experiments. The weak N–C and N–H bonds in these systems derive from the stability of the resulting borataverdazyl radical, in which the unpaired electron is delocalized over the four N atoms in the ligand backbone. We thus demonstrate the ability of this system to take up two electrons and an electrophile (E⁺ = Bn⁺, H⁺) in a process that takes place on the organic ligand. In addition, we show that the [2e^–/E⁺] stored on the ligand can be converted to E^• radicals, reactivity that has implications in energy storage applications such as hydrogen evolution.

A boron complex with a redox-active formazanate ligand in its two-electron-reduced state is shown to react with electrophiles (BnBr and H⁺). The resulting “borataleucoverdazyl” products have weak N−C and N−H bonds; homolytic cleavage reactions lead to stable ligand-based radicals. Thus, the accumulation of [2e⁻/E⁺] on the formazanate ligand and conversion to E^• radicals are demonstrated, and their potential relevance in energy-related electrocatalysis (e.g., proton reduction) is discussed.

Toward Quantitatively Accurate Calculation of the Redox-Associated Acid-Base and Ligand Binding Equilibria of Aquacobalamin

Redox processes in complex transition metal-containing species are often intimately associated with changes in ligand protonation states and metal coordination number. A major challenge is therefore to develop consistent computational approaches for computing pH-dependent redox and ligand dissociation properties of organometallic species. Reduction of the Co center in the vitamin B12 derivative aquacobalamin can be accompanied by ligand dissociation, protonation, or both, making these properties difficult to compute accurately. We examine this challenge here by using density functional theory and continuum solvation to compute Co-ligand binding equilibrium constants (Kon/off), pKas, and reduction potentials for models of aquacobalamin in aqueous solution. We consider two models for cobalamin ligand coordination: the first follows the hexa, penta, tetra coordination scheme for Co(III), Co(II), and Co(I) species, respectively, and the second model features saturation of each vacant axial coordination site on Co(II) and Co(I) species with a single, explicit water molecule to maintain six directly interacting ligands or water molecules in each oxidation state. Comparing these two coordination schemes in combination with five dispersion-corrected density functionals, we find that the accuracy of the computed properties is largely independent of the scheme used, but including only a continuum representation of the solvent yields marginally better results than saturating the first solvation shell around Co throughout. PBE performs best, displaying balanced accuracy and superior performance overall, with RMS errors of 80 mV for seven reduction potentials, 2.0 log units for five pKas and 2.3 log units for two log Kon/off values for the aquacobalamin system. Furthermore, we find that the BP86 functional commonly used in corrinoid studies suffers from erratic behavior and inaccurate descriptions of Co-axial ligand binding, leading to substantial errors in predicted pKas and Kon/off values. These findings demonstrate the effectiveness of the present approach for computing electrochemical and thermodynamic properties of a complex transition metal-containing cofactor.

A cobalt complex with a bioinspired molybdopterin-like ligand: a catalyst for hydrogen evolution

A cobalt complex using a bioinspired ligand, that mimics the molybdopterin cofactor, displays very good activity for electrochemical proton reduction in terms of turnover frequency, faradic yields and stability.

DFT study of the mechanism of hydrogen evolution catalysed by molecular Ni, Co and Fe catalysts containing a diamine–tripyridine ligand

Electrolysis of water to obtain hydrogen is a practical way to transform surplus electrical power into clean and sustainable hydrogen fuels.