Dimitrios A. Pantazis

All-Electron Scalar Relativistic Basis Sets for Third-Row Transition Metal Atoms.

A family of segmented all-electron relativistically contracted (SARC) basis sets for the elements Hf-Hg is constructed for use in conjunction with the Douglas-Kroll-Hess (DKH) and zeroth-order regular approximation (ZORA) scalar relativistic Hamiltonians. The SARC basis sets are loosely contracted and thus offer computational advantages compared to generally contracted relativistic basis sets, while their sufficiently small size allows them to be used in place of effective core potentials (ECPs) for routine studies of molecules. Practical assessments of the SARC basis sets in DFT calculations of atomic (ionization energies) as well as molecular properties (geometries and bond dissociation energies for MHn complexes) confirm that the basis sets yield accurate and reliable results, providing a balanced description of core and valence electron densities. CCSD(T) calculations on a series of gold diatomic compounds also demonstrate the applicability of the basis sets to correlated methods. The SARC basis sets will be of most utility in calculating molecular properties for which the core electrons cannot be neglected, such as studies of electron paramagnetic resonance, Mössbauer and X-ray absorption spectra, and topological analysis of electron densities.

Biological Water Oxidation

Photosystem II (PSII), a multisubunit pigment-protein supercomplex found in cyanobacteria, algae, and plants, catalyzes a unique reaction in nature: the light-driven oxidation of water. Remarkable recent advances in the structural analysis of PSII now give a detailed picture of the static supercomplex on the molecular level. These data provide a solid foundation for future functional studies, in particular the mechanism of water oxidation and oxygen release. The catalytic core of the PSII is a tetramanganese-calcium cluster (Mn₄O₅Ca), commonly referred to as the oxygen-evolving complex (OEC). The function of the OEC rests on its ability to cycle through five metastable states (Si, i = 0-4), transiently storing four oxidizing equivalents, and in so doing, facilitates the four electron water splitting reaction. While the latest crystallographic model of PSII gives an atomic picture of the OEC, the exact connectivity within the inorganic core and the S-state(s) that the X-ray model represents remain uncertain. In this Account, we describe our joint experimental and theoretical efforts to eliminate these ambiguities by combining the X-ray data with spectroscopic constraints and introducing computational modeling. We are developing quantum chemical methods to predict electron paramagnetic resonance (EPR) parameters for transition metal clusters, especially focusing on spin-projection approaches combined with density functional theory (DFT) calculations. We aim to resolve the geometric and electronic structures of all S-states, correlating their structural features with spectroscopic observations to elucidate reactivity. The sequence of manganese oxidations and concomitant charge compensation events via proton transfer allow us to rationalize the multielectron S-state cycle. EPR spectroscopy combined with theoretical calculations provides a unique window into the tetramangenese complex, in particular its protonation states and metal ligand sphere evolution, far beyond the scope of static techniques such as X-ray crystallography. This approach has led, for example, to a detailed understanding of the EPR signals in the S₂-state of the OEC in terms of two interconvertible, isoenergetic structures. These two structures differ in their valence distribution and spin multiplicity, which has important consequences for substrate binding and may explain its low barrier exchange with solvent water. New experimental techniques and innovative sample preparations are beginning to unravel the complex sequence of substrate uptake/inclusion, which is coupled to proton release. The introduction of specific site perturbations, such as replacing Ca²⁺ with Sr²⁺, provides discrete information about the ligand environment of the individual Mn ions. In this way, we have identified a potential open coordination site for one Mn center, which may serve as a substrate binding site in the higher S-states, such as S₃ and S₄. In addition, we can now monitor the binding of the substrate water in the lower S-states (S₁ and S₂) using new EPR-detected NMR spectroscopies. These studies provided the first evidence that one of the substrates is subsumed into the complex itself and forms an oxo-bridge between two Mn ions. This result places important new restrictions on the mechanism of O-O bond formation. These new insights from natures water splitting catalyst provide important criteria for the rational design of bioinspired synthetic catalysts.

Electronic structure of the oxygen-evolving complex in photosystem II prior to O-O bond formation

