Assessing the performance of ab initio classical valence bond methods for hydrogen transfer reactions

Abstract Ab initio classical valence bond theory in terms of localized orbitals has several advantages over electronic structure methods based on canonical delocalized Hartree-Fock molecular orbitals, with two key distinctions being greater applicability to inherently multiconfigurational systems (also called multireference systems or strongly correlated systems) and great interpretability in terms of covalent and ionic valence bond configurations (also called valence bond structures). This can be especially advantageous for applications to chemical reactions. However, until now only limited work tested the quantitative accuracy for the energetics of chemical reactions. The present study provides such tests and validations for a representative test set of six barrier heights corresponding to forward and reverse barriers of three hydrogen transfer reactions. In particular, we test the valence bond self-consistent-field theory (VBSCF) and three post-VBSCF methods that use VBSCF as a reference function for adding dynamic correlation, in particular valence bond configuration interaction (VBCI), breathing orbital valence bond (BOVB), and valence bond second-order perturbation theory (VBPT2). The VBSCF method itself is, as expected, not quantitatively accurate, with a mean unsigned error (MUE) for the six barrier heights of ∼17 kcal/mol. But the post-VBSCF methods are found to be quite successful. Depending on the basis set, and valence bond structure selection we obtain MUEs (in kcal/mol) as low as 3.7 for VBCI, 4.5 for BOVB, and (using a bigger basis-set) 1.3 for VBPT2. These compare well, on the same data set, with 1.6 for coupled clusters with singles and doubles (CCSD) and 1.5 for multireference second order perturbation theory based on a complete active space self-consistent field reference function (MRPT2). We discuss the results in terms of the pros and cons of each of these methods.

[1]  M. Head‐Gordon,et al.  A fifth-order perturbation comparison of electron correlation theories , 1989 .

[2]  Donald G. Truhlar,et al.  Effectiveness of Diffuse Basis Functions for Calculating Relative Energies by Density Functional Theory , 2003 .

[3]  Sason Shaik,et al.  Classical valence bond approach by modern methods. , 2011, Chemical reviews.

[4]  D. Truhlar,et al.  The DBH24/08 Database and Its Use to Assess Electronic Structure Model Chemistries for Chemical Reaction Barrier Heights. , 2009, Journal of chemical theory and computation.

[5]  Philippe C. Hiberty,et al.  Compact and accurate valence bond functions with different orbitals for different configurations: application to the two-configuration description of F2 , 1992 .

[6]  Y. Mo,et al.  Theoretical analysis of electronic delocalization , 1998 .

[7]  Peifeng Su,et al.  Ab initio nonorthogonal valence bond methods , 2013 .

[8]  Yuchun Lin,et al.  Block-localized wavefunction (BLW) method at the density functional theory (DFT) level. , 2007, The journal of physical chemistry. A.

[9]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[10]  Wei Wu,et al.  XMVB : A program for ab initio nonorthogonal valence bond computations , 2005, J. Comput. Chem..

[11]  Donald G. Truhlar,et al.  Robust and Affordable Multicoefficient Methods for Thermochemistry and Thermochemical Kinetics: The MCCM/3 Suite and SAC/3 , 2003 .

[12]  P. Hiberty,et al.  Explicit solvation effects on the conventional resonance model for protonated imine, carbonyl, and thiocarbonyl compounds , 2010 .

[13]  F. London,et al.  Wechselwirkung neutraler Atome und homöopolare Bindung nach der Quantenmechanik , 1927 .

[14]  Y. Mo.  The resonance energy of benzene: a revisit. , 2009, The journal of physical chemistry. A.

[15]  Arieh Warshel,et al.  The Empirical Valence Bond (EVB) Method , 2004 .

[16]  P. Hiberty.  Vibronic coupling: Cage-breaking cascade. , 2011, Nature chemistry.

[17]  C. Umrigar,et al.  Quantum Monte Carlo with Jastrow-valence-bond wave functions. , 2011, The Journal of chemical physics.

[18]  Coleen T Nemes,et al.  The effect of diffuse basis functions on valence bond structural weights , 2014 .

[19]  P. Taylor,et al.  A diagnostic for determining the quality of single‐reference electron correlation methods , 2009 .

