Are DFT Methods Accurate in Mononuclear Ruthenium-Catalyzed Water Oxidation? An ab Initio Assessment

Coupled Cluster Benchmark of New DFT and Local Correlation Methods: Mechanisms of Hydroarylation and Oxidative Coupling Catalyzed by Ru(II, III) Chloride Carbonyls

We have evaluated a set of accurate canonical CCSD(T) energies for stationary points on the potential energy surface for Ru(II, III) chloride carbonyl catalysis of two competing reactions between benzene and methyl acrylate (MA), namely, hydroarylation and oxidative coupling. We have then applied this set to evaluate the performance of localized orbital coupled-cluster methods and several new and common density functionals. We find that (a) DLPNO-CCSD(T) with TightPNO cutoffs is an acceptable substitute for full canonical CCSD(T) calculations on this system; (b) for the closed-shell systems where it could be applied, LNO-CCSD(T) with tight convergence criteria is very close to the canonical results; (c) the recent ωB97X-V and ωB97M-V functionals exhibit superior performance to commonly used DFT functionals in both closed- and open-shell calculations; (d) the revDSD-PBEP86 revision of the DSD-PBEP86 double hybrid represents an improvement over the original, even though transition metals were not involved in its parametrization; and (e) DSD-SCAN and DOD-SCAN show comparable efficiency. Most tested (meta)-GGA and hybrid density functionals perform better for open-shell than for closed-shell complexes; this is not the case for the double hybrids considered.

Isolation of a Ru(IV) side-on peroxo intermediate in the water oxidation reaction

The electrons that nature uses to reduce CO₂ during photosynthesis come from water oxidation at the oxygen-evolving complex of photosystem II. Molecular catalysts have served as models to understand its mechanism, in particular the O-O bond-forming reaction, which is still not fully understood. Here we report a Ru(IV) side-on peroxo complex that serves as a 'missing link' for the species that form after the rate-determining O-O bond-forming step. The Ru(IV) side-on peroxo complex (η²-1^iv-OO) is generated from the isolated Ru(IV) oxo complex (1^iv=O) in the presence of an excess of oxidant. The oxidation (IV) and spin state (singlet) of η²-1^iv-OO were determined by a combination of experimental and theoretical studies. ¹⁸O- and ²H-labelling studies evidence the direct evolution of O₂ through the nucleophilic attack of a H₂O molecule on the highly electrophilic metal-oxo species via the formation of η²-1^iv-OO. These studies demonstrate water nucleophilic attack as a viable mechanism for O-O bond formation, as previously proposed based on indirect evidence.

Ruthenium Water Oxidation Catalysts based on Pentapyridyl Ligands

Ruthenium complexes containing the pentapyridyl ligand 6,6''-(methoxy(pyridin-2-yl)methylene)di-2,2'-bipyridine (L-OMe) of general formula trans-[Ru^II (X)(L-OMe-κ-N⁵ )]ⁿ⁺ (X=Cl, n=1, trans-1⁺ ; X=H₂ O, n=2, trans-2²⁺ ) have been isolated and characterized in solution (by NMR and UV/Vis spectroscopy) and in the solid state by XRD. Both complexes undergo a series of substitution reactions at oxidation state Ru^II and Ru^III when dissolved in aqueous triflic acid-trifluoroethanol solutions as monitored by UV/Vis spectroscopy, and the corresponding rate constants were determined. In particular, aqueous solutions of the Ru^III -Cl complex trans-[Ru^III (Cl)(L-OMe-κ-N⁵ )]²⁺ (trans-1²⁺ ) generates a family of Ru aquo complexes, namely trans-[Ru^III (H₂ O)(L-OMe-κ-N⁵ )]³⁺ (trans-2³⁺ ), [Ru^III (H₂ O)₂ (L-OMe-κ-N⁴ )]³⁺ (trans-3³⁺ ), and [Ru^III (Cl)(H₂ O)(L-OMe-κ-N⁴ )]²⁺ (trans-4²⁺ ). Although complex trans-4²⁺ is a powerful water oxidation catalyst, complex trans-2³⁺ has only a moderate activity and trans-3³⁺ shows no activity. A parallel study with related complexes containing the methyl-substituted ligand 6,6''-(1-pyridin-2-yl)ethane-1,1-diyl)di-2,2'-bipyridine (L-Me) was carried out. The behavior of all of these catalysts has been rationalized based on substitution kinetics, oxygen evolution kinetics, electrochemical properties, and density functional theory calculations. The best catalyst, trans-4²⁺ , reaches turnover frequencies of 0.71 s^-1 using Ce^IV as a sacrificial oxidant, with oxidative efficiencies above 95 %.

Modeling σ-Bond Activations by Nickel(0) Beyond Common Approximations: How Accurately Can We Describe Closed-Shell Oxidative Addition Reactions Mediated by Low-Valent Late 3d Transition Metal?

