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

First-Principles Investigations of Magnetic Anisotropy and Spin-Crossover Behavior of Fe(III)-TBP Complexes

With the ongoing effort to obtain mononuclear 3d-transition-metal complexes that manifest slow relaxation of magnetization and, hence, can behave as single-molecule magnets (SMMs), we have modeled 14 Fe(III) complexes based on an experimentally synthesized (PMe₃)₂FeCl₃ complex [J. Am. Chem. Soc. 2017, 139 (46), 16474-16477], by varying the axial ligands with group XV elements (N, P, and As) and equatorial halide ligands from F, Cl, Br, and I. Out of these, nine complexes possess large zero field splitting (ZFS) parameter D in the range of -40 to -60 cm^-1. The first-principles investigation of the ground-spin state applying density functional theory (DFT) and wave function-based multiconfigurations methods, e.g., SA-CASSCF/NEVPT2, are found to be quite consistent except for few delicate cases with near-degenerate spin states. In such cases, the hybrid B3LYP functional is found to be biased toward high-spin (HS) state. Altering the percentage of exact exchange admixed in the B3LYP functional leads to intermediate-spin (IS) ground state consistent with the multireference calculations. The origin of large zero field splitting (ZFS) in the Fe(III)-based trigonal bipyramidal (TBP) complexes is investigated. Furthermore, a number of complexes are identified with very small ΔG_HS-IS^adia. values indicating the possible spin-crossover phenomenon between the bistable spin states.

Zeolites at the Molecular Level: What Can Be Learned from Molecular Modeling

This review puts the development of molecular modeling methods in the context of their applications to zeolitic active sites. We attempt to highlight the utmost necessity of close cooperation between theory and experiment, resulting both in advances in computational methods and in progress in experimental techniques.

Theoretical modeling of the singlet-triplet spin transition in different Ni(II)-diketo-pyrphyrin-based metal-ligand octahedral complexes

The structural stability, charge transfer effects and strength of the spin-orbit couplings in different Ni(ii)-ligand complexes have been studied at the DFT (B3LYP and CAM-B3LYP) and coupled cluster (DLPNO-CCSD(T)) levels of theory. Accordingly, two different, porphyrin- and diketo-pyrphyrin-based four-coordination macrocycles as planar ligands as well as pyridine (or pyrrole) and mesylate anion molecular groups as vertical ligands were considered in order to build metal-organic complexes with octahedral coordination configurations. For each molecular system, the identification of equilibrium geometries and the intersystem crossing (the minimum energy crossing) points between the potential energy surfaces of the singlet and triplet spin states is followed by computing the spin-orbit couplings between the two spin states. Structures, based on the diketo-pyrphyrin macrocycle as the planar ligand, show stronger six-coordination metal-organic complexes due to the extra electrostatic interaction between the positively charged central metal cation and the negatively charged vertical ligands. The results also show that the magnitude of the spin-orbit coupling is influenced by the atomic positions of deprotonations of the ligands, and implicitly the direction of the charge transfer between the ligand and the central metal ion.

Scalable Electron Correlation Methods. 8. Explicitly Correlated Open-Shell Coupled-Cluster with Pair Natural Orbitals PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12

