Quantum Chemical Modeling of Homogeneous Water Oxidation Catalysis

Mechanistic elucidation of O2 production from tBuOOH in water using the Mn(II) catalyst [Mn2(mcbpen)2(H2O)2]2+: a DFT study

This study employs density functional theory at the SMD/B3LYP-D3/6-311+G(2d,p),def2-TZVPP//SMD/B3LYP-D3/6-31G(d),SDD level of theory to explore the mechanistic details of O2 generation from tBuOOH, using H218O as the solvent, in the presence of the Mn(II) catalyst [Mn2(mcbpen)2(H2O)2]2+. Since this chemistry was reported to occur through the reaction of Mn(III)(μ-O)Mn(IV)-O˙ with water, we first revaluated this proposal and found that it occurs with an activation barrier greater than 36 kcal mol-1, ruling out the functioning of such a dimer as the active catalyst. Experimental evidence has shown that the oxidation of [Mn2(mcbpen)2(H2O)2]2+ by tBuOOH in H218O produces the Mn(IV) species [Mn(18O)(mcbpen)]+. Our investigations revealed a plausible mechanism for this observation in which [Mn (18O)(mcbpen)]+ acts as the active catalyst, generating the tert-butyl peroxyl radical (tBuOO˙) through its reaction with tBuOOH. In this proposed mechanism, the O-O bond is formed through the interaction of tBuOO˙ with another [Mn(18O)(mcbpen)]+, finally leading to the formation of the 16O18O product. Our findings underscore the pivotal role of [Mn(18O)(mcbpen)]+ in both generating the active species tBuOO˙ and consuming it to produce 16O18O. With activation barriers as low as about 9 kcal mol-1, these elementary steps highlight the feasibility of our proposed mechanism. Moreover, this mechanism elucidates why, experimentally, one of the oxygen atoms in the released O2 comes from water, while the other originates from tBuOOH. This research broadens our understanding of high oxidation state manganese chemistry, setting the stage for the development of more efficient Mn-based catalysts, aimed at improving processes in both renewable energy and synthetic chemistry.

Mechanisms of Mn(V)-Oxo to Mn(IV)-Oxyl Conversion: From Closed-Cubane Photosystem II to Mn(V) Catalysts and the Role of the Entering Ligands

The low activation barrier for O-O coupling in the closed-cubane Oxygen-Evolving Centre (OEC) of Photosystem II (PSII) requires water coordination with the Mn4 'dangler' ion in the Mn(V)-oxo fragment. This coordination transforms the Mn(V)-oxo complex into a more reactive Mn4(IV)-oxyl species, enhancing O-O coupling. This study explains the mechanism behind the coordination and indicates that in the most stable form of the OEC, the Mn4 fragment adopts a trigonal bipyramidal geometry but needs to transition to a square pyramidal form to be activated for O-O coupling. This transition stabilizes the Mn4 dxy orbital, enabling electron transfer from the oxo ligand to the dxy orbital, converting the oxo ligand into an oxyl species. The role of the water is to coordinate with the square pyramidal structure, reducing the energy gap between the oxo and oxyl forms, thereby lowering the activation energy for O-O coupling. This mechanism applies not only to the OEC system but also to other Mn(V)-based catalysts. For other catalysts, ligands such as OH- stabilize the Mn(IV)-oxyl species better than water, improving catalyst activation for reactions like C-H bond activation. This study is the first to explain the Mn(V)-oxo to Mn(IV)-oxyl conversion, providing a new foundation for Mn-based catalyst design.

Elucidating the catalytic mechanisms of O2 generation by [Mn2(μ-O)2(terpy)2(OH2)2]3+ using DFT calculations: a focus on ClO- as oxidant

The experimentally reported Mn(IV)Mn(III) complex [Mn2(μ-O)2(terpy)2(OH2)2]3+ has been observed catalyzing O2 generation with oxidants like ClO- and HSO5-. Previous mechanistic studies primarily focused on O2 generation with HSO5-, concluding that Mn(IV)Mn(III) acts as a catalyst, generating a Mn(IV)Mn(IV)-oxyl species as a key intermediate responsible for O-O bond formation. This computational study employs DFT calculations to investigate whether the catalytic generation of O2 using ClO- follows the same mechanism previously identified with HSO5- as the oxidant, or if it proceeds through an alternate pathway. To this end, we explored multiple pathways using ClO- as the oxidant. Interestingly, our findings confirm that in the case of ClO- as the oxidant, similar to what was observed with HSO5-, the Mn(IV)Mn(IV)-oxyl species indeed plays a crucial role in driving the catalytic evolution of O2 with the potential formation of the binuclear complexes Mn(IV)Mn(IV)-oxy and Mn(IV)Mn(IV)-OH during the reaction. These complexes are reactive in producing O2, with activation free energies of 15.9 and 14.3 kcal mol-1, respectively. However, our calculations revealed that the Mn(IV)Mn(IV)-oxyl complex is significantly more reactive in producing O2 than Mn(IV)Mn(IV)-oxy and Mn(IV)Mn(IV)-OH, with a lower free energy barrier of 8.1 kcal mol-1. Consequently, even though Mn(IV)Mn(IV)-oxyl is predicted to be present in much lower concentrations than Mn(IV)Mn(IV)-oxy and Mn(IV)Mn(IV)-OH, it emerges as the species acting as the active catalyst for catalytic O2 generation. This study enhances our knowledge of high oxidation state (+3 and +4) manganese chemistry, highlighting its key role in catalysis and paving the way for more efficient Mn-based catalysts with broad applications.

