Relative stability between the manganese hydroxide- and oxo-models for water oxidation by CCSD, DMRG CASCI, CASSCF, CASPT2 and CASDFT methods; Importance of static and dynamical electron correlation effects for OEC of PSII

Theoretical elucidation of the structure, bonding, and reactivity of the CaMn4Ox clusters in the whole Kok cycle for water oxidation embedded in the oxygen evolving center of photosystem II. New molecular and quantum insights into the mechanism of the O-O bond formation

This paper reviews our historical developments of broken-symmetry (BS) and beyond BS methods that are applicable for theoretical investigations of metalloenzymes such as OEC in PSII. The BS hybrid DFT (HDFT) calculations starting from high-resolution (HR) XRD structure in the most stable S1 state have been performed to elucidate structure and bonding of whole possible intermediates of the CaMn4Ox cluster (1) in the Si (i = 0 ~ 4) states of the Kok cycle. The large-scale HDFT/MM computations starting from HR XRD have been performed to elucidate biomolecular system structures which are crucial for examination of possible water inlet and proton release pathways for water oxidation in OEC of PSII. DLPNO CCSD(T0) computations have been performed for elucidation of scope and reliability of relative energies among the intermediates by HDFT. These computations combined with EXAFS, XRD, XFEL, and EPR experimental results have elucidated the structure, bonding, and reactivity of the key intermediates, which are indispensable for understanding and explanation of the mechanism of water oxidation in OEC of PSII. Interplay between theory and experiments have elucidated important roles of four degrees of freedom, spin, charge, orbital, and nuclear motion for understanding and explanation of the chemical reactivity of 1 embedded in protein matrix, indicating the participations of the Ca(H2O)n ion and tyrosine(Yz)-O radical as a one-electron acceptor for the O-O bond formation. The Ca-assisted Yz-coupled O-O bond formation mechanisms for water oxidation are consistent with recent XES and very recent time-resolved SFX XFEL and FTIR results.

Reconciling Local Coupled Cluster with Multireference Approaches for Transition Metal Spin-State Energetics

Spin-state energetics of transition metal complexes remain one of the most challenging targets for electronic structure methods. Among single-reference wave function approaches, local correlation approximations to coupled cluster theory, most notably the domain-based local pair natural orbital (DLPNO) approach, hold the promise of bringing the accuracy of coupled cluster theory with single, double, and perturbative triple excitations, CCSD(T), to molecular systems of realistic size with acceptable computational cost. However, recent studies on spin-state energetics of iron-containing systems raised doubts about the ability of the DLPNO approach to adequately and systematically approximate energetics obtained by the reference-quality complete active space second-order perturbation theory with coupled-cluster semicore correlation, CASPT2/CC. Here, we revisit this problem using a diverse set of iron complexes and examine several aspects of the application of the DLPNO approach. We show that DLPNO-CCSD(T) can accurately reproduce both CASPT2/CC and canonical CCSD(T) results if two basic principles are followed. These include the consistent use of the improved iterative (T₁) versus the semicanonical perturbative triple corrections and, most importantly, a simple two-point extrapolation to the PNO space limit. The latter practically eliminates errors arising from the default truncation of electron-pair correlation spaces and should be viewed as standard practice in applications of the method to transition metal spin-state energetics. Our results show that reference-quality results can be readily achieved with DLPNO-CCSD(T) if these principles are followed. This is important also in view of the applicability of the method to larger single-reference systems and multinuclear clusters, whose treatment of dynamic correlation would be challenging for multireference-based approaches.

Are lanthanide-transition metal direct bonds a route to achieving new generation {3d-4f} SMMs?

