Probing the limits of accuracy in electronic structure calculations: is theory capable of results uniformly better than "chemical accuracy"?

Current limitations in electronic structure methods are discussed from the perspective of their potential to contribute to inherent uncertainties in predictions of molecular properties, with an emphasis on atomization energies (or heats of formation). The practical difficulties arising from attempts to achieve high accuracy are illustrated via two case studies: the carbon dimer (C2) and the hydroperoxyl radical (HO2). While the HO2 wave function is dominated by a single configuration, the carbon dimer involves considerable multiconfigurational character. In addition to these two molecules, statistical results will be presented for a much larger sample of molecules drawn from the Computational Results Database. The goal of this analysis will be to determine if a combination of coupled cluster theory with large 1-particle basis sets and careful incorporation of several computationally expensive smaller corrections can yield uniform agreement with experiment to better than "chemical accuracy" (+/-1 kcal/mol). In the case of HO2, the best current theoretical estimate of the zero-point-inclusive, spin-orbit corrected atomization energy (SigmaD0=166.0+/-0.3 kcal/mol) and the most recent Active Thermochemical Table (ATcT) value (165.97+/-0.06 kcal/mol) are in excellent agreement. For C2 the agreement is only slightly poorer, with theory (D0=143.7+/-0.3 kcal/mol) almost encompassing the most recent ATcT value (144.03+/-0.13 kcal/mol). For a larger collection of 68 molecules, a mean absolute deviation of 0.3 kcal/mol was found. The same high level of theory that produces good agreement for atomization energies also appears capable of predicting bond lengths to an accuracy of +/-0.001 A.

[1]  D. Dixon,et al.  Heats of Formation of Simple Boron Compounds , 1998 .

[2]  D. Dixon,et al.  Predicting the heats of formation of model hydrocarbons up to benzene , 2000 .

[3]  R. Bartlett,et al.  The full CCSDT model for molecular electronic structure , 1987 .

[4]  Trygve Helgaker,et al.  Basis-set convergence in correlated calculations on Ne, N2, and H2O , 1998 .

[5]  D. Dixon,et al.  Heats of Formation of Simple Perfluorinated Carbon Compounds , 1999 .

[6]  D. Dixon,et al.  Coupled Cluster Theory Determination of the Heats of Formation of Combustion-Related Compounds: CO, HCO, CO2, HCO2, HOCO, HC(O)OH, and HC(O)OOH , 2003 .

[7]  J. L. Dunham.  The Wentzel-Brillouin-Kramers Method of Solving the Wave Equation , 1932 .

[8]  M. Litorja,et al.  A PHOTOIONIZATION STUDY OF THE HYDROPEROXYL RADICAL, HO2, AND HYDROGEN PEROXIDE, H2O2 , 1998 .

[9]  D. Dixon,et al.  Extended benchmark studies of coupled cluster theory through triple excitations , 2001 .

[10]  K. Peterson,et al.  An examination of intrinsic errors in electronic structure methods using the Environmental Molecular Sciences Laboratory computational results database and the Gaussian-2 set , 1998 .

[11]  M. Gordon,et al.  Enthalpies of formation of gas-phase N3, N3-, N5+, and N5- from Ab initio molecular orbital theory, stability predictions for N5(+)N3(-) and N5(+)N5(-), and experimental evidence for the instability of N5(+)N3(-). , 2004, Journal of the American Chemical Society.

[12]  B. Shepler,et al.  Chemically accurate thermochemistry of cadmium: An ab initio study of Cd + XY (X = H, O, Cl, Br; Y = Cl, Br). , 2006, The journal of physical chemistry. A.

[13]  L. Curtiss,et al.  Gaussian-3X (G3X) theory : use of improved geometries, zero-point energies, and Hartree-Fock basis sets. , 2001 .

[14]  Peter J. Knowles,et al.  Perturbative corrections to account for triple excitations in closed and open shell coupled cluster theories , 1994 .

[15]  Curtis L. Janssen,et al.  New diagnostics for coupled-cluster and Møller–Plesset perturbation theory , 1998 .

[16]  David Feller,et al.  Application of systematic sequences of wave functions to the water dimer , 1992 .

[17]  W. C. Lineberger,et al.  Heat of Formation of the Hydroperoxyl Radical HOO Via Negative Ion Studies , 2002 .

[18]  Juana Vázquez,et al.  HEAT: High accuracy extrapolated ab initio thermochemistry. , 2004, The Journal of chemical physics.

[19]  Theresa L Windus,et al.  Thermodynamic properties of the C5, C6, and C8 n-alkanes from ab initio electronic structure theory. , 2005, The journal of physical chemistry. A.

[20]  Marvin Douglas,et al.  Quantum electrodynamical corrections to the fine structure of helium , 1971 .

[21]  K. Peterson,et al.  Re-examination of atomization energies for the Gaussian-2 set of molecules , 1999 .

[22]  Vincenzo Barone,et al.  Vibrational zero-point energies and thermodynamic functions beyond the harmonic approximation. , 2004, The Journal of chemical physics.