Exploring the Cooperation of the Redox Non-Innocent Ligand and Di-Cobalt Center for the Water Oxidation Reaction Catalyzed by a Binuclear Complex

Water oxidation is a crucial reaction in the artificial photosynthesis system. In the present work, density functional calculations were employed to decipher the mechanism of water oxidation catalyzed by a binuclear cobalt complex, which was disclosed to be a homogeneous water oxidation catalyst in pH=7 phosphate buffer. The calculations showed that the catalytic cycle starts from the CoIII,III-OH2 species. Then, a proton-coupled electron transfer followed by a one-electron transfer process leads to the generation of the formal CoIV,IV-OH intermediate. The subsequent PCET produces the active species, namely the formal CoIV,V=O intermediate (4). The oxidation processes mainly occur on the ligand moiety, including the coordinated water moiety, implying a redox non-innocent behavior. Two cobalt centers keep their oxidation states and provide one catalytic center for water activation during the oxidation process. 4 triggers the O-O bond formation via the water nucleophilic attack pathway, in which the phosphate buffer ion functions as the proton acceptor. The O-O bond formation is the rate-limiting step with a calculated total barrier of 17.7 kcal/mol. The last electron oxidation process coupled with an intramolecular electron transfer results in the generation of O2.

Evolution in the Design of Water Oxidation Catalysts with Transition-Metals: A Perspective on Biological, Molecular, Supramolecular, and Hybrid Approaches

Increased demand for a carbon-neutral sustainable energy scheme augmented by climatic threats motivates the design and exploration of novel approaches that reserve intermittent solar energy in the form of chemical bonds in molecules and materials. In this context, inspired by biological processes, artificial photosynthesis has garnered significant attention as a promising solution to convert solar power into chemical fuels from abundantly found H2O. Among the two redox half-reactions in artificial photosynthesis, the four-electron oxidation of water according to 2H2O → O2 + 4H+ + 4e- comprises the major bottleneck and is a severe impediment toward sustainable energy production. As such, devising new catalytic platforms, with traditional concepts of molecular, materials and biological catalysis and capable of integrating the functional architectures of the natural oxygen-evolving complex in photosystem II would certainly be a value-addition toward this objective. In this review, we discuss the progress in construction of ideal water oxidation catalysts (WOCs), starting with the ingenuity of the biological design with earth-abundant transition metal ions, which then diverges into molecular, supramolecular and hybrid approaches, blurring any existing chemical or conceptual boundaries. We focus on the geometric, electronic, and mechanistic understanding of state-of-the-art homogeneous transition-metal containing molecular WOCs and summarize the limiting factors such as choice of ligands and predominance of environmentally unrewarding and expensive noble-metals, necessity of high-valency on metal, thermodynamic instability of intermediates, and reversibility of reactions that create challenges in construction of robust and efficient water oxidation catalyst. We highlight how judicious heterogenization of atom-efficient molecular WOCs in supramolecular and hybrid approaches put forth promising avenues to alleviate the existing problems in molecular catalysis, albeit retaining their fascinating intrinsic reactivities. Taken together, our overview is expected to provide guiding principles on opportunities, challenges, and crucial factors for designing novel water oxidation catalysts based on a synergy between conventional and contemporary methodologies that will incite the expansion of the domain of artificial photosynthesis.

Mechanistic Insights into Photochemical CO2 Reduction to CH4 by a Molecular Iron-Porphyrin Catalyst

Iron tetraphenylporphyrin complex modified with four trimethylammonium groups (Fe-p-TMA) is found to be capable of catalyzing the eight-electron eight-proton reduction of CO₂ to CH₄ photochemically in acetonitrile. In the present work, density functional theory (DFT) calculations have been performed to investigate the reaction mechanism and to rationalize the product selectivity. Our results revealed that the initial catalyst Fe-p-TMA ([Cl-Fe(III)-LR₄]⁴⁺, where L = tetraphenylporphyrin ligand with a total charge of -2, and R₄ = four trimethylammonium groups with a total charge of +4) undergoes three reduction steps, accompanied by the dissociation of the chloride ion to form [Fe(II)-L^••2-R₄]²⁺. [Fe(II)-L^••2-R₄]²⁺, bearing a Fe(II) center ferromagnetically coupled with a tetraphenylporphyrin diradical, performs a nucleophilic attack on CO₂ to produce the ¹η-CO₂ adduct [CO₂^•--Fe(II)-L^•-R₄]²⁺. Two intermolecular proton transfer steps then take place at the CO₂ moiety of [CO₂^•--Fe(II)-L^•-R₄]²⁺, resulting in the cleavage of the C-O bond and the formation of the critical intermediate [Fe(II)-CO]⁴⁺ after releasing a water molecule. Subsequently, [Fe(II)-CO]⁴⁺ accepts three electrons and one proton to generate [CHO-Fe(II)-L^•-R₄]²⁺, which finally undergoes a successive four-electron-five-proton reduction to produce methane without forming formaldehyde, methanol, or formate. Notably, the redox non-innocent tetraphenylporphyrin ligand was found to play an important role in CO₂ reduction since it could accept and transfer electron(s) during catalysis, thus keeping the ferrous ion at a relatively high oxidation state. Hydrogen evolution reaction via the formation of Fe-hydride ([Fe(II)-H]³⁺) turns out to endure a higher total barrier than the CO₂ reduction reaction, therefore providing a reasonable explanation for the origin of the product selectivity.

Copper Complex Catalyzed Two-Electron and Proton Shuttle Mechanism of O-O Bond Formation from DFT-Based Metadynamics Simulations

