Which density functional is close to CCSD accuracy to describe geometry and interaction energy of small non‐covalent dimers? A benchmark study using gaussian09. Issue 15 (1st March 2013)

• Record Type: Journal Article
• Authors: Remya, Karunakaran; Suresh, Cherumuttathu H.
• Is Part Of: Journal of computational chemistry. Volume 34:Issue 15(2013)
• Page Start: 1341
• Page End: 1353
• Publication Date: 2013-03-01
• Subjects: Chemistry -- Data processing -- Periodicals; 542.85
• DOI: 10.1002/jcc.23263
• Languages: English
• ISSNs: 0192-8651

Abstract

<abstract abstract-type="main" xml:lang="en"> <title>Abstract</title> <p>A benchmark study on all possible density functional theory (DFT) methods in <italic>Gaussian09</italic> is done to locate functionals that agree well with CCSD/aug‐cc‐pVTZ geometry and Ave‐CCSD(T)/(Q‐T) interaction energy (<italic>E</italic><sub>int</sub>) for small non‐covalently interacting molecular dimers in "dispersion‐dominated" (class 1), "dipole‐induced dipole" (class 2), and "dipole‐dipole" (class 3) classes. A DFT method is recommended acceptable if the geometry showed close agreement to CCSD result (RMSD < 0.045) and <italic>E</italic><sub>int</sub> was within 80–120% accuracy. Among 382 tested functionals, 1–46% gave good geometry, 13–44% gave good <italic>E</italic><sub>int</sub>, while 1–33% satisfied geometry and energy criteria. Further screening to locate the best performing functionals for all the three classes was made by counting the acceptable values of energy and geometry given by each functionals. The meta‐generalized gradient approximation (GGA) functional M06L was the best performer with total 14 hits; seven acceptable energies and seven acceptable geometries. This was the only functional "recommended" for at least two dimers in each class. The functionals M05, B2PLYPD, B971, mPW2PLYPD, PBEB95, and CAM‐B3LYP gave 11 hits while PBEhB95, PW91B95, Wb97x, BRxVP86, BRxP86, HSE2PBE, HSEh1PBE, PBE1PBE, PBEh1PBE, and PW91TPSS gave 10 hits. Among these, M05, B971, mPW2PLYPD, Wb97x, and PW91TPSS were among the "recommended" list of at least one dimer from each class. Long‐range correction (LC) of Hirao and coworkers to exchange‐correlation functionals showed massive improvement in geometry and <italic>E</italic><sub>int</sub>. The best performing LC‐functionals were LC‐G96KCIS and LC‐PKZBPKZB. Our results predict that M06L is the most trustworthy DFT method in <italic>Gaussian09</italic> to study small non‐covalently interacting systems. © 2013 Wiley Periodicals, Inc.</p> </abstract>