[23]  P. Mayer,et al.  Heats of formation of oxygen-containing organic free radicals from appearance energy measurements , 1991 .

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

[25]  Henry F. Schaefer,et al.  A new implementation of the full CCSDT model for molecular electronic structure , 1988 .

[26]  Gustavo E. Scuseria,et al.  The open-shell restricted Hartree—Fock singles and doubles coupled-cluster method including triple excitations CCSD (T): application to C+3 , 1991 .

[27]  D. C. Mckean,et al.  Equilibrium structures for butadiene and ethylene: compelling evidence for pi-electron delocalization in butadiene. , 2006, The journal of physical chemistry. A.

[28]  Trygve Helgaker,et al.  Basis-set convergence of correlated calculations on water , 1997 .

[29]  Robert J. Harrison,et al.  Parallel Douglas-Kroll Energy and Gradients in NWChem. Estimating Scalar Relativistic Effects Using Douglas-Kroll Contracted Basis Sets. , 2001 .

[30]  D. Dixon,et al.  Theoretical study of the heats of formation of small silicon-containing compounds , 1999 .

[31]  Mihaly Kallay,et al.  W3 theory: robust computational thermochemistry in the kJ/mol accuracy range. , 2003, Journal of Chemical Physics.

[32]  John A. Montgomery,et al.  A complete basis set model chemistry. V. Extensions to six or more heavy atoms , 1996 .

[33]  Trygve Helgaker,et al.  Highly accurate calculations of molecular electronic structure , 1999 .

[34]  Nevin Horace Oliphant,et al.  A multireference coupled-cluster method using a single-reference formalism. , 1991 .

[35]  K. Peterson.  Systematically convergent basis sets with relativistic pseudopotentials. I. Correlation consistent basis sets for the post-d group 13–15 elements , 2003 .

[36]  T. Crawford,et al.  Sources of error in electronic structure calculations on small chemical systems. , 2006, The Journal of chemical physics.

[37]  David Feller,et al.  The use of systematic sequences of wave functions for estimating the complete basis set, full configuration interaction limit in water , 1993 .

[38]  J. Olsen,et al.  Coupled-cluster connected quadruples and quintuples corrections to the harmonic vibrational frequencies and equilibrium bond distances of HF, N(2), F(2), and CO. , 2004, The Journal of chemical physics.

[39]  J. Bearden,et al.  Atomic energy levels , 1965 .

[40]  W. A. Jong,et al.  Performance of coupled cluster theory in thermochemical calculations of small halogenated compounds , 2003 .

[41]  W. A. Jong,et al.  Heats of Formation of CBr, CHBr, and CBr2 from Ab Initio Quantum Chemistry , 2002 .

[42]  H. Stoll,et al.  Systematically convergent basis sets with relativistic pseudopotentials. II. Small-core pseudopotentials and correlation consistent basis sets for the post-d group 16–18 elements , 2003 .

[43]  D. Schwenke.  The extrapolation of one-electron basis sets in electronic structure calculations: how it should work and how it can be made to work. , 2005, The Journal of chemical physics.

[44]  Krishnan Raghavachari,et al.  Gaussian‐1 theory of molecular energies for second‐row compounds , 1990 .

[45]  Timothy J. Lee,et al.  The atomization energy and proton affinity of NH3. An ab initio calibration study , 1996 .

[46]  L. Radom,et al.  G3-RAD and G3X-RAD: Modified Gaussian-3 (G3) and Gaussian-3X (G3X) procedures for radical thermochemistry , 2003 .

[47]  M. Heckert,et al.  Molecular equilibrium geometries based on coupled-cluster calculations including quadruple excitations , 2005 .

[48]  D. Feller,et al.  A CCSDT study of the effects of higher order correlation on spectroscopic constants. I. First row diatomic hydrides , 2000 .

[49]  Kirk A. Peterson,et al.  Benchmark calculations with correlated molecular wave functions. IV. The classical barrier height of the H+H2→H2+H reaction , 1994 .

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

[51]  H. Schaefer,et al.  Theoretical investigations of molecules composed only of fluorine, oxygen and nitrogen: determination of the equilibrium structures of FOOF, (NO)2 and FNNF and the transition state structure for FNNF cis-trans isomerization , 1989 .

[52]  Krishnan Raghavachari,et al.  Gaussian-3 theory using reduced Mo/ller-Plesset order , 1999 .

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

[54]  W. D. Allen,et al.  Toward subchemical accuracy in computational thermochemistry: focal point analysis of the heat of formation of NCO and [H,N,C,O] isomers. , 2004, The Journal of chemical physics.

[55]  J. Gauss,et al.  The re Structure of Cyclopropane , 2000 .

[56]  D. Z. Goodson.  Extrapolating the coupled-cluster sequence toward the full configuration-interaction limit , 2002 .

[57]  C. Zhan,et al.  Decomposition Pathways of Peroxynitrous Acid: Gas-Phase and Solution Energetics , 2002 .

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

[59]  Kaizar Amin,et al.  Introduction to Active Thermochemical Tables: Several “Key” Enthalpies of Formation Revisited† , 2004 .