We performed first-principles metadynamics simulations to explore the mechanistic pathway of oxygen-oxygen bond formation catalyzed by cis-bis(hydroxo) and cis-(hydroxo)oxo copper complexes. The ligands of considered complexes involve modified bipyridine ligands with oxo and hydroxo groups on 6, 6' positions. The study focuses on the kinetics and thermodynamics of the oxygen-oxygen bond formation. The individual migration of the proton to the hydroxyl group and hydroxide to the oxo and hydroxo moieties of the complexes was examined. The proton transfer requires more kinetic barrier than the hydroxide migration. The nature of the electronic density was analyzed with the help of spin population analysis. The molecular orbitals and natural orbital analysis were carried out to examine the nature of the orbitals involved in the oxygen-oxygen bond formation. The σ*(d_x²_-y²-p_x) molecular orbital of the Cu-O or Cu-OH bond overlaps with the p_z orbital of the hydroxide ion in forming the oxygen-oxygen bond. The two-electron two-centered (2e^--2C) bond is observed in the oxygen-oxygen bond formation. In the oxidation process, these ligands stabilize the electron density from the water or hydroxide ion. These redox-active ligands also help stabilize the formed hydrogen peroxide or peroxide complexes.

QM Calculations Revealed that Outer-Sphere Electron Transfer Boosted O-O Bond Cleavage in the Multiheme-Dependent Cytochrome bd Oxygen Reductase

The cytochrome bd oxygen reductase catalyzes the four-electron reduction of dioxygen to two water molecules. The structure of this enzyme reveals three heme molecules in the active site, which differs from that of heme-copper cytochrome c oxidase. The quantum chemical cluster approach was used to uncover the reaction mechanism of this intriguing metalloenzyme. The calculations suggested that a proton-coupled electron transfer reduction occurs first to generate a ferrous heme b₅₉₅. This is followed by the dioxygen binding at the heme d center coupled with an outer-sphere electron transfer from the ferrous heme b₅₉₅ to the dioxygen moiety, affording a ferric ion superoxide intermediate. A second proton-coupled electron transfer produces a heme d ferric hydroperoxide, which undergoes efficient O-O bond cleavage facilitated by an outer-sphere electron transfer from the ferrous heme b₅₉₅ to the O-O σ* orbital and an inner-sphere proton transfer from the heme d hydroxyl group to the leaving hydroxide. The synergistic benefits of the two types of hemes rationalize the highly efficient oxygen reduction repertoire for the multi-heme-dependent cytochrome bd oxygen reductase family.

Theoretical study on the mechanism of water oxidation catalyzed by a mononuclear copper complex: important roles of a redox non-innocent ligand and HPO4 2- anion

The water oxidation reaction is the bottleneck problem of the artificial photosynthetic system. In this work, the mechanism of water oxidation catalyzed by a mononuclear copper complex in alkaline conditions was studied by density functional calculations. Firstly, a water molecule coordinating with the copper center of the complex (Cuii, 1) generates Cuii-H2O (2). 2 undergoes two proton-coupled electron transfer processes to produce intermediate (4). The oxidation process occurs mainly on the ligand moiety, and 4 (˙L-Cuii-O˙) can be described as a Cuii center interacting with a ligand radical antiferromagnetically and an oxyl radical ferromagnetically. 4 is the active species that can trigger O-O bond formation via the water nucleophilic attack mechanism. This process occurs in a step-wise manner. The attacking water transfers one of the protons to the HPO4 2- coupled with an electron transfer to the ligand radical, which generates a transient OH˙ interacting with the oxyl radical and H2PO4 -. Then the O-O bond is formed through the direct coupling of the oxo radical and the OH radical. The triplet di-oxygen could be released after two oxidation processes. According to the Gibbs free energy diagram, the O-O bond formation was suggested to be the rate-limiting step with a calculated total barrier of 19.5 kcal mol-1. More importantly, the copper complex catalyzing water oxidation with the help of a redox non-innocent ligand and HPO4 2- was emphasized.

Three-Electron Two-Centered Bond and Single-Electron Transfer Mechanism of Water Splitting via a Copper-Bipyridine Complex

We report the atomistic and electronic details of the mechanistic pathway of the oxygen-oxygen bond formation catalyzed by a copper-2,2'-bipyridine complex. Density functional theory-based molecular dynamics simulations and enhanced sampling methods were employed for this study. The thermodynamics and electronic structure of the oxygen-oxygen bond formation are presented in this study by considering the cis-bishydroxo, [Cu^III(bpy)(OH)₂]⁺, and cis-(hydroxo)oxo, [Cu^IV(bpy)(OH)(═O)]⁺, complexes as active catalysts. In the cis-bishydroxo complex, the hydroxide transfer requires a higher kinetic barrier than the proton transfer process. In the case of [Cu^IV(bpy)(OH)(═O)]⁺, the proton transfer requires a higher free energy than the hydroxide one. The peroxide bond formation is thermodynamically favorable for the [Cu^IV(bpy)(OH)(═O)]⁺ complex compared with the other. The hydroxide ion is transferred to one of the Cu-OH moieties, and the proton is transferred to the solvent. The free energy barrier for this migration is higher than that for the former transfer. From the analysis of molecular orbitals, it is found that the electron density is primarily present on the water molecules near the active sites in the highest occupied molecular orbital (HOMO) state and lowest unoccupied molecular orbital (LUMO) of the ligands. Natural bond orbital (NBO) analysis reveals the electron transfer process during the oxygen-oxygen bond formation. The σ*Cu(d_xz)-O(p) orbitals are involved in the oxygen-oxygen bond formation. During the bond formation, three-electron two-centered (3e^--2C) bonds are observed in [Cu^III(bpy)(OH)₂]⁺ during the transfer of the hydroxide before the formation of the oxygen-oxygen bond.

