Which hybrid GGA DFT is suitable for Cu2O2 systems if the spin contamination error is removed?

Geometric and electronic structures of the synthetic Mn₄CaO₄ model compound mimicking the photosynthetic oxygen-evolving complex

Water oxidation by photosystem II (PSII) converts light energy into chemical energy with the concomitant production of molecular oxygen, both of which are indispensable for sustaining life on Earth. This reaction is catalyzed by an oxygen-evolving complex (OEC) embedded in the huge PSII complex, and its mechanism remains elusive in spite of the extensive studies of the geometric and electronic structures. In order to elucidate the water-splitting mechanism, synthetic approaches have been extensively employed to mimic the native OEC. Very recently, a synthetic complex [Mn4CaO4(Bu(t)COO)8(py)(Bu(t)COOH)2] (1) closely mimicking the structure of the native OEC was obtained. In this study, we extensively examined the geometric, electronic and spin structures of 1 using the density functional theory method. Our results showed that the geometric structure of 1 can be accurately reproduced by theoretical calculations, and revealed many similarities in the ground valence and spin states between 1 and the native OEC. We also revealed two different valence states in the one-electron oxidized state of 1 (corresponding to the S2 state), which lie in the lower and higher ground spin states (S = 1/2 and S = 5/2), respectively. One remarkable difference between 1 and the native OEC is the presence of a non-negligible antiferromagnetic interaction between the Mn1 and Mn4 sites, which slightly influenced their ground spin structures (spin alignments). The major reason causing the difference can be attributed to the short Mn1-O5 and Mn1-Mn4 distances in 1. The introduction of the missing O4 atom and the reorientation of the Ca coordinating ligands improved the Mn1-O5 and Mn1-Mn4 distances comparable to the native OEC. These modifications will therefore be important for the synthesis of further advanced model complexes more closely mimicking the native OEC beyond 1.

B2N2O4: Prediction of a Magnetic Ground State for a Light Main-Group Molecule

Cyclobutanetetrone, (CO)4, has a triplet ground state. Here we predict, based on electronic structure calculations, that the B2N2O4 molecule also has a triplet ground state and is therefore paramagnetic; the structure is an analogue of (CO)4 in which the carbon ring is replaced by a (BN)2 ring. Similar to (CO)4, the triplet ground-state structure of B2N2O4 is also thermodynamically unstable. Besides analysis of the molecular orbitals, we found that the partial atomic charges are good indicators for predicting magnetic ground states.

Computational Study of Catalytic Reaction of Quercetin 2,4-Dioxygenase

We present a quantum mechanics/molecular mechanics (QM/MM) and QM-only study on the oxidative ring-cleaving reaction of quercetin catalyzed by quercetin 2,4-dioxygenase (2,4-QD). 2,4-QD has a mononuclear type 2 copper center and incorporates two oxygen atoms at C2 and C4 positions of the substrate. It has not been clear whether dioxygen reacts with a copper ion or a substrate radical as the first step. We have found that dioxygen is more likely to bind to a Cu(2+) ion, involving the dissociation of the substrate from the copper ion. Then a Cu(2+)-alkylperoxo complex can be generated. Comparison of geometry and stability between QM-only and QM/MM results strongly indicates that steric effects of the protein environment contribute to maintain the orientation of the substrate dissociated from the copper center. The present QM/MM results also highlight that a prior rearrangement of the Cu(2+)-alkylperoxo complex and a subsequent hydrogen bond switching assisted by the movement of Glu73 can facilitate formation of an endoperoxide intermediate selectively.

Adaptive multiconfigurational wave functions

