Search

Strain effects on g and metal hyperfine coupling tensors, A , are often manifested in Electron Paramagnetic Resonance (EPR) spectra of transition metal complexes, as a result of their intrinsic and/or solvent-mediated structural variations. Although distributions of these tensors are quite common and well understood in continuous-wave (cw) EPR spectroscopy, reported strain effects on ligand hyperfine coupling constants are rather scarce. Here we explore the case of a conformationally flexible Cu(II) complex, [Cu{Ph 2 P(O)NP(O)Ph 2 -κ 2 O,O '} 2 ], bearing P atoms in its second coordination sphere and exhibiting two structurally distinct CuO 4 coordination spheres, namely a square planar and a tetrahedrally distorted one, as revealed by X-ray crystallography. The Hyperfine Sublevel Correlation (HYSCORE) spectra of this complex exhibit 31 P correlation ridges that have unusual inverse or so-called "boomerang" shapes and features that cannot be reproduced by standard simulation procedures assuming only one set of magnetic parameters. Our work shows that a distribution of isotropic hyperfine coupling constants (hfc) spanning a range between negative and positive values is necessary in order to describe in detail the unusual shapes of HYSCORE spectra. By employing DFT calculations we show that these hfc correspond to molecules showing variable distortions from square planar to tetrahedral geometry, and we demonstrate that line shape analysis of such HYSCORE spectra provides new insight into the conformation-dependent spectroscopic response of the spin system under investigation.

Calcium is an abundant, nontoxic metal that finds many roles in synthetic and biological systems including the oxygen-evolving complex (OEC) of photosystem II. Characterization methods for calcium centers, however, are underdeveloped compared to those available for transition metals. Valence-to-core X-ray emission spectroscopy (VtC XES) selectively probes the electronic structure of an element's chemical environment, providing insight that complements the geometric information available from other techniques. Here, the utility of calcium VtC XES is established using an in-house dispersive spectrometer in combination with density functional theory. Spectral trends are rationalized within a molecular orbital framework, and Kβ 2,5 transitions, derived from molecular orbitals with primarily ligand p character, are found to be a promising probe of the calcium coordination environment. In particular, it is shown that calcium VtC XES is sensitive to the electronic structure changes that accompany oxo protonation in Mn 3 CaO 4 -based molecular mimics of the OEC. Through correlation to calculations, the potential of calcium VtC XES to address unresolved questions regarding the mechanism of biological water oxidation is highlighted.

A new tool for the interpretation of multiconfigurational wave functions representing the spin states of exchange-coupled transition metal complexes is introduced. Based on orbital entanglement measures, herein derived from multiconfigurational density matrix renormalization group calculations, the complexity of the wave function is reduced, thus facilitating a connection with established concepts for the interpretation of magnetically coupled systems. We show that the entanglement of localized orbitals with a small basis set is a good representation of the magnetic coupling topology and that it is sensitive to chemical changes in homologous complexes. Furthermore, we introduce a measure for the magnetic relevance of orbitals in the active subspace and a concept for the quantitative comparison of different chemical species. The approach presented here will be easily applicable to higher nuclearity clusters, providing a direct insight into all states of the Heisenberg spin ladder for systems previously accessible only by single-configurational methods.

Heterogeneity in intermediate catalytic states of the oxygen-evolving complex (OEC) of Photosystem II is known from a wide range of experimental and theoretical data, but its potential implications for the mechanism of water oxidation remain unexplored. We delineate the consequences of structural heterogeneity for the final step of the catalytic cycle by tracing the evolution of three spectroscopically relevant and structurally distinct components of the last metastable S 3 state to the transient O 2 -evolving S 4 state of the OEC. Using quantum chemical calculations, we show that each S 3 isomer leads to a different electronic structure formulation for the active S 4 state. Crucially, in addition to previously hypothesized Mn(IV)-oxyl species, we establish for the first time, how a genuine Mn(V)-oxo can be obtained in the catalytically active S 4 state: this takes the form of a five-coordinate and locally high-spin (S Mn  = 1) Mn(V) site. This formulation for the S 4 state evolves naturally from a preceding S 3 -state structural intermediate that contains a quasi-trigonal-bipyramidal Mn(IV) ion. The results strongly suggest that water binding in the S 3 state is not prerequisite for reaching the oxygen-evolving S 4 state of the complex, supporting the notion that both substrates are preloaded at the beginning of the catalytic cycle. This scenario allows true four-electron metal-centered hole accumulation to precede OO bond formation and hence the latter can proceed via a genuine even-electron mechanism. This can occur as intramolecular nucleophilic coupling of two oxo units synchronously with the binding of a water substrate for the next catalytic cycle., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The variation in the singlet-triplet energy gap of diphenylcarbene (DPC) upon interaction with hydrogen (water and methanol) or halogen bond (XCF 3 , X = Cl, Br, I) donors to form van der Waals (vdW) complexes is investigated in relation to the electrostatic and dispersion components of such intermolecular interactions. The domain-based local pair natural orbital coupled cluster method, DLPNO-CCSD(T), is used for calculating accurate single-triplet energy gaps and interaction energies for both spin states. The local energy decomposition scheme is used to provide an accurate quantification to the various interaction energy components at the DLPNO-CCSD(T) level. It is shown that the formation of vdW adducts stabilizes the singlet state of DPC, and in the case of water, methanol, and ICF 3 , it reverses the ground state from triplet to singlet. Electrostatic interactions are significant in both spin states, but preferentially stabilize the singlet state. For methanol and ClCF 3 , London dispersion forces have the opposite effect, stabilizing preferentially the triplet state. The quantification of the energetic components of the interactions through the local energy decomposition analysis correlates well with experimental findings and provides the basis for more elaborate treatments of microsolvation in carbenes.

In biological water oxidation, a redox-active tyrosine residue (D1-Tyr161 or Y Z ) mediates electron transfer between the Mn 4 CaO 5 cluster of the oxygen-evolving complex and the charge-separation site of photosystem II (PSII), driving the cluster through progressively higher oxidation states S i ( i = 0-4). In contrast to lower S-states (S 0 , S 1 ), in higher S-states (S 2 , S 3 ) of the Mn 4 CaO 5 cluster, Y Z cannot be oxidized at cryogenic temperatures due to the accumulation of positive charge in the S 1 → S 2 transition. However, oxidation of Y Z by illumination of S 2 at 77-190 K followed by rapid freezing and charge recombination between Y Z • and the plastoquinone radical Q A •- allows trapping of an S 2 variant, the so-called S 2 trapped state (S 2 t ), that is capable of forming Y Z • at cryogenic temperature. To identify the differences between the S 2 and S 2 t states, we used the S 2 t Y Z • intermediate as a probe for the S 2 t state and followed the S 2 t Y Z • /Q A •- recombination kinetics at 10 K using time-resolved electron paramagnetic resonance spectroscopy in H 2 O and D 2 O. The results show that while S 2 t Y Z • /Q A •- recombination can be described as pure electron transfer occurring in the Marcus inverted region, the S 2 t → S 2 reversion depends on proton rearrangement and exhibits a strong kinetic isotope effect. This suggests that Y Z oxidation in the S 2 t state is facilitated by favorable proton redistribution in the vicinity of Y Z , most likely within the hydrogen-bonded Y Z -His190-Asn298 triad. Computational models show that tautomerization of Asn298 to its imidic acid form enables proton translocation to an adjacent asparagine-rich cavity of water molecules that functions as a proton reservoir and can further participate in proton egress to the lumen.

Photosystem II (PSII) of oxygenic photosynthesis captures sunlight to drive the catalytic oxidation of water and the reduction of plastoquinone. Among the several redox-active cofactors that participate in intricate electron transfer pathways there are two tyrosine residues, Y Z and Y D . They are situated in symmetry-related electron transfer branches but have different environments and play distinct roles. Y Z is the immediate oxidant of the oxygen-evolving Mn 4 CaO 5 cluster, whereas Y D serves regulatory and protective functions. The protonation states and hydrogen-bond network in the environment of Y D remain debated, while the role of microsolvation in stabilizing different redox states of Y D and facilitating oxidation or mediating deprotonation, as well the fate of the phenolic proton, is unclear. Here we present detailed structural models of Y D and its environment using large-scale quantum mechanical models and all-atom molecular dynamics of a complete PSII monomer. The energetics of water distribution within a hydrophobic cavity adjacent to Y D are shown to correlate directly with electron paramagnetic resonance (EPR) parameters such as the tyrosyl g-tensor, allowing us to map the correspondence between specific structural models and available experimental observations. EPR spectra obtained under different conditions are explained with respect to the mode of interaction of the proximal water with the tyrosyl radical and the position of the phenolic proton within the cavity. Our results revise previous models of the energetics and build a detailed view of the role of confined water in the oxidation and deprotonation of Y D . Finally, the model of microsolvation developed in the present work rationalizes in a straightforward way the biphasic oxidation kinetics of Y D , offering new structural insights regarding the function of the radical in biological photosynthesis.

Face-sharing octahedral dinuclear Cr(III) compounds with d 3 -d 3 electronic configurations represent nontrivial examples of electronic complexity, posing particular challenges for theoretical and computational studies. A tris-hydroxy-bridged Cr(III)-Cr(III) system has proven to be a richly rewarding target for studies of magnetism and electron paramagnetic resonance spectroscopy. It was also reported to be a peculiarly difficult system to treat with density functional theory (DFT). In this work the magnetic coupling problem for this dimer is approached with broken-symmetry (BS)-DFT and multireference calculations that utilize the density matrix renormalization group (DMRG) to handle full-valence active spaces. BS-DFT is shown to recover the correct ordering and energy spacing of Heisenberg spin states if used in conjunction with appropriate spin projection procedures, albeit with pronounced functional sensitivity. The contrasting conclusions of previous studies are traced to incorrect inclusion of electronically excited configurations. Analysis of the direct and differential overlap of corresponding orbital pairs from the BS-DFT solution indicates that metal-metal through-space interaction is the dominant contributor to antiferromagnetic coupling. At the DFT level a procedure that utilizes pseudopotential substitution is demonstrated that allows evaluation of the direct exchange vs superexchange contributions. A complementary description is obtained with DMRG-SCF calculations that enable state-averaged CASSCF calculations with both metal and bridge orbitals in the active space. A localized orbital subspace analysis supports the DFT conclusions that in contrast to doubly bridged isoelectronic analogues, antiferromagnetic coupling in the chromium dimer arises primarily from direct metal-metal interaction but is significantly enhanced by ligand-mediated superexchange.

Quinones play vital roles as electron carriers in fundamental biological processes; therefore, the ability to accurately predict their electron affinities is crucial for understanding their properties and function. The increasing availability of cost-effective implementations of correlated wave function methods for both closed-shell and open-shell systems offers an alternative to density functional theory approaches that have traditionally dominated the field despite their shortcomings. Here, we define a benchmark set of quinones with experimentally available electron affinities and evaluate a range of electronic structure methods, setting a target accuracy of 0.1 eV. Among wave function methods, we test various implementations of coupled cluster (CC) theory, including local pair natural orbital (LPNO) approaches to canonical and parameterized CCSD, the domain-based DLPNO approximation, and the equations-of-motion approach for electron affinities, EA-EOM-CCSD. In addition, several variants of canonical, spin-component-scaled, orbital-optimized, and explicitly correlated (F12) Møller-Plesset perturbation theory are benchmarked. Achieving systematically the target level of accuracy is challenging and a composite scheme that combines canonical CCSD(T) with large basis set LPNO-based extrapolation of correlation energy proves to be the most accurate approach. Methods that offer comparable performance are the parameterized LPNO-pCCSD, the DLPNO-CCSD(T 0 ), and the orbital optimized OO-SCS-MP2. Among DFT methods, viable practical alternatives are only the M06 and the double hybrids, but the latter should be employed with caution because of significant basis set sensitivity. A highly accurate yet cost-effective DLPNO-based coupled cluster approach is used to investigate the methoxy conformation effect on the electron affinities of ubiquinones found in photosynthetic bacterial reaction centers. © 2018 Wiley Periodicals, Inc., (© 2018 Wiley Periodicals, Inc.)

A test set of 12 aryl carbenes (AC12) is compiled with the purpose of establishing their adiabatic singlet-triplet energy splittings using correlated wave function based methods. The set covers both singlet and triplet ground state aryl carbenes, as well as a range of magnitudes for the ground state to excited state gap. The performance of coupled cluster methods is examined with respect to the reference wave function, the basis set, and a number of additional methodological parameters that enter the calculation. Inclusion of perturbative triples and basis set extrapolation with a combination of triple and quadruple-ζ basis sets are both required to ensure high accuracy. When canonical coupled cluster calculations become too expensive, the domain-based local pair natural orbital approach DLPNO-CCSD(T) can be used as a reliable method for larger systems, as it achieves a mean absolute error of only 0.2 kcal/mol for the singlet-triplet gaps in the present test set. Other first-principles wave function methods and selected density functional methods are also evaluated. Second-order Møller-Plesset perturbation theory approaches are only applicable in conjunction with orbital optimization (OO-MP2). Among the representative density functional methods tested, only double hybrid functionals perform sufficiently accurately to be considered useful for systems with small singlet-triplet gaps. On the basis of the reference coupled cluster results, projected gas-phase free energies are reported for all aryl carbenes. Finally, the treatment of singlet-triplet gaps in solution is discussed in terms of both implicit and explicit solvation.

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.

Nature's water-splitting catalyst, an oxygen-bridged tetramanganese calcium (Mn 4 O 5 Ca) complex, sequentially activates two substrate water molecules generating molecular O 2 . Its reaction cycle is composed of five intermediate (S i ) states, where the index i indicates the number of oxidizing equivalents stored by the cofactor. After formation of the S 4 state, the product dioxygen is released and the cofactor returns to its lowest oxidation state, S 0 . Membrane-inlet mass spectrometry measurements suggest that at least one substrate is bound throughout the catalytic cycle, as the rate of 18 O-labeled water incorporation into the product O 2 is slow, on a millisecond to second time scale depending on the S state. Here, we demonstrate that the Mn 4 O 5 Ca complex poised in the S 0 state contains an exchangeable hydroxo bridge. On the basis of a combination of magnetic multiresonance (EPR) spectroscopies, comparison to biochemical models and theoretical calculations we assign this bridge to O5, the same bridge identified in the S 2 state as an exchangeable fully deprotonated oxo bridge [Pérez Navarro, M.; et al. Proc. Natl. Acad. Sci. U.S.A. 2013, 110, 15561]. This oxygen species is the most probable candidate for the slowly exchanging substrate water in the S 0 state. Additional measurements provide new information on the Mn ions that constitute the catalyst. A structural model for the S 0 state is proposed that is consistent with available experimental data and explains the observed evolution of water exchange kinetics in the first three states of the catalytic cycle.

The site of biological water oxidation is highly conserved across photosynthetic organisms, but differences of unidentified structural and electronic origin exist between taxonomically discrete clades, revealed by distinct spectroscopic signatures of the oxygen-evolving Mn 4 CaO 5 cluster and variations in active-site accessibility. Comparison of atomistic models of a native cyanobacterial form (Thermosynechococcus vulcanus) and a chimeric spinach-like form of photosystem II allows us to identify the precise atomic-level differences between organisms in the vicinity of the manganese cluster. Substitution of cyanobacterial D1-Asn87 by higher-plant D1-Ala87 is the principal discriminating feature: it drastically rearranges a network of proximal hydrogen bonds, modifying the local architecture of a water channel and the interaction of second coordination shell residues with the manganese cluster. The two variants explain species-dependent differences in spectroscopic properties and in the interaction of substrate analogues with the oxygen-evolving complex, enabling assignment of a substrate delivery channel to the active site.

A recently reported synthetic complex with a Mn 4 CaO 4 core represents a remarkable structural mimic of the Mn 4 CaO 5 cluster in the oxygen-evolving complex (OEC) of photosystem II (Zhang et al., Science 2015, 348, 690). Oxidized samples of the complex show electron paramagnetic resonance (EPR) signals at g ≈ 4.9 and 2, similar to those associated with the OEC in its S 2 state (g ≈ 4.1 from an S = 5 / 2 form and g ≈ 2 from an S = 1 / 2 form), suggesting similarities in the electronic as well as geometric structure. We use quantum-chemical methods to characterize the synthetic complex in various oxidation states, to compute its magnetic and spectroscopic properties, and to establish connections with reported data. Only one energetically accessible form is found for the oxidized "S 2 state" of the complex. It has a ground spin state of S = 5 / 2 , and EPR simulations confirm it can be assigned to the g ≈ 4.9 signal. However, no valence isomer with an S = 1 / 2 ground state is energetically accessible, a conclusion supported by a wide range of methods, including density matrix renormalization group with full valence active space. Alternative candidates for the g ≈ 2 signal were explored, but no low-spin/low-energy structure was identified. Therefore, our results suggest that despite geometric similarities the synthetic model does not mimic the valence isomerism that is the hallmark of the OEC in its S 2 state, most probably because it lacks a coordinatively flexible oxo bridge. Only one of the observed EPR signals can be explained by a structurally intact high-spin one-electron-oxidized form, while the other originates from an as-yet-unidentified rearrangement product. Nevertheless, this model provides valuable information for understanding the high-spin EPR signals of both the S 1 and S 2 states of the OEC in terms of the coordination number and Jahn-Teller axis orientation of the Mn ions, with important consequences for the development of magnetic spectroscopic probes to study S-state intermediates immediately prior to O-O bond formation.

The pentanuclear Mn complex [Mn 5 (μ 3 -O)(bpp) 6 ] n+ can access six consecutive total oxidation levels from Mn(ii) 5 to Mn(iii) 5 . The electronic structure and redox behavior of this cluster are studied computationally and the results are compared with experimental data. The tunability of the redox potential span and of the absolute potential position is explored through systematic modifications of the 3,5-bis(pyridin-2-yl)-pyrazole (Hbpp) ligand. By substitutions with electron-donating and electron-withdrawing groups, the position of the redox events can be shifted by ca. 1 eV, however no ligand modification is predicted to alter the total redox span. By changing the pyridyl groups to benzimidazole groups, yielding the 3,5-bis(benzimidazol-2-yl)pyrazole (H 3 bbp) ligand, it is predicted that higher oxidation states can be stabilized, from Mn(iii) 4 Mn(iv) up to Mn(iv) 5 . In this system, the redox span for the same number of redox equivalents accumulated is slightly smaller than that in the original system. The manganese system is compared with its pentanuclear iron analogue that was recently reported to be catalytically active in oxygen evolution (Okamura et al., Nature, 2016, 530, 465). The electronic and structural requirements for utilization of the stored oxidizing equivalents in water oxidation are discussed.

Methanol has long being used as a substrate analogue to probe access pathways and investigate water delivery at the oxygen-evolving complex (OEC) of photosystem-II. In this contribution we study the interaction of methanol with the OEC by assembling available spectroscopic data into a quantum mechanical treatment that takes into account the local channel architecture of the active site. The effect on the magnetic energy levels of the Mn 4 Ca cluster in the S 2 state of the catalytic cycle can be explained equally well by two models that involve either methanol binding to the calcium ion of the cluster, or a second-sphere interaction in the vicinity of the "dangler" Mn4 ion. However, consideration of the latest 13 C hyperfine interaction data shows that only one model is fully consistent with experiment. In contrast to previous hypotheses, methanol is not a direct ligand to the OEC, but is situated at the end-point of a water channel associated with the O4 bridge. Its effect on magnetic properties of plant PS-II results from disruption of hydrogen bonding between O4 and proximal channel water molecules, thus enhancing superexchange (antiferromagnetic coupling) between the Mn3 and Mn4 ions. The same interaction mode applies to the dark-stable S 1 state and possibly to all other states of the complex. Comparison of protein sequences from cyanobacteria and plants reveals a channel-altering substitution (D1-Asn87 versus D1-Ala87) in the proximity of the methanol binding pocket, explaining the species-dependence of the methanol effect. The water channel established as the methanol access pathway is the same that delivers ammonia to the Mn4 ion, supporting the notion that this is the only directly solvent-accessible manganese site of the OEC. The results support the pivot mechanism for water binding at a component of the S 3 state and would be consistent with partial inhibition of water delivery by methanol. Mechanistic implications for enzymatic regulation and catalytic progression are discussed.

The calculation of redox potentials involves large energetic terms arising from gas phase ionization energies, thermodynamic contributions, and solvation energies of the reduced and oxidized species. In this work we study the performance of a wide range of wave function and density functional theory methods for the prediction of ionization energies and aqueous one-electron oxidation potentials of a set of 19 organic molecules. Emphasis is placed on evaluating methods that employ the computationally efficient local pair natural orbital (LPNO) approach, as well as several implementations of coupled cluster theory and explicitly correlated F12 methods. The electronic energies are combined with implicit solvation models for the solvation energies. With the exception of MP2 and its variants, which suffer from enormous errors arising at least partially from the poor Hartree-Fock reference, ionization energies can be systematically predicted with average errors below 0.1 eV for most of the correlated wave function based methods studies here, provided basis set extrapolation is performed. LPNO methods are the most efficient way to achieve this type of accuracy. DFT methods show in general larger errors and suffer from inconsistent behavior. The only exception is the M06-2X functional which is found to be competitive with the best LPNO-based approaches for ionization energies. Importantly, the limiting factor for the calculation of accurate redox potentials is the solvation energy. The errors in the predicted solvation energies by all continuum solvation models tested in this work dominate the final computed reduction potential, resulting in average errors typically in excess of 0.3 V and hence obscuring the gains that arise from choosing a more accurate electronic structure method.

The redox potential of synthetic oligonuclear transition metal complexes has been shown to correlate with the Lewis acidity of a redox-inactive cation connected to the redox-active transition metals of the cluster via oxo or hydroxo bridges. Such heterometallic clusters are important cofactors in many metalloenzymes, where it is speculated that the redox-inactive constituent ion of the cluster serves to optimize its redox potential for electron transfer or catalysis. A principal example is the oxygen-evolving complex in photosystem II of natural photosynthesis, a Mn4CaO5 cofactor that oxidizes water into dioxygen, protons and electrons. Calcium is critical for catalytic function, but its precise role is not yet established. In analogy to synthetic complexes it has been suggested that Ca(2+) fine-tunes the redox potential of the manganese cluster. Here we evaluate this hypothesis by computing the relative redox potentials of substituted derivatives of the oxygen-evolving complex with the cations Sr(2+), Gd(3+), Cd(2+), Zn(2+), Mg(2+), Sc(3+), Na(+) and Y(3+) for two sequential transitions of its catalytic cycle. The theoretical approach is validated with a series of experimentally well-characterized Mn3AO4 cubane complexes that are structural mimics of the enzymatic cluster. Our results reproduce perfectly the experimentally observed correlation between the redox potential and the Lewis acidities of redox-inactive cations for the synthetic complexes. However, it is conclusively demonstrated that this correlation does not hold for the oxygen evolving complex. In the enzyme the redox potential of the cluster only responds to the charge of the redox-inactive cations and remains otherwise insensitive to their precise identity, precluding redox-tuning of the metal cluster as a primary role for Ca(2+) in biological water oxidation.

Rapid progress has been made in the last five years towards resolution of the structure of nature's water splitting catalyst - a Mn4O5Ca cofactor embedded in Photosystem II - especially in the field of X-ray crystallography. In addition, recent magnetic resonance data have allowed the structure of the cofactor to be accessed in its last metastable intermediate state, prior to O-O bond formation. This activated form of the catalyst is geometrically similar to that seen by X-ray crystallography, which represents the resting state of the cofactor, but requires the coordination of an additional water molecule to the cofactor, rendering all Mn ions six coordinate. Importantly, it locates two water derived, Mn bound oxygen ligands in close proximity. It is these two oxygen ligands that likely form the product O2 molecule, as proposed earlier by quantum chemical modeling. Current views on the molecular level events that facilitate catalyst activation, that is, catalyst/substrate deprotonation, Mn oxidation and water molecule insertion are briefly described., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Improved versions of the segmented all-electron relativistically contracted (SARC) basis sets for the lanthanides are presented. The second-generation SARC2 basis sets maintain efficient construction of their predecessors and their individual adaptation to the DKH2 and ZORA Hamiltonians, but feature exponents optimized with a completely new orbital shape fitting procedure and a slightly expanded f space that results in sizable improvement in CASSCF energies and in significantly more accurate prediction of spin-orbit coupling parameters. Additionally, an extended set of polarization/correlation functions is constructed that is appropriate for multireference correlated calculations and new auxiliary basis sets for use in resolution-of-identity (density-fitting) approximations in combination with both DFT and wave function based treatments. Thus, the SARC2 basis sets extend the applicability of the first-generation DFT-oriented basis sets to routine all-electron wave function-based treatments of lanthanide complexes. The new basis sets are benchmarked with respect to excitation energies, radial distribution functions, optimized geometries, orbital eigenvalues, ionization potentials, and spin-orbit coupling parameters of lanthanide systems and are shown to be suitable for the description of magnetic and spectroscopic properties using both DFT and multireference wave function-based methods.

In transition-metal complexes, the geometric structure is intimately connected with the spin state arising from magnetic coupling between the paramagnetic ions. The tetramanganese-calcium cofactor that catalyzes biological water oxidation in photosystem II cycles through five catalytic intermediates, each of which adopts a specific geometric and electronic structure and is thus characterized by a specific spin state. Here, we review spin-structure correlations in Nature's water-splitting catalyst. The catalytic cycle of the Mn4O5Ca cofactor can be described in terms of spin-dependent reactivity. The lower "inactive" S states of the catalyst, S0 and S1, are characterized by low-spin ground states, SGS = 1/2 and SGS = 0. This is connected to the "open cubane" topology of the inorganic core in these states. The S2 state exhibits structural and spin heterogeneity in the form of two interconvertible isomers and is identified as the spin-switching point of the catalytic cycle. The first S2 state form is an open cubane structure with a low-spin SGS = 1/2 ground state, whereas the other represents the first appearance of a closed cubane topology in the catalytic cycle that is associated with a higher-spin ground state of SGS = 5/2. It is only this higher-spin form of the S2 state that progresses to the "activated" S3 state of the catalyst. The structure of this final metastable catalytic state was resolved in a recent report, showing that all manganese ions are six-coordinate. The magnetic coupling is dominantly ferromagnetic, leading to a high-spin ground state of SGS = 3. The ability of the Mn4O5Ca cofactor to adopt two distinct structural and spin-state forms in the S2 state is critical for water binding in the S3 state, allowing spin-state crossing from the inactive, low-spin configuration of the catalyst to the activated, high-spin configuration. Here we describe how an understanding of the magnetic properties of the catalyst in all S states has allowed conclusions on the catalyst function to be reached. A summary of recent literature results is provided that constrains the sequence of molecular level events: catalyst/substrate deprotonation, manganese oxidation, and water molecule insertion.

Nature relies on a unique and intricate biochemical setup to achieve sunlight-driven water splitting. Combined experimental and computational efforts have produced significant insights into the structural and functional principles governing the operation of the water-oxidizing enzyme Photosystem II in general, and of the oxygen-evolving manganese-calcium cluster at its active site in particular. Here we review the most important aspects of biological water oxidation, emphasizing current knowledge on the organization of the enzyme, the geometric and electronic structure of the catalyst, and the role of calcium and chloride cofactors. The combination of recent experimental work on the identification of possible substrate sites with computational modeling have considerably limited the possible mechanistic pathways for the critical O-O bond formation step. Taken together, the key features and principles of natural photosynthesis may serve as inspiration for the design, development, and implementation of artificial systems.

Among the four photo-driven transitions of the water-oxidizing tetramanganese-calcium cofactor of biological photosynthesis, the second-last step of the catalytic cycle, that is the S 2 to S 3 state transition, is the crucial step that poises the catalyst for the final O-O bond formation. This transition, whose intermediates are not yet fully understood, is a multi-step process that involves the redox-active tyrosine residue and includes oxidation and deprotonation of the catalytic cluster, as well as the binding of a water molecule. Spectroscopic data has the potential to shed light on the sequence of events that comprise this catalytic step, which still lacks a structural interpretation. In this work the S 2 -S 3 state transition is studied and a key intermediate species is characterized: it contains a Mn 3 O 4 Ca cubane subunit linked to a five-coordinate Mn(iv) ion that adopts an approximately trigonal bipyramidal ligand field. It is shown using high-level density functional and multireference wave function calculations that this species accounts for the near-infrared absorption and electron paramagnetic resonance observations on metastable S 2 -S 3 intermediates. The results confirm that deprotonation and Mn oxidation of the cofactor must precede the coordination of a water molecule, and lead to identification of a novel low-energy water binding mode that has important implications for the identity of the substrates in the mechanism of biological water oxidation.

Multifrequency pulsed EPR data are reported for a series of oxygen bridged (μ-oxo/μ-hydroxo) bimetallic manganese complexes where the oxygen is labeled with the magnetically active isotope (17)O (I = 5/2). Two synthetic complexes and two biological metallocofactors are examined: a planar bis-μ-oxo bridged complex and a bent, bis-μ-oxo-μ-carboxylato bridge complex; the dimanganese catalase, which catalyzes the dismutation of H2O2 to H2O and O2, and the recently identified manganese/iron cofactor of the R2lox protein, a homologue of the small subunit of the ribonuclotide reductase enzyme (class 1c). High field (W-band) hyperfine EPR spectroscopies are demonstrated to be ideal methods to characterize the (17)O magnetic interactions, allowing a magnetic fingerprint for the bridging oxygen ligand to be developed. It is shown that the μ-oxo bridge motif displays a small positive isotropic hyperfine coupling constant of about +5 to +7 MHz and an anisotropic/dipolar coupling of -9 MHz. In addition, protonation of the bridge is correlated with an increase of the hyperfine coupling constant. Broken symmetry density functional theory is evaluated as a predictive tool for estimating hyperfine coupling of bridging species. Experimental and theoretical results provide a framework for the characterization of the oxygen bridge in Mn metallocofactor systems, including the water oxidizing cofactor of photosystem II, allowing the substrate/solvent interface to be examined throughout its catalytic cycle.

First principle calculations of extended X-ray absorption fine structure (EXAFS) data have seen widespread use in bioinorganic chemistry, perhaps most notably for modeling the Mn4Ca site in the oxygen evolving complex (OEC) of photosystem II (PSII). The logic implied by the calculations rests on the assumption that it is possible to a priori predict an accurate EXAFS spectrum provided that the underlying geometric structure is correct. The present study investigates the extent to which this is possible using state of the art EXAFS theory. The FEFF program is used to evaluate the ability of a multiple scattering-based approach to directly calculate the EXAFS spectrum of crystallographically defined model complexes. The results of these parameter free predictions are compared with the more traditional approach of fitting FEFF calculated spectra to experimental data. A series of seven crystallographically characterized Mn monomers and dimers is used as a test set. The largest deviations between the FEFF calculated EXAFS spectra and the experimental EXAFS spectra arise from the amplitudes. The amplitude errors result from a combination of errors in calculated S0(2) and Debye-Waller values as well as uncertainties in background subtraction. Additional errors may be attributed to structural parameters, particularly in cases where reliable high-resolution crystal structures are not available. Based on these investigations, the strengths and weaknesses of using first-principle EXAFS calculations as a predictive tool are discussed. We demonstrate that a range of DFT optimized structures of the OEC may all be considered consistent with experimental EXAFS data and that caution must be exercised when using EXAFS data to obtain topological arrangements of complex clusters.

Carboxylate-bridged Mn(II)-Ca(II) complexes are potentially relevant for mimicking the first stages of the Oxygen-Evolving Complex (OEC) assembly process. Here, we report on new homonuclear Ca(II) and heteronuclear Mn(II)–Ca(II) complexes with carboxylate-functionalized tripodal tris(2-pyridylmethyl)amine ligands, the heptadentate H3tpaa, previously reported, and the new hexadentate H2tpada, containing respectively three and two carboxylate units. The mononuclear [Ca(Htpaa)(OH2)] (Ca1) and dinuclear [Ca(tpada)(OH2)2]2 (Ca2) calcium complexes, as well as the tetranuclear [{Mn(tpaa)}2{Ca(OH2)5(μ-OH2)}2][Mn(tpaa)]2 (Mn2Ca2·2Mn) and dinuclear [Mn(tpada)ClCa(OH2)2.67(MeOH)2.33]Cl (MnCa) heterometallic species have been structurally characterized; the syntheses of Ca1 and Mn2Ca2·2Mn being previously reported by us (Inorg. Chem., 2015, 54, 1283). The Mn(II) and Ca(II) are linked by two μ1,1-bridging carboxylates in MnCa, while only one μ1,3-carboxylate bridge connects each Ca2+ ion to each Mn(II) in Mn2Ca2. A variable number of water molecules (n = 1 to 7) are coordinated to Ca in all complexes, most of them being involved in hydrogen-bond networks, in analogy to what occurs in the photosystem II. All donor atoms of the tpaa3- and tpada2- ligands are coordinated to the Mn2+ ions, despite the unusually long distance between the Mn2+ ion and the tertiary amine imposed by the constraining nature of the ligands, as supported by theoretical calculations. Solid state EPR spectroscopy, in combination with DFT calculations, has also shown that the Ca2+ ion has an effect on the electronic parameters (zero field splitting) of the linked Mn(II) in the case of MnCa (μ1,1-carboxylate bridges). In Mn2Ca2 (μ1,3-carboxylate bridge) the Ca2+ ion induces only slight structural changes in the Mn coordination sphere.

Herein, we describe an uncommon example of a manganese-thiolate complex, which is capable of activating dioxygen and catalyzing its two-electron reduction to generate H2O2. The structurally characterized dimercapto-bridged Mn(II) dimer [Mn(II)2(LS)(LSH)]ClO4 (Mn(II)2SH) is formed by reaction of the LS ligand (2,2'-(2,2'-bipyridine-6,6'-diyl)bis(1,1-diphenylethanethiolate)) with Mn(II). The unusual presence of a pendant thiol group bound to one Mn(II) ion in Mn(II)2SH is evidenced both in the solid state and in solution. The Mn(II)2SH complex reacts with dioxygen in CH3CN, leading to the formation of a rare mono-μ-hydroxo dinuclear Mn(III) complex, [(Mn(III)2(LS)2(OH)]ClO4 (Mn(III)2OH), which has also been structurally characterized. When Mn(II)2SH reacts with O2 in the presence of a proton source, 2,6-lutidinium tetrafluoroborate (up to 50 equiv), the formation of a new Mn species is observed, assigned to a bis-μ-thiolato dinuclear Mn(III) complex with two terminal thiolate groups (Mn(III)2), with the concomitant production of H2O2 up to ∼40% vs Mn(II)2SH. The addition of a catalytic amount of Mn(II)2SH to an air-saturated solution of MenFc (n = 8 or 10) and 2,6-lutidinium tetrafluoroborate results in the quantitative and efficient oxidation of MenFc by O2 to afford the respective ferrocenium derivatives (MenFc(+), with n = 8 or 10). Hydrogen peroxide is mainly produced during the catalytic reduction of dioxygen with 80-84% selectivity, making the Mn(II)2SH complex a rare Mn-based active catalyst for two-electron O2 reduction.

In the context of a global artificial photosynthesis (GAP) project, we review our current work on nature's water splitting catalyst. In a recent report (Cox et al. 2014 Science 345, 804-808 (doi:10.1126/science.1254910)), we showed that the catalyst-a Mn4O5Ca cofactor-converts into an 'activated' form immediately prior to the O-O bond formation step. This activated state, which represents an all Mn(IV) complex, is similar to the structure observed by X-ray crystallography but requires the coordination of an additional water molecule. Such a structure locates two oxygens, both derived from water, in close proximity, which probably come together to form the product O2 molecule. We speculate that formation of the activated catalyst state requires inherent structural flexibility. These features represent new design criteria for the development of biomimetic and bioinspired model systems for water splitting catalysts using first-row transition metals with the aim of delivering globally deployable artificial photosynthesis technologies.

A central question in biological water splitting concerns the oxidation states of the manganese ions that comprise the oxygen-evolving complex of photosystem II. Understanding the nature and order of oxidation events that occur during the catalytic cycle of five S i states ( i = 0-4) is of fundamental importance both for the natural system and for artificial water oxidation catalysts. Despite the widespread adoption of the so-called "high-valent scheme"-where, for example, the Mn oxidation states in the S 2 state are assigned as III, IV, IV, IV-the competing "low-valent scheme" that differs by a total of two metal unpaired electrons ( i.e. III, III, III, IV in the S 2 state) is favored by several recent studies for the biological catalyst. The question of the correct oxidation state assignment is addressed here by a detailed computational comparison of the two schemes using a common structural platform and theoretical approach. Models based on crystallographic constraints were constructed for all conceivable oxidation state assignments in the four (semi)stable S states of the oxygen evolving complex, sampling various protonation levels and patterns to ensure comprehensive coverage. The models are evaluated with respect to their geometric, energetic, electronic, and spectroscopic properties against available experimental EXAFS, XFEL-XRD, EPR, ENDOR and Mn K pre-edge XANES data. New 2.5 K 55 Mn ENDOR data of the S 2 state are also reported. Our results conclusively show that the entire S state phenomenology can only be accommodated within the high-valent scheme by adopting a single motif and protonation pattern that progresses smoothly from S 0 (III, III, III, IV) to S 3 (IV, IV, IV, IV), satisfying all experimental constraints and reproducing all observables. By contrast, it was impossible to construct a consistent cycle based on the low-valent scheme for all S states. Instead, the low-valent models developed here may provide new insight into the over-reduced S states and the states involved in the assembly of the catalytically active water oxidizing cluster.

The interpretation of electron paramagnetic resonance spectra of polynuclear transition metal complexes in terms of individual contributions from each paramagnetic center can be greatly facilitated by the availability of theoretical methods that enable the reliable prediction of local spectroscopic parameters. In this work we report an approach that enables the application of multireference ab initio methods for the calculation of local zero field splitting tensors, one of the leading terms in the spin Hamiltonian for exchange-coupled systems of high nuclearity. The method referred to as local complete active space configuration interaction (L-CASCI) represents a multireference calculation with an active space composed of local orbitals of the center of interest. By successive permutation of the active space to include the localized orbitals corresponding to a particular center of the complex, all on-site parameters can be easily obtained at a high-level of theory with a corresponding low computational cost. Benchmark calculations on synthetic complexes confirm the validity of the approach. As an example of the applicability of the L-CASCI method to large systems, we determine the local anisotropy of the Mn(III) ion of the tetranuclear manganese cluster of photosystem II in both structural forms of its S2 state.

The electronic structure of the Mn/Fe cofactor identified in a new class of oxidases (R2lox) described by Andersson and Högbom [Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 5633] is reported. The R2lox protein is homologous to the small subunit of class Ic ribonucleotide reductase (R2c) but has a completely different in vivo function. Using multifrequency EPR and related pulse techniques, it is shown that the cofactor of R2lox represents an antiferromagnetically coupled Mn(III)/Fe(III) dimer linked by a μ-hydroxo/bis-μ-carboxylato bridging network. The Mn(III) ion is coordinated by a single water ligand. The R2lox cofactor is photoactive, converting into a second form (R2loxPhoto) upon visible illumination at cryogenic temperatures (77 K) that completely decays upon warming. This second, unstable form of the cofactor more closely resembles the Mn(III)/Fe(III) cofactor seen in R2c. It is shown that the two forms of the R2lox cofactor differ primarily in terms of the local site geometry and electronic state of the Mn(III) ion, as best evidenced by a reorientation of its unique (55)Mn hyperfine axis. Analysis of the metal hyperfine tensors in combination with density functional theory (DFT) calculations suggests that this change is triggered by deprotonation of the μ-hydroxo bridge. These results have important consequences for the mixed-metal R2c cofactor and the divergent chemistry R2lox and R2c perform.

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., (Copyright © 2014, American Association for the Advancement of Science.)

The EPR "split signals" represent key intermediates of the S-state cycle where the redox active D1-Tyr161 (YZ) has been oxidized by the reaction center of the photosystem II enzyme to its tyrosyl radical form, but the successive oxidation of the Mn4CaO5 cluster has not yet occurred (SiYZ˙). Here we focus on the S2YZ˙ state, which is formed en route to the final metastable state of the catalyst, the S3 state, the state which immediately precedes O-O bond formation. Quantum chemical calculations demonstrate that both isomeric forms of the S2 state, the open and closed cubane isomers, can form states with an oxidized YZ˙ residue without prior deprotonation of the Mn4CaO5 cluster. The two forms are expected to lie close in energy and retain the electronic structure and magnetic topology of the corresponding S2 state of the inorganic core. As expected, tyrosine oxidation results in a proton shift towards His190. Analysis of the electronic rearrangements that occur upon formation of the tyrosyl radical suggests that a likely next step in the catalytic cycle is the deprotonation of a terminal water ligand (W1) of the Mn4CaO5 cluster. Diamagnetic metal ion substitution is used in our calculations to obtain the molecular g-tensor of YZ˙. It is known that the gx value is a sensitive probe not only of the extent of the proton shift between the tyrosine-histidine pair, but also of the polarization environment of the tyrosine, especially about the phenolic oxygen. It is shown for PSII that this environment is determined by the Ca(2+) ion, which locates two water molecules about the phenoxyl oxygen, indirectly modulating the oxidation potential of YZ.

The S2 state of the oxygen-evolving complex of photosystem II, which consists of a Mn4O5Ca cofactor, is EPR-active, typically displaying a multiline signal, which arises from a ground spin state of total spin ST = 1/2. The precise appearance of the signal varies amongst different photosynthetic species, preparation and solvent conditions/compositions. Over the past five years, using the model species Thermosynechococcus elongatus, we have examined modifications that induce changes in the multiline signal, i.e. Ca(2+)/Sr(2+)-substitution and the binding of ammonia, to ascertain how structural perturbations of the cluster are reflected in its magnetic/electronic properties. This refined analysis, which now includes high-field (W-band) data, demonstrates that the electronic structure of the S2 state is essentially invariant to these modifications. This assessment is based on spectroscopies that examine the metal centres themselves (EPR, (55)Mn-ENDOR) and their first coordination sphere ligands ((14)N/(15)N- and (17)O-ESEEM, -HYSCORE and -EDNMR). In addition, extended quantum mechanical models from broken-symmetry DFT now reproduce all EPR, (55)Mn and (14)N experimental magnetic observables, with the inclusion of second coordination sphere ligands being crucial for accurately describing the interaction of NH3 with the Mn tetramer. These results support a mechanism of multiline heterogeneity reported for species differences and the effect of methanol [Biochim. Biophys. Acta, Bioenerg., 2011, 1807, 829], involving small changes in the magnetic connectivity of the solvent accessible outer MnA4 to the cuboidal unit Mn3O3Ca, resulting in predictable changes of the measured effective (55)Mn hyperfine tensors. Sr(2+) and NH3 replacement both affect the observed (17)O-EDNMR signal envelope supporting the assignment of O5 as the exchangeable μ-oxo bridge and it acting as the first site of substrate inclusion.

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.

The latest crystal structure of photosystem II at 1.9 Å resolution, which resolves the topology of the Mn4CaO5 oxygen evolving complex (OEC) at atomistic detail, enables a better correlation between structural features and spectroscopic properties than ever before. Building on the refined crystallographic model of the OEC and the protein, we present combined quantum mechanical/molecular mechanical (QM/MM) studies of the spectroscopic properties of the natural catalyst embedded in the protein matrix. Focusing on the S2 state of the catalytic cycle, we examine the convergence of not only structural parameters but also of the intracluster magnetic interactions in terms of exchange coupling constants and of experimentally relevant (55)Mn, (17)O, and (14)N hyperfine coupling constants with respect to QM/MM partitioning using five QM regions of increasing size. This enables us to assess the performance of the method and to probe second sphere effects by identifying amino acid residues that principally affect the spectroscopic properties of the OEC. Comparison between QM-only and QM/MM treatments reveals that whereas QM/MM models converge quickly to stable values, the QM cluster models need to incorporate significantly larger parts of the second coordination sphere and surrounding water molecules to achieve convergence for certain properties. This is mainly due to the sensitivity of the QM-only models to fluctuations in the hydrogen bonding network and ligand acidity. Additionally, a hydrogen bond that is typically omitted in QM-only treatments is shown to determine the hyperfine coupling tensor of the unique Mn(III) ion by regulating the rotation plane of the ligated D1-His332 imidazole ring, the only N-donor ligand of the OEC.

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 nature's water splitting catalyst provide important criteria for the rational design of bioinspired synthetic catalysts.

The Mn(IV)3CaO4 cubane is a structural motif present in the oxygen-evolving complex (OEC) of photosystem II and in water-oxidizing Mn/Ca layered oxides. This work investigates the magnetic and spectroscopic properties of two recently synthesized complexes and a series of idealized models that incorporate this structural unit. Magnetic interactions, accessible spin states, and (55)Mn isotropic hyperfine couplings are computed with quantum chemical methods and form the basis for structure-property correlations. Additionally, the effects of oxo-bridge protonation and one-electron reduction are examined. The calculated properties are found to be in excellent agreement with available experimental data. It is established that all synthetic and model Mn(IV)3CaO4 cubane complexes have the same high-spin S = (9)/2 ground state. The magnetic coupling conditions under which different ground spin states can be accessed are determined. Substitution of Mn(IV) magnetic centers by diamagnetic ions [e.g., Ge(IV)] allows one to "switch off" specific spin sites in order to examine the magnetic orbitals along individual Mn-Mn exchange pathways, which confirms the predominance of ferromagnetic interactions within the cubane framework. The span of the Heisenberg spin ladder is found to correlate inversely with the number of protonated oxo bridges. Energetic comparisons for protonated models show that the tris-μ-oxo bridge connecting only Mn ions in the cubane has the lowest proton affinity and that the average relaxation energy per additional proton is on the order of 18 kcal·mol(-1), thus making access to ground states other than the high-spin S = (9)/2 state in these cubanes unlikely. The relevance of these cubanes for the OEC and synthetic oxides is discussed.

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)., (Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The performance of different wave-function-based and density functional theory (DFT) methods was evaluated with respect to the prediction of relative energies for gas-phase monosaccharide isomers. A test set of 58 structures was employed, representing all forms of isomerism encountered in d-aldohexoses. The set was built from eight hexopyranose epimers by deriving subsets of isomers that include hydroxymethyl rotamers, anomers, ring conformers, furanose, and open-chain forms. Each subset of isomers spans a different energy range and involves various stereoelectronic effects. Reference energy values were obtained with coupled-cluster calculations extrapolated to the complete basis set limit, CCSD(T)/CBS. The tested CBS-extrapolated ab initio methods include various types of Møller-Plesset (MP) perturbation theory and the localized paired natural orbital coupled electron pair approach (LPNO-CEPA). Extensive benchmarking of DFT methods was carried out with 31 functionals. The results allow us to establish a hierarchy of methods that forms a reference guide for further computational studies. Among wave-function-based methods, LPNO-CEPA proved indistinguishable from CCSD(T), offering a promising alternative for a reference method that can be applied to larger systems. MP2 and SCS-MP2 follow closely, surpassing SOS-MP2 and MP3. The mPW2PLYP-D double hybrid and the Minnesota M06-2X hybrid meta-GGA are the best performing density functionals and are directly competitive with wave-function-based ab initio methods. Among the remaining functionals, B3PW91, TPSSh, mPW1PW91, and PBE0 yield the best results on average, while PBE is the best general-purpose GGA functional, surpassing meta-GGAs and several hybrids such as B3LYP. The choice of method strongly depends on the type of isomerism that needs to be considered, since many DFT methods perform well for purely conformational isomerism, but most of them fail to describe ring versus open-chain isomerism, where LYP-based GGA functionals perform particularly poorly.

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., (© 2011 American Chemical Society)

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., (© 2011 American Chemical Society)

The electronic properties of the Mn(4)O(x)Ca cluster in the S(2) state of the oxygen-evolving complex (OEC) were studied using X- and Q-band EPR and Q-band (55)Mn-ENDOR using photosystem II preparations isolated from the thermophilic cyanobacterium T. elongatus and higher plants (spinach). The data presented here show that there is very little difference between the two species. Specifically it is shown that: (i) only small changes are seen in the fitted isotropic hyperfine values, suggesting that there is no significant difference in the overall spin distribution (electronic coupling scheme) between the two species; (ii) the inferred fine-structure tensor of the only Mn(III) ion in the cluster is of the same magnitude and geometry for both species types, suggesting that the Mn(III) ion has the same coordination sphere in both sample preparations; and (iii) the data from both species are consistent with only one structural model available in the literature, namely the Siegbahn structure [Siegbahn, P. E. M. Accounts Chem. Res.2009, 42, 1871-1880, Pantazis, D. A. et al., Phys. Chem. Chem. Phys.2009, 11, 6788-6798]. These measurements were made in the presence of methanol because it confers favorable magnetic relaxation properties to the cluster that facilitate pulse-EPR techniques. In the absence of methanol the separation of the ground state and the first excited state of the spin system is smaller. For cyanobacteria this effect is minor but in plant PS II it leads to a break-down of the S(T)=½ spin model of the S(2) state. This suggests that the methanol-OEC interaction is species dependent. It is proposed that the effect of small organic solvents on the electronic structure of the cluster is to change the coupling between the outer Mn (Mn(A)) and the other three Mn ions that form the trimeric part of the cluster (Mn(B), Mn(C), Mn(D)), by perturbing the linking bis-μ-oxo bridge. The flexibility of this bridging unit is discussed with regard to the mechanism of O-O bond formation., (Copyright © 2011 Elsevier B.V. All rights reserved.)

The electronic structures of the native Mn(4)O(x)Ca cluster and the biosynthetically substituted Mn(4)O(x)Sr cluster of the oxygen evolving complex (OEC) of photosystem II (PSII) core complexes isolated from Thermosynechococcus elongatus, poised in the S(2) state, were studied by X- and Q-band CW-EPR and by pulsed Q-band (55)Mn-ENDOR spectroscopy. Both wild type and tyrosine D less mutants grown photoautotrophically in either CaCl(2) or SrCl(2) containing media were measured. The obtained CW-EPR spectra of the S(2) state displayed the characteristic, clearly noticeable differences in the hyperfine pattern of the multiline EPR signal [Boussac et al. J. Biol. Chem.2004, 279, 22809-22819]. In sharp contrast, the manganese ((55)Mn) ENDOR spectra of the Ca and Sr forms of the OEC were remarkably similar. Multifrequency simulations of the X- and Q-band CW-EPR and (55)Mn-pulsed ENDOR spectra using the Spin Hamiltonian formalism were performed to investigate this surprising result. It is shown that (i) all four manganese ions contribute to the (55)Mn-ENDOR spectra; (ii) only small changes are seen in the fitted isotropic hyperfine values for the Ca(2+) and Sr(2+) containing OEC, suggesting that there is no change in the overall spin distribution (electronic coupling scheme) upon Ca(2+)/Sr(2+) substitution; (iii) the changes in the CW-EPR hyperfine pattern can be explained by a small decrease in the anisotropy of at least two hyperfine tensors. It is proposed that modifications at the Ca(2+) site may modulate the fine structure tensor of the Mn(III) ion. DFT calculations support the above conclusions. Our data analysis also provides strong support for the notion that in the S(2) state the coordination of the Mn(III) ion is square-pyramidal (5-coordinate) or octahedral (6-coordinate) with tetragonal elongation. In addition, it is shown that only one of the currently published OEC models, the Siegbahn structure [Siegbahn, P. E. M. Acc. Chem. Res.2009, 42, 1871-1880, Pantazis, D. A. et al. Phys. Chem. Chem. Phys.2009, 11, 6788-6798], is consistent with all data presented here. These results provide important information for the structure of the OEC and the water-splitting mechanism. In particular, the 5-coordinate Mn(III) is a potential site for substrate 'water' (H(2)O, OH(-)) binding. Its location within the cuboidal structural unit, as opposed to the external 'dangler' position, may have important consequences for the mechanism of O-O bond formation.

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.

We report the synthesis of a niobium cyclopropyl complex, Tp(Me2)NbMe(c-C(3)H(5))(MeCCMe), and show that thermal loss of methane from this compound generates an intermediate that is capable of activating both aliphatic and aromatic C-H bonds. Isotopic labeling, trapping studies, a detailed kinetic analysis, and density functional theory all suggest that the active intermediate is an η(2)-cyclopropene complex formed via β-hydrogen abstraction rather than an isomeric cyclopropylidene species. C-H activation chemistry of this type represents a rather unusual reactivity pattern for η(2)-alkene complexes but is favored in this case by the strain in the C(3) ring which prevents the decomposition of the key intermediate via loss of cyclopropene.

Activation of a water molecule by the electrochemical oxidation of a Mn-aquo complex accompanied by the loss of protons is reported. The sequential (2 × 1 electron/1 proton) and direct (2 electron/2 proton) proton-coupled electrochemical oxidation of a non-porphyrinic six-coordinated Mn(II)OH 2 complex into a mononuclear Mn(O) complex is described. The intermediate Mn(III)OH 2 and Mn(III)OH complexes are electrochemically prepared and analysed. Complete deprotonation of the coordinated water molecule in the Mn(O) complex is confirmed by electrochemical data while the analysis of EXAFS data reveals a gradual shortening of an Mn-O bond upon oxidation from Mn(II)OH 2 to Mn(III)OH and Mn(O). Reactivity experiments, DFT calculations and XANES pre-edge features provide strong evidence that the bonding in Mn(O) is best characterized by a Mn(III)-oxyl description. Such oxyl species could play a crucial role in natural and artificial water splitting reactions. We provide here a synthetic example for such species, obtained by electrochemical activation of a water ligand.

A systematic study of the magnetic coupling in homovalent (III-III and IV-IV) and heterovalent (III-IV) manganese dimers as a function of the chemical identity and coordination mode of the bridging ligands is carried out with the aim of establishing a reference library of magnetostructural correlations. Emphasis is placed on rationalising the results through analysis of the superexchange pathways in terms of corresponding orbitals. Additionally, the influence of specific structural distortions on magnetic properties is explored. Consistent with chemical intuition and decades of experience, oxo bridges are shown to be efficient mediators of superexchange, primarily through pi-type pathways, whereas the introduction of bridging carboxylates inhibits the magnetic coupling of the metal centres by introducing structural distortions in the core and by reducing the antiferromagnetic contribution to the exchange. Protonation of oxo bridges is shown to induce a predictably systematic reduction in the magnitude of antiferromagnetic coupling by switching off the dominant antiferromagnetic exchange pathways. In the case of weakly coupled dimers, this can even induce a reversal of the coupling from antiferromagnetic to ferromagnetic.

The trinuclear oxo bridged manganese cluster, [Mn(IV)(3)O(4)(terpy)(terpyO(2))(2)(H(2)O)](S(2)O(8))(2) (5) (terpy = 2,2':2'',6'-terpyridine and terpyO(2) = 2,2':2'',6'-terpyridine 1,1''-dioxide), was isolated in an acidic aqueous medium from the reaction of MnSO(4), terpy, and oxone as chemical oxidant. The terpyO(2) ligands were generated in situ during the synthesis by partial oxidation of terpy. The complex crystallizes in the monoclinic space group P21/n with a = 14.251(5) A, b = 15.245(5) A, c = 24.672(5) A, alpha = 90.000(5) degrees, beta = 92.045(5) degrees, gamma = 90.000(5) degrees, and Z = 4. The triangular {Mn(IV)(3)O(4)}(4+) core observed in this complex is built up of a basal Mn(mu-O)(2)Mn unit where each Mn ion is linked to an apical Mn ion via mono(mu-O) bridges. The facial coordination of the two tridentate terpyO(2) ligands to the Mn(mu-O)(2)Mn unit allows the formation of the triangular core. 5 is also the first structurally characterized Mn complex with polypyridinyl N-oxide ligands. The variable-temperature magnetic susceptibility data for this complex, in the range of 10-300 K, are consistent with an S = 1/2 ground state and were fit using the spin Hamiltonian H(eff) with S(1) = S(2) = S(3) = 3/2, J(a) = -37 (+/-0.5) and J(b) = -53 (+/-1) cm(-1), where J(a) and J(b) are exchange constants through the mono-mu-oxo and the di-mu-oxo bridges, respectively. The doublet ground spin state of 5 is confirmed by EPR spectroscopic measurements. Density functional theory (DFT) calculations based on the broken symmetry approach reproduce the magnetic properties of 5 very well (calculated values: J(a) = -39.4 and J(b) = -55.9 cm(-1)), thus confirming the capability of this quantum chemical method for predicting the magnetic behavior of clusters involving more than two metal ions. The nature of the ground spin state of the magnetic {Mn(IV)(3)O(4)}(4+) core and the role of ancillary ligands on the magnitude of J are also discussed.