Hydrogen peroxide decomposition by a non-heme iron(III) catalase mimic: a DFT study

Non-heme iron(III) complexes of 14-membered tetraaza macrocycles have previously been found to catalytically decompose hydrogen peroxide to water and molecular oxygen, like the native enzyme catalase. Here the mechanism of this reaction is theoretically investigated by DFT calculations at the (U)B3LYP/6-31G* level, with focus on the reactivity of the possible spin states of the FeIII complexes. The computations suggest that H2O2 decomposition follows a homolytic route with intermediate formation of an iron(IV) oxo radical cation species (L.+FeIV==O) that resembles Compound I of natural iron porphyrin systems. Along the whole catalytic cycle, no significant energetic differences were found for the reaction proceeding on the doublet (S=1/2) or on the quartet (S=3/2) hypersurface, with the single exception of the rate-determining O--O bond cleavage of the first associated hydrogen peroxide molecule, for which reaction via the doublet state is preferred. The sextet (S=5/2) state of the FeIII complexes appears to be unreactive in catalase-like reactions.     

103Rh NMR chemical shifts in organometallic complexes: a combined experimental and density functional study

Experimental 103Rh NMR chemical shifts of mono- and binuclear rhodium(I) complexes containing s- or as-hydroindacenide and indacenediide bridging ligands with different ancillary ligands (1,5-cyclooctadiene, ethylene, carbonyl) are presented. A protocol, based on density functional theory calculations, was established to determine 103Rh NMR shielding constants in order to rationalise the effects of electronic and structural variations on the spectroscopic signal, and to gain insight into the efficiency of this computational method when applied to organometallic systems. Scalar and spin-orbit relativistic effects based on the ZORA (zeroth order regular approximation) level have been taken into account and discussed. A good agreement was found for model compounds over a wide range of chemical shifts of rhodium (approximately 10,000 ppm). This allowed us to discuss the experimental and calculated delta(103Rh) in larger complexes and to relate it to their electronic structure.     

Characterization of three members of the electron-transfer series [Fe(pda)2]n (n=2-, 1-, 0) by spectroscopy and density functional theoretical calculations [pda=redox non-innocent derivatives of N,N'-bis(pentafluorophenyl)-o-phenylenediamide(2-, 1.-, 0)

The four-coordinate iron complexes, [Fe(III)(pda(2-))(pda(.-))] (1) and [AsPh(4)](2)[Fe(II)(pda(2-))(2)] (2) were synthesized and fully characterized; pda(2-) is the closed-shell ligand N,N'-bis(pentafluorophenyl)-o-phenylenediamido(2-), and pda(.-) represents its one-electron-oxidized pi-radical anion. Single-crystal X-ray diffraction studies of 1 and 2 performed at 100(2) K reveal a distorted tetrahedral coordination environment at the iron centers, as a result of the intramolecular pi-pi interactions between C(6)F(5) rings. The electronic structures of 1 and 2 were unambiguously determined by a combination of (57)Fe M?ssbauer and electronic spectroscopy, magnetic susceptibility measurements, X-ray crystallography, and DFT calculations. Compound 1 contains an intermediate-spin Fe(III) ion (S(Fe)=3/2) strongly antiferromagnetically coupled to a pi-ligand radical (S(R)=1/2) yielding an S(t)=1 ground state. Complex 2 possesses a high-spin Fe(II) center (S(Fe)=2) with two closed-shell dianionic ligands. Complexes 1 and 2 are members of the redox series [Fe(pda)(2)](n) with n=0 for 1 and n=2- for 2. The anion n=1- has been reported previously in the coordination salt [Fe(dad)(3)][Fe(pda)(2)] (3; dad=N,N'-bis(phenyl)-2,3-dimethyl-1,4-diaza-1,3-butadiene). A complicated temperature-dependent electronic structure has been observed for this salt. Here, DFT calculations performed on 3 confirm the previous assignments of spin- and oxidation-states. Thus, [Fe(pda)(2)](n) (n=0, 1-, 2-) constitutes an electron-transfer series, which has also been established by cyclic voltammetry; the mono- and dications (n=1+ and 2+) are also accessible in solution, but have not been further investigated. The (57)Fe M?ssbauer spectra of [Fe(pda)(2)](n) species in 1 and 3 show extremely large quadrupole splitting constants due to addition of the valence and covalence contributions that have been confirmed by DFT calculations.     