A method is suggested to build simple multiconfigurational wave functions specified uniquely by an energy cutoff Λ. These are constructed from a model space containing determinants with energy relative to that of the most stable determinant no greater than Λ. The resulting Λ-CI wave function is adaptive, being able to represent both single-reference and multireference electronic states. We also consider a more compact wave function parameterization (Λ+SD-CI), which is based on a small Λ-CI reference and adds a selection of all the singly and doubly excited determinants generated from it. We report two heuristic algorithms to build Λ-CI wave functions. The first is based on an approximate prescreening of the full configuration interaction space, while the second performs a breadth-first search coupled with pruning. The Λ-CI and Λ+SD-CI approaches are used to compute the dissociation curve of N2 and the potential energy curves for the first three singlet states of C2. Special attention is paid to the issue of energy discontinuities caused by changes in the size of the Λ-CI wave function along the potential energy curve. This problem is shown to be solvable by smoothing the matrix elements of the Hamiltonian. Our last example, involving the Cu2O2(2+) core, illustrates an alternative use of the Λ-CI method: as a tool to both estimate the multireference character of a wave function and to create a compact model space to be used in subsequent high-level multireference coupled cluster computations.

Quantum mechanics/molecular mechanics study of oxygen binding in hemocyanin

We report a combined quantum mechanics/molecular mechanics (QM/MM) study on the mechanism of reversible dioxygen binding in the active site of hemocyanin (Hc). The QM region is treated by broken-symmetry density functional theory (DFT) with spin projection corrections. The X-ray structures of deoxygenated (deoxyHc) and oxygenated (oxyHc) hemocyanin are well reproduced by QM/MM geometry optimizations. The computed relative energies strongly depend on the chosen density functional. They are consistent with the available thermodynamic data for oxygen binding in hemocyanin and in synthetic model complexes when the BH&HLYP hybrid functional with 50% Hartree-Fock exchange is used. According to the QM(BH&HLYP)/MM results, the reaction proceeds stepwise with two sequential electron transfer (ET) processes in the triplet state followed by an intersystem crossing to the singlet product. The first ET step leads to a nonbridged superoxo CuB(II)-O2(•-) intermediate via a low-barrier transition state. The second ET step is even more facile and yields a side-on oxyHc complex with the characteristic Cu2O2 butterfly core, accompanied by triplet-singlet intersystem crossing. The computed barriers are very small so that the two ET processes are expected to very rapid and nearly simultaneous.

Unrestricted prescriptions for open-shell singlet diradicals: using economical ab initio and density functional theory to calculate singlet-triplet gaps and bond dissociation curves

Here we present and test several computational prescriptions for calculating singlet-triplet (ST) gap energies and bond dissociation curves for open-shell singlet diradicals using economical unrestricted single reference type calculations. For ST gap energies from Slipchenko and Krylov's atom and molecule test set (C, O, Si, NH, NF, OH(+), O(2), CH(2), and NH(2)(+)) spin unrestricted Hartree-Fock and MP2 energies result in errors greater than 15 kcal/mol. However, spin-projected (SP) Hartree-Fock theory in combination with spin-component-scaled (SCS) or scaled-opposite-spin (SOS) second-order perturbation theory gives ST gap energies with a mean unsigned error (MUE) of less than 2 kcal/mol. Density functionals generally give poor results for unrestricted energies and only the ωB97X-D, the M06, and the M06-2X functionals provide reasonable accuracy after spin-projection with MUE values of 4.7, 4.3, and 3.0 kcal/mol, respectively, with the 6-311++G(2d,2p) basis set. We also present a new one parameter hybrid density functional, diradical-1 (DR-1), based on Adamo and Barone's modified PW exchange functional with the PW91 correlation functional. This DR-1 method gives a mean error (ME) of 0.0 kcal/mol and a MUE value of 1.3 kcal/mol for ST gap energies. As another test of unrestricted methods the bond dissociation curves for methane (CH(4)) and hydrofluoric acid (H-F) were calculated with the M06-2X, DR-1, and ωB97X-D density functionals. All three of these functionals give reasonable results for the methane C-H bond but result in errors greater than 50 kcal/mol for the H-F bond dissociation. Spin-projection is found to significantly degrade bond dissociation curves past ~2.2 Å. Although unrestricted Hartree-Fock theory provides a very poor description of H-F bond dissociation, unrestricted SCS-MP2 and SOS-MP2 methods give accurate results.