Theoretical Study on the Electro-Reduction of Carbon Dioxide to Methanol Catalyzed by Cobalt Phthalocyanine

Density functional theory (DFT) calculations have been conducted to investigate the mechanism of cobalt(II) tetraamino phthalocyanine (CoPc-NH₂) catalyzed electro-reduction of CO₂. Computational results show that the catalytically active species 1 (⁴[Co^II(H₄L)]⁰) is formed by a four-electron-four-proton reduction of the initial catalyst CoPc-NH₂. Complex 1 can attack CO₂ after a one-electron reduction to give a [Co^III-CO₂^2-]^- intermediate, followed by a protonation and a one-electron reduction to give intermediate [Co^II-COOH]^- (4). Complex 4 is then protonated on its hydroxyl group by a carbonic acid to generate the critical species 6 (Co^III-L^•--CO), which can release the carbon monoxide as an intermediate (and also as a product). In parallel, complex 6 can go through a successive four-electron-four-proton reduction to produce the targeted product methanol without forming formaldehyde as an intermediate product. The high-lying π orbital and the low-lying π* orbital of the phthalocyanine endow the redox noninnocent nature of the ligand, which could be a dianion, a radical monoanion, or a radical trianion during the catalysis. The calculated results for the hydrogen evolution reaction indicate a higher energy barrier than the carbon dioxide reduction. This is consistent with the product distribution in the experiments. Additionally, the amino group on the phthalocyanine ligand was found to have a minor effect on the barriers of critical steps, and this accounts for the experimentally observed similar activity for these two catalysts, namely, CoPc-NH₂ and CoPc.

Water Oxidation Catalyzed by a Bioinspired Tetranuclear Manganese Complex: Mechanistic Study and Prediction

Density functional theory calculations were utilized to elucidate the water oxidation mechanism catalyzed by polyanionic tetramanganese complex a [Mn^III ₃ Mn^IV O₃ (CH₃ COO)₃ (A-α-SiW₉ O₃₄ )]^6- . Theoretical results indicated that catalytic active species 1 (Mn₄ ^{III,III,IV,IV} ) was formed after O₂ formation in the first turnover. From 1, three sequential proton-coupled electron transfer (PCET) oxidations led to the Mn^IV -oxyl radical 4 (Mn₄ ^IV,IV,IV,IV -O⋅). Importantly, 4 had an unusual butterfly-shaped Mn₂ O₂ core for the two substrate-coordinated Mn sites, which facilitated O-O bond formation via direct coupling of the oxyl radical and the adjacent Mn^IV -coordinated hydroxide to produce the hydroperoxide intermediate Int1 (Mn₄ ^III,IV,IV,IV -OOH). This step had an overall energy barrier of 24.9 kcal mol^-1 . Subsequent PCET oxidation of Int1 to Int2 (Mn₄ ^III,IV,IV,IV -O₂ ⋅) enabled the O₂ release in a facile process. Furthermore, apart from the Si-centered complex, computational study suggested that tetramanganese polyoxometalates with Ge, P, and S could also catalyze the water oxidation process, where those bearing P and S likely present higher activities.

Non-heme oxoiron complexes as active intermediates in the water oxidation process with hydrogen/oxygen atom transfer reactions

In this study, we explore the water oxidation process with the help of density functional theory. The formation of an oxygen-oxygen bond is crucial in the water oxidation process. Here, we report the formation of the oxygen-oxygen bond by the N5-coordinate oxoiron species with a higher oxidation state of Fe^IV and Fe^V. This bond formation is studied through the nucleophilic addition of water molecules and the transfer of the oxygen atom from meta-chloroperbenzoic acid (mCPBA). Our study reveals that the oxygen-oxygen bond formation by reacting mCPBA with Fe^VO requires less activation barrier (13.7 kcal mol^-1) than the other three pathways. This bond formation by the oxygen atom transfer (OAT) pathway is more favorable than that achieved by the hydrogen atom transfer (HAT) pathway. In both cases, the oxygen-oxygen bond formation occurs by interacting the σ*d_z²-2p_z molecular orbital of the iron-oxo intermediate with the 2p_x orbital of the oxygen atom. From this study, we understand that the oxygen-oxygen bond formation by Fe^IVO with the OAT process is also feasible (16 kcal mol^-1), suggesting that Fe^VO may not always be required for the water oxidation process by non-heme N5-oxoiron. After the oxygen-oxygen bond formation, the release of the dioxygen molecule occurs with the addition of the water molecule. The release of dioxygen requires a barrier of 7.0 kcal mol^-1. The oxygen-oxygen bond formation is found to be the rate-determining step.

Mechanistic Insights into Cobalt-Based Water Oxidation Catalysis by DFT-Based Molecular Dynamics Simulations