[20]  Sason Shaik,et al.  Valence bond configuration interaction: A practical ab initio valence bond method that incorporates dynamic correlation , 2002 .

[21]  R. D. Harcourt,et al.  Ab initio valence bond calculations and the spin-paired diradical character of S42+ , 1985 .

[22]  Piotr Piecuch,et al.  Single-reference, size-extensive, non-iterative coupled-cluster approaches to bond breaking and biradicals , 2006 .

[23]  Mark S. Gordon,et al.  Chapter 41 – Advances in electronic structure theory: GAMESS a decade later , 2005 .

[24]  W. Goddard,et al.  The Self-Consistent Field Equations for Generalized Valence Bond and Open-Shell Hartree—Fock Wave Functions , 1977 .

[25]  C. Landis,et al.  Valence bond concepts applied to the molecular mechanics description of molecular shapes. 4. Transition metals with pi-bonds. , 2001, Journal of the American Chemical Society.

[26]  Rosendo Valero,et al.  Perspective on Diabatic Models of Chemical Reactivity as Illustrated by the Gas-Phase S(N)2 Reaction of Acetate Ion with 1,2-Dichloroethane. , 2009, Journal of chemical theory and computation.

[27]  Wei Wu,et al.  A practical valence bond method: A configuration interaction method approach with perturbation theoretic facility , 2004, J. Comput. Chem..

[28]  Jiali Gao,et al.  Energy decomposition analysis based on a block-localized wavefunction and multistate density functional theory. , 2011, Physical chemistry chemical physics : PCCP.

[29]  A. Faessler,et al.  Quasi-elastic electroproduction of charged $ \rho$ -mesons on nucleons , 2010 .

[30]  Y. Mo,et al.  An ab initio molecular orbital-valence bond (MOVB) method for simulating chemical reactions in solution , 2000 .

[31]  Lewis-based valence bond scheme: application to the allyl cation. , 2006, The journal of physical chemistry. A.

[32]  A. Warshel.  Computer simulations of enzyme catalysis: methods, progress, and insights. , 2003, Annual review of biophysics and biomolecular structure.

[33]  P. Hiberty,et al.  Charge-shift bonding and its manifestations in chemistry. , 2009, Nature chemistry.

[34]  Kimihiko Hirao,et al.  Multireference Møller-Plesset method , 1992 .

[35]  Wei Wu,et al.  DFVB: A Density-Functional-Based Valence Bond Method. , 2012, Journal of chemical theory and computation.

[36]  P. Hiberty,et al.  Can aromaticity coexist with diradical character? An ab initio valence bond study of S2N2 and related 6π-electron four-membered rings E2N2 and E4(2+) (E=S, Se, Te). , 2012, Chemphyschem : a European journal of chemical physics and physical chemistry.

[37]  Sason Shaik,et al.  Valence bond perturbation theory. A valence bond method that incorporates perturbation theory. , 2009, The journal of physical chemistry. A.

[38]  Y. Mo,et al.  Ab initio QM/MM simulations with a molecular orbital‐valence bond (MOVB) method: application to an SN2 reaction in water , 2000 .

[39]  Y. Mo.  Two-state model based on the block-localized wave function method. , 2007, The Journal of chemical physics.

[40]  Angela K. Wilson,et al.  Gaussian basis sets for use in correlated molecular calculations. IX. The atoms gallium through krypton , 1993 .

[41]  Philippe C. Hiberty,et al.  Compact valence bond functions with breathing orbitals: Application to the bond dissociation energies of F2 and FH , 1994 .

[42]  Donald G. Truhlar,et al.  Valence bond theory for chemical dynamics , 2007, J. Comput. Chem..

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

[44]  David L. Cooper,et al.  Applications of spin-coupled valence bond theory , 1991 .

[45]  Sason Shaik,et al.  A Chemist's Guide to Valence Bond Theory , 2007 .

[46]  J. H. van Lenthe,et al.  The valence‐bond self‐consistent field method (VB–SCF): Theory and test calculations , 1983 .

[47]  Arieh Warshel,et al.  The EVB as a quantitative tool for formulating simulations and analyzing biological and chemical reactions. , 2010, Faraday discussions.