[60]  P. Armentrout,et al.  Neutral and ionic metal-hydrogen and metal-carbon bond energies: Reactions of Co+, Ni+, and Cu+ with ethane, propane, methylpropane, and dimethylpropane , 1989 .

[61]  Jan M. L. Martin.  Basis set convergence study of the atomization energy, geometry, and anharmonic force field of SO2: The importance of inner polarization functions , 1998 .

[62]  Branko Ruscic,et al.  Evidence for a Lower Enthalpy of Formation of Hydroxyl Radical and a Lower Gas-Phase Bond Dissociation Energy of Water , 2001 .

[63]  Branko Ruscic,et al.  On the Enthalpy of Formation of Hydroxyl Radical and Gas-Phase Bond Dissociation Energies of Water and Hydroxyl , 2002 .

[64]  L. Curtiss,et al.  Gaussian‐1 theory: A general procedure for prediction of molecular energies , 1989 .

[65]  S. Benson,et al.  The oxidation of HI at low temperatures and the heat of formation of HO2 , 1983 .

[66]  R. Bartlett,et al.  The coupled‐cluster single, double, and triple excitation model for open‐shell single reference functions , 1990 .

[67]  Jan M. L. Martin.  Ab initio total atomization energies of small molecules — towards the basis set limit , 1996 .

[68]  Mihály Kállay,et al.  Basis-set extrapolation techniques for the accurate calculation of molecular equilibrium geometries using coupled-cluster theory. , 2006, The Journal of chemical physics.

[69]  Heat of formation of the hydroperoxyl radical HO2. A direct determination from guided ion beam studies of oxygen-methane [O2+(2.PI.g, v = 0) + CH4] reaction , 1990 .

[70]  D. Feller,et al.  Heats of Formation of CF2, FCO, and CF2O , 1998 .

[71]  Kirk A. Peterson,et al.  Accurate correlation consistent basis sets for molecular core–valence correlation effects: The second row atoms Al–Ar, and the first row atoms B–Ne revisited , 2002 .

[72]  D. Feller,et al.  Performance of CCSDT for diatomic dissociation energies , 2000 .

[73]  Juana Vázquez,et al.  High-accuracy extrapolated ab initio thermochemistry. II. Minor improvements to the protocol and a vital simplification. , 2006, The Journal of chemical physics.

[74]  Curtis L. Janssen,et al.  Concerning zero‐point vibrational energy corrections to electronic energies , 1991 .

[75]  B. Ruscic,et al.  W4 theory for computational thermochemistry: In pursuit of confident sub-kJ/mol predictions. , 2006, The Journal of chemical physics.

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

[77]  K. Peterson,et al.  Mercury and Reactive Halogens: The Thermochemistry of Hg + {Cl2, Br2, BrCl, ClO, and BrO} , 2003 .

[78]  R. Bartlett,et al.  Recursive intermediate factorization and complete computational linearization of the coupled-cluster single, double, triple, and quadruple excitation equations , 1991 .

[79]  P. Knowles,et al.  Erratum: “Coupled cluster theory for high spin, open shell reference wave functions” [ J. Chem. Phys. 99, 5219 (1993)] , 2000 .

[80]  Hans-Joachim Werner,et al.  Coupled cluster theory for high spin, open shell reference wave functions , 1993 .

[81]  Hans-Joachim Werner,et al.  A comparison of the efficiency and accuracy of the quadratic configuration interaction (QCISD), coupled cluster (CCSD), and Brueckner coupled cluster (BCCD) methods , 1992 .

[82]  D. Z. Goodson.  Improving the accuracy of ab initio methods with summation approximants and singularity analysis , 2003 .

[83]  Kirk A Peterson,et al.  Systematically convergent basis sets for transition metals. I. All-electron correlation consistent basis sets for the 3d elements Sc-Zn. , 2005, The Journal of chemical physics.

[84]  Gregor von Laszewski,et al.  Active Thermochemical Tables: thermochemistry for the 21st century , 2005 .

[85]  Benchmark Thermochemistry of the Hydroperoxyl Radical , 2004 .

[86]  R. Bartlett,et al.  The coupled‐cluster single, double, triple, and quadruple excitation method , 1992 .

[87]  Branko Ruscic,et al.  Active Thermochemical Tables: accurate enthalpy of formation of hydroperoxyl radical, HO2. , 2006, The journal of physical chemistry. A.

[88]  Bernd A. Hess,et al.  Revision of the Douglas-Kroll transformation. , 1989, Physical review. A, General physics.

[89]  C. W. Bauschlicher,et al.  A modification of the Gaussian‐2 approach using density functional theory , 1995 .

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

[91]  John A. Montgomery,et al.  A complete basis set model chemistry. IV. An improved atomic pair natural orbital method , 1994 .

[92]  Krishnan Raghavachari,et al.  Gaussian-2 theory for molecular energies of first- and second-row compounds , 1991 .

[93]  R. Urdahl,et al.  An experimental determination of the heat of formation of C2 and the CH bond dissociation energy in C2H , 1991 .

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

[95]  Review of the heat of formation of the hydroperoxyl radical , 1983 .