We report the mechanistic details of the water oxidation process by the complex, [Co^II(bpbH₂)Cl₂], where bpbH₂ = N, N'-bis(2'-pyridinecarboxamide)-1,2-benzene. An experimental study reported the complex as the efficient catalyst for the water oxidation process. We performed density functional theory calculations at the M06-L level and first-principles molecular dynamics simulations to study the catalytic nature of the complex. We investigated the energetics of the total catalytic cycle, which combines the oxygen-oxygen bond formation, proton-coupled electron transfer (PCET) steps, and release of oxygen molecule. The formed peroxide and superoxide intermediates in the catalytic cycle were characterized with the help of the Mulliken spin density parameters. Mulliken spin densities of the metal-oxo bond reveal that the triplet state of Co^V═O has a double-bond nature, but the quintet state of the complex has a radical nature (Co^IV-O^•-). In an alternative way, the deprotonation of the amide groups of the ligand is also considered. The deprotonation and formation of higher oxidation metal-oxo intermediates are also possible. In addition to this, we have considered the effect of phosphate buffer on water nucleophilic addition. The oxygen-oxygen bond formation is favorable by the catalyst with the deprotonated form of the ligand, with the addition of water as the nucleophile. In the oxidation process, the C═O bonds of the ligand transfer the electron density to nitrogen atoms, stabilizing the higher order oxo, peroxide, and superoxide bonds. The oxygen-oxygen bond formation is the rate-determining step in the overall water oxidation process. This bond was further investigated using first-principles molecular dynamics at the PBE-D2 level. The dynamics of proton, hydroxide ion, and the nature of the ligand structure on the oxygen-oxygen bond were examined. We find that the oxygen molecule is released from the superoxide complex with the addition of water molecules.

Molecular Mechanism of the Mononuclear Copper Complex-Catalyzed Water Oxidation from Cluster-Continuum Model Calculations

Cluster-continuum model calculations were conducted to decipher the mechanism of water oxidation catalyzed by a mononuclear copper complex. Among various O-O bond formation mechanisms investigated in this study, the most favorable pathway involved the nucleophilic attack of OH^- onto the ^.+ L-Cu^II -OH^- intermediate. During such process, the initial binding of OH^- to the proximity of ^.+ L-Cu^II -OH^- would result in the spontaneous oxidation of OH^- , leading to OH⋅ radical and Cu^II -OH^- species. The further O-O coupling between OH⋅ radical and Cu^II -OH^- was associated with a barrier of 14.8 kcal mol^-1 , leading to the formation of H₂ O₂ intermediate. Notably, the formation of "Cu^III -O^.- " species, a widely proposed active species for O-O bond formation, was found to be thermodynamically unfavorable and could be bypassed during the catalytic reactions. On the basis the present calculations, a catalytic cycle of the mononuclear copper complex-catalyzed water oxidation was proposed.

Mechanism and selectivity of photocatalyzed CO2 reduction by a function-integrated Ru catalyst

The phosphine-substituted Ru(II) polypyridyl complex, [Ru^II-(tpy)(pqn)(MeCN)]²⁺ (RuP), was disclosed to be an efficient photocatalyst for the reduction of CO₂ to CO with excellent selectivity. In this work, density functional calculations were performed to elucidate the reaction mechanism and understand the origin of selectivity. The calculations showed that RuP was first excited to the singlet excited state, followed by intersystem crossing to produce a triplet species (3RuIII(L˙-)-S), which was then reduced by the sacrificial electron donor BIH to generate a Ru^II(L˙^-) intermediate. The ligand of Ru^II(L˙^-) was further reduced to produce a Ru^II(L^2-) intermediate. The redox non-innocent nature of the tpy and pqn ligands endows the Ru center with an oxidation state of +2 after two one-electron reductions. Ru^II(L^2-) nucleophilically attacks CO₂, in which two electrons are delivered from the ligands to CO₂, affording a Ru^II-COOH species after protonation. This is followed by the protonation of the hydroxyl moiety of Ru^II-COOH, coupled with the C-O bond cleavage, resulting in the formation of Ru^II-CO. Ultimately, CO is dissociated after two one-electron reductions. Protonation of Ru^II(L^2-) to generate a Ru^II-hydride, a critical intermediate for the production of formate and H₂, turns out to be kinetically less favorable, even though it is thermodynamically more favorable. This fact is due to the presence of a Ru²⁺ ion in the reduced catalyst, which disfavors its protonation.

Computational investigation into intramolecular hydrogen bonding controlling the isomer formation and pKa of octahedral nickel (II) proton reduction catalysts

This work demonstrates the impact of intramolecular hydrogen bonding (H-bonding) on the calculated pKa of octahedral tris-(pyridinethiolato)nickel (II), [Ni(PyS)3]-, proton reduction catalysts. Density Functional Theory (DFT) calculations on a [Ni(PyS)3]-...

Theoretical Investigation into Selective Benzene Hydroxylation by Ruthenium-Substituted Keggin-Type Polyoxometalates

Benzene hydroxylation catalyzed by ruthenium-substituted Keggin-type polyoxometalates [Ru^V(O)XW₁₁O₃₉]^n- (Ru^VOX; X = Al, Ga, Si, Ge, P, As, S; heteroatoms; 3 ≤ n ≤ 6) is investigated using the density functional theory approach. As a possible side reaction, the water oxidation reaction is also considered. We found that the rate-determining step for water oxidation by Ru^VOX requires a higher activation free energy than the benzene hydroxylation reaction, suggesting that all of the Ru^VOX catalysts show high chemoselectivity toward benzene hydroxylation. Additionally, the heteroatom effect in benzene hydroxylation by Ru^VOX is discussed. The replacement of Si by X induces changes in the bond length of μ₄O-X, resulting in a change in the activation free energy for benzene hydroxylation by Ru^VOX. Consequentially, Ru^VOS is expected to be the most effective catalyst among the (Ru^VOX) catalysts for the benzene hydroxylation reaction.

Mechanistic Insight into the O2 Evolution Catalyzed by Copper Complexes with Tetra- and Pentadentate Ligands