Setting the stage for release of oxygen Plants transform water into the oxygen we breathe using a protein-bound cluster of four manganese (Mn) ions and a calcium ion. Cox et al. now establish the precise electronic structure in that cluster immediately before formation of the O-O bond (see the Perspective by Britt and Oyala). Using the technique of electron paramagnetic resonance spectroscopy, they confirm a hypothesis that all four Mn ions are octahedrally coordinated and in the 4+ oxidation state. Such clues to the efficiency of the photosynthetic process, so essential to life on Earth, may also facilitate the development of artificial waters-plitting catalysts. Science, this issue p. 804; see also p. 736 Electron paramagnetic resonance spectroscopy reveals that photosynthetic production of O2 proceeds from four MnIV ions. [Also see Perspective by Britt and Oyala] The photosynthetic protein complex photosystem II oxidizes water to molecular oxygen at an embedded tetramanganese-calcium cluster. Resolving the geometric and electronic structure of this cluster in its highest metastable catalytic state (designated S3) is a prerequisite for understanding the mechanism of O-O bond formation. Here, multifrequency, multidimensional magnetic resonance spectroscopy reveals that all four manganese ions of the catalyst are structurally and electronically similar immediately before the final oxygen evolution step; they all exhibit a 4+ formal oxidation state and octahedral local geometry. Only one structural model derived from quantum chemical modeling is consistent with all magnetic resonance data; its formation requires the binding of an additional water molecule. O-O bond formation would then proceed by the coupling of two proximal manganese-bound oxygens in the transition state of the cofactor.

Two Interconvertible Structures that Explain the Spectroscopic Properties of the Oxygen-Evolving Complex of Photosystem II in the S2 State†

Using models derived from the X-ray structure of photosystem II, it is shown that the oxygen evolving complex in the S(2) state exists in two energetically similar and interconvertible forms. A longstanding question regarding the spectroscopy of the catalyst is thus answered: one form corresponds to the multiline g=2.0 EPR signal (see picture, right; O red, Mn purple, Ca yellow), and the other to the g≥4.1 signals (left).

Theoretical evaluation of structural models of the S2 state in the oxygen evolving complex of Photosystem II: protonation states and magnetic interactions.

Protonation states of water ligands and oxo bridges are intimately involved in tuning the electronic structures and oxidation potentials of the oxygen evolving complex (OEC) in Photosystem II, steering the mechanistic pathway, which involves at least five redox state intermediates S(n) (n = 0-4) resulting in the oxidation of water to molecular oxygen. Although protons are practically invisible in protein crystallography, their effects on the electronic structure and magnetic properties of metal active sites can be probed using spectroscopy. With the twin purpose of aiding the interpretation of the complex electron paramagnetic resonance (EPR) spectroscopic data of the OEC and of improving the view of the cluster at the atomic level, a complete set of protonation configurations for the S(2) state of the OEC were investigated, and their distinctive effects on magnetic properties of the cluster were evaluated. The most recent X-ray structure of Photosystem II at 1.9 Å resolution was used and refined to obtain the optimum structure for the Mn(4)O(5)Ca core within the protein pocket. Employing this model, a set of 26 structures was constructed that tested various protonation scenarios of the water ligands and oxo bridges. Our results suggest that one of the two water molecules that are proposed to coordinate the outer Mn ion (Mn(A)) of the cluster is deprotonated in the S(2) state, as this leads to optimal experimental agreement, reproducing the correct ground state spin multiplicity (S = 1/2), spin expectation values, and EXAFS-derived metal-metal distances. Deprotonation of Ca(2+)-bound water molecules is strongly disfavored in the S(2) state, but dissociation of one of the two water ligands appears to be facile. The computed isotropic hyperfine couplings presented here allow distinctions between models to be made and call into question the assumption that the largest coupling is always attributable to Mn(III). The present results impose limits for the total charge and the proton configuration of the OEC in the S(2) state, with implications for the cascade of events in the Kok cycle and for the water splitting mechanism.

What is not required to make a single molecule magnet

