Electronic Structure and Formation of Simple Ferryloxo Complexes: Mechanism of the Fenton Reaction

Approaching the Complete Basis Set Limit for Spin-State Energetics of Mononuclear First-Row Transition Metal Complexes

Convergence to the complete basis set (CBS) limit is analyzed for the problem of spin-state energetics in mononuclear first-row transition metal (TM) complexes by taking under scrutiny a benchmark set of 18 energy differences between spin states for 13 chemically diverse TM complexes. The performance of conventional CCSD(T) and explicitly correlated CCSD(T)-F12a/b calculations in approaching the CCSD(T)/CBS limits is systematically studied. An economic computational protocol is developed based on the CCSD-F12a approximation and (here proposed) modified scaling of the perturbative triples term (T#). This computational protocol recovers the relative spin-state energetics of the benchmark set in excellent agreement with the reference CCSD(T)/CBS limits (mean absolute deviation of 0.4, mean signed deviation of 0.2, and maximum deviation of 0.8 kcal/mol) and enables performing canonical CCSD(T) calculations for mononuclear TM complexes sized up to ca. 50 atoms, which is illustrated by application to heme-related metalloporphyrins. Furthermore, a good transferability of the basis set incompleteness error (BSIE) is demonstrated for spin-state energetics computed using CCSD(T) and other wave function methods (MP2, CASPT2, CASPT2/CC, NEVPT2, and MRCI + Q), which justifies efficient focal-point approximations and simplifies the construction of multimethod benchmark studies.

On The Nature of FeIV =Oaq in Aqueous Media: A DFT analysis

FeIV =Oaq is a key intermediate in many advanced oxidation processes and probably in biological systems. It is usually referred to as FeIV =O2+ . The pKa's of FeIV =Oaq as derived by DFT are: pKa1=2.37 M06 L/6-311++G(d,p) (SDD for Fe) and pKa2=7.79 M06 L/6-311++G(d,p) (SDD for Fe). This means that in neutral solutions, FeIV =Oaq is a mixture of (H2 O)4 (OH)FeIV =O+ and (H2 O)2 (OH)2 FeIV =O. The oxidation potential of FeIV =Oaq in an acidic solution, E0 {(H2 O)5 FeIV =O2+ /FeIII (H2 O)6 3+ , pH 0.0} is calculated with and without a second solvation sphere and the recommended value is between 2.86 V (B3LYP/Def2-TZVP, with a second solvation sphere) and 2.23 V (M06 L/Def2-TZVP without a second solvation sphere). This means that FeIV =Oaq is the strongest oxidizing agent formed in systems involving FeVI O4 2- even in neutral media.

Comprehensive Theoretical View of the [Cu2O2] Side-on-Peroxo-/Bis-μ-Oxo Equilibria

Coupled binuclear copper (CBC) sites are employed by many metalloenzymes to catalyze a broad set of biochemical transformations. Typically, the CBC catalytic sites are activated by the O2 molecule to form various [Cu2O2] reactive species. This has also inspired synthesis and development of various biomimetic inorganic complexes featuring the CBC core. From theoretical perspective, the [Cu2O2] reactivity often hinges on the side-on-peroxo-dicopper(II) (P) vs. bis-μ-oxo-dicopper(III) (O) isomerism - an equilibrium that has become almost iconic in theoretical bioinorganic chemistry. Herein, we present a comprehensive calibration and evaluation of the performance of various composite computational protocols available in contemporary computational chemistry, involving coupled-cluster and multi-reference (relativistic) wave function methods, popular density functionals and solvation models. Starting with the well-studied reference [Cu2O2(NH3)6]2+ system, we compared the performance of electronic structure methods and discussed the relativistic effects. This allowed us to select several 'calibrated' DFT functionals that can be conveniently employed to study ten experimentally well-characterized [Cu2O2] inorganic systems. We mostly correctly predicted the lowest-energy structures (P vs. O) of the studied systems. In addition, we present calibration of the used electronic structure methods for prediction of the spectroscopic features of the [Cu2O2] core.

Heterometal-Organic Cages as Photo-Fenton-like Catalysts

Metal-organic cages, a class of supramolecular containers constructed by the self-assembly of metal ions and organic ligands, show great promise as catalytic agents. In this work, we designed and synthesized a series of rhombic dodecahedral Ni-Cu heterometal imidazolate cages (Ni₈Cu₆L₂₄) that can act as highly active photo-Fenton-like catalysts. These cages possess a high ability to generate hydroxyl radicals (^•OH) under visible light in the presence of H₂O₂, which can rapidly degrade organic pollutants (e.g., rhodamine B, methylene blue, and methyl orange) into CO₂ and H₂O. Besides, they are robust catalysts, with high catalytic activity and reusability under conditions in high H₂O₂ concentration, providing potentially advanced materials for degrading persistent organic pollutants.

DFT mechanistic study on the formation of 8-oxoguanine and spiroiminodihydantoin mediated by iron Fenton reactions

Fenton reactions unavoidably take place in the human body and have been demonstrated to cause oxidative DNA damage. However, the molecular-level understanding of DNA damage mediated by Fenton reactions is limited. Herein, density functional theory (DFT) calculations were made to investigate the counterion effects on aqueous Fenton reactions and the detailed mechanisms of chemical modifications to guanine induced by Fenton reactions. Our calculations show that the activation energy of the Fenton reaction catalyzed by a pure aquo complex [FeII(H2O)6]2+ is too high to agree with experiments, whereas complexation with counteranions reduces the activation energy to a reasonable range. This result suggests that FeII-counteranion complexes are the real catalyst for fast aqueous Fenton reactions. In addition, we found that the Fenton oxidation mediated by FeII bonded to the N7 atom of guanine can result in the formation of 8-oxoguanine and spiroiminodihydantoin through multiple reaction pathways, including the electrophilic addition of ˙OH, H-abstraction by ˙OH, and oxygen atom transfer of oxoiron(iv) species. The activation of hydrogen peroxide by ferrous iron is the rate-determining step. The guanine N7-bound iron ion and the coordinated counteranion were found to play an important role in the Fenton oxidation of guanine.

Overlooked enhancement of chloride ion on the transformation of reactive species in peroxymonosulfate/Fe(II)/NH2OH system

Though hydroxylamine (NH₂OH) is effective for accelerating pollutants degradation in Fenton and Fenton-like systems, the effect of anions simultaneously introduced by the hydroxylamine salts have always been ignored. Herein, effect of two commonly used hydroxylamine salts, hydroxylamine hydrochloride (NH₂OH·HCl) and hydroxylamine sulfate [(NH₂OH)₂·H₂SO₄], for the degradation of dimethyl phthalate (DMP) in peroxymonosulfate (PMS)/Fe(II) system was comparatively investigated. Degradation efficiency of DMP with NH₂OH·HCl was 1.6 times of that with same dosages of (NH₂OH)₂·H₂SO₄. SO₄^·-, Fe(IV) and ·OH formed in the PMS/Fe(II)/NH₂OH system, but ·OH was the major species for DMP degradation. Addition of Cl^- significantly improved the production of ·OH and Cl·, and the exposure dose of ·OH (CT_·OH) was more than 10 times that of CT_Cl· as the concentration of Cl^- increased to 1 mM. Calculations based on branching ratios of Cl· and ·OH indicated that the reactions of Cl^- with SO₄^·- and Cl· with H₂O were not the only production sources of ·OH in the system. Further experiments with methyl phenyl sulfoxide (PMSO) as the probe indicated that Cl^- would facilitate the shift of reactive species from Fe(IV) to radicals (SO₄^·- or ·OH) in the system. Both hydroxylation and nitration intermediate products were detected in the oxidation of DMP. Cl^- promoted the formation of hydroxylation intermediates and reduced the formation of nitration intermediates. This study revealed for the first time that Cl^- could shift reactive species from Fe(IV) to radicals in PMS/Fe(II) system, raising attention to the influence of the coexisting anions (especially Cl^-) for pollutants oxidation in iron-related oxidation processes.

Reaction pathways and deactivation mechanisms of isostructural Cr and Fe MIL-101 during liquid-phase styrene oxidation by hydrogen peroxide

Isostructural MIL-101(Cr, Fe) is investigated as a modular platform to quantify differences in reactivity, selectivity, and deactivation as functions of intrinsic material properties for styrene oxidation by hydrogen peroxide at mild conditions.

Spin-state energetics of metallocenes: How do best wave function and density functional theory results compare with the experimental data?

We benchmark the accuracy of quantum-chemical methods, including wave function theory methods [coupled cluster theory at the CCSD(T) level, multiconfigurational perturbation-theory (CASPT2, NEVPT2) and internally contracted multireference configuration interaction (MRCI)] and 30 density functional theory (DFT) approximations, in reproducing the spin-state splittings of metallocenes. The reference values of the electronic energy differences are derived from the experimental spin-crossover enthalpy for manganocene and the spectral data of singlet-triplet transitions for ruthenocene, ferrocene, and cobaltocenium. For ferrocene and cobaltocenium we revise the previous experimental interpretations regarding the lowest triplet energy; our argument is based on the comparison with the lowest singlet excitation energy and herein reported, carefully determined absorption spectrum of ferrocene. When deriving vertical energies from the experimental band maxima, we go beyond the routine vertical energy approximation by introducing vibronic corrections based on simulated vibrational envelopes. The benchmarking result confirms the high accuracy of the CCSD(T) method (in particular, for UCCSD(T) based on Hartree-Fock orbitals we find for our dataset: maximum error 0.12 eV, weighted mean absolute error 0.07 eV, weighted mean signed error 0.01 eV). The high accuracy of the single-reference method is corroborated by the analysis of a multiconfigurational character of the complete active space wave function for the triplet state of ferrocene. On the DFT side, our results confirm the non-universality problem with approximate functionals. The present study is an important step toward establishing an extensive and representative benchmark set of experiment-derived spin-state energetics for transition metal complexes.

Mechanistic understanding of methane-to-methanol conversion on graphene-stabilized single-atom iron centers

DFT calculations and kinetic modeling elucidate solvent effects and complex mechanisms for the room-temperature methane-to-methanol conversion on an FeN4/graphene catalyst.

Metal Complexes Containing Redox-Active Ligands in Oxidation of Hydrocarbons and Alcohols: A Review

Ligands are innocent when they allow oxidation states of the central atoms to be defined. A noninnocent (or redox) ligand is a ligand in a metal complex where the oxidation state is not clear. Dioxygen can be a noninnocent species, since it exists in two oxidation states, i.e., superoxide (O2−) and peroxide (O22−). This review is devoted to oxidations of C–H compounds (saturated and aromatic hydrocarbons) and alcohols with peroxides (hydrogen peroxide, tert-butyl hydroperoxide) catalyzed by complexes of transition and nontransition metals containing innocent and noninnocent ligands. In many cases, the oxidation is induced by hydroxyl radicals. The mechanisms of the formation of hydroxyl radicals from H2O2 under the action of transition (iron, copper, vanadium, rhenium, etc.) and nontransition (aluminum, gallium, bismuth, etc.) metal ions are discussed. It has been demonstrated that the participation of the second hydrogen peroxide molecule leads to the rapture of O–O bond, and, as a result, to the facilitation of hydroxyl radical generation. The oxidation of alkanes induced by hydroxyl radicals leads to the formation of relatively unstable alkyl hydroperoxides. The data on regioselectivity in alkane oxidation allowed us to identify an oxidizing species generated in the decomposition of hydrogen peroxide: (hydroxyl radical or another species). The values of the ratio-of-rate constants of the interaction between an oxidizing species and solvent acetonitrile or alkane gives either the kinetic support for the nature of the oxidizing species or establishes the mechanism of the induction of oxidation catalyzed by a concrete compound. In the case of a bulky catalyst molecule, the ratio of hydroxyl radical attack rates upon the acetonitrile molecule and alkane becomes higher. This can be expanded if we assume that the reactions of hydroxyl radicals occur in a cavity inside a voluminous catalyst molecule, where the ratio of the local concentrations of acetonitrile and alkane is higher than in the whole reaction volume. The works of the authors of this review in this field are described in more detail herein.

Why the Reactive Oxygen Species of the Fenton Reaction Switches from Oxoiron(IV) Species to Hydroxyl Radical in Phosphate Buffer Solutions? A Computational Rationale

It has been shown that the major reactive oxygen species (ROS) generated by the aqueous reaction of Fe(II) and H2O2 (i.e., the Fenton reaction) are high-valent oxoiron(IV) species, whereas the hydroxyl radical plays a role only in very acidic conditions. Nevertheless, when the Fenton reaction is conducted in phosphate buffer solutions, the resulting ROS turns into hydroxyl radical even in neutral pH conditions. The present density functional theory (DFT) study discloses the underlying principle for this phenomenon. Static and dynamic DFT calculations indicate that in phosphate buffer solutions, the iron ion is highly coordinated by phosphoric acid anions. Such a coordination environment substantially raises the pK a of coordinated water on Fe(III). As a consequence, the Fe(III)-OH intermediate, resulting from the reductive decomposition of H2O2 by ferrous ion is relatively unstable and will be readily protonated by phosphoric acid ligand or by free proton in solution. These proton-transfer reactions, which become energetically favorable when the number of phosphate coordination goes up to three, prevent the Fe(III)-OH from hydrogen abstraction by nascent •OH to form Fe(IV)=O species. On the basis of this finding, a ligand design strategy toward controlling the nature of ROS produced in the Fenton reaction is put forth. In addition, it is found that while phosphate buffers facilitate •OH radical generation in the Fenton reaction, phosphoric acid anions can act as •OH radical scavengers through hydrogen atom transfer reactions.

Advancing Fenton and photo-Fenton water treatment through the catalyst design

The review is devoted to modern Fenton, photo-Fenton, as well as Fenton-like and photo-Fenton-like reactions with participation of iron species in liquid phase and as heterogeneous catalysts. Mechanisms of these reactions were considered that include hydroxyl radical and oxoferryl species as the reactive intermediates. The barriers in the way of application of these reactions to wastewater treatment were discussed. The following fundamental problems need further research efforts: inclusion of more mechanism steps and quantum calculations of all rate constants lacking in the literature, checking the outer sphere electron transfer contribution, determination of the causes for the key changes in the homogeneous Fenton reaction mechanism with a change in the reagents concentration. The key advances for Fenton reactions implementation for the water treatment are related to tremendous hydrodynamical effects on the catalytic activity, design of ligands for high rate and completeness of mineralization in short time, and design of highly active heterogeneous catalysts. While both homogeneous and heterogeneous Fenton and photo-Fenton systems are open for further improvements, heterogeneous photo-Fenton systems are most promising for practical applications because of the inherent higher catalyst stability. Modern methods of quantum chemistry are expected to play a continuously increasing role in development of such catalysts.

Benchmarking quantum chemistry methods for spin-state energetics of iron complexes against quantitative experimental data

The accuracy of relative spin-state energetics predicted by selected quantum chemistry methods: coupled cluster theory at the CCSD(T) level, multiconfigurational perturbation theory (CASPT2, NEVPT2), multireference configuration interaction at the MRCISD+Q level, and a number of DFT methods, is quantitatively evaluated by comparison with the experimental data of four octahedral iron complexes. The available experimental data, either spin-forbidden transition energies or spin crossover enthalpies, are corrected for relevant environmental effects in order to derive the quantitative benchmark set of iron spin-state energetics. Comparison of theory predictions with the resulting reference data: (1) validates the high accuracy of the CCSD(T) method, particularly when based on Kohn-Sham orbitals, giving the maximum error below 2 kcal mol-1 and the mean absolute error (MAE) below 1 kcal mol-1; (2) corroborates the tendency of CASPT2 to systematically overstabilize higher-spin states by up to 5.5 kcal mol-1; (3) confirms that the latter problem is partly remedied by the recently proposed CASPT2/CC approach [Phung et al., J. Chem. Theory Comput., 2018, 14, 2446-2455]; (4) demonstrates that NEVPT2 performs worse than CASPT2, by giving errors up to 7 kcal mol-1; (5) shows that the accuracy of MRCISD+Q spin-state energetics strongly depends on the size-consistency correction: the Davidson-Silver and Pople corrections perform best (MAE < 3 kcal mol-1), whereas the standard Davidson correction is not recommended (MAE of 7 kcal mol-1). Only a few DFT methods (including the best performing ones identified in this study: B2PLYP-D3 and OPBE) are able to provide a balanced description of the spin-state energetics for all four studied iron complexes simultaneously, corroborating the non-universality problem of approximate density functionals.

A computational study of the Fenton reaction in different pH ranges

The mechanism of the Fenton reaction is pH dependent and four distinct reactive species have been identified and found to display quite different oxidation reactivities.

New oxidovanadium(iv) complex with a BIAN ligand: synthesis, structure, redox properties and catalytic activity

The combination of a new oxidovanadium(iv) complex1with pyrazine-2-carboxylic acid (PCA; a cocatalyst) affords a catalytic system for the efficient oxidation of saturated hydrocarbons.

Mechanism of O2 Activation by α-Ketoglutarate Dependent Oxygenases Revisited. A Quantum Chemical Study

Four mechanisms previously proposed for dioxygen activation catalyzed by α-keto acid dependent oxygenases (α-KAO) were studied with dispersion-corrected DFT methods employing B3LYP and TPSSh functionals in combination with triple-ζ basis set (cc-pVTZ). The aim of this study was to revisit mechanisms suggested in the past decade and resolve remaining issues related to dioxygen activation. Mechanism A, which runs on the quintet potential energy surface (PES) and includes formation of an Fe(III)-superoxide radical anion complex, subsequent oxidative decarboxylation, and O-O bond cleavage, was found to be most likely. However, mechanism B taking place on the septet PES involves a rate limiting barrier comparable to the one found for mechanism A, and thus it cannot be excluded, though two other mechanisms (C and D) were ruled out. Mechanism C is a minor variation of mechanism A, whereas mechanism D proceeds through formation of a triplet Fe(IV)-alkyl peroxo bridged intermediate. The study covered also full optimization of relevant minimum energy crossing points (MECPs). The relative energy of critical intermediates was also studied with the CCSD(T) method in order to benchmark TPSSh and B3LYP functionals with respect to their credibility in predicting relative energies of septet and triplet spin states of the ternary enzyme-Fe-α-keto glutarate (α-KG)-O2 complex.

Ammonia-modified Co(II) sites in zeolites: spin and electron density redistribution through the Co(II)-NO bond

Electronic factors essential for the bonding of a non-innocent NO ligand to ammonia-modified Co(2+) sites in cobalt-exchanged zeolites are examined for small cluster models using DFT and advanced correlated wave function calculations. The analysis of charge transfer processes between the NO ligand and the cobalt center involves two protocols: valence-bond expansion of the multiconfiguration CASSCF wave function (in terms of fragment-localized active orbitals) and spin-resolved natural orbitals for chemical valence (SR-NOCV). Applicability of SR-NOCV analysis to transition metal complexes involving non-innocent fragments is critically assessed and the approach based on the CASSCF wave function turns out to be much more robust and systematic for all studied models. It is shown that the character and direction of electron density redistribution through the Co-N-O bond, quantified by relative share of the Co(II)-NO(0), Co(III)-NO(-), and Co(I)-NO(+) resonance structures in the total wave function, fully rationalize the activation of the N-O bond upon NH3 co-ligation (evidenced by calculated and measured red-shift of the NO stretching frequency and commonly ascribed to enhanced backdonation). The huge red-shift of νN-O is attributed to an effective electron transfer between the ammonia-modified Co(ii) centers and the NO antibonding π*-orbitals (related to the increased share of the Co(III)-NO(-) form). Unexpectedly, the effect is stronger for the singlet complex with three NH3 ligands than for that with five NH3 ligands bound to the cobalt center. Our results also indicate that high-efficiency electron transfers between the Co(ii) center and the NO ligand may be enabled for the selected spin state and disabled for the other spin state of the adduct. This illustrates how the cobalt center may serve to fine-tune the electronic communication between the NO ligand and its binding site.

Role of Spin States in Nitric Oxide Binding to Cobalt(II) and Manganese(II) Porphyrins. Is Tighter Binding Always Stronger?

Binding of nitric oxide (NO) to metalloporphyrins and heme groups is important in biochemistry while challenging to describe accurately by density functional theory (DFT) calculations. Here, the structural and thermochemical aspect of NO binding to Co(II) and Mn(II) porphyrins is investigated by DFT and DFT-D (dispersion-corrected) calculations, supported by reliable coupled-cluster methodology (CCSD(T)), and critically correlated with the experimental data. It is argued that whereas the bonding of NO to Co(II) porphyrin is a simple radical recombination, the bonding of NO to Mn(II) porphyrin is accompanied by a crossing of spin states. For this reason, the spin-state conversion energy contributes to the Mn-NO bond energy, and the paradigmatic correlation between bond length and bond energy is violated for the considered nitrosyl complexes: the Mn-NO bond is (structurally) shorter by ∼0.2 Å, albeit (energetically) weaker by a few kcal/mol, compared with the Co-NO bond. Moreover, none of the many tested DFT methods can reproduce the Co-NO and Mn-NO bond energies simultaneously, except for calculations with B3LYP*-D3, TPSSh-D3, and M06-D3 methods supplemented with the proposed spin-state energy correction (to compensate for an error on the calculated spin-state conversion energy). The results of this study are important to appreciate the role of spin-state changes in ligand binding properties of heme-related models. They also highlight the need for accurate calculations for correct interpretation of experimental data, including the qualitative structure-energy relationship.