Validation of density functional modeling protocols on experimental bis(μ-oxo)/μ-η2:η2-peroxo dicopper equilibria

The bis(μ-oxo)/μ-η2:η2-peroxo equilibria for seven supported Cu2O2 cores were studied with different hybrid and nonhybrid density functional theory models, namely, BLYP, mPWPW, TPSS, TPSSh, B3LYP, mPW1PW, and MPW1K. Supporting ligands 3,3′-iminobis(N,N-dimethylpropylamine), N,N,N′,N′,N″-pentamethyldipropylenetriamine, N-[2-(pyridin-2-yl)ethyl]-N,N,N′-trimethylpropane-1,3-diamine, bis[2-(2-pyridin-2-yl)ethyl]methylamine, bis[2-(4-methoxy-2-pyridin-2-yl)ethyl]methylamine, bis[2-(4-N,N-dimethylamino-2-pyridin-2-yl)ethyl]methylamine, and 1,4,7-triisopropyl-1,4,7-triazacyclononane were chosen on the basis of the availability of experimental data for comparison. Density functionals were examined with respect to their ability accurately to reproduce experimental properties, including, in particular, geometries and relative energies for the bis(μ-oxo) and side-on peroxo forms. While geometries from both hybrid and nonhybrid functionals were in good agreement with experiment, the incorporation of Hartree–Fock (HF) exchange in hybrid density functionals was found to have a large, degrading effect on predicted relative isomer energies. Specifically, hybrid functionals predicted the μ-η2:η2-peroxo isomer to be too stable by roughly 5–10 kcal mol−1 for each 10% of HF exchange incorporated into the model. Continuum solvation calculations predict electrostatic effects to favor bis(μ-oxo) isomers by 1–4 kcal mol−1 depending on ligand size, with larger ligands having smaller differential solvation effects. Analysis of computed molecular partition functions suggests that nonzero measured entropies of isomerization are likely to be primarily associated with interactions between molecular solutes and their first solvation shell.

[1]  J. Tomasi,et al.  Ab initio study of ionic solutions by a polarizable continuum dielectric model , 1998 .

[2]  Jacopo Tomasi,et al.  A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics , 1997 .

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

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

[5]  Thom H. Dunning,et al.  Gaussian basis sets for use in correlated molecular calculations. V. Core-valence basis sets for boron through neon , 1995 .

[6]  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 .

[7]  K. Karlin,et al.  Resonance raman investigation of equatorial ligand donor effects on the Cu(2)O(2)(2+) core in end-on and side-on mu-peroxo-dicopper(II) and bis-mu-oxo-dicopper(III) complexes. , 2003, Journal of the American Chemical Society.

[8]  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.

[9]  M. Musiał,et al.  Method of moments of coupled-cluster equations: a new formalism for designing accurate electronic structure methods for ground and excited states , 2004 .

[10]  Patrick L. Holland,et al.  Resonance Raman spectroscopy as a probe of the bis(μ-oxo)dicopper core , 2000 .

[11]  Frank Jensen,et al.  A spin correction procedure for unrestricted Hartree-Fock and Møller-Plesset wavefunctions for singlet diradicals and polyradicals , 1988 .

[12]  Jan B. F. N. Engberts,et al.  STRUCTURE AND REACTIVITY IN AQUEOUS SOLUTION , 1994 .

[13]  K. Karlin,et al.  Solvent Effects on the Conversion of Dicopper(II) μ-η2:η2-Peroxo to Bis-μ-oxo Dicopper(III) Complexes: Direct Probing of the Solvent Interaction , 2004 .

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

[15]  Michael J. McGuire,et al.  Recent advances in electronic structure theory: Method of moments of coupled-cluster equations and renormalized coupled-cluster approaches , 2002 .

[16]  J. Klinman.  Mechanisms Whereby Mononuclear Copper Proteins Functionalize Organic Substrates. , 1996, Chemical reviews.

[17]  A. Bérces.  Ligand Effects in the Models and Mimics of Oxyhemocyanin and Oxytyrosinase. A Density Functional Study of Reversible Dioxygen Binding and Reversible O-O Bond Cleavage. , 1997, Inorganic chemistry.

