No abstract

In this work, an implementation of an approach to calculate the zero-field splitting (ZFS) constants in the framework of ab initio methods such as complete active space self-consistent field, multireference configuration interaction, or spectroscopy oriented configuration interaction is reported. The spin-orbit coupling (SOC) contribution to ZFSs is computed using an accurate multicenter mean-field approximation for the Breit-Pauli Hamiltonian. The SOC parts of ZFS constants are obtained directly after diagonalization of the SOC operator in the basis of a preselected number of roots of the spin-free Hamiltonian. This corresponds to an infinite order treatment of the SOC in terms of perturbation theory. The spin-spin (SS) part is presently estimated in a mean-field fashion and appears to yield results close to the more complete treatments available in the literature. Test calculations for the first- and second-row atoms as well as first-row transition metal atoms and a set of diatomic molecules show accurate results for the SOC part of ZFSs. SS contributions have been found to be relatively small but not negligible (exceeding 1 cm(-1) for oxygen molecule). At least for the systems studied in this work, it is demonstrated that the presented method provides much more accurate estimations for the SOC part of ZFS constants than the emerging density functional theory approaches.

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

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.