Multireference Model Chemistries for Thermochemical Kinetics.

By combining the generalized valence bond ansatz of correlated participating orbitals (CPO) with the complete-active-space prescription for selecting configurations and with the use of multireference second order perturbation theory (MRMP2) for including dynamical correlation, we define three levels of multireference (MR) theoretical model chemistries for electronic structure calculations of chemical reaction energies and barrier heights. The three levels differ in their choice of which orbitals are considered to be participating; the choices are called nominal (nom-CPO), moderate (mod-CPO), and extended (ext-CPO). Combining any of these three choices with a method for treatment of dynamical correlation energy and a one-electron basis set yields a theoretical model chemistry. Unlike the full-valence choice of active orbitals, the CPO choices lead to active spaces that contain the orbitals needed to include important static correlation effects on chemical reactions but do not increase with the size of the nonparticipating portion of the system, and hence they remain viable computational options even for many large and complex reacting systems. The accuracies of the new levels, combined with the MG3S basis set (a partially augmented, multiply polarized valence triple-ζ basis with appropriately tight d functions for 3p-block elements) and with the fully augmented correlation-consistent aug-cc-pVTZ basis set, are assessed against a previously presented database of barrier heights for diverse reaction types. We find that nom-CPO level captures the bulk of the static correlation energy, and MRMP2/nom-CPO calculations have an average error of only 1.4 kcal/mol in barrier heights, which may be compared to 5.0 kcal/mol for single-reference MP2 theory, 2.5 kcal/mol for CCSD, and 4.1 and 1.0 kcal/mol for the B3LYP and M06-2X density functionals, respectively. The accuracy of MRMP2/CPO for transition structure bond lengths and donor-acceptor distances is excellent, with a mean unsigned error of only 0.007 Å as compared to 0.018 Å for CCSD, 0.019 Å for M06-2X, and 0.039 Å for MP2 and B3LYP. We also introduce a new multireference diagnostic, called the M diagnostic, that allows one to measure the importance of static correlation in a given reagent or transition state.

[1]  Hans-Joachim Werner,et al.  Matrix-formulated direct multiconfiguration self-consistent field and multiconfiguration reference configuration-interaction methods , 2007 .

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

[3]  Donald G. Truhlar,et al.  Dual-Level Direct Dynamics Calculations of Deuterium and Carbon-13 Kinetic Isotope Effects for the Reaction Cl + CH4 , 1998 .

[4]  Josef Paldus,et al.  General-model-space state-universal coupled-cluster theory: Connectivity conditions and explicit equations , 2003 .

[5]  D. Truhlar,et al.  Energetic and structural features of the CH4+O(3P)→CH3+OH abstraction reaction: Does perturbation theory from a multiconfiguration reference state (finally) provide a balanced treatment of transition states? , 1999 .

[6]  S. Chattopadhyay,et al.  State-specific multi-reference coupled electron-pair approximation like methods: formulation and molecular applications☆ , 2002 .

[7]  Jiří Pittner,et al.  Continuous transition between Brillouin-Wigner and Rayleigh-Schrödinger perturbation theory, generalized Bloch equation, and Hilbert space multireference coupled cluster , 2003 .

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

[9]  Wolfram Koch,et al.  Stabilities and nature of the attractive interactions in HeBeO, NeBeO, and ArBeO and a comparison with analogs NGLiF, NGBN, and NGLiH (NG = He, Ar). A theoretical investigation , 1988 .

[10]  Peter Botschwina,et al.  Stationary points of the potential surface for the reaction F− + CH3Cl → FCH3 + Cl−: Results of large‐scale coupled cluster calculations , 1997 .

[11]  L. Radom,et al.  Towards Multireference Equivalents of the G2 and G3 Methods , 2001 .

[12]  B. Roos.  The Complete Active Space Self‐Consistent Field Method and its Applications in Electronic Structure Calculations , 2007 .

[13]  Benjamin A. Ellingson,et al.  Reactions of hydrogen atom with hydrogen peroxide. , 2007, The journal of physical chemistry. A.

[14]  Angela K. Wilson,et al.  Gaussian basis sets for use in correlated molecular calculations. X. The atoms aluminum through argon revisited , 2001 .

[15]  Angela K. Wilson,et al.  Effects of Basis Set Choice upon the Atomization Energy of the Second-Row Compounds SO2, CCl, and ClO2 for B3LYP and B3PW91 , 2003 .

[16]  Donald G Truhlar,et al.  Assessment of Model Chemistries for Noncovalent Interactions. , 2006, Journal of chemical theory and computation.