Generating and dimerizing the transient 16-electron phosphinidene complex [Cp*Ir=PAr]: a theoretical and experimental study

The properties of the 16-electron phosphinidene complex [CpRIr=PR] were investigated experimentally and theoretically. Density functional theory calculations show a preferred bent geometry for the model complex [CpIr=PH], in contrast to the linear structure of [CpIr=NH]. Dimerization to give [[CpIr=PH]2] and ligand addition to afford [Cp(L)Ir=PH] (L=PH3, CO) were calculated to give compounds that were energetically highly favorable, but which differed from the related imido complexes. Transient 16-electron phosphinidene complex [Cp*Ir=PAr] could not be detected experimentally. Dehydrohalogenation of [Cp*IrCl2(PH2Ar)] in CH2Cl2 at low temperatures resulted in the novel fused-ring systems 17 (Ar=Mes*) and 20 (Ar=Mes), with dimeric [[Cp*Ir=PAr]2] being the likely intermediate. Intramolecular C-H bond activation induced by steric factors is considered to be the driving force for the irreversible formation of 17 and 20. ONIOM calculations suggest this arises because of the large steric congestion in [[Cp*Ir=PAr]2], which forces it toward a more reactive planar structure that is apt to rearrange.     

The mechanism of the (bispidine)copper(II)-catalyzed aziridination of styrene: a combined experimental and theoretical study

Experimental and DFT-based computational results on the aziridination mechanism and the catalytic activity of (bispidine)copper(I) and -copper(II) complexes are reported and discussed (bispidine=tetra- or pentadentate 3,7-diazabicyclo[3.1.1]nonane derivative with two or three aromatic N donors in addition to the two tertiary amines). There is a correlation between the redox potential of the copper(II/I) couple and the activity of the catalyst. The most active catalyst studied, which has the most positive redox potential among all (bispidine)copper(II) complexes, performs 180 turnovers in 30 min. A detailed hybrid density functional theory (DFT) study provides insight into the structure, spin state, and stability of reactive intermediates and transition states, the oxidation state of the copper center, and the denticity of the nitrene source. Among the possible pathways for the formation of the aziridine product, the stepwise formation of the two N-C bonds is shown to be preferred, which also follows from experimental results. Although the triplet state of the catalytically active copper nitrene is lowest in energy, the two possible spin states of the radical intermediate are practically degenerate, and there is a spin crossover at this stage because the triplet energy barrier to the singlet product is exceedingly high.     

Experimental and DFT Evidence for the Fractional Non‐Innocence of a β‐Diketonate Ligand

New compounds [Ru(pap)₂(L)](ClO₄), [Ru(pap)(L)₂], and [Ru(acac)₂(L)] (pap=2‐phenylazopyridine, L^?=9‐oxidophenalenone, acac^?=2,4‐pentanedionate) have been prepared and studied regarding their electron‐transfer behavior, both experimentally and by using DFT calculations. [Ru(pap)₂(L)](ClO₄) and [Ru(acac)₂(L)] were characterized by crystal‐structure analysis. Spectroelectrochemistry (EPR, UV/Vis/NIR), in conjunction with cyclic voltammetry, showed a wide range of about 2 V for the potential of the Ru^III/II couple, which was in agreement with the very different characteristics of the strongly π‐accepting pap ligand and the σ‐donating acac^? ligand. At the rather high potential of +1.35 V versus SCE, the oxidation of L^? into L^. could be deduced from the near‐IR absorption of [Ru^III(pap)(L^.)(L^?)]²⁺. Other intense long‐wavelength transitions, including LMCT (L^?→Ru^III) and LL/CT (pap^.?→L^?) processes, were confirmed by TD‐DFT results. DFT calculations and EPR data for the paramagnetic intermediates allowed us to assess the spin densities, which revealed two cases with considerable contributions from L‐radical‐involving forms, that is, [Ru^III(pap⁰)₂(L^?)]²⁺?[Ru^II(pap⁰)₂(L^.)]²⁺ and [Ru^III(pap⁰)(L^?)₂]⁺?[Ru^II(pap⁰)(L^?)(L^?)]⁺. Calculations of electrogenerated complex [Ru^II(pap^.?)(pap⁰)(L^?)] displayed considerable negative spin density (?0.188) at the bridging metal.     

