FeS clusters are a universal biological motif. They carry out electron transfer, redox chemistry, and even oxygen sensing, in diverse processes including nitrogen fixation, respiration, and photosynthesis. The low-lying electronic states are key to their remarkable reactivity, but cannot be directly observed.Here we present the first ever quantum calculation of the electronic levels of [2Fe-2S] and [4Fe-4S] clusters free from any model assumptions. Our results highlight limitations of long-standing models of their electronic structure. In particular, we demonstrate that the widely used Heisenberg-Double-Exchange model underestimates the number of states by 1-2 orders of magnitude, which can conclusively be traced to the absence of Fe d→d excitations, thought to be important in these clusters. Further, the electronic energy levels of even the same spin are dense on the scale of vibrational fluctuations, and this provides a natural explanation for the ubiquity of these clusters in nature for catalyzing reactions.

A theoretical, computational, and conceptual framework for the interpretation and prediction of the magnetic anisotropy of transition metal complexes with orbitally degenerate or orbitally nearly degenerate ground states is explored. The treatment is based on complete active space self-consistent field (CASSCF) wave functions in conjunction with N-electron valence perturbation theory (NEVPT2) and quasidegenerate perturbation theory (QDPT) for treatment of magnetic field- and spin-dependent relativistic effects. The methodology is applied to a series of Fe(II) complexes in ligand fields of almost trigonal pyramidal symmetry as provided by several variants of the tris-pyrrolylmethyl amine ligand (tpa). These systems have recently attracted much attention as mononuclear single-molecule magnet (SMM) complexes. This study aims to establish how the ligand field can be fine tuned in order to maximize the magnetic anisotropy barrier. In trigonal ligand fields high-spin Fe(II) complexes adopt an orbitally degenerate 5 E ground state with strong in-state spin–orbit coupling (SOC). We study the competing effects of SOC and the 5 E⊗ε multimode Jahn–Teller effect as a function of the peripheral substituents on the tpa ligand. These subtle distortions were found to have a significant effect on the magnetic anisotropy. Using a rigorous treatment of all spin multiplets arising from the triplet and quintet states in the d6 configuration the parameters of the effective spin-Hamiltonian (SH) approach were predicted from first principles. Being based on a nonperturbative approach we investigate under which conditions the SH approach is valid and what terms need to be retained. It is demonstrated that already tiny geometric distortions observed in the crystal structures of four structurally and magnetically well-documented systems, reported recently, i.e., [Fe(tpaR)]− (R = tert-butyl, Tbu (1), mesityl, Mes (2), phenyl, Ph (3), and 2,6-difluorophenyl, Dfp (4), are enough to lead to five lowest and thermally accessible spin sublevels described sufficiently well by S = 2 SH provided that it is extended with one fourth order anisotropy term. Using this most elementary parametrization that is consistent with the actual physics, the reported magnetization data for the target systems were reinterpreted and found to be in good agreement with the ab initio results. The multiplet energies from the ab initio calculations have been fitted with remarkable consistency using a ligand field (angular overlap) model (ab initio ligand field, AILFT). This allows for determination of bonding parameters and quantitatively demonstrates the correlation between increasingly negative D values and changes in the σ-bond strength induced by the peripheral ligands. In fact, the sigma-bonding capacity (and hence the Lewis basicity) of the ligand decreases along the series 1 > 2 > 3 > 4.

No abstract

The multireference n-electron Valence State Perturbation Theory is applied to a benchmark set of 28 organic molecules compiled by Schreiber et al. J. Chem. Phys. (2008) 128, 13. Different types of low-lying vertical excitation energies are computed using the same geometries and TZVP basis set as in the original work. The previously published coupled cluster CC3 results are used as a reference. The complete active space second order perturbation theory (CASPT2) results, as well as the results of second order N-electron valence perturbation theory (NEVPT2) (both in their single-state variants) are evaluated against this reference set, which includes 153 singlet and 72 triplet vertical transition energies. NEVPT2 calculations are carried out in two variants: the partially contracted (PC) and the strongly contracted (SC) scheme. The statistical evaluation with respect to CC3 is found to be similar for both: the mean unsigned deviations is 0.28 eV for singlets and 0.16 eV for triplets for PC-NEVPT2, while it is 0.23 and 0.17 eV for SC-NEVPT2, respectively. Further analysis has shown that deficiencies in the zeroth-order wave functions, in particular for the subset of π → π* singlet excitations, are responsible for the largest deviations from CC3. Those states have either a charge transfer or an ionic character. For the remaining singlet and all triplet excitations the general trend was established that NEVPT2 tends to slightly overestimate excitation energies while CASPT2 slightly underestimates them. However, overall, both methods are of very similar accuracy provided that the IPEA shift is used in the CASPT2 method. Interestingly, the conclusions reached in this study are independent of the orbital canonicalization scheme used in the NEVPT2 calculation.

This paper presents a detailed evaluation of the performance of density functional theory (DFT) as well as complete active space self-consistent field (CASSCF)-based methods (CASSCF and second-order N-electron valence state perturbation theory, NEVPT2) to predict the zero-field splitting (zfs) parameters for a series of coordination complexes containing the Mn(III) ion. The physical origin of the experimentally determined zfs's was investigated by studying the different contributions to these parameters. To this end, a series of mononuclear Mn(III) complexes was chosen for which the structures have been resolved by X-ray diffraction and the zfs parameters have been accurately determined by high-field EPR spectroscopy. In a second step, small models have been constructed to allow for a systematic assessment of the factors that dominate the variations in the observed zfs parameters and to establish magnetostructural correlations. Among the tested functionals, the best predictions have been obtained with B3LYP, followed by the nonhybrid BP86 functional, which in turn is more successful than the meta-hybrid GGA functional TPSSh. For the estimation of the spin-orbit coupling (SOC) part of the zfs, it was found that the coupled perturbed SOC approach CP is more successful than the Pederson-Khanna method. Concerning the spin-spin interaction (SS), the restricted open-shell Kohn-Sham (ROKS) approach led to a slightly better agreement with the experiment than the unrestricted KS (UKS) approach. The ab initio state-averaged CASSCF (SA-CASSCF) method with a minimal active space and the most recent implementation that treats the SOC and SS contributions on an equal footing provides the best predictions for the zfs. The analysis demonstrates that the major contribution to the axial zfs parameter (D) originates from the SOC interaction but that the SS part is far from being negligible (between 10 and 20% of D). Importantly, the various excited triplet ligand field states account for roughly half of the value of D, contrary to popular ligand field models. Despite covering dynamic correlation contributions to the transition energies, NEVPT2 does not lead to large improvements in the results as the excitation energies of the Mn(III) d-d transitions are already fairly accurate at the SA-CASSCF level. For a given type of coordination sphere (e.g., elongated or compressed octahedron), the magnetic anisotropy of the Mn(III) ion, D, does not appear to be highly sensitive to the nature of the ligands, while the E/D ratio is notably affected by all octahedral distortions. Furthermore, the introduction of different halides into the coordination sphere of Mn(III) only leads to small effects on D. Nevertheless, it appears that oxygen-based ligands afford larger D values than nitrogen-based ligands.