[48]  F. Weinhold,et al.  High Bond Orders in Metal-Metal Bonding , 2007, Science.

[49]  R. D. Harcourt,et al.  Ab initio valence bond calculations and the spin-paired diradical character of sulfur nitride dimer (S2N2) , 1984 .

[50]  Wei Wu,et al.  Nonorthogonal orbital based N-body reduced density matrices and their applications to valence bond theory. II. An efficient algorithm for matrix elements and analytical energy gradients in VBSCF method. , 2013, The Journal of chemical physics.

[51]  William A. Goddard,et al.  Improved Quantum Theory of Many‐Electron Systems. V. The Spin‐Coupling Optimized GI Method , 1969 .

[52]  David L. Cooper,et al.  Spin-coupled valence bond theory , 1988 .

[53]  Wei Wu,et al.  Nonorthogonal orbital based N-body reduced density matrices and their applications to valence bond theory. I. Hamiltonian matrix elements between internally contracted excited valence bond wave functions. , 2013, The Journal of chemical physics.

[54]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[55]  Björn O. Roos,et al.  Multiconfigurational Quantum Chemistry , 2016 .

[56]  Michel Dupuis,et al.  A complete active space valence bond (CASVB) method , 1996 .

[57]  D. Truhlar,et al.  What are the best affordable multi-coefficient strategies for calculating transition state geometries and barrier heights? , 2002 .

[58]  Arieh Warshel,et al.  Computer modeling of enzyme catalysis and its relationship to concepts in physical organic chemistry , 2005 .

[59]  P. Hiberty,et al.  Quadruple bonding in C2 and analogous eight-valence electron species. , 2012, Nature chemistry.

[60]  Donald G Truhlar,et al.  Representative Benchmark Suites for Barrier Heights of Diverse Reaction Types and Assessment of Electronic Structure Methods for Thermochemical Kinetics. , 2007, Journal of chemical theory and computation.

[61]  William A. Goddard,et al.  Self‐Consistent Procedures for Generalized Valence Bond Wavefunctions. Applications H3, BH, H2O, C2H6, and O2 , 1972 .

[62]  Sason Shaik,et al.  An accurate barrier for the hydrogen exchange reaction from valence bond theory: is this theory coming of age? , 2003, Chemistry.

[63]  Donald G. Truhlar,et al.  Small Representative Benchmarks for Thermochemical Calculations , 2003 .

[64]  The nature of resonance in allyl ions and radical. , 2008, The journal of physical chemistry. A.

[65]  Huaiyu Zhang,et al.  Nonorthogonal orbital based n-body reduced density matrices and their applications to valence bond theory. III. Second-order perturbation theory using valence bond self-consistent field function as reference. , 2014, The Journal of chemical physics.

[66]  William A. Goddard,et al.  The Description of Chemical Bonding From AB Initio Calculations , 1978 .

[67]  S. Levchenko,et al.  Breaking the curse of the non-dynamical correlation problem : The spin-flip method , 2007 .

[68]  Wei Wu,et al.  Nonorthogonal orbital based N-body reduced density matrices and their applications to valence bond theory. IV. The automatic implementation of the Hessian based VBSCF method. , 2014, The Journal of chemical physics.

[69]  F. London.  Zur Quantentheorie der homöopolaren Valenzzahlen , 1928 .

[70]  Yirong Mo,et al.  XMVB 2.0: A new version of Xiamen valence bond program , 2015 .

[71]  Y. Mo.  Computational evidence that hyperconjugative interactions are not responsible for the anomeric effect. , 2010, Nature chemistry.

[72]  Arieh Warshel,et al.  Computer Modeling of Chemical Reactions in Enzymes and Solutions , 1991 .

[73]  Wei Wu,et al.  Hamiltonian Matrix Correction Based Density Functional Valence Bond Method. , 2017, Journal of chemical theory and computation.

[74]  R. Bartlett,et al.  A full coupled‐cluster singles and doubles model: The inclusion of disconnected triples , 1982 .

[75]  Donald G Truhlar,et al.  Multireference Model Chemistries for Thermochemical Kinetics. , 2008, Journal of chemical theory and computation.