Accurate modelings of reactions involving 3d transition metals (TMs) are very challenging to both ab initio and DFT approaches. To gain more knowledge in this field, we herein explored typical σ-bond activations of H-H, C-H, C-Cl, and C-C bonds promoted by nickel(0), a low-valent late 3d TM. For the key parameters of activation energy (ΔE^‡) and reaction energy (ΔE_R) for these reactions, various issues related to the computational accuracy were systematically investigated. From the scrutiny of convergence issue with one-electron basis set, augmented (A) basis functions are found to be important, and the CCSD(T)/CBS level with complete basis set (CBS) limit extrapolation based on augmented double-ζ and triple-ζ basis pair (ADZ and ATZ), which produces deviations below 1 kcal/mol from the reference, is recommended for larger systems. As an alternative, the explicitly correlated F12 method can accelerate the basis set convergence further, especially after its CBS extrapolations. Thus, the CCSD(T)-F12/CBS(ADZ-ATZ) level with computational cost comparable to the conventional CCSD(T)/CBS(ADZ-ATZ) level, is found to reach the accuracy of the conventional CCSD(T)/A5Z level, which produces deviations below 0.5 kcal/mol from the reference, and is also highly recommendable. Scalar relativistic effects and 3s3p core-valence correlation are non-negligible for achieving chemical accuracy of around 1 kcal/mol. From the scrutiny of convergence issue with the N-electron basis set, in comparison with the reference CCSDTQ result, CCSD(T) is found to be able to calculate ΔE^‡ quite accurately, which is not true for the ΔE_R calculations. Using highest-level CCSD(T) results of ΔE^‡ in this work as references, we tested 18 DFT methods and found that PBE0 and CAM-B3LYP are among the three best performing functionals, irrespective of DFT empirical dispersion correction. With empirical dispersion correction included, ωB97XD is also recommendable due to its improved performance.

DFT Study on the Interaction of Tris(benzene-1,2-dithiolato)molybdenum Complex with Water. A Hydrolysis Mechanism Involving a Feasible Seven-Coordinate Aquomolybdenum Intermediate

In the present work, the reactivity of the tris(benzene-1,2-dithiolato)molybdenum complex ([Mo(bdt)₃]) toward water is studied by means of the density functional theory (DFT). DFT calculations were performed using the M06, B3P86, and B3PW91 hybrid functionals for comparison purposes. The M06 method was employed to elucidate the reaction pathway, relative stability of the intermediate products, nature of the Mo-S bond cleavage, and electronic structure of the involved molybdenum species. This functional was also used to study the transference of electrons from the molybdenum center toward the ligands. The reaction pathway confirms that [Mo(bdt)₃] undergoes hydrolysis, yielding dihydroxo-bis(benzene-1,2-dithiolato)molybdenum complex ([Mo(OH)₂(bdt)₂]) and benzenedithiol. The reaction takes place through seven transition structures, one of them involving an aquo seven-coordinate molybdenum intermediate stabilized by a lone pair (LP) LP_O→LP_Mo hyperconjugative interaction. This heptacoordinate species allows understanding of the observed oxygen atom exchange between water and tertiary phosphines mediated by these complexes. Calculations also show that [Mo(C₂H₄S₂)₃] and [Mo(OH)₂(C₂H₄S₂)₂] have d² and d⁰ electronic configuration, and hence an electron pair must be transferred during the course of the hydrolysis. The frontier molecular orbital (FMO) analysis concludes that the electron pair is transferred in the rupture of the second Mo-S bond, from the occupied donating Mo d_x²-y² orbital to the unoccupied C₂H₄(SH)₂ S-C σ* ligand orbital. This result is supported by the bond dissociation energy calculations, which demonstrate that the neutral dissociation of the second Mo-S bond is energetically the more favorable.

A mononuclear non-heme-iron dioxygen-carrying protein?

The ability of mononuclear non-heme iron complexes to function as molecular oxygen transporters is investigated by density functional theory. The factors governing the efficiency of the reversible binding of dioxygen at the active site of the dinuclear non-heme iron enzyme hemerythrin, including antiferromagnetic coupling and the conversion of dioxygen to hydroperoxo by a proton coupled 2-electron transfer mechanism, are revisited and considered as possible tools in mononuclear non-heme complexes. Several mononuclear non-heme model complexes, including active sites of enzymes already known to interact with dioxgenic ligands, are constructed and the molecular oxygen transportation capabilities of these complexes are examined computationally. The high-spin nature of the ground state of these complexes implies an intrinsic kinetic lability of the oxy structures, as also evident from potential energy surface calculations towards iron-dioxygen cleavage. Proton affinities as calibrated with reference compounds showed that these complexes are highly unlikely to undergo protonation to form hydroperoxo-like adducts. Mixed superoxo descriptions of the dissociated dioxygenic ligands in all complexes add to the overall conclusion that these model structures are significantly disadvantaged in any attempt to be employed for molecular oxygen transportation.

Spectroscopic, Electrochemical and Computational Characterisation of Ru Species Involved in Catalytic Water Oxidation: Evidence for a [Ru(V) (O)(Py2 (Me) tacn)] Intermediate