[17]  John A. Pople,et al.  Nobel Lecture: Quantum chemical models , 1999 .

[18]  O. Tishchenko,et al.  Oxidation of Alkali-Metal Atoms with Nitrous Oxide: Molecular Mechanisms from First Principles Calculations , 2004 .

[19]  F. B. Brown,et al.  A new semi-empirical method of correcting large-scale configuration interaction calculations for incomplete dynamic correlation of electrons , 1985 .

[20]  Hans-Joachim Werner,et al.  Global ab initio potential energy surfaces for the ClH2 reactive system , 2000 .

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

[22]  Donald G Truhlar,et al.  A comparative assessment of the perturbative and renormalized coupled cluster theories with a noniterative treatment of triple excitations for thermochemical kinetics, including a study of basis set and core correlation effects. , 2008, The Journal of chemical physics.

[23]  Hisao Nakamura,et al.  An efficient molecular orbital approach for self-consistent calculations of molecular junctions. , 2006, The Journal of chemical physics.

[24]  A. Wilson,et al.  SO3 revisited: Impact of tight d augmented correlation consistent basis sets on atomization energy and structure , 2004 .

[25]  Michael W. Schmidt,et al.  The construction and interpretation of MCSCF wavefunctions. , 1998, Annual review of physical chemistry.

[26]  L. Curtiss,et al.  Gaussian-3 theory using scaled energies , 2000 .

[27]  Ernest R. Davidson,et al.  Considerations in constructing a multireference second‐order perturbation theory , 1994 .

[28]  K. Hirao,et al.  Energies and Dipole Moments of Excited States of Ozone and Ozone Radical Cation Using Fock Space Multireference Coupled-Cluster Analytical Response Approach , 2003 .

[29]  Donald G. Truhlar,et al.  How Well Can Hybrid Density Functional Methods Predict Transition State Geometries and Barrier Heights , 2001 .

[30]  F. B. Brown,et al.  Multireference configuration interaction treatment of potential energy surfaces: symmetric dissociation of H2O in a double-zeta basis , 1984 .

[31]  W. D. Allen,et al.  Definitive ab initio studies of model SN2 reactions CH(3)X+F- (X=F, Cl, CN, OH, SH, NH(2), PH(2)). , 2003, Chemistry.

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

[33]  Leo Radom,et al.  Molecular orbital theory of the electronic structure of organic compounds. VII. Systematic study of energies, conformations, and bond interactions , 1971 .

[34]  P. Piecuch,et al.  Improved computational strategy for the state‐selective coupled‐cluster theory with semi‐internal triexcited clusters: Potential energy surface of the HF molecule , 1995 .

[35]  K. Hirao,et al.  The hydrogen abstraction reactions: a multireference Møller–Plesset perturbation (MRMP) theory study , 2000 .

[36]  Donald G Truhlar,et al.  The 6-31B(d) basis set and the BMC-QCISD and BMC-CCSD multicoefficient correlation methods. , 2005, The journal of physical chemistry. A.

[37]  Jan M. L. Martin,et al.  TOWARDS STANDARD METHODS FOR BENCHMARK QUALITY AB INITIO THERMOCHEMISTRY :W1 AND W2 THEORY , 1999, physics/9904038.

[38]  G. Herzberg,et al.  Molecular spectra and molecular structure. Vol.3: Electronic spectra and electronic structure of polyatomic molecules , 1966 .

[39]  G. Herzberg,et al.  Spectra of diatomic molecules , 1950 .

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

[41]  Josef Paldus,et al.  Orthogonally spin-adapted multi-reference Hilbert space coupled-cluster formalism: diagrammatic formulation , 1992 .

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

[43]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[44]  Peter Pulay,et al.  Consistent generalization of the Møller-Plesset partitioning to open-shell and multiconfigurational SCF reference states in many-body perturbation theory , 1987 .

[45]  O. Tishchenko,et al.  Theoretical study on the group 2 atoms + N2O reactions. , 2005, The journal of physical chemistry. A.

[46]  D. Truhlar,et al.  The Gaussian-2 method with proper dissociation, improved accuracy, and less cost , 1999 .

[47]  T. Dunning,et al.  Benchmark Calculations with Correlated Molecular Wave Functions. 11. Energetics of the Elementary Reactions F + H2, O + H2, and H‘ + HCl , 1997 .

[48]  P. Knowles,et al.  An efficient second-order MC SCF method for long configuration expansions , 1985 .

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

[50]  C. Rohlfing,et al.  Theoretical characterization of the minimum energy path for the reaction H+O2→HO2*→HO+O , 1988 .

