thermocalc — A poor man's approach to computational thermochemistry

We present thermocalc, a Perl module to perform the automated calculation of atomization energies and heats of formation for lists of molecules. The methods used are based on density functional theory and second‐order perturbation theory to ensure that data sets of medium sized to large molecules can be run at reasonable throughput rates. The quantum chemical calculations are performed using the program package TURBOMOLE in a three‐step protocol. In a first step, a pre‐optimization of the structure and a zero‐point energy calculation are performed. As second step, a geometry optimization is being carried out, and the last step is a single point energy calculation. The level of theory used in the different steps can be modified by the user to allow for customized protocols. The performance of example protocols is investigated on different test sets of molecules. In the course of this work, a simple, but efficient one‐parameter correction term based on the shared electron numbers has been developed, which reduces the error of calculated heats of formation significantly. © 2012 Wiley Periodicals, Inc.

[1]  Martin Head-Gordon,et al.  Scaled opposite-spin second order Møller-Plesset correlation energy: an economical electronic structure method. , 2004, The Journal of chemical physics.

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

[3]  Christof Hättig,et al.  Optimization of auxiliary basis sets for RI-MP2 and RI-CC2 calculations: Core–valence and quintuple-ζ basis sets for H to Ar and QZVPP basis sets for Li to Kr , 2005 .

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

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

[6]  F. Weigend.  Accurate Coulomb-fitting basis sets for H to Rn. , 2006, Physical chemistry chemical physics : PCCP.

[7]  J. Stewart.  Optimization of parameters for semiempirical methods II. Applications , 1989 .

[8]  Clémence Corminboeuf,et al.  Systematic errors in computed alkane energies using B3LYP and other popular DFT functionals. , 2006, Organic letters.

[9]  Jan M. L. Martin.  On the performance of large Gaussian basis sets for the computation of total atomization energies , 1992 .

[10]  Holger Patzelt,et al.  RI-MP2: optimized auxiliary basis sets and demonstration of efficiency , 1998 .

[11]  W. M. Haynes.  CRC Handbook of Chemistry and Physics , 1990 .

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

[13]  Donald G Truhlar,et al.  Computational Thermochemistry: Scale Factor Databases and Scale Factors for Vibrational Frequencies Obtained from Electronic Model Chemistries. , 2010, Journal of chemical theory and computation.

[14]  Stefan Grimme,et al.  n-Alkane isodesmic reaction energy errors in density functional theory are due to electron correlation effects. , 2010, Organic letters.

[15]  J. Perdew,et al.  Density-functional approximation for the correlation energy of the inhomogeneous electron gas. , 1986, Physical review. B, Condensed matter.

[16]  L. Curtiss,et al.  Assessment of Gaussian-3 and density functional theories for a larger experimental test set , 2000 .

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

[18]  Claus Ehrhardt,et al.  Population analysis based on occupation numbers II. Relationship between shared electron numbers and bond energies and characterization of hypervalent contributions , 1985 .

[19]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[20]  S. Grimme.  Accurate calculation of the heats of formation for large main group compounds with spin-component scaled MP2 methods. , 2005, The journal of physical chemistry. A.

[21]  L. Curtiss,et al.  Assessment of Gaussian-2 and density functional theories for the computation of enthalpies of formation , 1997 .

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

[23]  S. Grimme,et al.  Theoretical thermodynamics for large molecules: walking the thin line between accuracy and computational cost. , 2008, Accounts of chemical research.

[24]  Martin W. Feyereisen,et al.  Use of approximate integrals in ab initio theory. An application in MP2 energy calculations , 1993 .

[25]  Florian Weigend,et al.  Auxiliary basis sets for main row atoms and transition metals and their use to approximate Coulomb potentials , 1997 .

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

[27]  Keith J. Laidler,et al.  A SYSTEM OF MOLECULAR THERMOCHEMISTRY FOR ORGANIC GASES AND LIQUIDS , 1956 .

[28]  A. Schäfer,et al.  Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr , 1994 .

[29]  Branko Ruscic,et al.  High-accuracy extrapolated ab initio thermochemistry. III. Additional improvements and overview. , 2008, The Journal of chemical physics.

[30]  K. Joback,et al.  ESTIMATION OF PURE-COMPONENT PROPERTIES FROM GROUP-CONTRIBUTIONS , 1987 .

[31]  Marek Sierka,et al.  Fast evaluation of the Coulomb potential for electron densities using multipole accelerated resolution of identity approximation , 2003 .

[32]  J. Almlöf,et al.  Integral approximations for LCAO-SCF calculations , 1993 .

[33]  S. Benson,et al.  Estimation of heats of formation of organic compounds by additivity methods , 1993 .

[34]  F. Weigend,et al.  Gaussian basis sets of quadruple zeta valence quality for atoms H–Kr , 2003 .

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

[36]  D. Tew,et al.  Accurate computational thermochemistry from explicitly correlated coupled-cluster theory , 2010 .

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

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

[39]  Stefan Grimme,et al.  Semiempirical GGA‐type density functional constructed with a long‐range dispersion correction , 2006, J. Comput. Chem..

[40]  Hans W. Horn,et al.  Fully optimized contracted Gaussian basis sets for atoms Li to Kr , 1992 .