Theoretical study on the mechanism of water oxidation catalyzed by a mononuclear copper complex: important roles of a redox non-innocent ligand and HPO4 2- anion

The water oxidation reaction is the bottleneck problem of the artificial photosynthetic system. In this work, the mechanism of water oxidation catalyzed by a mononuclear copper complex in alkaline conditions was studied by density functional calculations. Firstly, a water molecule coordinating with the copper center of the complex (Cuii, 1) generates Cuii-H2O (2). 2 undergoes two proton-coupled electron transfer processes to produce intermediate (4). The oxidation process occurs mainly on the ligand moiety, and 4 (˙L-Cuii-O˙) can be described as a Cuii center interacting with a ligand radical antiferromagnetically and an oxyl radical ferromagnetically. 4 is the active species that can trigger O-O bond formation via the water nucleophilic attack mechanism. This process occurs in a step-wise manner. The attacking water transfers one of the protons to the HPO4 2- coupled with an electron transfer to the ligand radical, which generates a transient OH˙ interacting with the oxyl radical and H2PO4 -. Then the O-O bond is formed through the direct coupling of the oxo radical and the OH radical. The triplet di-oxygen could be released after two oxidation processes. According to the Gibbs free energy diagram, the O-O bond formation was suggested to be the rate-limiting step with a calculated total barrier of 19.5 kcal mol-1. More importantly, the copper complex catalyzing water oxidation with the help of a redox non-innocent ligand and HPO4 2- was emphasized.

Excited-State DMRG Made Simple with FEAST

We introduce DMRG[FEAST], a new method for optimizing excited-state many-body wave functions with the density matrix renormalization group (DMRG) algorithm. Our approach applies the FEAST algorithm, originally designed for large-scale diagonalization problems, to matrix product state wave functions. We show that DMRG[FEAST] enables the stable optimization of both low- and high-energy eigenstates, therefore overcoming the limitations of state-of-the-art excited-state DMRG algorithms. We demonstrate the reliability of DMRG[FEAST] by calculating anharmonic vibrational excitation energies of molecules with up to 30 fully coupled degrees of freedom.

Modern multireference methods and their application in transition metal chemistry

Transition metal chemistry is a challenging playground for quantum chemical methods owing to the simultaneous presence of static and dynamic electron correlation effects in many systems. Wavefunction based multireference (MR) methods constitute a physically sound and systematically improvable Ansatz to deal with this complexity but suffer from some conceptual difficulties and high computational costs. The latter problem partially arises from the unfavorable scaling of the Full Configuration Interaction (Full-CI) problem which in the majority of MR methods is solved for a subset of the molecular orbital space, the so-called active space. In the last years multiple methods such as modern variants of selected CI, Full-CI Quantum Monte Carlo (FCIQMC) and the density matrix renormalization group (DMRG) have been developed that solve the Full-CI problem approximately for a fraction of the computational cost required by conventional techniques thereby significantly extending the range of applicability of modern MR methods. This perspective review outlines recent advancements in the field of MR electronic structure methods together with the resulting chances and challenges for theoretical studies in the field of transition metal chemistry. In light of its emerging importance a special focus is put on the selection of adequate active spaces and the concomitant development of numerous selection aides in recent years.

Computational Modeling of Cobalt-Based Water Oxidation: Current Status and Future Challenges

A lot of effort is nowadays put into the development of novel water oxidation catalysts. In this context, mechanistic studies are crucial in order to elucidate the reaction mechanisms governing this complex process, new design paradigms and strategies how to improve the stability and efficiency of those catalysts. This review is focused on recent theoretical mechanistic studies in the field of homogeneous cobalt-based water oxidation catalysts. In the first part, computational methodologies and protocols are summarized and evaluated on the basis of their applicability toward real catalytic or smaller model systems, whereby special emphasis is laid on the choice of an appropriate model system. In the second part, an overview of mechanistic studies is presented, from which conceptual guidelines are drawn on how to approach novel studies of catalysts and how to further develop the field of computational modeling of water oxidation reactions.

Exchange Coupling Interactions from the Density Matrix Renormalization Group and N-Electron Valence Perturbation Theory: Application to a Biomimetic Mixed-Valence Manganese Complex

The accurate description of magnetic level energetics in oligonuclear exchange-coupled transition-metal complexes remains a formidable challenge for quantum chemistry. The density matrix renormalization group (DMRG) brings such systems for the first time easily within reach of multireference wave function methods by enabling the use of unprecedentedly large active spaces. But does this guarantee systematic improvement in predictive ability and, if so, under which conditions? We identify operational parameters in the use of DMRG using as a test system an experimentally characterized mixed-valence bis-μ-oxo/μ-acetato Mn(III,IV) dimer, a model for the oxygen-evolving complex of photosystem II. A complete active space of all metal 3d and bridge 2p orbitals proved to be the smallest meaningful starting point; this is readily accessible with DMRG and greatly improves on the unrealistic metal-only configuration interaction or complete active space self-consistent field (CASSCF) values. Orbital optimization is critical for stabilizing the antiferromagnetic state, while a state-averaged approach over all spin states involved is required to avoid artificial deviations from isotropic behavior that are associated with state-specific calculations. Selective inclusion of localized orbital subspaces enables probing the relative contributions of different ligands and distinct superexchange pathways. Overall, however, full-valence DMRG-CASSCF calculations fall short of providing a quantitative description of the exchange coupling owing to insufficient recovery of dynamic correlation. Quantitatively accurate results can be achieved through a DMRG implementation of second order N-electron valence perturbation theory (NEVPT2) in conjunction with a full-valence metal and ligand active space. Perspectives for future applications of DMRG-CASSCF/NEVPT2 to exchange coupling in oligonuclear clusters are discussed.

Communication: DMRG-SCF study of the singlet, triplet, and quintet states of oxo-Mn(Salen)

We use CheMPS2, our free open-source spin-adapted implementation of the density matrix renormalization group (DMRG) [S. Wouters, W. Poelmans, P. W. Ayers, and D. Van Neck, Comput. Phys. Commun. 185, 1501 (2014)], to study the lowest singlet, triplet, and quintet states of the oxo-Mn(Salen) complex. We describe how an initial approximate DMRG calculation in a large active space around the Fermi level can be used to obtain a good set of starting orbitals for subsequent complete-active-space or DMRG self-consistent field calculations. This procedure mitigates the need for a localization procedure, followed by a manual selection of the active space. Per multiplicity, the same active space of 28 electrons in 22 orbitals (28e, 22o) is obtained with the 6-31G(*), cc-pVDZ, and ANO-RCC-VDZP basis sets (the latter with DKH2 scalar relativistic corrections). Our calculations provide new insight into the electronic structure of the quintet.