The mononuclear complexes ([(bztpen)Cu] (BF₄)₂ (bztpen = N-benzyl-N,N',N'-tris (pyridin-2-yl methyl ethylenediamine))) and ([(dbzbpen)Cu(OH₂)] (BF₄)₂ (dbzbpen = N,N'-dibenzyl-N,N'-bis(pyridin-2-ylmethyl) ethylenediamine)) have been reported as water oxidation catalysts in basic medium (pH = 11.5). We explore the O₂ evolution process catalyzed by these copper catalysts with various ligands (L) by applying the first-principles molecular dynamics simulations. First, the oxidation of catalysts to the metal-oxo intermediates [LCu(O)]²⁺ occurs through the proton-coupled electron transfer (PCET) process. These intermediates are involved in the oxygen-oxygen bond formation through the water-nucleophilic addition process. Here, we have considered two types of oxygen-oxygen bond formation. The first one is the transfer of the hydroxide of the water molecule to the Cu═O moiety; the proton transfer to the solvent leads to the formation of the peroxide complex ([LCu(OOH)]⁺). The other is the formation of the hydrogen peroxide complex ([LCu(HOOH)]²⁺) by the transfer of proton and hydroxide of the water molecule to the metal-oxo intermediate. The formation of the peroxide complex requires less activation free energy than hydrogen peroxide formation for both catalysts. We found two transition states in the well-tempered metadynamics simulations: one for proton transfer and another for hydroxide transfer. In both cases, the proton transfer requires higher free energy. Following the formation of the oxygen-oxygen bond, we study the release of the dioxygen molecule. The formed peroxide and hydrogen peroxide complexes are converted into the superoxide complex ([LCu(OO)]²⁺) through the transfer of proton, electron, and PCET processes. The superoxide complex releases an oxygen molecule upon the addition of a water molecule. The free energy of activation for the release of the dioxygen molecule is lesser than that of the oxygen-oxygen bond formation. When we observe the entire water oxidation process, the oxygen-oxygen bond formation is the rate-determining step. We calculated the rates of reaction by using the Eyring equation and found them to be close to the experimental values.

Electrocatalytic Isoxazoline-Nanocarbon Metal Complexes

We report the synthesis of new carbon-nanomaterial-based metal chelates that enable effective electronic coupling to electrocatalytic transition metals. In particular, multiwalled carbon nanotubes (MWCNTs) and few-layered graphene (FLG) were covalently functionalized by a microwave-assisted cycloaddition with nitrile oxides to form metal-binding isoxazoline functional groups with high densities. The covalent attachment was evidenced by Raman spectroscopy, and the chemical identity of the surface functional groups was confirmed by X-ray photoelectron spectroscopy (XPS) and time-of-flight secondary ion mass spectrometry (ToF-SIMS). The functional carbon nanomaterials effectively chelate precious metals Ir(III), Pt(II), and Ru(III), as well as earth-abundant metals such as Ni(II), to afford materials with metal contents as high as 3.0 atom %. The molecularly dispersed nature of the catalysts was confirmed by X-ray absorption spectroscopy (XAS) and energy-dispersive X-ray spectroscopy (STEM-EDS) elemental mapping. The interplay between the chelate structure on the graphene surface and its metal binding ability has also been investigated by a combination of experimental and computational studies. The defined ligands on the graphene surfaces enable the formation of structurally precise heterogeneous molecular catalysts. The direct attachment of the isoxazoline functional group on the graphene surfaces provides strong electronic coupling between the chelated metal species and the conductive carbon nanomaterial support. We demonstrate that the metal-chelated carbon nanomaterials are effective heterogeneous catalysts in the oxygen evolution reaction with low overpotentials and tunable catalytic activity.

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.

Unveiling the High Catalytic Activity of a Dinuclear Iron Complex for the Oxygen Evolution Reaction