A new family of ruthenium complexes based on the N-pentadentate ligand Py2 (Me) tacn (N-methyl-N',N''-bis(2-picolyl)-1,4,7-triazacyclononane) has been synthesised and its catalytic activity has been studied in the water-oxidation (WO) reaction. We have used chemical oxidants (ceric ammonium nitrate and NaIO4 ) to generate the WO intermediates [Ru(II) (OH2 )(Py2 (Me) tacn)](2+) , [Ru(III) (OH2 )(Py2 (Me) tacn)](3+) , [Ru(III) (OH)(Py2 (Me) tacn)](2+) and [Ru(IV) (O)(Py2 (Me) tacn)](2+) , which have been characterised spectroscopically. Their relative redox and pH stability in water has been studied by using UV/Vis and NMR spectroscopies, HRMS and spectroelectrochemistry. [Ru(IV) (O)(Py2 (Me) tacn)](2+) has a long half-life (>48 h) in water. The catalytic cycle of WO has been elucidated by using kinetic, spectroscopic, (18) O-labelling and theoretical studies, and the conclusion is that the rate-determining step is a single-site water nucleophilic attack on a metal-oxo species. Moreover, [Ru(IV) (O)(Py2 (Me) tacn)](2+) is proposed to be the resting state under catalytic conditions. By monitoring Ce(IV) consumption, we found that the O2 evolution rate is redox-controlled and independent of the initial concentration of Ce(IV) . Based on these facts, we propose herein that [Ru(IV) (O)(Py2 (Me) tacn)](2+) is oxidised to [Ru(V) (O)(Py2 (Me) tacn)](2+) prior to attack by a water molecule to give [Ru(III) (OOH)(Py2 (Me) tacn)](2+) . Finally, it is shown that the difference in WO reactivity between the homologous iron and ruthenium [M(OH2 )(Py2 (Me) tacn)](2+) (M=Ru, Fe) complexes is due to the difference in the redox stability of the key M(V) (O) intermediate. These results contribute to a better understanding of the WO mechanism and the differences between iron and ruthenium complexes in WO reactions.

Comparative Assessment of DFT Performances in Ru- and Rh-Promoted σ-Bond Activations

In this work, the performances of 19 density functional theory (DFT) methods are calibrated comparatively on Ru- and Rh-promoted σ-bond (C-H, O-H, and H-H) activations. DFT calibration reference is generated from explicitly correlated coupled cluster CCSD(T)-F12 calculations, and the 4s4p core-valence correlation effect of the two 4d platinum group transition metals is also included. Generally, the errors of DFT methods for calculating energetics of Ru-/Rh-mediated reactions appear to correlate more with the magnitude of energetics itself than other factors such as metal identity. For activation energy calculations, the best performing functionals for both Ru and Rh systems are MN12SX < CAM-B3LYP < M06-L < MN12L < M06 < ωB97X < B3LYP < LC-ωPBE (in the order of increasing mean unsigned deviations, MUDs, of less than 2 kcal/mol). For reaction energy calculations, best functionals with MUDs less than 2 kcal/mol are PBE0 < CAM-B3LYP ≈ N12SX. The effect of the DFT empirical dispersion correction on the performance of the DFT methods is beneficial for most density functionals tested in this work, reducing their MUDs to different extents. After including empirical dispersion correction, ωB97XD, B3LYP-D3, and CAM-B3LYP-D3 (PBE0-D3, B3LYP-D3, and ωB97XD) are the three best performing DFs for activation energy (reaction energy) calculations, from which B3LYP-D3 and ωB97XD can notably be recommended uniformly for both the reaction energy and reaction barrier calculations. The good performance of B3LYP-D3 in quantitative description of the energetic trends further adds value to B3LYP-D3 and singles this functional out as a reasonable choice in the Ru/Rh-promoted σ-bond activation processes.

An ab Initio Benchmark and DFT Validation Study on Gold(I)-Catalyzed Hydroamination of Alkynes

High level ab initio calculations have been carried out on an archetypal gold(I)-catalyzed reaction: hydroamination of ethyne. We studied up to 12 structures of possible gold(I)-coordinated species modeling different intermediates potentially present in a catalytic cycle for the addition of a protic nucleophile to an alkyne. The benchmark is used to evaluate the performances of some popular density functionals for describing geometries and relative energies of stationary points along the reaction profile. Most functionals (including hybrid or meta-hybrid) give accurate structures but large nonsystematic errors (4-12 kcal/mol) along the reaction energy profile. The double hybrid functional B2PLYP outperforms all considered functionals and compares very nicely with our reference ab initio benchmark energies. Moreover, we present an assessment of the accuracy of commonly used approaches to include relativistic effects, such as relativistic effective potentials and a scalar ZORA Hamiltonian, by a comparison with the results obtained using a relativistic all-electron four-component Dirac-Kohn-Sham method. The contribution of nonscalar relativistic effects in gold(I)-catalyzed reactions, as we investigated here, is expected to be on the order of 1 kcal/mol.

Mechanistic analysis of water oxidation catalyzed by mononuclear copper in aqueous bicarbonate solutions

We characterize a mechanism for a monomeric copper catalyst reported to oxidize water in bicarbonate solution when subject to sufficiently high external potentials at near neutral pH values.