Lanthanide based single-molecule magnets are gaining wide attention due to their potential applications in emerging technologies. One of the main challenges in this area is quenching quantum tunnelling of magnetisation (QTM), which often undercuts the magnetisation reversal barrier. Among the several strategies employed, enhancing exchange coupling has been studied in detail, with large exchanges resulting in stronger quenching of QTM effects. Lanthanides, however, suffer from weak exchanges offered by the deeply buried 4f orbitals and the numerous attempts to enhance the exchange coupling in the {3d-4f} pairs have not exceeded values larger than 30 cm^-1. In this work, using a combination of DFT and the ab initio CASSCF/RASSI-SO method, we have explored lanthanide-transition metal direct bonds as a tool to quench QTM effects. In this direction, we have modelled [PyCp₂LnMCp(CO)₂] (Ln = Gd(III), Dy(III), and Er(III) and M = V(0), Mn(0), Co(0) and Fe(I) and here PyCp₂ = [2,6-(CH₂C₅H₃)₂C₅H₃N]₂^- using [PyCp₂DyFeCp(CO)₂] as an example as reported by Nippe et al. (C. P. Burns, X. Yang, J. D. Wofford, N. S. Bhuvanesh, M. B. Hall and M. Nippe, Angew. Chem., Int. Ed. 2018, 57, 8144). Bonding analysis reveals a dative Ln-TM bond with a donation of π(V/Mnd_xy-π*CO) to 5d_z² (Gd) in the case of Gd-V and Gd-Mn and 4s(Co) to 5d_xy/5d_yz (Gd) for Gd-Co with the transition metal ion being found in the low-spin S = ½ configurations in all the cases. B3LYP/TZV (Gd;CSDZ) calculations on [PyCp₂GdMCp(CO)₂] yield J_Gd-V = -46.1 cm^-1, J_Gd-Mn = -57.1 cm^-1, J_Gd-Co = +55.3 cm^-1, J_Gd-Fe⁺ = +13.9 cm^-1, J_Gd-Vhs = -162.1 cm^-1 and J_Gd-Mnhs = -343.9 cm^-1 and unveiling record-high J values for {3d-4f} complexes. The mechanism of magnetic coupling is developed, which discloses the dominating presence of strong 3d-4f orbital overlaps in most of the cases studied, leading to antiferromagnetic exchange. When these overlaps are weaker and 3d to Gd(5d_z²), charge transfer dominates, yielding a ferromagnetic coupling for the Gd-Co/Gd-Fe⁺ complexes. Calculations performed on the anisotropic Dy(III) and Er(III) complexes reveal that the ground state g_zz axis lies along the Cp-Ln-Cp axis and the Ln-TM bonds, respectively. Thus the Ln-TM bond hinders the single-ion anisotropy of Dy(III) by offering equatorial ligation and lowering the m_J = ±½ state energy, and at the same time, helping in enhancing the axiality of Er(III). When strong {3d-4f} exchange couplings are introduced, record-high barrier heights as high as 229 cm^-1 were accomplished. Furthermore, the exchange coupling annihilates the QTM effects and suggests the lanthanide-transition metal direct bond as a viable alternative to enhance exchange coupling to bring {3d-4f} complexes back in the race for high-blocking SMMs.

Spin, Valence, and Structural Isomerism in the S3 State of the Oxygen-Evolving Complex of Photosystem II as a Manifestation of Multimetallic Cooperativity

Photosynthetic water oxidation is catalyzed by a Mn₄CaO₅-cluster in photosystem II through an S-state cycle. Understanding the roles of heterogeneity in each S-state, as identified recently by the EPR spectroscopy, is very important to gain a complete description of the catalytic mechanism. We performed herein hybrid DFT calculations within the broken-symmetry formalism and associated analyses of Heisenberg spin models to study the electronic and spin structures of various isomeric structural motifs (hydroxo-oxo, oxyl-oxo, peroxo, and superoxo species) in the S₃ state. Our extensive study reveals several factors that affect the spin ground state: (1) (formal) Mn oxidation state; (2) metal-ligand covalency; (3) coordination geometry; and (4) structural change of the Mn cluster induced by alternations in Mn···Mn distances. Some combination of these effects could selectively stabilize/destabilize some spin states. We found that the high spin state ( S_total = 6) of the oxyl-oxo species can be causative for catalytic function, which manifests through mixing of the metal-ligand character in magnetic orbitals at relatively short O5···O6 distances (<2.0 Å) and long Mn_A···O5 distances (>2.0 Å). These results will serve as a basis to conceptually identify and rationalize the physicochemical synergisms that can be evoked by the unique "distorted chair" topology of the cluster through cooperative Jahn-Teller effects on multimetallic centers.