Accurate description of van der Waals complexes by density functional theory including empirical corrections

Journal of Computational Chemistry - Tập 25 Số 12 - Trang 1463-1473 - 2004
Stefan Grimme1
1Organisch-Chemisches Institut der Universität Münster, Corrensstraße 40, D-48149 Münster, Germany

Tóm tắt

AbstractAn empirical method to account for van der Waals interactions in practical calculations with the density functional theory (termed DFT‐D) is tested for a wide variety of molecular complexes. As in previous schemes, the dispersive energy is described by damped interatomic potentials of the form C6R−6. The use of pure, gradient‐corrected density functionals (BLYP and PBE), together with the resolution‐of‐the‐identity (RI) approximation for the Coulomb operator, allows very efficient computations for large systems. Opposed to previous work, extended AO basis sets of polarized TZV or QZV quality are employed, which reduces the basis set superposition error to a negligible extend. By using a global scaling factor for the atomic C6 coefficients, the functional dependence of the results could be strongly reduced. The “double counting” of correlation effects for strongly bound complexes is found to be insignificant if steep damping functions are employed. The method is applied to a total of 29 complexes of atoms and small molecules (Ne, CH4, NH3, H2O, CH3F, N2, F2, formic acid, ethene, and ethine) with each other and with benzene, to benzene, naphthalene, pyrene, and coronene dimers, the naphthalene trimer, coronene · H2O and four H‐bonded and stacked DNA base pairs (AT and GC). In almost all cases, very good agreement with reliable theoretical or experimental results for binding energies and intermolecular distances is obtained. For stacked aromatic systems and the important base pairs, the DFT‐D‐BLYP model seems to be even superior to standard MP2 treatments that systematically overbind. The good results obtained suggest the approach as a practical tool to describe the properties of many important van der Waals systems in chemistry. Furthermore, the DFT‐D data may either be used to calibrate much simpler (e.g., force‐field) potentials or the optimized structures can be used as input for more accurate ab initio calculations of the interaction energies. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1463–1473, 2004

Từ khóa


Tài liệu tham khảo

Kutzelnigg W., 1978, Einführung in die Theoretische Chemie: Band 2, Die chemische Bindung

10.1002/9781119019572

Cremer D., 1998, 1706

10.1063/1.1630953

Grimme S., Chem Eur J

Parr R. G., 1989, Density‐Functional Theory of Atoms and Molecules

10.1002/3527600043

10.1063/1.1522715

10.1063/1.1637034

10.1063/1.1344891

10.1002/(SICI)1096-987X(19990115)20:1<12::AID-JCC4>3.0.CO;2-U

10.1063/1.1424928

10.1016/S0009-2614(03)00083-6

10.1073/pnas.0308730100

10.1063/1.1412004

10.1063/1.1329889

10.1021/ja011490d

10.1016/0301-0104(77)85124-0

10.1021/jp011509z

10.1063/1.1480873

10.1021/jp030587e

Jeziorski B., 1998, 1376

The composing rule used has also been tested for the interactions between second‐row hydrides ((PH3)2 (H2S)2 and (HCl)2) where theC6coefficients of the atoms differ by almost two orders of magnitude. Taking MP2/aug‐cc‐pVTZ results as reference the DFT‐D models provide an accuracy comparable to that obtained for systems composed of first‐row atoms only.

10.1103/PhysRevA.38.3098

10.1103/PhysRevB.37.785

10.1103/PhysRevB.33.8822

Perdew J. P., 1986, Phys Rev B, 34, 7046

10.1103/PhysRevLett.77.3865

Ahlrichs R.;Bär M.;Baron H.‐P.;Bauernschmitt R.;Böcker S.;Ehrig M.;Eichkorn K.;Elliott S.;Furche F.;Haase F.;Häser M.;Horn H.;Huber C.;Huniar U.;Kattannek M.;Kölmel C.;Kollwitz M.;May K.;Ochsenfeld C.;Öhm H.;Schäfer A.;Schneider U.;Treutler O.;von Arnim M.;Weigend F.;Weis P.;Weiss H.TURBOMOLE (ver. 5.6); Universität Karlsruhe 2003. See also:http://www.turbomole.com.

10.1016/0009-2614(89)85118-8

10.1016/0009-2614(69)80118-1

10.1080/00268977000101561

10.1063/1.472902

10.1063/1.467146

10.1063/1.456153

10.1063/1.1627293

10.1016/0009-2614(95)00621-A

10.1016/0301-0104(73)80059-X

10.1063/1.1679012

10.1063/1.438728

10.1007/s002140050244

The basis sets are available from the TURBOMOLE homepagehttp://www.turbomole.comvia FTP Server Button (in the subdirectories basen jbasen and cbasen).

10.1021/cr9800255

Hyla‐Krispin I., Chem Eur J

10.1063/1.475105

10.1063/1.1569242

10.1021/ja036611j

10.1021/jp012457g

10.1063/1.456207

10.1021/jp991932w

10.1021/ja025896h

Gerenkamp M., Chem Phys Lett

10.1063/1.465035

10.1021/jp022605i

10.1063/1.1630953

10.1021/jp034810j

10.1021/jp984636c

10.1021/jp001766o

Birks J. B., 1970, Photophysics of Aromatic Molecules

10.1021/jp027217v

The high accuracy of the DFTB approach for stacking interactions is frequently claimed in ref. 17. According to the standards of the present work however these statements seem to be too optimistic. The DFTB potential curves in ref. 17 (Figs. 7–10) show relatively large deviations from the MP2 reference data especially with respect to shape and the position of the minima.

10.1063/1.1626543

10.1063/1.480575

10.1063/1.462569

10.1016/0022-2860(92)87094-C

10.1016/0022-2860(93)85053-W

10.1063/1.467939

Thông tin
Thông tin xuất bản

Journal of Computational Chemistry

Tập 25 Số 12

1463-1473

Thông tin tác giả