The Cu2O2 torture track for a real‐life system: [Cu2(btmgp)2O2]2+ oxo and peroxo species in density functional calculations†

Density functional theory (DFT) calculations of the equilibrium geometry, vibrational modes, ionization energies, electron affinities, and optical response of [Cu2(btmgp)2(μ‐O)2]2+ (oxo) and [Cu2(btmgp)2(μ‐η2:η2‐O2)]2+ (peroxo) are presented. Comprehensive benchmarking shows that the description of the oxo–peroxo energetics is still a torture track for DFT, but finds the molecular geometry to be comparatively robust with respect to changes in the exchange‐correlation functionals and basis sets. Pure functionals favor the oxo core found experimentally, whereas hybrid functionals shift the bias toward the peroxo core. Further stabilization of peroxo core results from relaxing the spin degrees of freedom using the broken‐symmetry (BS) approach. Dispersion effects, conversely, tend to favor the oxo configuration. Triple‐zeta basis sets are found to represent a sensible compromise between numerical accuracy and computational effort. Particular attention is paid to the modification of the electronic structure, optical transitions, and excited‐state energies along the transition path between the oxo and peroxo species. The excited‐state potential energy surface calculations indicate that two triplet states are involved in the transition that stabilize the BS solution. Charge decomposition and natural transition orbital analyses are used for obtaining microscopic insight into the molecular orbital interactions. Here, the crucial role of guanidine π‐interactions is highlighted for the stabilization of the Cu2O2 core. © 2015 Wiley Periodicals, Inc.

[1]  Christopher J. Cramer,et al.  Validation of density functional modeling protocols on experimental bis(μ-oxo)/μ-η2:η2-peroxo dicopper equilibria , 2007, JBIC Journal of Biological Inorganic Chemistry.

[2]  Y. Moro-oka,et al.  Copper-Dioxygen Complexes. Inorganic and Bioinorganic Perspectives , 1994 .

[3]  Per E. M. Siegbahn,et al.  Theoretical study of the catalytic mechanism of catechol oxidase , 2007, JBIC Journal of Biological Inorganic Chemistry.

[4]  Hidemi Nagao,et al.  Effective exchange integrals for open-shell species by density functional methods , 1994 .

[5]  H. Gray,et al.  Magnetic susceptibility studies of laccase and oxyhemocyanin. , 1978, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Gernot Frenking,et al.  Y-Conjugated compounds: the equilibrium geometries and electronic structures of guanidine, guanidinium cation, urea, and 1,1-diaminoethylene , 1993 .

[7]  S. Grimme,et al.  A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. , 2010, The Journal of chemical physics.

[8]  Sonja Herres-Pawlis,et al.  Bidentate guanidine ligands with ethylene spacer in copper-dioxygen chemistry: Structural characterization of bis(μ-hydroxo) dicopper complexes , 2011 .

[9]  Piero Fariselli,et al.  AB INITIO STUDY OF THE MECHANISM OF THE BINDING OF TRIPLET O2 TO HEMOCYANIN , 1996 .

[10]  P. Siegbahn,et al.  Is the bis-mu-oxo Cu2(III,III) state an intermediate in tyrosinase? , 2001, Journal of the American Chemical Society.

[11]  A. Becke Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[12]  Jackson,et al.  Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation. , 1992, Physical review. B, Condensed matter.

[13]  William B Tolman,et al.  Reactivity of dioxygen-copper systems. , 2004, Chemical reviews.

[14]  Friedhelm Bechstedt,et al.  Molecular electronic excitations calculated from a solid-state approach: Methodology and numerics , 2005 .

[15]  C. T. Lyons,et al.  Self-assembly of the oxy-tyrosinase core and the fundamental components of phenolic hydroxylation. , 2012, Nature chemistry.

[16]  Gernot Frenking,et al.  Investigation of Donor-Acceptor Interactions: A Charge Decomposition Analysis Using Fragment Molecular Orbitals , 1995 .

[17]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[18]  Richard L. Martin NATURAL TRANSITION ORBITALS , 2003 .

[19]  Christopher J. Cramer,et al.  Quantum chemical studies of molecules incorporating a Cu2O22+ core , 2009 .

[20]  Heinz Decker,et al.  Wie funktioniert die Tyrosinase? Neue Einblicke aus Modellchemie und Strukturbiologie , 2000 .

[21]  Haruki Nakamura,et al.  Theoretical studies of the magnetic couplings and the chemical indices of the biomimetic models of oxyhemocyanin and oxytyrosinase , 2009 .

[22]  W. Schmidt,et al.  Resolving the optical spectrum of water: coordination and electrostatic effects. , 2008, Physical review letters.

