Hyper Open-Shell Excited Spin States of Transition-Metal Compounds: FeF2, FeF2···Ethane, and FeF2···Ethylene

Local Oxidation States in {FeNO}6-8 Porphyrins: Insights from DMRG/CASSCF-CASPT2 Calculations

A first DMRG/CASSCF-CASPT2 study of a series of paradigmatic {FeNO}6, {FeNO}7, and {FeNO}8 heme-nitrosyl complexes has led to substantial new insight as well as uncovered key shortcomings of the DFT approach. By virtue of its balanced treatment of static and dynamic correlation, the calculations have provided some of the most authoritative information available to date on the energetics of low- versus high-spin states of different classes of heme-nitrosyl complexes. Thus, the calculations indicate low doublet-quartet gaps of 1-4 kcal/mol for {FeNO}7 complexes and high singlet-triplet gaps of ≳20 kcal/mol for both {FeNO}6 and {FeNO}8 complexes. In contrast, DFT calculations yield widely divergent spin state gaps as a function of the exchange-correlation functional. DMRG-CASSCF calculations also help calibrate DFT spin densities for {FeNO}7 complexes, pointing to those obtained from classic pure functionals as the most accurate. The general picture appears to be that nearly all the spin density of Fe[P](NO) is localized on the Fe, while the axial ligand imidazole (ImH) in Fe[P](NO)(ImH) pushes a part of the spin density onto the NO moiety. An analysis of the DMRG-CASSCF wave function in terms of localized orbitals and of the resulting configuration state functions in terms of resonance forms with varying NO(π*) occupancies has allowed us to address the longstanding question of local oxidation states in heme-nitrosyl complexes. The analysis indicates NO(neutral) resonance forms [i.e., Fe(II)-NO0 and Fe(III)-NO0] as the major contributors to both {FeNO}6 and {FeNO}7 complexes. This finding is at variance with the common formulation of {FeNO}6 hemes as Fe(II)-NO+ species but is consonant with an Fe L-edge XAS analysis by Solomon and co-workers. For the {FeNO}8 complex {Fe[P](NO)}-, our analysis suggests a resonance hybrid description: Fe(I)-NO0 ↔ Fe(II)-NO-, in agreement with earlier DFT studies. Vibrational analyses of the compounds studied indicate an imperfect but fair correlation between the NO stretching frequency and NO(π*) occupancy, highlighting the usefulness of vibrational data as a preliminary indicator of the NO oxidation state.

Benchmarking CASPT3 Vertical Excitation Energies

Based on 280 reference vertical transition energies of various natures (singlet, triplet, valence, Rydberg, n → π∗, π → π∗, and double excitations) extracted from the QUEST database, we assess the accuracy of third-order multireference perturbation theory, CASPT3, in the context of molecular excited states. When one applies the disputable ionization- potential-electron-affinity (IPEA) shift, we show that CASPT3 provides a similar accuracy as its second-order counterpart, CASPT2, with the same mean absolute error of 0.11 eV. However, as already reported, we also observe that the accuracy of CASPT3 is almost insensitive to the IPEA shift, irrespective of the transition type and system size, with a small reduction of the mean absolute error to 0.09 eV when the IPEA shift is switched off.

Spin-Orbit Coupling Changes the Identity of the Hyper-Open-Shell Ground State of Ce+, and the Bond Dissociation Energy of CeH+ Proves to Be Challenging for Theory

Cerium (Ce) plays important roles in catalysis. Its position in the sixth period of the periodic table leads to spin-orbit coupling (SOC) and other open-shell effects that make the quantum mechanical calculation of cerium compounds challenging. In this work, we investigated the low-lying spin states of Ce⁺ and the bond energy of CeH⁺, both by multiconfigurational methods, in particular, SA-CASSCF, MC-PDFT, CASPT2, XMS-PDFT, and XMS-CASPT2, and by single-configurational methods, namely, Hartree-Fock theory and unrestricted Kohn-Sham density functional theory with 34 choices of the exchange-correlation functional. We found that only CASPT2, XMS-CASPT2, and SA-CASSCF (among the five multiconfigurational methods) and GAM, HCTH, SOGGA11, and OreLYP (among the 35 single-configuration methods) successfully predict that the SOC-free ground spin state of Ce⁺ is a doublet state, and CASPT2 and GAM give the most accurate multireference and single-reference calculations, respectively, of the excitation energy of the first SOC-free excited state for Ce⁺. We calculated that the ground doublet state of Ce⁺ is an intra-atomic hyper-open-shell state. We calculated the spin-orbit energy (E^SO) of Ce⁺ by the five multiconfigurational methods and found that E^SO calculated by CASPT2 is the closest to the experimental value. Taking advantage of the availability of an experimental D₀ for CeH⁺ as a way to provide a unique test of theory, we showed that all the multiconfigurational methods overestimate D₀ by at least 246 meV (5.7 kcal/mol), and only three functionals, namely, SOGGA, MN15, and GAM, have an error of D₀ that is less than 200 meV (5 kcal/mol).

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.

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.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Is the Inversion of Phosphorus Trihalides (PF3, PCl3, PBr3, and PI3) a Diradical Process?

This work explores possible reaction paths for the inversion of a series of trigonal pyramidal phosphorus trihalides, PF₃, PCl₃, PBr₃, and PI₃, and it especially addresses the question of whether and when the bonding of the lowest-energy species along the inversion paths should be described as a hyper-open-shell diradical. The various paths for inversion are calculated using a single-reference method within the framework of Kohn-Sham density functional theory and also with multireference wave function methods. Our calculated results using both kinds of methods show that, for all the halogens studied (F, Cl, Br, and I), the lowest-energy singlet path for the inversion occurs by the formation of a C_{2 v} transition structure rather than a D_{3 h} transition structure. This geometrical preference agrees with what has been inferred previously based on closed-shell singlet calculations. But in the present study, we examined not only closed-shell singlet transition states but also open-shell singlet states and triplet states for calculating stationary points and inversion paths, and for some of the phosphorus trihalides, we found that paths involving open-shell configurations are lower in energy than those restricted to closed-shell configurations. We analyzed the changes along the paths in terms of hybridization and orientation of the frontier orbitals and in terms of locally avoided crossings, and the extent of the diradical character was quantified by calculating the effective number of unpaired electrons. Even for the singlet inversion path that goes via a D_{3 h} structure, the barrier for PF₃, PCl₃, and PBr₃ is higher for a closed-shell singlet spin state than for the open-shell singlet configuration. Furthermore, the energy of the triplet D_{3 h} structure is below even that of the open-shell D_{3 h} singlet for PCl₃, PBr₃, and PI₃. This necessitates rethinking the role of open-shell states in nominally closed-shell processes.