Impact of the Ligand Flexibility and Solvent on the O-O Bond Formation Step in a Highly Active Ruthenium Water Oxidation Catalyst

Cholinergic Inhibition and Antioxidant Potential of Gongronema latifolium Benth Leaf in Neurodegeneration: Experimental and In Silico Study

The use of Gongronema latifolium for the management of various forms of neurological disorders has generated a lot of interest in the need to further investigate its neurotherapeutic constituents. This work, therefore, focused on assessing the inhibitory potential of selected bioactive components derived from G. latifolium against key neurotherapeutic targets and oxidant species associated with neurodegeneration using in vitro analysis and biomolecular modelling. G. latifolium methanol extract (GLME), solvent partition, chromatographic fractions (A-F) of GLME and pregnane compounds (Iloneoside and marsectohexol) derived from fraction-B with the highest activity were investigated for in vitro acetylcholinesterase (AChE), butyrylcholinesterase (BChE), and monoamine oxidase (MAO) inhibition in addition to their in vitro antioxidant activities. The interactions of iloneoside, marsectohexol, and reference drugs with human acetylcholinesterase, butyrylcholinesterase, and β-secretase (BACE-1) were further assessed using molecular docking, binding free energy calculations, cluster analysis, and molecular dynamics simulations. The GLME and fractions inhibited the activities of both acetylcholinesterase and butyrylcholinesterase in a dose-dependent manner. Iloneoside and marsectohexol exhibited in vitro concentration-dependent inhibitory activities against acetylcholinesterase (IC50 = 19.28, 184.9 µM, respectively) and butyrylcholinesterase (IC50 = 30.75, 43.4 µM, respectively). These compounds also possess ferric ion-reducing, hydroxyl, and superoxide radical-scavenging activities. Iloneoside had the highest docking scores of -9.8, -9.9 -9.4 Kcal for AChE, BChE, and BACE1, respectively. The stability of the interaction of the bioactive compounds with the catalytic residues of the protein targets was preserved in a 100 ns molecular dynamics simulation. Iloneoside, a rare pregnane glycoside, was identified as a neurotherapeutic constituent of G. latifolium leaf. Further studies are suggested to investigate the neurotherapeutic potential in animal models.

Effect of Anchoring Dynamics on Proton-Coupled Electron Transfer in the Ru(bda) Coordination Oligomer on a Graphitic Surface

