Progressive Systematic Underestimation of Reaction Energies by the B3LYP Model as the Number of C−C Bonds Increases: Why Organic Chemists Should Use Multiple DFT Models for Calculations Involving Polycarbon Hydrocarbons
- 21 October 2005
- journal article
- research article
- Published by American Chemical Society (ACS) in The Journal of Organic Chemistry
- Vol. 70 (24) , 9828-9834
- https://doi.org/10.1021/jo051545k
Abstract
Computational studies of three different reaction types involving hydrocarbons (homolytic C-C bond breaking of alkanes, progressive insertions of triplet methylene into C 14 bonds of ethane, and [2+2] cyclizations of methyl-substituted alkenes to form polymethyleyclobutanes) show that the B3LYP model consistently underestimates the reaction energy, even when extremely large basis sets are employed. The error is systematic and cumulative, such that the reaction energies of reactions involving hydrocarbons with more than 4-6 C-C bonds are predicted quite poorly. Energies are underestimated for slightly and highly methyl-substitutcd cyclic and acyclic hydrocarbons, so the errors do not arise from structural issues such as steric repulsion or ring strain energy. We trace the error associated with the B3LYP approach to its consistent underestimation of the C-C bond energy. Other DFT models show this problem to lesser extents, while the MP2 method avoids it. As a consequence, we discourage the use of the B3LYP model for reaction energy calculations for organic compounds containing more than four carbon atoms. We advocate use of a collection of pure and hybrid DFT models (and ab initio models where possible) to provide computational "error bars".This publication has 37 references indexed in Scilit:
- Test of a nonempirical density functional: Short-range part of the van der Waals interaction in rare-gas dimersThe Journal of Chemical Physics, 2005
- Experimental and Computational Studies of Trialkylaluminum and Alkylaluminum Chloride Reactions with SilicaThe Journal of Physical Chemistry B, 2005
- Climbing the Density Functional Ladder: Nonempirical Meta–Generalized Gradient Approximation Designed for Molecules and SolidsPhysical Review Letters, 2003
- Assertion and validation of the performance of the B3LYP⋆ functional for the first transition metal row and the G2 test setThe Journal of Chemical Physics, 2002
- Can Fulvenes Form from Enediynes? A Systematic High-Level Computational Study on Parent and Benzannelated Enediyne and Enyne−Allene CyclizationsThe Journal of Physical Chemistry A, 2001
- A new parametrization of exchange–correlation generalized gradient approximation functionalsThe Journal of Chemical Physics, 2001
- A Critical Validation of Density Functional and Coupled-Cluster Approaches for the Calculation of EPR Hyperfine Coupling Constants in Transition Metal ComplexesThe Journal of Physical Chemistry A, 1999
- Exchange functionals with improved long-range behavior and adiabatic connection methods without adjustable parameters: The mPW and mPW1PW modelsThe Journal of Chemical Physics, 1998
- The group equivalent reaction: An improved method for determining ring strain energyJournal of Chemical Education, 1990
- Note on an Approximation Treatment for Many-Electron SystemsPhysical Review B, 1934