The widely accepted assumption that the development of more efficient single molecule magnets must involve ever higher total spin values has so far driven synthetic efforts towards molecular clusters of increasing nuclearity. In the present paper it is suggested that it might be worthwhile to reconsider this approach. There is evidence from theory and experiment to suggest that the race for multinuclear complexes with higher total spin might not necessarily be fruitful as a strategy for maximizing the magnetic relaxation barrier. Instead, we propose that more effort should be directed in understanding the parameters involved in maximizing the anisotropy of small, perhaps even mononuclear, molecules. Using multi-reference ab initio calculations we demonstrate the theory that can be applied and the principles of the computational approach for representative mononuclear complexes. Such small units may subsequently be employed as building blocks for the controlled assembly of larger and maximally anisotropic single molecule magnets.

All-Electron Scalar Relativistic Basis Sets for the Lanthanides.

Segmented all-electron relativistically contracted (SARC) basis sets are constructed for the elements 57La-71Lu and optimized for density functional theory (DFT) applications. The basis sets are intended for use in combination with the DKH2 or ZORA scalar relativistic Hamiltonians for which individually optimized contractions are provided. Significant computational advantages can be realized owing to the loose contraction of the SARC basis sets compared to generally contracted basis sets, while their compact size allows them to replace effective core potentials for routine studies of lanthanide complexes. The new basis sets are evaluated in DFT calculations of the first four ionization energies of the lanthanides. They yield results that accurately reproduce the experimental trends, confirming a balanced treatment of different electronic configurations. The performance of the basis sets is further assessed in molecular systems with a comprehensive study of the lanthanide trihalides. Despite their compact size, the SARC basis sets demonstrate consistent, efficient, and reliable performance and will be especially useful in calculations of molecular properties that require explicit treatment of the core electrons.

Density functional theory

Density functional theory (DFT) finds increasing use in applications related to biological systems. Advancements in methodology and implementations have reached a point where predicted properties of reasonable to high quality can be obtained. Thus, DFT studies can complement experimental investigations, or even venture with some confidence into experimentally unexplored territory. In the present contribution, we provide an overview of the properties that can be calculated with DFT, such as geometries, energies, reaction mechanisms, and spectroscopic properties. A wide range of spectroscopic parameters is nowadays accessible with DFT, including quantities related to infrared and optical spectra, X-ray absorption and Mössbauer, as well as all of the magnetic properties connected with electron paramagnetic resonance spectroscopy except relaxation times. We highlight each of these fields of application with selected examples from the recent literature and comment on the capabilities and limitations of current methods.

Ammonia binding to the oxygen-evolving complex of photosystem II identifies the solvent-exchangeable oxygen bridge (μ-oxo) of the manganese tetramer

The assignment of the two substrate water sites of the tetra-manganese penta-oxygen calcium (Mn4O5Ca) cluster of photosystem II is essential for the elucidation of the mechanism of biological O-O bond formation and the subsequent design of bio-inspired water-splitting catalysts. We recently demonstrated using pulsed EPR spectroscopy that one of the five oxygen bridges (μ-oxo) exchanges unusually rapidly with bulk water and is thus a likely candidate for one of the substrates. Ammonia, a water analog, was previously shown to bind to the Mn4O5Ca cluster, potentially displacing a water/substrate ligand [Britt RD, et al. (1989) J Am Chem Soc 111(10):3522–3532]. Here we show by a combination of EPR and time-resolved membrane inlet mass spectrometry that the binding of ammonia perturbs the exchangeable μ-oxo bridge without drastically altering the binding/exchange kinetics of the two substrates. In combination with broken-symmetry density functional theory, our results show that (i) the exchangable μ-oxo bridge is O5 {using the labeling of the current crystal structure [Umena Y, et al. (2011) Nature 473(7345):55–60]}; (ii) ammonia displaces a water ligand to the outer manganese (MnA4-W1); and (iii) as W1 is trans to O5, ammonia binding elongates the MnA4-O5 bond, leading to the perturbation of the μ-oxo bridge resonance and to a small change in the water exchange rates. These experimental results support O-O bond formation between O5 and possibly an oxyl radical as proposed by Siegbahn and exclude W1 as the second substrate water.

Detailed Ab Initio First-Principles Study of the Magnetic Anisotropy in a Family of Trigonal Pyramidal Iron(II) Pyrrolide Complexes

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 d(6) 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(tpa(R))](-) (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.