3 years ago

Heats of Formation of Medium-Sized Organic Compounds from Contemporary Electronic Structure Methods

Heats of Formation of Medium-Sized Organic Compounds from Contemporary Electronic Structure Methods
S. Mani Sarathy, Heng Wang, Yury Minenkov, Luigi Cavallo, Zhandong Wang
Computational electronic structure calculations are routinely undertaken to predict thermodynamic properties of various species. However, the application of highly accurate wave function theory methods, such as the “gold standard” coupled cluster approach including single, double, and partly triple excitations in perturbative fashion, CCSD(T), to large molecules is limited due to high computational cost. In this work, the promising domain based local pair natural orbital coupled cluster approach, DLPNO–CCSD(T), has been tested to reproduce 113 accurate formation enthalpies of medium-sized molecules (few dozens heavy atoms) important for bio- and combustion chemistry via the reaction based Feller–Peterson–Dixon approach. For comparison, eight density functional theory (B3LYP, B3LYP-D3, PBE0, PBE0-D3, M06, M06-2X, ωB97X-D3, and ωB97M-V) and MP2-based (B2PLYP-D3, PWPB95-D3, B2T-PLYP, B2T-PLYP-D, B2GP-PLYP, DSD-PBEP86-D3, SCS-MP2, and OO-SCS-MP2) methods have been tested. The worst performance has been obtained for the standard hybrid DFT functionals, PBE0 (mean unsigned error (MUE)/mean signed error (MSE) = 9.1/6.0 kcal/mol) and B3LYP (MUE/MSE = 13.5/–13.3 kcal/mol). The influence of an empirical dispersion correction term on these functionals’ performance is not homogeneous: B3LYP performance is improved (B3LYP-D3 (MUE/MSE = 6.0/0.8 kcal/mol)); meanwhile PBE0 performance is worse (PBE0-D3 (MUE/MSE = 14.1/13.6 kcal/mol)). The Minnesota functionals, M06 (MUE/MSE = 3.8/–2.0 kcal/mol) and M06-2X (MUE/MSE = 3.5/3.0 kcal/mol), and recently developed ωB97X-D3 (MUE/MSE = 3.2/0.2 kcal/mol) and ωB97M-V (MUE/MSE = 2.2/1.3 kcal/mol) methods provided significantly better formation enthalpies. Enthalpies of similar quality can also be obtained from some double hybrid methods (B2PLYP-D3 (MUE/MSE = 4.7/2.0 kcal/mol), PWPB95-D3 (MUE/MSE = 4.3/3.2 kcal/mol), B2T-PLYP (MUE/MSE = 4.1/–3.0 kcal/mol), and B2T-PLYP-D (MUE/MSE = 3.3/1.7 kcal/mol)). The two spin component scaled (SCS) MP2 methods resulted in even smaller errors (SCS-MP2 (MUE/MSE = 1.9/1.2 kcal/mol) and OO-SCS-MP2 (MUE/MSE = 1.6/0.1 kcal/mol)). The best performance was found for the frozen core (FC) DLPNO–CCSD(T) method with a MUE/MSE of 1.6/–1.2 kcal/mol. The performance of the DLPNO–CCSD(T) method can be further improved by running the post-SCF calculations on the B3LYP orbitals: the MUE/MSE for the DLPNO–CCSD(T,B3LYP) approximation is 1.2/–0.4 kcal/mol. We recommend the DLPNO–CCSD(T,B3LYP) method for black box applications in thermodynamics of medium-sized organic molecules when the canonical CCSD(T) calculations with basis sets of reasonable quality are prohibitively expensive.

Publisher URL: http://dx.doi.org/10.1021/acs.jctc.7b00335

DOI: 10.1021/acs.jctc.7b00335

You might also like
Discover & Discuss Important Research

Keeping up-to-date with research can feel impossible, with papers being published faster than you'll ever be able to read them. That's where Researcher comes in: we're simplifying discovery and making important discussions happen. With over 19,000 sources, including peer-reviewed journals, preprints, blogs, universities, podcasts and Live events across 10 research areas, you'll never miss what's important to you. It's like social media, but better. Oh, and we should mention - it's free.

  • Download from Google Play
  • Download from App Store
  • Download from AppInChina

Researcher displays publicly available abstracts and doesn’t host any full article content. If the content is open access, we will direct clicks from the abstracts to the publisher website and display the PDF copy on our platform. Clicks to view the full text will be directed to the publisher website, where only users with subscriptions or access through their institution are able to view the full article.