[23]  Yasuyuki Nakanishi,et al.  Theoretical studies on chemical bonding between Cu(II) and oxygen molecule in type 3 copper proteins , 2009 .

[24]  Axel D. Becke,et al.  Exploring the limits of gradient corrections in density functional theory , 1999, J. Comput. Chem..

[25]  Sonja Herres-Pawlis,et al.  Phenolate hydroxylation in a bis(mu-oxo)dicopper(III) complex: lessons from the guanidine/amine series. , 2009, Journal of the American Chemical Society.

[26]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[27]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[28]  Laura Gagliardi,et al.  Copper(I)-alpha-ketocarboxylate complexes: characterization and O2 reactions that yield copper-oxygen intermediates capable of hydroxylating arenes. , 2007, Journal of the American Chemical Society.

[29]  William B. Tolman,et al.  Copper-Dioxygen and Copper-Oxo Species Relevant to Copper Oxygenases and Oxidases , 2000 .

[30]  G. Kresse,et al.  From ultrasoft pseudopotentials to the projector augmented-wave method , 1999 .

[31]  Per E. M. Siegbahn,et al.  The performance of hybrid DFT for mechanisms involving transition metal complexes in enzymes , 2006, JBIC Journal of Biological Inorganic Chemistry.

[32]  Edward I. Solomon,et al.  An electronic structural comparison of copper-peroxide complexes of relevance to hemocyanin and tyrosinase active sites , 1991 .

[33]  Per E M Siegbahn,et al.  Significant van der Waals Effects in Transition Metal Complexes. , 2010, Journal of chemical theory and computation.

[34]  Per E. M. Siegbahn,et al.  A comparison of the thermodynamics of O–O bond cleavage for dicopper complexes in enzymes and synthetic systems , 2003, JBIC Journal of Biological Inorganic Chemistry.

[35]  G. Scuseria,et al.  Climbing the density functional ladder: nonempirical meta-generalized gradient approximation designed for molecules and solids. , 2003, Physical review letters.

[36]  Yasuyuki Nakanishi,et al.  Which hybrid GGA DFT is suitable for Cu2O2 systems if the spin contamination error is removed , 2010 .

[37]  Laura Gagliardi,et al.  The restricted active space followed by second-order perturbation theory method: theory and application to the study of CuO2 and Cu2O2 systems. , 2008, The Journal of chemical physics.

[38]  Luigi Cavallo,et al.  Probing the mechanism of O2 activation by a copper(I) biomimetic complex of a C-H hydroxylating copper monooxygenase. , 2009, Inorganic chemistry.

[39]  S. Grimme,et al.  A thorough benchmark of density functional methods for general main group thermochemistry, kinetics, and noncovalent interactions. , 2011, Physical chemistry chemical physics : PCCP.

[40]  Luigi Cavallo,et al.  Dancing multiplicity states supported by a carboxylated group in dicopper structures bonded to O2 , 2013, Theoretical Chemistry Accounts.

[41]  Sonja Herres-Pawlis,et al.  Bis‐μ‐oxo and μ‐η2:η2‐peroxo dicopper complexes studied within (time‐dependent) density‐functional and many‐body perturbation theory , 2013, J. Comput. Chem..

[42]  K. Pierloot,et al.  Theoretical Study of the Interconversion of O2-Binding Dicopper Complexes , 1997 .

[43]  Edward I. Solomon,et al.  Spectroscopic and Electronic Structural Studies of the Cu(III)2 Bis-μ-oxo Core and Its Relation to the Side-On Peroxo-Bridged Dimer , 1999 .

[44]  Arvi Rauk,et al.  On the calculation of multiplet energies by the hartree-fock-slater method , 1977 .

[45]  Britt Hedman,et al.  Bis(μ-oxo) dicopper(III) species of the simplest peralkylated diamine: enhanced reactivity toward exogenous substrates. , 2010, Inorganic chemistry.

[46]  Cristina Puzzarini,et al.  Theoretical models on the Cu2O2 torture track: mechanistic implications for oxytyrosinase and small-molecule analogues. , 2006, The journal of physical chemistry. A.

[47]  Xavi Ribas,et al.  Tyrosinase-like reactivity in a Cu(III)2(mu-O)2 species. , 2008, Chemistry.

[48]  Hermann Stoll,et al.  Results obtained with the correlation energy density functionals of becke and Lee, Yang and Parr , 1989 .

[49]  T. D. Stack,et al.  Structure and spectroscopy of copper-dioxygen complexes. , 2004, Chemical reviews.