The dinuclear iron complex [(H₂O)-Fe^III-(ppq)-O-(ppq)-Fe^III-Cl]³⁺ (Fe^III(ppq), ppq = 2-(pyrid-2'-yl)-8-(1″,10″-phenanthrolin-2″-yl)-quinoline) demonstrates a catalytic activity about one order of magnitude higher than the mononuclear iron complex [Cl-Fe^III(dpa)-Cl]⁺ (Fe^III(dpa), dpa = N,N-di(1,10-phenanthrolin-2-yl)-N-isopentylamine) for the oxygen evolution reaction (OER). However, the mechanism behind such an unusually high activity has remained largely unclear. To solve this puzzle, a decomposition-and-reaction mechanism is proposed for the OER with the dinuclear Fe^III(ppq) complex as the initial state of the catalytic agent. In this mechanism, the high-valent dinuclear iron complex first dissociates into two mononuclear moieties, and the oxidized mononuclear iron complexes directly catalyze the formation of an O-O bond through a nitrate attack pathway with nitrate functioning as a cocatalyst. Density functional theory calculations reveal that it is the electron-deficient microenvironment around the iron center that gives rise to the remarkable catalytic activity observed experimentally. Therefore, the outstanding performance of the Fe^III(ppq) catalyst can be ascribed to the high reactivity of its mononuclear moieties in a high oxidation state, which is concomitant with the structural stability of the low-valent dinuclear complex. The theoretical insights provided by this study could be useful for the optimization and design of novel iron-based water oxidation catalysts.

Electronic Effect on Bimetallic Catalysts: Cleavage of Phosphodiester Mediated by Fe(III)-Zn(II) Purple Acid Phosphatase Mimics

The bimetallic system is an important strategy for the catalytic hydrolysis of phosphodiester. The purple acid phosphatase (PAPs) enzyme is a typical bimetallic catalyst in this field. Mechanistic details for the hydrolysis cleavage of the DNA dinucleotide analogue BNPP^- (BNPP^- = bis(p-nitrophenyl) phosphate) by hetero-binuclear [Fe^III(μ-OH)Zn^IIL]²⁺ complexes (L = 2-[N-bis(2-pyridylmethyl)-aminomethyl]-4-methyl-6-[N'-(2-pyridylmethyl)(2-hydroxybenzyl) aminomethyl] phenol) were investigated using density functional theory calculations. The catalysts with single-bridged hydroxyl and double-bridged hydroxyl groups were compared. The calculation results show that the doubly hydroxide-bridged complex could better bind to substrates. For the BNPP^- hydrolysis, the doubly hydroxide-bridged reactant isomerizes into a single hydroxide-bridged complex, and then the attack is initiated by the hydroxyl group on the iron center. In addition, the catalyst with the electron-donating group (Me) was determined to take precedence over electron-withdrawing groups (Br and NO₂ groups) in the hydrolysis reaction. This is because the substituents affect the high-lying occupied molecular orbitals, tuning the Lewis acidity of iron and pK_a values of the metal-bonded water. These factors influence the hydroxyl nucleophilicity, leading to changes in catalytic activity. To further examine substituent effects, the occupied orbital energies were calculated with several different substituent groups (-CF₃, -OMe, -OH, -NH₂, and -N(Me)₂). It was found that the HOMO or HOMO-1 energy decreases with the increase of the σ_p value. Further, the catalyst activity of the [Fe^III(μ-OH)Zn^IIL]²⁺ complexes was found to be mainly affected by the phenolate ligand (B) coordinated to the iron and zinc centers. These fundamental aspects of the hydrolysis reactions of BNPP^- catalyzed by [Fe^III(μ-OH)Zn^IIL]²⁺ complexes should contribute to improved understanding of the mechanism and to catalyst design involving hetero-binuclear metals complexes.

Understanding the performance of a bisphosphonate Ru water oxidation catalyst

Water oxidation catalysts (WOCs) are a key part of generating H₂ from water and sunlight, consequently, it is a promising process for the production of clean energy.

Theoretical Insight into the Performance of MnII/III-Monosubstituted Heteropolytungstates as Water Oxidation Catalysts

The performance of Mn^II/III-monosubstituted heteropolytungstates [Mn^III(H₂O)GeW₁₁O₃₉]^5- ([GT-Mn^III-OH₂]^5-, where GT = GeW₁₁O₃₉) and [Mn^II(H₂O)GeW₁₁O₃₉]^6- ([GT-Mn^II-OH₂]^6-) as water oxidation catalysts at pH 9 was explored using density functional theory calculations. The counterion effect was fully considered, in which five and six Na⁺ ions were included in the calculations for water oxidation catalyzed by [GT-Mn^III-OH₂]^5- and [GT-Mn^II-OH₂]^6-, respectively. The process of water oxidation catalysis was divided into three elemental stages: (i) oxidative activation, (ii) O-O bond formation, and (iii) O₂ evolution. In the oxidative activation stage, two electrons and two protons are removed from [Na₅-GT-Mn^III-OH₂] and three electrons and two protons are removed from [Na₆-GT-Mn^II-OH₂]. Therefore, the Mn^IV-O^• species [Na₅-GT-Mn^IV-O^•] is obtained. Two mechanisms, (i) water nucleophilic attack and (ii) oxo-oxo coupling, were demonstrated to be competitive in O-O bond formation triggered from [Na₅-GT-Mn^IV-O^•]. In the last stage, the O₂ molecule could be readily evolved from the peroxo or dinuclear species and the catalyst returns to the ground state after the coordination of a water molecule(s).

O-O Bond Formation and Liberation of Dioxygen Mediated by N5 -Coordinate Non-Heme Iron(IV) Complexes

Formation of the O-O bond is considered the critical step in oxidative water cleavage to produce dioxygen. High-valent metal complexes with terminal oxo (oxido) ligands are commonly regarded as instrumental for oxygen evolution, but direct experimental evidence is lacking. Herein, we describe the formation of the O-O bond in solution, from non-heme, N5 -coordinate oxoiron(IV) species. Oxygen evolution from oxoiron(IV) is instantaneous once meta-chloroperbenzoic acid is administered in excess. Oxygen-isotope labeling reveals two sources of dioxygen, pointing to mechanistic branching between HAT (hydrogen atom transfer)-initiated free-radical pathways of the peroxides, which are typical of catalase-like reactivity, and iron-borne O-O coupling, which is unprecedented for non-heme/peroxide systems. Interpretation in terms of [FeIV (O)] and [FeV (O)] being the resting and active principles of the O-O coupling, respectively, concurs with fundamental mechanistic ideas of (electro-) chemical O-O coupling in water oxidation catalysis (WOC), indicating that central mechanistic motifs of WOC can be mimicked in a catalase/peroxidase setting.

How Metal Ion Lewis Acidity and Steric Properties Influence the Barrier to Dioxygen Binding, Peroxo O-O Bond Cleavage, and Reactivity

Herein we quantitatively investigate how metal ion Lewis acidity and steric properties influence the kinetics and thermodynamics of dioxygen binding versus release from structurally analogous Mn-O₂ complexes, as well as the barrier to Mn peroxo O-O bond cleavage, and the reactivity of Mn oxo intermediates. Previously we demonstrated that the steric and electronic properties of Mn^III-OOR complexes containing N-heterocyclic (N^Ar) ligand scaffolds can have a dramatic influence on alkylperoxo O-O bond lengths and the barrier to alkylperoxo O-O bond cleavage. Herein, we examine the dioxygen reactivity of a new Mn^II complex containing a more electron-rich, less sterically demanding N^Ar ligand scaffold, and compare it with previously reported Mn^II complexes. Dioxygen binding is shown to be reversible with complexes containing the more electron-rich metal ions. The kinetic barrier to O₂ binding and peroxo O-O bond cleavage is shown to correlate with redox potentials, as well as the steric properties of the supporting N^Ar ligands. The reaction landscape for the dioxygen chemistry of the more electron-rich complexes is shown to be relatively flat. A total of four intermediates, including a superoxo and peroxo species, are observed with the most electron-rich complex. Two new intermediates are shown to form following the peroxo, which are capable of cleaving strong X-H bonds. In the absence of a sacrificial H atom donor, solvent, or ligand, serves as a source of H atoms. With TEMPOH as sacrificial H atom donor, a deuterium isotope effect is observed (k_H/k_D = 3.5), implicating a hydrogen atom transfer (HAT) mechanism. With 1,4-cyclohexadiene, 0.5 equiv of benzene is produced prior to the formation of an EPR detected Mn^IIIMn^IV bimetallic species, and 0.5 equiv after its formation.

Computational Approaches to Photoelectrode Design through Molecular Functionalization for Enhanced Photoelectrochemical Water Splitting

Photoelectrochemical water splitting is a promising carbon-free approach to produce hydrogen from water. A photoelectrochemical cell consists of a semiconductor photoelectrode in contact with an aqueous electrolyte. Its performance is sensitive to properties of the photoelectrode/electrolyte interface, which may be tuned through functionalization of the photoelectrode surface with organic molecules. This can lead to improvements in the photoelectrode's properties. This Minireview summarizes key computational investigations on using molecular functionalization to modify photoelectrode stability, barrier height, and catalytic activity. It is discussed how first-principles density functional theory, first-principles molecular dynamics, and device modeling simulations can provide predictive insights and complement experimental investigations of functionalized photoelectrodes. Challenges and future directions in the computational modeling of functionalized photoelectrode/electrolyte interfaces within the context of experimental studies are also highlighted.

Quantum Chemical Study of the Mechanism of Water Oxidation Catalyzed by a Heterotrinuclear Ru2 Mn Complex

The heterotrinuclear complex A {[Ru^II (H₂ O)(tpy)]₂ (μ-[Mn^II (H₂ O)₂ (bpp)₂ ])}⁴⁺ [tpy=2,2':6',2''-terpyridine, bpp=3,5-bis(2-pyridyl)pyrazolate] was found to catalyze water oxidation both electrochemically and photochemically with [Ru(bpy)₃ ]³⁺ (bpy=2,2'-bipyridine) as the photosensitizer and Na₂ S₂ O₈ as the electron acceptor in neutral phosphate buffer. The mechanism of water oxidation catalyzed by this unprecedented trinuclear complex was studied by density functional calculations. The calculations showed that a series of oxidation and deprotonation events take place from A, leading to the formation of complex 1 (formal oxidation state of Ru1^IV Mn^III Ru2^III ), which is the starting species for the catalytic cycle. Three sequential oxidations of 1 result in the generation of the catalytically competing species 4 (formal oxidation state of Ru1^IV Mn^V Ru2^IV ), which triggers the O-O bond formation. The direct coupling of two adjacent oxo ligands bound to Ru and Mn leads to the production of a superoxide intermediate Int1. This step was calculated to have a barrier of 7.2 kcal mol^-1 at the B3LYP*-D3 level. Subsequent O₂ release from Int1 turns out to be quite facile. Other possible pathways were found to be much less favorable, including water nucleophilic attack, the coupling of an oxo and a hydroxide, and the direct coupling pathway at a lower oxidation state (Ru^IV Mn^IV Ru^IV ).

Computational Characterization of Single-Electron Transfer Steps in Water Oxidation

The presence of single-electron transfer (SET) steps in water oxidation processes catalyzed by first-row transition metal complexes has been recently recognized, but the computational characterization of this type of process is not trivial. We report a systematic theoretical study based on density functional theory (DFT) calculations on the reactivity of a specific copper complex active in water oxidation that reacts through two consecutive single-electron transfers. Both inner-sphere (through transition state location) and outer-sphere (through Marcus theory) mechanisms are analyzed. The first electron transfer is found to operate through outer-sphere, and the second one through inner-sphere. The current work proposes a scheme for the systematic study of single-electron transfer in water oxidation catalysis and beyond.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

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.

Energetic Effects of a Closed System Approach Including Explicit Proton and Electron Acceptors as Demonstrated by a Mononuclear Ruthenium Water Oxidation Catalyst

When considering water oxidation catalysis theoretically, accounting for the transfer of protons and electrons from one catalytic intermediate to the next remains challenging: correction factors are usually employed to approximate the energetics of electron and proton transfer. Here these energetics were investigated using a closed system approach, which places the catalytic intermediate in a simulation box including proton and electron acceptors, as well as explicit solvent. As a proof of principle, the first two catalytic steps of the mononuclear ruthenium-based water oxidation catalyst [Ru(cy)(bpy)(H2O)]2+ were examined using Car-Parrinello Molecular Dynamics. This investigation shows that this approach offers added insight, not only into the free energy profile between two stable intermediates, but also into how the solvent environment impacts this dynamic evolution.