Domain-based local pair natural orbital CCSD(T) calculations of strongly correlated electron systems: Examination of dynamic equilibrium models based on multiple intermediates in S1 state of photosystem II

Redox Isomerism in the S3 State of the Oxygen-Evolving Complex Resolved by Coupled Cluster Theory

The electronic and geometric structures of the water-oxidizing complex of photosystem II in the steps of the catalytic cycle that precede dioxygen evolution remain hotly debated. Recent structural and spectroscopic investigations support contradictory redox formulations for the active-site Mn4 CaOx cofactor in the final metastable S3 state. These range from the widely accepted MnIV 4 oxo-hydroxo model, which presumes that O-O bond formation occurs in the ultimate transient intermediate (S4 ) of the catalytic cycle, to a MnIII 2 MnIV 2 peroxo model representative of the contrasting "early-onset" O-O bond formation hypothesis. Density functional theory energetics of suggested S3 redox isomers are inconclusive because of extreme functional dependence. Here, we use the power of the domain-based local pair natural orbital approach to coupled cluster theory, DLPNO-CCSD(T), to present the first correlated wave function theory calculations of relative stabilities for distinct redox-isomeric forms of the S3 state. Our results enabled us to evaluate conflicting models for the S3 state of the oxygen-evolving complex (OEC) and to quantify the accuracy of lower-level theoretical approaches. Our assessment of the relevance of distinct redox-isomeric forms for the mechanism of biological water oxidation strongly disfavors the scenario of early-onset O-O formation advanced by literal interpretations of certain crystallographic models. This work serves as a case study in the application of modern coupled cluster implementations to redox isomerism problems in oligonuclear transition metal systems.

Electronic and spin structures of CaMn4Ox clusters in the S0 state of the oxygen evolving complex of photosystem II. Domain-based local pair natural orbital (DLPNO) coupled-cluster (CC) calculations using optimized geometries and natural orbitals (UNO) by hybrid density functional theory (HDFT) calculations

Domain-based local pair natural orbital (DLPNO) coupled cluster single and double (CCSD) with triple perturbation (T) correction methods were performed to elucidate the relative stabilities of ten different intermediate structures of the CaMn4Ox cluster in the S0 state of the oxygen evolving complex (OEC) of photosystem II (PSII). Full geometry optimizations of all the S0 intermediates were performed by the UB3LYP-D3/Def2-TZVP methods, providing the assumed geometrical structures and starting natural orbitals (UNO) for DLPNO-CCSD(T)/Def2TZVP calculations. The effective exchange integrals (J) for the spin Hamiltonian models for the ten intermediates were obtained by the UB3LYP/Def2-TZVP calculations followed by the general spin projections. DLPNO-CCSD(T) calculations followed by the CBS extrapolation procedure elucidated that the (II, III, IV, IV) and (III, III, III, IV) valence states in the CaMn4O5 cluster of the OEC of the PS II were nearly degenerated in energy in the S0 state, indicating an important role of dynamical electron correlation effects for the valence and spin fluctuations in strongly correlated electron systems (SCESs) consisting of 3d transition metals.