We present explicitly correlated open-shell pair natural orbital local coupled-cluster methods, PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. The methods are extensions of our previously reported PNO-R/UCCSD methods (J. Chem. Theory Comput., 2020, 16, 3135-3151, https://pubs.acs.org/doi/10.1021/acs.jctc.0c00192) with additions of explicit correlation and perturbative triples corrections. The explicit correlation treatment follows the spin-orbital CCSD-F12b theory using Ansatz 3*A, which is found to yield comparable or better basis set convergence than the more rigorous Ansatz 3C in computed ionization potentials and reaction energies using double- to quaduple-ζ basis sets. The perturbative triples correction is adapted from the spin-orbital (T) theory to use triples natural orbitals (TNOs). To address the coupling due to off-diagonal Fock matrix elements, the local triples amplitudes are iteratively solved using small domains of TNOs, and a semicanonical (T0) domain correction with larger domains is applied to reduce the domain errors. The performance of the methods is demonstrated through benchmark calculations on ionization potentials, radical stabilization energies, reaction energies of fragmentations and rearrangements in radical cations, and spin-state energy differences of iron complexes. For a few test sets where canonical calculations are feasible, PNO-RCCSD(T)-F12 results agree with the canonical ones to within 0.4 kcal mol^-1, and this maximum error is reduced to below 0.2 kcal mol^-1 when large local domains are used. For larger systems, results using different thresholds for the local approximations are compared to demonstrate that 1 kcal mol^-1 level of accuracy can be achieved using our default settings. For a couple of difficult cases, it is demonstrated that the errors from individual approximations are only a fraction of 1 kcal mol^-1, and the overall accuracy of the method does not rely on error compensations. In contrast to canonical calculations, the use of spin-orbitals does not lead to a significant increase of computational time and memory usage in the most expensive steps of PNO-R/UCCSD(T)-F12 calculations. The only exception is the iterative solution of the (T) amplitudes, which can be avoided without significant errors by using a perturbative treatment of the off-diagonal coupling, known as (T1) approximation. For most systems, even the semicanonical approximation (T0) leads only to small errors in relative energies. Our program is well parallelized and capable of computing accurate correlation energies for molecules with 100-200 atoms using augmented triple-ζ basis sets in less than a day of elapsed time on a small computer cluster.

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.

Spin-state dependence of exchange-correlation holes

Applications of density-functional theory (DFT) in computational chemistry rely on an approximate exchange-correlation (xc) functional. However, existing approximations can fail dramatically for open-shell molecules, in particular for transition-metal complexes or radicals. Most importantly, predicting energy differences between different spin-states with approximate exchange-correlation functionals remains extremely challenging. Formally, it is known that the exact xc functional should be spin-state dependent, but none of the available approximations feature such an explicit spin-state dependence [C. R. Jacob and M. Reiher, Int. J. Quantum Chem., 2012, 112, 3661-3684]. Thus, to find novel approximations for the xc functional for open-shell systems, the development of spin-state dependent xc functionals appears to be a promising avenue. Here, we set out to shed light on the spin-state dependence of the xc functional by investigating the underlying xc holes, which we extract from configuration interaction calculations for model systems. We analyze the similarities and differences between the xc holes of the lowest-energy singlet and triplet states of the dihydrogen molecule, the helium atom, and the lithium dimer. To shed further light on the spin-state dependence of these xc holes we also discuss exact conditions that can be derived from the spin structure of the reduced two-electron density matrix. Altogether, our results suggest several possible routes towards the construction of explicitly spin-state dependent approximations for the xc functional.

A Review of Density Functional Models for the Description of Fe(II) Spin-Crossover Complexes

Spin-crossover (SCO) materials have for more than 30 years stood out for their vast application potential in memory, sensing and display devices. To reach magnetic multistability conditions, the high-spin (HS) and low-spin (LS) states have to be carefully balanced by ligand field stabilization and spin-pairing energies. Both effects could be effectively modelled by electronic structure theory, if the description would be accurate enough to describe these concurrent influences to within a few kJ/mol. Such a milestone would allow for the in silico-driven development of SCO complexes. However, so far, the ab initio simulation of such systems has been dominated by general gradient approximation density functional calculations. The latter can only provide the right answer for the wrong reasons, given that the LS states are grossly over-stabilized. In this contribution, we explore different venues for the parameterization of hybrid functionals. A fitting set is provided on the basis of explicitly correlated coupled cluster calculations, with single- and multi-dimensional fitting approaches being tested to selected classes of hybrid functionals (hybrid, range-separated, and local hybrid). Promising agreement to benchmark data is found for a rescaled PBE0 hybrid functional and a local version thereof, with a discussion of different atomic exchange factors.

Valence and Structure Isomerism of Al2FeO4+: Synergy of Spectroscopy and Quantum Chemistry

We provide spectroscopic and computational evidence for a substantial change in structure and gas phase reactivity of Al₃O₄⁺ upon Fe-substitution, which is correctly predicted by multireference (MR) wave function calculations. Al₃O₄⁺ exhibits a cone-like structure with a central trivalent O atom (C_3v symmetry). The replacement of the Al- by an Fe atom leads to a planar bicyclic frame with a terminal Al-O^•- radical site, accompanied by a change from the Fe^+III/O^-II to the Fe^+II/O^-I valence state. The gas phase vibrational spectrum of Al₂FeO₄⁺ is exclusively reproduced by the latter structure, which MR wave function calculations correctly identify as the most stable isomer. This isomer of Al₂FeO₄⁺ is predicted to be highly reactive with respect to C-H bond activation, very similar to Al₈O₁₂⁺ which also features the terminal Al-O^•- radical site. Density functional theory, in contrast, predicts a less reactive Al₃O₄⁺-like "isomorphous substitution" structure of Al₂FeO₄⁺ to be the most stable one, except for functionals with very high admixture of Fock exchange (50%, BHLYP).

Dioxygen Binding to all 3d, 4d, and 5d Transition Metals from Coupled-Cluster Theory

Understanding how transition metals bind and activate dioxygen (O₂ ) is limited by experimental and theoretical uncertainties, making accurate quantum mechanical descriptors of interest. Here we report coupled-cluster CCSD(T) energies with large basis sets and vibrational and relativistic corrections for 160 3d, 4d, and 5d metal-O₂ systems. We define four reaction energies (120 in total for the 30 metals) that quantify O-O activation and reveal linear relationships between metal-oxygen and O-O binding energies. The CCSD(T) data can be combined with thermochemical cycles to estimate chemisorption and physisorption energies for each metal from metal oxide embedding energies, in good correlation with atomization enthalpies (R² =0.75). Spin-geometry variations can break the linearities, of interest to circumventing the Sabatier principle. Pt, Pd, Co, and Fe form a distinct group with the weakest O₂ binding. R² up to 0.84 between surface adsorption energies and our energies for MO₂ systems indicate relevance also to real catalytic systems.

A DFT Study on FeI/FeII/FeIII Mechanism of the Cross-Coupling between Haloalkane and Aryl Grignard Reagent Catalyzed by Iron-SciOPP Complexes

To explore plausible reaction pathways of the cross-coupling reaction between a haloalkane and an aryl metal reagent catalyzed by an iron–phosphine complex, we examine the reaction of FeBrPh(SciOPP) 1 and bromocycloheptane employing density functional theory (DFT) calculations. Besides the cross-coupling, we also examined the competitive pathways of β-hydrogen elimination to give the corresponding alkene byproduct. The DFT study on the reaction pathways explains the cross-coupling selectivity over the elimination in terms of Fe^I/Fe^II/Fe^III mechanism which involves the generation of alkyl radical intermediates and their propagation in a chain reaction manner. The present study gives insight into the detailed molecular mechanic of the cross-coupling reaction and revises the Fe^II/Fe^II mechanisms previously proposed by us and others.

Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe²⁺, Fe³⁺, and Co³⁺, including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. We find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method-in terms of computational cost and in terms of accuracy on monatomic systems-and therefore we chose this method for molecular benchmarks. Then Kohn-Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF₂, FeCl₂, and CoF₂. We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.

Assessment of the SCAN Functional for Spin-State Energies in Spin-Crossover Systems

The strongly constrained and appropriately normed (SCAN) functional has been tested toward the calculation of spin-state energy differences in a data set of 20 spin-crossover (SCO) systems, ranging from d⁴ to d⁷. Results show that the SCAN functional is able to correctly predict the low-spin state as the ground state for all systems, and the energy window provided by the calculations falls in the approximate range of energies that will allow for SCO to occur. Moreover, the SCAN functional can be used in periodic boundary condition calculations, accounting for the effect of collective crystal vibrations and counterions in the thermochemistry of the spin transition. Our results validate this functional as a potential method for in silico screening of new SCO systems at both, molecular and crystal-packed levels.

Scalable Electron Correlation Methods. 7. Local Open-Shell Coupled-Cluster Methods Using Pair Natural Orbitals: PNO-RCCSD and PNO-UCCSD

We present well-parallelized local implementations of high-spin open-shell coupled cluster methods with single and double excitations (CCSD) using pair natural orbitals (PNOs). The methods are based on the spin-orbital coupled cluster theory using restricted open-shell Hartree-Fock (ROHF) reference functions. Two variants, namely, PNO-UCCSD and PNO-RCCSD are implemented and compared. In PNO-UCCSD, the coupled cluster amplitudes are spin-unrestricted, while in PNO-RCCSD the linear terms are spin-adapted by a spin-projection approach as described in J. Chem. Phys. 1993, 99, 5219-5227. Near linear scaling of the computational cost with the number of correlated electrons is achieved by applying domain and pair approximations. The PNOs are spin-independent and obtained using a semicanonical spin-restricted MP2 approximation with large domains of projected atomic orbitals (PAOs). The pair approximations of our previously described closed-shell PNO-LCCSD method are carefully revised so that they are compatible to the UCCSD theory, and PNO-UCCSD or PNO-RCCSD calculations for closed-shell molecules yield exactly the same results as corresponding spin-free closed-shell PNO-LCCSD calculations. The convergence of the results with respect to the thresholds and options that control the domain and pair approximations is demonstrated. It is found that large domains are required for the single excitations in open-shell calculations in order to obtain converged results. In general, the errors of relative energies caused by the local approximations can be reduced to below 1 kcal mol^-1, even for difficult cases. Presently, PNO-RCCSD and PNO-UCCSD calculations for molecules with 100-200 atoms and augmented triple-ζ basis sets can be carried out in a few hours of elapsed time using ∼100 CPU cores. In addition, the program is also capable of performing distinguishable cluster (PNO-RDCSD and PNO-UDCSC) calculations. The present work is a critical step in developing fully local open-shell PNO-RCCSD(T)-F12 methods.

Anharmonic Frequencies of (MO)2 and Related Hydrides for M = Mg, Al, Si, P, S, Ca, and Ti and Heuristics for Predicting Anharmonic Corrections of Inorganic Oxides

The low-frequency vibrational fundamentals of D_2h inorganic oxides are readily modeled by heuristic scaling factors at fractions of the computational cost compared to explicit anharmonic frequency computations. Oxygen and the other elements in the present study are abundant in geochemical environments and have the potential to aggregate into minerals in planet-forming regions or in the remnants of supernovae. Explicit quartic force field computations at the CCSD(T)-F12b/cc-pVTZ-F12 level of theory generate scaling factors that accurately predict the anharmonic frequencies with an average error of less than 1.0 cm^-1 for both the metal-oxygen stretching frequencies and the torsion and antisymmetric stretching frequencies. Inclusion of hydrogen motions is less absolutely accurate but is similarly relatively predictive. The fundamental vibrational frequencies for the seven tetra-atomic inorganic oxides examined presently fall below 876 cm^-1 and most of the hydrogenated species do as well. Additionally, ν₆ for the SiO dimer is shown to have an intensity of 562 km mol^-1, with each of the other molecules having one or more frequencies with intensities greater than 80 km mol^-1, again with most in the low-frequency infrared range. These intensities and the frequencies computed in the present study should assist in laboratory characterization and potential interstellar or circumstellar observation.

The Molpro quantum chemistry package

Molpro is a general purpose quantum chemistry software package with a long development history. It was originally focused on accurate wavefunction calculations for small molecules but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding, and anharmonic vibrational spectra. In addition to conventional input-file specification of calculations, Molpro calculations can now be specified and analyzed via a new graphical user interface and through a Python framework.

Thermochemistry of FeOmHnz Species: Assessment of Some DFT Functionals

Thermochemical data for 20 anionic, cationic, and neutral gas-phase species, including Fe^0/+, FeO^-/0/+/2+, FeOH^0/+/2+, FeO₂^-/0/+, OFeOH^0/+, Fe(OH)₂^0/+, Fe(H₂O)^+/2+, and Fe(H₂O)₂^+/2+ with oxidation states between +I and +IV for Fe and -I and -II for O, compiled by Schröder [ J. Phys. Chem. A 2008, 112, 13215], are used to assess the performance of the "Jacob's ladder" functionals PBE, TPSS, PBE0, and TPSSh for the SVP, TZVP, and QZVP basis sets. In addition, the BP86 and B3LYP functionals are considered. The TPSSh functional performs best. With the TZVP basis set (recommended), the mean absolute and the maximum errors are 24 and 63 kJ/mol, respectively. With 32 and 78 kJ/mol, respectively, BP86 is second best, better than PBE.

Transition Metal Spin-State Energetics by MC-PDFT with High Local Exchange

The energetics of the spin states of transition metal complexes have been explored with a variety of electronic structure methods, but the calculations require a compromise between accuracy and affordability. In this work, the spin splittings of several iron complexes are studied with multiconfiguration pair-density functional theory (MC-PDFT). The results are compared to previously published results obtained by complete active space second-order perturbation theory (CASPT2) and CASPT2 with coupled-cluster semicore correlation (CASPT2/CC). In contrast to CASPT2's systematic overstabilization of high-spin states with respect to the CASPT2/CC reference, MC-PDFT with the tPBE on-top functional understabilizes high-spin states. This systematic understabilization is largely corrected by revising the exchange and correlation contributions to the on-top functional using the high local-exchange approximation (tPBE-HLE). Moreover, tPBE-HLE correctly predicts the spin of the ground state in most cases, while CASPT2 incorrectly predicts high-spin ground states in all cases. This is encouraging for practical work because tPBE and tPBE-HLE are faster than CASPT2 by a factor of 50 even in a moderately sized example.

Spin modification of iron(ii) complexes via covalent (dative) and dispersion guided non-covalent bonding with N-heterocyclic carbenes: DFT, DLPNO-CCSD(T) and MCSCF studies

Spin crossover from high spin Fe(ii)-phthalocyanine to low or intermediate spin via either dative covalent or non-covalent interaction by just varying the substituent using the same core ligand.

Accurate computed spin-state energetics for Co(iii) complexes: implications for modelling homogeneous catalysis

DLPNO-CCSD(T) calculations provide accurate spin state energetics for a range of Co(iii) complexes and so represent a promising approach to modelling homogeneous catalysis based on Co(iii) species.

Thermal spin crossover in Fe(ii) and Fe(iii). Accurate spin state energetics at the solid state

Parametrization of PBE+U under the D3 and D3-BJ dispersion corrections to study Fe^II and Fe^III-based Spin Crossover complexes.

Approximate Analytical Gradients and Nonadiabatic Couplings for the State-Average Density Matrix Renormalization Group Self-Consistent-Field Method

We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wave function. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) ansatz, which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy.