[50]  Alexander Hoffmann,et al.  Hiking on the potential energy surface of a functional tyrosinase model--implications of singlet, broken-symmetry and triplet description. , 2014, Chemical communications.

[51]  Hans-Joachim Werner,et al.  Ab initio study of the O2 binding in dicopper complexes , 2005 .

[52]  Y. Matoba,et al.  Crystallographic Evidence That the Dinuclear Copper Center of Tyrosinase Is Flexible during Catalysis* , 2006, Journal of Biological Chemistry.

[53]  Matthias Bauer,et al.  BOX Ligands in Biomimetic Copper‐Mediated Dioxygen Activation: A Hemocyanin Model , 2014 .

[54]  Serge I. Gorelsky,et al.  Electronic structure and spectra of ruthenium diimine complexes by density functional theory and INDO/S. Comparison of the two methods , 2001 .

[55]  Jason A. Halfen,et al.  Reversible Cleavage and Formation of the Dioxygen O-O Bond Within a Dicopper Complex , 1996, Science.

[56]  G. Kresse,et al.  Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set , 1996 .

[57]  Per E M Siegbahn,et al.  Proton pumping mechanism in cytochrome c oxidase. , 2008, The journal of physical chemistry. A.

[58]  Blöchl,et al.  Projector augmented-wave method. , 1994, Physical review. B, Condensed matter.

[59]  Shane R. Yost,et al.  Assessment of the ΔSCF density functional theory approach for electronic excitations in organic dyes. , 2011, The Journal of chemical physics.

[60]  G L Romani,et al.  Letter: Susceptibility studies of laccase and oxyhemocyanin using an ultrasensitive magnetometer. Antiferromagnetic behavior of the type 3 copper in Rhus laccase. , 1976, Journal of the American Chemical Society.

[61]  N. Handy,et al.  A new hybrid exchange–correlation functional using the Coulomb-attenuating method (CAM-B3LYP) , 2004 .

[62]  T. D. Stack,et al.  Reaction coordinate of a functional model of tyrosinase: spectroscopic and computational characterization. , 2009, Journal of the American Chemical Society.

[63]  Richard Grunzke,et al.  Insights into the influence of dispersion correction in the theoretical treatment of guanidine‐quinoline copper(I) complexes , 2014, J. Comput. Chem..

[64]  Li Tian,et al.  Copper active sites in biology. , 2014, Chemical reviews.

[65]  Per E. M. Siegbahn,et al.  Theoretical study of the hydroxylation of phenolates by the Cu2O2(N,N′-dimethylethylenediamine)22+ complex , 2009, JBIC Journal of Biological Inorganic Chemistry.

[66]  Stefan Grimme,et al.  Effect of the damping function in dispersion corrected density functional theory , 2011, J. Comput. Chem..

[67]  Sonja Herres-Pawlis,et al.  Tuning of Copper(I)–Dioxygen Reactivity by Bis(guanidine) Ligands , 2005 .

[68]  Jens Eberhard,et al.  2-Aminopyrimidine-silver(I) based organic semiconductors: Electronic structure and optical response , 2012 .

[69]  F. Matthias Bickelhaupt,et al.  Chemistry with ADF , 2001, J. Comput. Chem..

[70]  Piotr Piecuch,et al.  Theoretical characterization of end-on and side-on peroxide coordination in ligated Cu2O2 models. , 2006, The journal of physical chemistry. A.

[71]  M. Reiher,et al.  Density matrix renormalization group calculations on relative energies of transition metal complexes and clusters. , 2008, The Journal of chemical physics.

[72]  Matthias Tamm,et al.  Synthesis and reactivity of copper(I) complexes with an ethylene-bridged bis(imidazolin-2-imine) ligand. , 2008, Dalton transactions.

[73]  Michael Vance,et al.  Tyrosinase Reactivity in a Model Complex: An Alternative Hydroxylation Mechanism , 2005, Science.

[74]  Frank Neese,et al.  A critical evaluation of DFT, including time-dependent DFT, applied to bioinorganic chemistry , 2006, JBIC Journal of Biological Inorganic Chemistry.

[75]  Sonja Herres-Pawlis,et al.  A Library of Peralkylated Bis‐guanidine Ligands for Use in Biomimetic Coordination Chemistry , 2005 .

[76]  Sonja Herres-Pawlis,et al.  Geometrical and optical benchmarking of copper guanidine–quinoline complexes: Insights from TD‐DFT and many‐body perturbation theory† , 2014, J. Comput. Chem..

[77]  Lucia Reining,et al.  Quasiparticle excitations in the photoemission spectrum of CuO from first principles: AGWstudy , 2015 .