Four‐Center Oxidation State Combinations and Near‐Infrared Absorption in [Ru(pap)(Q)2]n (Q=3,5‐Di‐tert‐butyl‐N‐aryl‐1,2‐benzoquinonemonoimine,pap=2‐Phenylazopyridine)

The complex series [Ru(pap)(Q)₂]ⁿ ([ 1 ]ⁿ–[ 4 ]ⁿ; n=+2, +1, 0, ?1, ?2) contains four redox non‐innocent entities: one ruthenium ion, 2‐phenylazopyridine (pap), and two o‐iminoquinone moieties, Q=3,5‐di‐tert‐butyl‐N‐aryl‐1,2‐benzoquinonemonoimine (aryl=C₆H₅ ( 1⁺ ); m‐(Cl)₂C₆H₃ ( 2⁺ ); m‐(OCH₃)₂C₆H₃ ( 3⁺ ); m‐(tBu)₂C₆H₃ ( 4 ⁺)). A crystal structure determination of the representative compound, [ 1 ]ClO₄, established the crystallization of the ctt‐isomeric form, that is, cis and trans with respect to the mutual orientations of O and N donors of two Q ligands, and the coordinating azo N atom trans to the O donor of Q. The sensitive C? O (average: 1.299(3) Å), C? N (average: 1.346(4) Å) and intra‐ring C? C (meta; average: 1.373(4) Å) bond lengths of the coordinated iminoquinone moieties in corroboration with the N?N length (1.292(3) Å) of pap in 1 ⁺ establish [Ru^III(pap⁰)(Q^.?)₂]⁺ as the most appropriate electronic structural form. The coupling of three spins from one low‐spin ruthenium(III) (t_2g⁵) and two Q^.? radicals in 1 ⁺– 4 ⁺ gives a ground state with one unpaired electron on Q^.?, as evident from g=1.995 radical‐type EPR signals for 1 ⁺– 4 ⁺. Accordingly, the DFT‐calculated Mulliken spin densities of 1 ⁺ (1.152 for two Q, Ru: ?0.179, pap: 0.031) confirm Q‐based spin. Complex ions 1 ⁺– 4 ⁺ exhibit two near‐IR absorption bands at about λ=2000 and 920 nm in addition to intense multiple transitions covering the visible to UV regions; compounds [ 1 ]ClO₄–[ 4 ]ClO₄ undergo one oxidation and three separate reduction processes within ±2.0 V versus SCE. The crystal structure of the neutral (one‐electron reduced) state ( 2 ) was determined to show metal‐based reduction and an EPR signal at g=1.996. The electronic transitions of the complexes 1 ⁿ– 4 ⁿ (n=+2, +1, 0, ?1, ?2) in the UV, visible, and NIR regions, as determined by using spectroelectrochemistry, have been analyzed by TD‐DFT calculations and reveal significant low‐energy absorbance (λ_max>1000 nm) for cations, anions, and neutral forms. The experimental studies in combination with DFT calculations suggest the dominant valence configurations of 1 ⁿ– 4 ⁿ in the accessible redox states to be [Ru^III(pap⁰)(Q^.?)(Q⁰)]²⁺ ( 1 ²⁺– 4 ²⁺)→[Ru^III(pap⁰)(Q^.?)₂]⁺ ( 1 ⁺– 4 ⁺)→[Ru^II(pap⁰)(Q^.?)₂] ( 1 – 4 )→[Ru^II(pap^.?)(Q^.?)₂]^? ( 1 ^?– 4 ^?)→[Ru^III(pap^.?)(Q^2?)₂]^2? ( 1 ^2?– 4 ^2?).     

A Unified Treatment of the Relationship Between Ligand Substituents and Spin State in a Family of Iron(II) Complexes

The influence of ligands on the spin state of a metal ion is of central importance for bioinorganic chemistry, and the production of base‐metal catalysts for synthesis applications. Complexes derived from [Fe(bpp)₂]²⁺ (bpp=2,6‐di{pyrazol‐1‐yl}pyridine) can be high‐spin, low‐spin, or spin‐crossover (SCO) active depending on the ligand substituents. Plots of the SCO midpoint temperature (T ) in solution vs. the relevant Hammett parameter show that the low‐spin state of the complex is stabilized by electron‐withdrawing pyridyl (“X”) substituents, but also by electron‐donating pyrazolyl (“Y”) substituents. Moreover, when a subset of complexes with halogeno X or Y substituents is considered, the two sets of compounds instead show identical trends of a small reduction in T for increasing substituent electronegativity. DFT calculations reproduce these disparate trends, which arise from competing influences of pyridyl and pyrazolyl ligand substituents on Fe‐L σ and π bonding.     