Calculating spin crossover temperatures by a first-principles LDA+U scheme with parameter U evaluated from GW

The prediction of spin crossover (SCO) temperatures (T_1/2) depends sensitively on the description of local Coulomb correlation. Due to its balance between accuracy and computational cost, local density approximation combined with Hubbard U model (LDA+U) is an appealing tool for this purpose. Despite its accurate performance on energetic properties, such as spin adiabatic energy difference, it is well-known that the LDA+U approach would lose its predictive power if U is tuned to achieve close agreement with experiment for a certain property. On the other hand, a static U value cannot account for changes in the electronic structure. Here, we propose a framework to derive dynamical U (U_dyn) values for iron(ii) complexes from the many-body GW calculations. By performing model calculations on a series of compounds with varying ligand fields, we show that the U values determined in this way are local environment dependent, and the resulting LDA+U_dyn method could reproduce their experimental ground spin states. We present applications to selected SCO complexes illustrating that U_dyn considerably overcomes some of the drawbacks of employing a constant U in the calculation of thermochemical quantities. Using the described calculation procedure, the T_1/2 values are predicted with a small mean absolute error of 176 K with respect to experiment.

Computational Study on the Mechanisms and Rate Constants for the O(3P,1D) + OCS Reactions

The mechanisms and kinetics of O(³P,¹D) + OCS(X¹Σ⁺) reactions have been studied by the high-level G2M(CC2) and CCSD(T)/6-311+G(3df)//B3LYP/6-311+G(3df) methods in conjunction with the transition-state theory and variational Rice-Ramsperger-Kassel-Marcus theory calculations. The result shows that the triplet surface proceeds directly by abstraction and substitution channels to produce SO(³P) + CO(X¹Σ⁺) and S(³P) + CO₂(X¹ Σ_g⁺) by passing the barriers of 7.6 and 9.1 kcal·mol^-1 at the G2M(CC2)//B3LYP/6-311+G(3df) level, respectively, while two stable intermediates, LM1 (OSCO¹) and LM2 (SC(O)O¹), are formed barrierlessly from O(¹D) + OCS(X¹Σ⁺) in the singlet surface, which lie at -40.5 and -50.1 kcal·mol^-1 relative to O(³P) + OCS(X¹Σ⁺) reactants and decompose to CO(X¹Σ⁺) + SO(a¹Δ) and S(¹D) + CO₂(X¹Σ_g⁺). LM1 and LM2 may also be produced by singlet-triplet surface crossings via MSX1 and MSX2; the predicted total rate constant for the O(³P) + OCS(X¹Σ⁺) reaction including the crossings, 9.2 × 10^-11 exp(-5.18 kcal·mol^-1/RT) cm³ molecule^-1 s^-1, is in good agreement with available experimental data. The branching ratio of the CO₂ product channel, 0.22-0.32, between 1200 and 1600 K, is also in excellent agreement with the value of 0.2-0.3 measured by Isshiki et al. (J. Phys. Chem. A. 2003, 107, 2464).