The oligomeric ruthenium-based water oxidation catalyst, Ru(bda), is known to be experimentally anchored on graphitic surfaces through CH-π stacking interactions between the auxiliary bda ([2,2'-bipyridine]-6,6'-dicarboxylate) ligand bonded to ruthenium and the hexagonal rings of the surface. This anchoring provides control over their molecular coverage and enables efficient catalysis of water oxidation to dioxygen. The oligomeric nature of the molecule offers multiple anchoring sites at the surface, greatly enhancing the overall stability of the hybrid catalyst-graphitic surface anode through dynamic bonding. However, the impact of this dynamic anchoring on the overall catalytic mechanism is still a topic of debate. In this study, a crucial proton-coupled electron transfer event in the catalytic cycle is investigated using DFT-based molecular dynamics simulations plus metadynamics. The CH-π stacking anchoring plays a critical role not only in stabilizing this hybrid system but also in facilitating the proton-coupled electron transfer event with possible vibronic couplings between the anchoring bonds motion and charge fluctuations at the catalyst - graphitic surface interface. Furthermore, this computational investigation displays the presence of a quartet spin state intermediate that can lead to the experimentally observed and thermodynamically more stable doublet spin state.

Elemental evaluation, nutritional analysis, GC-MS analysis and ameliorative effects of Artocarpus communis J.R.Forst. & G.Forst. seeds' phytoconstituents on metabolic syndrome via in silico approach

The nutritional as well as beneficial effects of the Artocarpus communis seed on metabolic syndrome complications have not been studied. In this research, the aim was to investigate the nutritional composition and beneficial effects of Artocarpus communis seeds' phytoconstituents on the p53 core, fat mass and obesity-associated (FTO) protein and cytochrome P450 CYP11A1 domains. The elements and phytochemicals in the seed were determined through atomic absorption spectroscopy assay and gas chromatography-mass spectrometry (GC-MS) analysis, respectively. Also, the compounds detected were docked to the p53 core, FTO protein and cytochrome P450 CYP11A1 domains protein. Artocarpus communis seed contains sodium (7.824 ± 0.0134 ppm), magnesium (10.187 ± 0.0239 ppm) and iron (1.924 ± 0.0017), while zinc and cadmium were undetected. Phenolics and flavonoids were the most abundant phytochemicals in the seed. Phytoconstituents, such as pentadecanoic acid, hexadecanoic acid and methyl ester, possessing different therapeutic effects were identified via GC-MS analysis. In A. communis seed, 3-methyl-4-nitro-5-(1-pyrazolyl) pyrazole and phenanthrene were able to bind more peculiarly and specifically to the p53 core, FTO protein and cytochrome P450 CYP11A1 domains. One of the important processes that were hypothesized for the recovery of metabolic syndrome in affected victims is shown by the molecular dynamics analysis, which shows that the binding of these chemicals to the targeted structure stabilized the proteins. Therefore, Artocarpus communis seeds could be a new strategy for the management of metabolic syndrome.Communicated by Ramaswamy H. Sarma.

Therapeutic Study of Cinnamic Acid Derivative for Oxidative Stress Ablation: The Computational and Experimental Answers

This study aimed to examine the therapeutic activity of the cinnamic acid derivative KAD-7 (N'-(2,4-dichlorobenzylidene)-3-(4-methoxyphenyl) acrylohydrazide) on Fe2+-induced oxidative hepatic injury via experimental and computational models. In addition, the role of ATPase and ectonucleoside triphosphate diphosphohydrolase (ENTPDase) in the coordination of cellular signals is speculated upon to proffer suitable therapeutics for metabolic stress disorder upon their inhibition. While we know little about therapeutics with flexible dual inhibitors for these protein targets, this study was designed to screen KAD-7's (N'-(2,4-dichlorobenzylidene)-3-(4-methoxyphenyl) acrylohydrazide) inhibitory potential for both protein targets. We induced oxidative hepatic damage via the incubation of hepatic tissue supernatant with 0.1 mM FeSO4 for 30 min at 37 °C. We achieved the treatment by incubating the hepatic tissues with KAD-7 under the same conditions. The catalase (CAT), glutathione (GSH), malondialdehyde (MDA), ATPase, and ENTPDase activity were all measured in the tissues. We predicted how the drug candidate would work against ATPase and ENTPDase targets using molecular methods. When hepatic injury was induced, there was a significant decrease in the levels of the GSH, CAT, and ENTPDase (p < 0.05) activities. In contrast, we found a noticeable rise in the MDA levels and ATPase activity. KAD-7 therapy resulted in lower levels of these activities overall (p < 0.05), as compared to the control levels. We found the compound to have a strong affinity for ATPase (-7.1 kcal/mol) and ENTPDase (-7.4 kcal/mol), and a better chemical reactivity than quercetin. It also met all drug-likeness parameters. Our study shows that KAD-7 can protect the liver from damage caused by FeSO4 by reducing oxidative stress and purinergic actions. Our studies indicate that KAD-7 could be developed as a therapeutic option since it can flexibly inhibit both ATPase and ENTPDase.

The Role of Chloride ion in the Silicate Condensation Reaction from ab Initio Molecular Dynamics Simulations

The comprehension of silicate oligomer formation during the initial stage of zeolite synthesis is of significant importance. In this study, we investigated the effect of chloride ions (Cl-) on silicate oligomerization using ab initio molecular dynamics simulations with explicit water molecules. The results show that the presence of Cl- increases the free energy barriers of all reactions compared to the case without the anion. The formation of the 4-ring structure has the lowest free energy barrier (73 kJ/mol), while the formation of the 3-ring structure has the highest barrier (98 kJ/mol) in the presence of Cl-. These findings suggest that Cl- suppresses the formation of 3-rings and favors the formation of larger oligomers in the process of zeolite synthesis. Our study provides important insights into the directing role of Cl- in silicate oligomerization by regulating thermodynamic and kinetic parameters. An important point to consider is the impact of the anion on aqueous reactions, particularly in altering the hydrogen bond network around reactive species. These results also provide a basis for further studies of the formations of larger silicate oligomers in solution.

Insight into the role of excess hydroxide ions in silicate condensation reactions

The formation of silicate oligomers in the early stages is key to zeolite synthesis. The pH and the presence of hydroxide ions are important in regulating the reaction rate and the dominant species in solutions. This paper describes the formation of silicate species, from dimers to 4-membered rings, using ab initio molecular dynamics simulations in explicit water molecules with an excess hydroxide ion. The thermodynamic integration method was used to calculate the free energy profile of the condensation reactions. The hydroxide group's role is not only to control the pH of the environment, but also to actively participate in the condensation reaction. The results show that the most favorable reactions are linear-tetramer and 4-membered-ring formation, with overall barriers of 71 kJ mol^-1 and 73 kJ mol^-1, respectively. The formation of trimeric silicate, with the largest free-energy barrier of 102 kJ mol^-1, is the rate-limiting step under these conditions. The excess hydroxide ion aids in the stabilization of the 4-membered-ring structure over the 3-membered-ring structure. Due to a relatively high free-energy barrier, the 4-membered ring is the most difficult of the small silicate structures to dissolve in the backward reaction. This study is consistent with the experimental observation that silicate growth in zeolite synthesis is slower in a very-high-pH environment.

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.

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.

Modulation of Metal Carbonyl Stretching Frequencies in the Second Coordination Sphere through the Internal Stark Effect

Spectroscopic and computational examination of a homologous series of rhodium(I) pybox carbonyl complexes has revealed a correlation between the conformation of the flanking aryl-substituted oxazoline donors and the carbonyl stretching frequency. This relationship is also observed experimentally for octahedral rhodium(III) and ruthenium(II) variants and cannot be explained through the classical, Dewar-Chatt-Duncanson, interpretation of metal-carbonyl bonding. Instead, these findings are reconciled by local changes in the magnitude of the electric field that is projected along the metal-carbonyl vector: the internal Stark effect.

Water electrolysis: from textbook knowledge to the latest scientific strategies and industrial developments

Replacing fossil fuels with energy sources and carriers that are sustainable, environmentally benign, and affordable is amongst the most pressing challenges for future socio-economic development. To that goal, hydrogen is presumed to be the most promising energy carrier. Electrocatalytic water splitting, if driven by green electricity, would provide hydrogen with minimal CO2 footprint. The viability of water electrolysis still hinges on the availability of durable earth-abundant electrocatalyst materials and the overall process efficiency. This review spans from the fundamentals of electrocatalytically initiated water splitting to the very latest scientific findings from university and institutional research, also covering specifications and special features of the current industrial processes and those processes currently being tested in large-scale applications. Recently developed strategies are described for the optimisation and discovery of active and durable materials for electrodes that ever-increasingly harness first-principles calculations and machine learning. In addition, a technoeconomic analysis of water electrolysis is included that allows an assessment of the extent to which a large-scale implementation of water splitting can help to combat climate change. This review article is intended to cross-pollinate and strengthen efforts from fundamental understanding to technical implementation and to improve the 'junctions' between the field's physical chemists, materials scientists and engineers, as well as stimulate much-needed exchange among these groups on challenges encountered in the different domains.

Computational Scaling Relationships Predict Experimental Activity and Rate-Limiting Behavior in Homogeneous Water Oxidation

While computational screening with first-principles density functional theory (DFT) is essential for evaluating candidate catalysts, limitations in accuracy typically prevent the prediction of experimentally relevant activities. Exemplary of these challenges are homogeneous water oxidation catalysts (WOCs) where differences in experimental conditions or small changes in ligand structure can alter rate constants by over an order of magnitude. Here, we compute mechanistically relevant electronic and energetic properties for 19 mononuclear Ru transition-metal complexes (TMCs) from three experimental water oxidation catalysis studies. We discover that 15 of these TMCs have experimental activities that correlate with a single property, the ionization potential of the Ru(II)-O₂ catalytic intermediate. This scaling parameter allows the quantitative understanding of activity trends and provides insight into the rate-limiting behavior. We use this approach to rationalize differences in activity with different experimental conditions, and we qualitatively analyze the source of distinct behavior for different electronic states in the other four catalysts. Comparison to closely related single-atom catalysts and modified WOCs enables rationalization of the source of rate enhancement in these WOCs.

Beyond Continuum Solvent Models in Computational Homogeneous Catalysis

In homogeneous catalysis solvent is an inherent part of the catalytic system. As such, it must be considered in the computational modeling. The most common approach to include solvent effects in quantum mechanical calculations is by means of continuum solvent models. When they are properly used, average solvent effects are efficiently captured, mainly those related with solvent polarity. However, neglecting atomistic description of solvent molecules has its limitations, and continuum solvent models all alone cannot be applied to whatever situation. In many cases, inclusion of explicit solvent molecules in the quantum mechanical description of the system is mandatory. The purpose of this article is to highlight through selected examples what are the reasons that urge to go beyond the continuum models to the employment of micro-solvated (cluster-continuum) of fully explicit solvent models, in this way setting the limits of continuum solvent models in computational homogeneous catalysis. These examples showcase that inclusion of solvent molecules in the calculation not only can improve the description of already known mechanisms but can yield new mechanistic views of a reaction. With the aim of systematizing the use of explicit solvent models, after discussing the success and limitations of continuum solvent models, issues related with solvent coordination and solvent dynamics, solvent effects in reactions involving small, charged species, as well as reactions in protic solvents and the role of solvent as reagent itself are successively considered.

Efficient workflow for the investigation of the catalytic cycle of water oxidation catalysts: Combining GFN-xTB and density functional theory

Photocatalytic water oxidation remains the bottleneck in many artificial photosynthesis devices. The efficiency of this challenging process is inherently linked to the thermodynamic and electronic properties of the chromophore and the water oxidation catalyst (WOC). Computational investigations can facilitate the search for favorable chromophore‐catalyst combinations. However, this remains a demanding task due to the requirements on the computational method that should be able to correctly describe different spin and oxidation states of the transition metal, the influence of solvation and the different rates of the charge transfer and water oxidation processes. To determine a suitable method with favorable cost/accuracy ratios, the full catalytic cycle of a molecular ruthenium based WOC is investigated using different computational methods, including density functional theory (DFT) with different functionals (GGA, Hybrid, Double Hybrid) as well as the semi‐empirical tight binding approach GFN‐xTB. A workflow with low computational cost is proposed that combines GFN‐xTB and DFT and provides reliable results. GFN‐xTB geometries and frequencies combined with single‐point DFT energies give free energy changes along the catalytic cycle that closely follow the full DFT results and show satisfactory agreement with experiment, while significantly decreasing the computational cost. This workflow allows for cost efficient determination of energetic, thermodynamic and dynamic properties of WOCs.

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.

Comparative first principles-based molecular dynamics study of catalytic mechanism and reaction energetics of water oxidation reaction on 2D-surface

The study of the water-splitting process, which can proceed in 2e^- as well as 4e^- pathway, reveals that the process is entirely an uphill process, and the third step, that is, the oxooxo bond formation is the rate-determining step. The kinetic barrier of the oxygen evolution reaction (OER) on the 2D material catalysts in the presence of explicit solvents is scarcely studied. Here, we investigate the dynamics of the OER on the undoped graphene and the activation energy barrier of each step using first principles molecular dynamics simulations. Here we provide a detailed analysis of the kinetics of all the 4e^- transfer steps of OER on the graphene surface. We also compare the accuracy of one of the density functional theory (DFT) functionals and density functional based tight binding (DFTB) method in explaining the OER steps. The comparative study reveals that DFTB can be used for performing metadynamics simulations quipped with much less computational cost than DFT functionals. By both Perdew-Burke-Ernzerhof and DFTB methods, the third step is revealed to be the rate-determining step with an energy barrier of 21.19 ± 0.51 and 20.23 ± 0.20 kcal mol^-1 , respectively. DFTB gives an impression of being successful in predicting the energy barriers of OER in 4e- transfer pathway and comparable to the DFT method, and we would like to extend the use of DFTB for further studies with a sizable and complex system.

Mechanism and Dynamics of Formation of Bisoxo Intermediates and O-O Bond in the Catalytic Water Oxidation Process

This work elucidates the reactivity of water molecules toward the tridentate nitrogen-containing iron complex in the water oxidation process. Here, we consider the Fe^V-bisoxo complex {[Fe^V(Me₃tacn)(OH₂)(═O)₂]⁺} to be responsible for the oxygen-oxygen bond formation. This O-O bond formation happens through the addition of water as a nucleophile. The transition state was determined by the synchronous transit-guided quasi-Newton method using reactants and products and verified by intrinsic reaction coordinates (IRCs). From the IRC calculations, we observe that the Fe^V═O moiety is attacked by water and assisted by the H-bonded interaction with the oxygen atom of the bisoxo complex. The hydrogen atom is transferred to the oxygen atom of the bisoxo complex through the transition state, and subsequently, the hydroxide is transferred to another oxygen of the bisoxo complex, resulting in the formation of the oxygen-oxygen bond. This work also explains the effect of explicit water molecules on the oxygen-oxygen bond formation. Our results also show how the formation of superoxide plays an essential role in O₂ evolution. We used the potential energy scan method to compute the transition state in the oxygen evolution step. In the present work, we study the effect of chlorine on the formation of the oxygen-oxygen bond formation. In this study, the changes in the oxidation state, spin density, and spin multiplicity of the complexes are investigated for each successive step. Apart from these static theoretical calculations, we also studied the oxygen-oxygen bond formation through first-principles molecular dynamics with the aid of the well-tempered metadynamics sampling technique. From the observation of the free energy surfaces from metadynamics simulations, it is evident that the hydroxide transfer has a lesser free energetic reaction as compared to the proton transfer. This complete mechanistic study may give an idea to design a suitable water oxidation catalyst.

Elucidating the Role of Tetraethylammonium in the Silicate Condensation Reaction from Ab Initio Molecular Dynamics Simulations

The understanding of the formation of silicate oligomers in the initial stage of zeolite synthesis is important. The use of organic structure-directing agents (OSDAs) is known to be a key factor in the formation of different silicate species and the final zeolite structure. For example, tetraethylammonium ion (TEA⁺) is a commonly used organic template for zeolite synthesis. In this study, ab initio molecular dynamics (AIMD) simulation is used to provide an understanding of the role of TEA⁺ in the formation of various silicate oligomers, ranging from dimer to 4-ring. Calculated free-energy profiles of the reaction pathways show that the formation of a 4-ring structure has the highest energy barrier (97 kJ/mol). The formation of smaller oligomers such as dimer, trimer, and 3-ring has lower activation barriers. The TEA⁺ ion plays an important role in regulating the predominant species in solution via its coordination with silicate structures during the condensation process. The kinetics and thermodynamics of the oligomerization reaction indicate a more favorable formation of the 3-ring over the 4-ring structure. The results from AIMD simulations are in line with the experimental observation that TEA⁺ favors the 3-ring and double 3-ring in solution. The results of this study imply that the role of OSDAs is not only important for the host-guest interaction but also crucial for controlling the reactivity of different silicate oligomers during the initial stage of zeolite formation.

In need of a second-hand? The second coordination sphere of ruthenium complexes enables water oxidation with improved catalytic activity

A ligand dangling arm, acting as an intramolecular proton acceptor, drastically increasing the catalytic activity of Ru-complexes for water oxidation.

Proton Acceptor near the Active Site Lowers Dramatically the O-O Bond Formation Energy Barrier in Photocatalytic Water Splitting

The O-O bond formation process via water nucleophilic attack represents a thermodynamic and kinetic bottleneck in photocatalytic water oxidation because of the considerably high activation free energy barrier. It is therefore of fundamental significance and yet challenging to find strategies to facilitate this reaction. The microscopic details of the photocatalytic water oxidation step involving the O-O bond formation in a catalyst-dye supramolecular complex are here elucidated by density functional theory-based Car-Parrinello molecular dynamics simulations in the presence of an extra proton acceptor. Introducing a proton acceptor group (OH-) in the hydration shell near the catalytic active site accelerates the rate-limiting O-O bond formation by inducing a cooperative event proceeding via a concerted proton-coupled electron-transfer mechanism and thus significantly lowering the activation free energy barrier. The in-depth insight provides a strategy for facilitating the photocatalytic water oxidation and for improving the efficiency of dye-sensitized photoelectrochemical cells.

Universal scaling relations for the rational design of molecular water oxidation catalysts with near-zero overpotential

A major roadblock in realizing large-scale production of hydrogen via electrochemical water splitting is the cost and inefficiency of current catalysts for the oxygen evolution reaction (OER). Computational research has driven important developments in understanding and designing heterogeneous OER catalysts using linear scaling relationships derived from computed binding energies. Herein, we interrogate 17 of the most active molecular OER catalysts, based on different transition metals (Ru, Mn, Fe, Co, Ni, and Cu), and show they obey similar scaling relations to those established for heterogeneous systems. However, we find that the conventional OER descriptor underestimates the activity for very active OER complexes as the standard approach neglects a crucial one-electron oxidation that many molecular catalysts undergo prior to O-O bond formation. Importantly, this additional step allows certain molecular catalysts to circumvent the "overpotential wall", leading to enhanced performance. With this knowledge, we establish fundamental principles for the design of ideal molecular OER catalysts.

Modeling the Catalyst Activation Step in a Metal–Ligand Radical Mechanism Based Water Oxidation System

Designing catalysts for water oxidation (WOCs) that operate at low overpotentials plays an important role in developing sustainable energy conversion schemes. Recently, a mononuclear ruthenium WOC that operates via metal–ligand radical coupling pathway was reported, with a very low barrier for O–O bond formation, that is usually the rate-determining step in most WOCs. A detailed mechanistic understanding of this mechanism is crucial to design highly active oxygen evolution catalysts. Here, we use density functional theory based molecular dynamics (DFT-MD) with an explicit description of the solvent to investigate the catalyst activation step for the [Ru(bpy) 2 (bpy–NO)] 2 + complex, that is considered to be the rate-limiting step in the metal–ligand radical coupling pathway. We find that a realistic description of the solvent environment, including explicit solvent molecules and thermal motion, is crucial for an accurate description of the catalyst activation step, and for the estimation of the activation barriers.

High-performance water oxidation catalysts based on the spontaneous deposition of ruthenium on electrochemically exfoliated graphene oxide

Sustainable highly active water oxidation reaction in acid over ruthenium loaded electrochemically exfoliated graphene oxide.