[18]  Michael Dolg,et al.  Energy‐adjusted ab initio pseudopotentials for the first row transition elements , 1987 .

[19]  Vincenzo Barone,et al.  Exchange functionals with improved long-range behavior and adiabatic connection methods without adjustable parameters: The mPW and mPW1PW models , 1998 .

[20]  E. C. Wilkinson,et al.  A New Intermediate in Copper Dioxygen Chemistry: Breaking the O-O Bond To Form a {Cu2(.mu.-O)2}2+ Core , 1995 .

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

[22]  A. Palmer,et al.  Oxygen Binding, Activation, and Reduction to Water by Copper Proteins. , 2001, Angewandte Chemie.

[23]  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.

[24]  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.

[25]  Karol Kowalski,et al.  Renormalized CCSD(T) and CCSD(TQ) approaches: Dissociation of the N2 triple bond , 2000 .

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

[27]  Brian B. Laird,et al.  Chemical Applications of Density-Functional Theory , 1996 .

[28]  Patrick L. Holland,et al.  Experimental Studies of the Interconversion of μ-η2:η2-Peroxo- and Bis(μ-oxo)dicopper Complexes , 1999 .

[29]  Bradley A. Smith,et al.  Ab Initio Characterization of the Isomerism between the μ-η2:η2-Peroxo- and Bis(μ-oxo)dicopper Cores , 1996 .

[30]  Adam P. Cole,et al.  Bis(mu-oxo)dicopper(III) complexes of a homologous series of simple peralkylated 1,2-diamines: steric modulation of structure, stability, and reactivity. , 2005, Inorganic chemistry.

[31]  G. Scuseria,et al.  Comparative assessment of a new nonempirical density functional: Molecules and hydrogen-bonded complexes , 2003 .

[32]  Wang,et al.  Accurate and simple analytic representation of the electron-gas correlation energy. , 1992, Physical review. B, Condensed matter.

[33]  C. Cramer,et al.  Effects of thioether substituents on the O2 reactivity of beta-diketiminate-Cu(I) complexes: probing the role of the methionine ligand in copper monooxygenases. , 2006, Journal of the American Chemical Society.

[34]  T. D. Stack,et al.  Differential Reactivity between Interconvertible Side-On Peroxo and Bis-μ-oxodicopper Isomers Using Peralkylated Diamine Ligands , 2000 .

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

[36]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

[37]  Karol Kowalski,et al.  The method of moments of coupled-cluster equations and the renormalized CCSD[T], CCSD(T), CCSD(TQ), and CCSDT(Q) approaches , 2000 .

[38]  L. Que,et al.  O2 chemistry of dicopper complexes with alkyltriamine ligands. Comparing synergistic effects on O2 binding. , 2006, Inorganic chemistry.

[39]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[40]  J. Reedijk,et al.  Bio-mimicking galactose oxidase and hemocyanin, two dioxygen-processing copper proteins. , 2004, Dalton transactions.

[41]  K. Yoshizawa,et al.  A Theoretical Study of Dioxygen Activation on the Dicopper Enzyme Model. , 1997 .

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

[43]  Warren J. Hehre,et al.  AB INITIO Molecular Orbital Theory , 1986 .

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

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

[46]  Donald G. Truhlar,et al.  Adiabatic connection for kinetics , 2000 .

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

[48]  E. Solomon,et al.  Multicopper Oxidases and Oxygenases. , 1996, Chemical reviews.

[49]  J. Tomasi,et al.  Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects , 1981 .

[50]  K. Houk,et al.  Extended Hartree—Fock (EHF) theory of chemical reactions VI: hybrid DFT and post-Hartree—Fock approaches for concerted and non-concerted transition structures of the Diels—Alder reaction , 2002 .

[51]  Edward I. Solomon,et al.  ELECTRONIC STRUCTURES OF ACTIVE SITES IN COPPER PROTEINS : CONTRIBUTIONS TO REACTIVITY , 1992 .

[52]  Jürgen Gräfenstein,et al.  On the diagnostic value of (Ŝ2) in Kohn-Sham density functional theory , 2001 .