Electronic Structure of N-Bridged High-Valent Diiron-Oxo

Density functional theory (DFT) and an advanced ab initio technique based on density matrix renormalization group (DMRG-CASPT2) were employed to investigate a reactive N-bridged high-valent diiron-oxo species involved in H-abstraction reactions. We studied in detail two important doublet states, the ground state with two iron(IV) centers and a mixed valence Fe^V -Fe^IV excited state. We found that the latter state is low-lying. Furthermore, its electronic structure and spin density imply that it has significantly higher H-abstraction reactivity than the ground state. This low-lying excited state might be the reason behind the high oxidation reactivity of this diiron-oxo species towards methane.

Effect of Exocyclic Substituents and π-System Length on the Electronic Structure of Chichibabin Diradical(oid)s

The ground state (GS) of Chichibabin's polycyclic hydrocarbons (CPHs) can be singlet [open- or closed-shell (OSS or CS)] or triplet (T), depending on the elongation of the π-system and the exocyclic substituents. CPHs with either a small singlet-triplet energy gap (ΔE ST) or even a triplet GS have potential applications in optoelectronics. To analyze the effect of the size and exocyclic substituents on the nature of the GS of CPHs, we have selected a number of them with different substituents in the exocyclic carbon atoms and different ring chain lengths. The OPBE/cc-pVTZ level of theory was used for the optimization of the systems. The aromaticity of the resulting electronic structures was evaluated with HOMA, NICS, FLU, PDI, Iring, and MCI aromaticity indices. Our results show that the shortest π-systems (one or two rings) have a singlet GS. However, systems with three to five rings favor OSS GSs. Electron-withdrawing groups (EWGs) and aromatic substituents in the exocyclic carbons tend to stabilize the OSS and T states, whereas electron-donating groups slightly destabilize them. For CS, OSS, and T states, aromaticity measures indicate a gain of aromaticity of the 6-membered rings of the CPHs with the increase in their size and when CPHs incorporate EWGs or aromatic substituents. In general, the CPHs analyzed present small singlet-triplet energy gaps, and in particular, the ones containing EWGs or aromatic substituents present the smallest singlet-triplet energy gaps.