Homolysis of N-alkoxyamines: a computational study.

During nitroxide-mediated polymerization (NMP) in the presence of a nitroxide R2(R1)NO*, the reversible formation of N-alkoxyamines [P-ON(R1)R2] reduces significantly the concentration of polymer radicals (P*) and their involvement in termination reactions. The control of the livingness and polydispersity of the resulting polymer depends strongly on the magnitude of the bond dissociation energy (BDE) of the C-ON(R1)R2 bond. In this study, theoretical BDEs of a large series of model N-alkoxyamines are calculated with the PM3 method. In order to provide a predictive tool, correlations between the calculated BDEs and the cleavage temperature (Tc), and the dissociation rate constant (k(d)), of the N-alkoxyamines are established. The homolytic cleavage of the N-OC bond is also investigated at the B3P86/6-311++G(d,p)//B3LYP/6-31G(d), level. Furthermore, a natural bond orbital analysis is carried out for some N-alkoxyamines with a O-C-ON(R1)R2 fragment, and the strengthening of their C-ON(R1)R2 bond is interpreted in terms of stabilizing anomeric interactions.     

Growing graphene sheets from reactions with methyl radicals: a quantum chemical study.

Hydrogen abstraction reactions by methyl radicals on the zigzag and armchair edges of perylene are studied by density functional theory (DFT) to explore various growth pathways that seem to be in line with experimental observations. The DFT approach is validated by comparing the results obtained from calculations with six different functionals with those obtained from correlated ab initio methods, thereby emphasizing the calculation of reaction barriers. A useful compromise between accuracy and computational cost is provided by DFT, and possible pathways are studied in detail at this level of calculation. Our computational study is carried out by ordering, as a first step, all of the isomers that arise from the abstraction of one or two H atoms from 1,12-dimethyl-1,12-dihydroperylene and 3,4-dimethyl-3,4-dihydroperylene with respect to their energies. Subsequently, only those pathways that connect low-energy isomers are investigated. The calculations reveal that the selected pathways are favored thermodynamically, and also that the reaction barriers are somewhat higher than the energy locally available for the respective reaction. Notably, in the case of 3,4-dimethyl-3,4-dihydroperylene, the first two reaction steps have no or only a very low reaction barrier. The final conclusion of our study is that a cascade of reactions is possible that leads to the growth of a graphene sheet on a graphite surface.     

Wolves in Sheep's Clothing: μ‐Oxo‐Diiron Corroles Revisited

For well over 20 years, μ‐oxo‐diiron corroles, first reported by Vogel and co‐workers in the form of μ‐oxo‐bis[(octaethylcorrolato)iron] (Mössbauer δ 0.02 mm s^?1, ΔE_Q 2.35 mm s^?1), have been thought of as comprising a pair antiferromagnetically coupled low‐spin Fe^IV centers. The remarkable stability of these complexes, which can be handled at room temperature and crystallographically analyzed, present a sharp contrast to the fleeting nature of enzymatic, iron(IV)‐oxo intermediates. An array of experimental and theoretical methods have now shown that the iron centers in these complexes are not Fe^IV but intermediate‐spin Fe^III coupled to a corrole^.2?. The intramolecular spin couplings in {Fe[TPC]}₂(μ‐O) were analyzed via DFT(B3LYP) calculations in terms of the Heisenberg–Dirac–van Vleck spin Hamiltonian H=J_Fe–corrole(S_Fe?S_corrole)+J_Fe–Fe′(S_Fe?S_Fe′)+J_{Fe′–corrole}(S_Fe′?S_corrole′), which yielded J_Fe–corrole=J_{Fe′–corrole′}=0.355 eV (2860 cm^?1) and J_Fe–Fe′=0.068 eV (548 cm^?1). The unexpected stability of μ‐oxo‐diiron corroles thus appears to be attributable to charge delocalization via ligand noninnocence.