[51]  B. Roos,et al.  A complete active space SCF method (CASSCF) using a density matrix formulated super-CI approach , 1980 .

[52]  H. Monkhorst,et al.  Coupled-cluster method for multideterminantal reference states , 1981 .

[53]  M. Plesset,et al.  Note on an Approximation Treatment for Many-Electron Systems , 1934 .

[54]  Aleksandar Sabljić,et al.  Theoretical Study of the Mechanism and Kinetics of Gas-Phase Ozone Additions to Ethene, Fluoroethene, and Chloroethene: A Multireference Approach , 2002 .

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

[56]  L. Curtiss,et al.  Gaussian-3 (G3) theory for molecules containing first and second-row atoms , 1998 .

[57]  George S. Hammond,et al.  A Correlation of Reaction Rates , 1955 .

[58]  Martin Head-Gordon,et al.  Quantum chemistry and molecular processes , 1996 .

[59]  B. Roos,et al.  The complete active space SCF (CASSCF) method in a Newton–Raphson formulation with application to the HNO molecule , 1981 .

[60]  Michael W. Schmidt,et al.  Are atoms intrinsic to molecular electronic wavefunctions? I. The FORS model , 1982 .

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

[62]  Péter G. Szalay,et al.  New Versions of Approximately Extensive Corrected Multireference Configuration Interaction Methods , 1996 .

[63]  T. Dunning,et al.  The HSO−SOH Isomers Revisited: The Effect of Tight d Functions† , 2004 .

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

[65]  Kerstin Andersson,et al.  Second-order perturbation theory with a CASSCF reference function , 1990 .

[66]  Michael J. Frisch,et al.  Self‐consistent molecular orbital methods 25. Supplementary functions for Gaussian basis sets , 1984 .

[67]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

[68]  Rodney J. Bartlett,et al.  The multireference coupled‐cluster method in Hilbert space: An incomplete model space application to the LiH molecule , 1991 .

[69]  J. Paldus,et al.  A truncated version of reduced multireference coupled-cluster method with singles and doubles and noniterative triples: application to F2 and Ni(CO)n (n=1, 2, and 4). , 2006, The Journal of chemical physics.

[70]  Jan M.L. Martin,et al.  Benchmark ab Initio Energy Profiles for the Gas-Phase SN2 Reactions Y- + CH3X → CH3Y + X- (X,Y = F,Cl,Br). Validation of Hybrid DFT Methods , 2000 .

[71]  Jonathan Tennyson,et al.  Ab initio global potential, dipole, adiabatic, and relativistic correction surfaces for the HCN-HNC system , 2001 .

[72]  Kimihiko Hirao,et al.  Multireference Møller–Plesset perturbation treatment of potential energy curve of N2 , 1992 .

[73]  J. Pople,et al.  Effect of electron correlation of theoretical equilibrium geometries. 2. Comparison of third-order perturbation and configuration interaction results with experiment , 1982 .

[74]  Arthur Greenberg,et al.  Structure and reactivity , 1988 .

[75]  Klaus Ruedenberg,et al.  Electronic rearrangements during chemical reactions. II. Planar dissociation of ethylene , 1979 .

[76]  P. Knowles,et al.  A second order multiconfiguration SCF procedure with optimum convergence , 1985 .

[77]  D. Truhlar,et al.  Multi-coefficient Gaussian-3 method for calculating potential energy surfaces , 1999 .

[78]  Francesco A Evangelista,et al.  High-order excitations in state-universal and state-specific multireference coupled cluster theories: model systems. , 2006, The Journal of chemical physics.

[79]  G. Herzberg.  Molecular Spectra and Molecular Structure IV. Constants of Diatomic Molecules , 1939 .

[80]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[81]  Hans Lischka,et al.  Automerization reaction of cyclobutadiene and its barrier height: an ab initio benchmark multireference average-quadratic coupled cluster study. , 2006, Journal of Chemical Physics.

[82]  Henry F. Schaefer,et al.  Large multiconfiguration self-consistent-field wave functions for the ozone molecule , 1981 .

[83]  Rodney J. Bartlett,et al.  A multireference coupled‐cluster study of the ground state and lowest excited states of cyclobutadiene , 1994 .

[84]  Josef Paldus,et al.  N-reference, M-state coupled-cluster method: Merging the state-universal and reduced multireference coupled-cluster theories , 2003 .

[85]  Ludwik Adamowicz,et al.  A state-selective multireference coupled-cluster theory employing the single-reference formalism , 1993 .

[86]  Mark S. Gordon,et al.  The Parallel Implementation of a Full Configuration Interaction Program , 2003 .