Low-Lying Electromeric States in Chloro-Ligated Iron(IV)-Oxo Porphyrin as a Model for Compound I, Studied with Second-Order Perturbation Theory Based on Density Matrix Renormalization Group

Employing second-order perturbation theory based on the density matrix renormalization group (DMRG-CASPT2), this work aims at providing a quantitative description of the spin state energetics of a chloro-ligated iron(IV)-oxo porphyrin as a model for the cytochromes P450 active species, also known as compound I (Cpd I). We explored DMRG-CASPT2 to its full extent with an extensive active space (up to 31 active orbitals) as well as a large number of renormalized states m (up to 10000). Different flavors of DMRG-CASPT2, using either the costly exact 4-particle reduced density matrix (4-RDM) or the cheaper cumulant approximated 4-RDM (cu(4)), were analyzed. All flavors essentially converge to similar relative energies between different spin states. Including a correction for the protein environment, we found a quartet Fe^IVO ground state and, more importantly, a thermally accessible doublet Fe^VO excited state that might directly contribute to the reactivity of this iron-oxo species. Our results also showed that cheaper approaches, such as CASPT2 based on a smaller active space or the cumulant approximation DMRG-cu(4)-CASPT2, are capable of accurately describing the spin state energetics of this species.

The S3 State of the Oxygen-Evolving Complex: Overview of Spectroscopy and XFEL Crystallography with a Critical Evaluation of Early-Onset Models for O–O Bond Formation