[53]  P. Holland,et al.  Dioxygen activation by copper sites: relative stability and reactivity of (μ-η2:η2-peroxo)- and bis(μ-oxo)dicopper cores , 1999 .

[54]  K. Karlin,et al.  Contrasting copper-dioxygen chemistry arising from alike tridentate alkyltriamine copper(I) complexes. , 2002, Journal of the American Chemical Society.

[55]  W. Tolman,et al.  Bis(μ‐oxo)dimetal “Diamond” Cores in Copper and Iron Complexes Relevant to Biocatalysis , 2002 .

[56]  P. Piecuch,et al.  In Search of the Relationship between Multiple Solutions Characterizing Coupled-Cluster Theories , 2000 .

[57]  C. Cramer,et al.  Aqueous solvation free energies of ions and ion-water clusters based on an accurate value for the absolute aqueous solvation free energy of the proton. , 2006, The journal of physical chemistry. B.

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

[59]  Patrick L. Holland,et al.  Ligand macrocycle structural effects on copper-dioxygen reactivity. , 2000, Inorganic chemistry.

[60]  E. C. Wilkinson,et al.  Structural, Spectroscopic, and Theoretical Characterization of Bis(μ-oxo)dicopper Complexes, Novel Intermediates in Copper-Mediated Dioxygen Activation , 1996 .

[61]  T. D. Stack,et al.  Complexity with simplicity: a steric continuum of chelating diamines with copper(I) and dioxygen , 2003 .

[62]  Giovanni Scalmani,et al.  New developments in the polarizable continuum model for quantum mechanical and classical calculations on molecules in solution , 2002 .

[63]  W. Tolman.  MAKING AND BREAKING THE DIOXYGEN 0-0 BOND : NEW INSIGHTS FROM STUDIES OF SYNTHETIC COPPER COMPLEXES , 1997 .

[64]  H. Hayashi,et al.  Ligand effect on reversible conversion between copper(I) and bis(mu-oxo)dicopper(III) complex with a sterically hindered tetradentate tripodal ligand and monooxygenase activity of bis(mu-oxo)dicopper(III) complex. , 2003, Inorganic chemistry.

[65]  S. Mann.  Biological Inorganic Chemistry: Structure and Reactivity , 2001 .

[66]  W. Kaim,et al.  Copper—A “Modern” Bioelement , 1996 .

[67]  Jacopo Tomasi,et al.  A new definition of cavities for the computation of solvation free energies by the polarizable continuum model , 1997 .

[68]  Montserrat Corbella,et al.  Copper(II) hexaaza macrocyclic binuclear complexes obtained from the reaction of their copper(I) derivates and molecular dioxygen. , 2006, Inorganic chemistry.

[69]  K. Karlin,et al.  Copper-dioxygen adducts and the side-on peroxo dicopper(II)/bis(mu-oxo) dicopper(III) equilibrium: Significant ligand electronic effects. , 2006, Inorganic chemistry.

[70]  Jacopo Tomasi,et al.  Continuum solvation models: A new approach to the problem of solute’s charge distribution and cavity boundaries , 1997 .

[71]  D. Root,et al.  Spectroscopy of Binuclear Dioxygen Complexes , 1994 .

[72]  Akira Nakamura,et al.  A new model for dioxygen binding in hemocyanin. Synthesis, characterization, and molecular structure of the .mu.-.eta.2:.eta.2 peroxo dinuclear copper(II) complexes, [Cu(HB(3,5-R2pz)3)]2(O2) (R = isopropyl and Ph) , 1992 .

[73]  J. Perdew,et al.  Electronic structure of solids , 2021, Modern Physics with Modern Computational Methods.

[74]  K. Karlin,et al.  Oxidant types in copper–dioxygen chemistry: the ligand coordination defines the Cun-O2 structure and subsequent reactivity , 2004, JBIC Journal of Biological Inorganic Chemistry.

[75]  K. Hodgson,et al.  A Study of Solid [{Cu(MePY2)}2O2]2+ Using Resonance Raman and X-ray Absorption Spectroscopies: An Intermediate Cu2O2 Core Structure or a Solid Solution? , 1999 .