The catalytic cycle of the oxygen-evolving complex (OEC) of photosystem II (PSII) comprises five intermediate states Si (i = 0–4), from the most reduced S0 state to the most oxidized S4, which spontaneously evolves dioxygen. The precise geometric and electronic structure of the Si states, and hence the mechanism of O–O bond formation in the OEC, remain under investigation, particularly for the final steps of the catalytic cycle. Recent advances in protein crystallography based on X-ray free-electron lasers (XFELs) have produced new structural models for the S3 state, which indicate that two of the oxygen atoms of the inorganic Mn4CaO6 core of the OEC are in very close proximity. This has been interpreted as possible evidence for “early-onset” O–O bond formation in the S3 state, as opposed to the more widely accepted view that the O–O bond is formed in the final state of the cycle, S4. Peroxo or superoxo formation in S3 has received partial support from computational studies. Here, a brief overview is provided of spectroscopic information, recent crystallographic results, and computational models for the S3 state. Emphasis is placed on computational S3 models that involve O–O formation, which are discussed with respect to their agreement with structural information, experimental evidence from various spectroscopic studies, and substrate exchange kinetics. Despite seemingly better agreement with some of the available crystallographic interpretations for the S3 state, models that implicate early-onset O–O bond formation are hard to reconcile with the complete line of experimental evidence, especially with X-ray absorption, X-ray emission, and magnetic resonance spectroscopic observations. Specifically with respect to quantum chemical studies, the inconclusive energetics for the possible isoforms of S3 is an acute problem that is probably beyond the capabilities of standard density functional theory.