Density functional theory augmented with an empirical dispersion term. Interaction energies and geometries of 80 noncovalent complexes compared with ab initio quantum mechanics calculations

Standard density functional theory (DFT) is augmented with a damped empirical dispersion term. The damping function is optimized on a small, well balanced set of 22 van der Waals (vdW) complexes and verified on a validation set of 58 vdW complexes. Both sets contain biologically relevant molecules such as nucleic acid bases. Results are in remarkable agreement with reference high‐level wave function data based on the CCSD(T) method. The geometries obtained by full gradient optimization are in very good agreement with the best available theoretical reference. In terms of the standard deviation and average errors, results including the empirical dispersion term are clearly superior to all pure density functionals investigated—B‐LYP, B3‐LYP, PBE, TPSS, TPSSh, and BH‐LYP—and even surpass the MP2/cc‐pVTZ method. The combination of empirical dispersion with the TPSS functional performs remarkably well. The most critical part of the empirical dispersion approach is the damping function. The damping parameters should be optimized for each density functional/basis set combination separately. To keep the method simple, we optimized mainly a single factor, sR, scaling globally the vdW radii. For good results, a basis set of at least triple‐ζ quality is required and diffuse functions are recommended, since the basis set superposition error seriously deteriorates the results. On average, the dispersion contribution to the interaction energy missing in the DFT functionals examined here is about 15 and 100% for the hydrogen‐bonded and stacked complexes considered, respectively. © 2006 Wiley Periodicals, Inc. J Comput Chem 28: 555–569, 2007

[1]  Krzysztof Szalewicz,et al.  Dispersion energy from density-functional theory description of monomers. , 2003, Physical review letters.

[2]  P. Hobza,et al.  Quantum chemical and statistical thermodynamic investigations of anesthetic activity. 1. The interaction between chloroform, fluoroform, cyclopropane and an O-H.cntdot..cntdot..cntdot.O hydrogen bond , 1981 .

[3]  G. Scuseria,et al.  Comparative assessment of a new nonempirical density functional: Molecules and hydrogen-bonded complexes , 2003 .

[4]  Pavel Hobza,et al.  Density functional theory and molecular clusters , 1995, J. Comput. Chem..

[5]  Stefan Grimme,et al.  Accurate description of van der Waals complexes by density functional theory including empirical corrections , 2004, J. Comput. Chem..

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

[7]  Marco Häser,et al.  Auxiliary basis sets to approximate Coulomb potentials , 1995 .

[8]  Bouke P. van Eijck,et al.  Transferable ab Initio Intermolecular Potentials. 1. Derivation from Methanol Dimer and Trimer Calculations , 1999 .

[9]  D. Lacks,et al.  Pair interactions of rare-gas atoms as a test of exchange-energy-density functionals in regions of large density gradients. , 1993, Physical review. A, Atomic, molecular, and optical physics.

[10]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[11]  Noriyuki Kurita,et al.  Adjustment of Perdew-Wang exchange functional for describing van der Waals and DNA base-stacking interactions , 2003 .

[12]  Hans W. Horn,et al.  ELECTRONIC STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS: THE PROGRAM SYSTEM TURBOMOLE , 1989 .

[13]  Qin Wu,et al.  Empirical correction to density functional theory for van der Waals interactions , 2002 .

[14]  Vincenzo Barone,et al.  Exchange functionals with improved long-range behavior and adiabatic connection methods without adjustable parameters: The mPW and mPW1PW models , 1998 .

[15]  F. Weigend,et al.  Efficient use of the correlation consistent basis sets in resolution of the identity MP2 calculations , 2002 .

[16]  Serguei Patchkovskii,et al.  An Efficient a Posteriori Treatment for Dispersion Interaction in Density-Functional-Based Tight Binding. , 2005, Journal of chemical theory and computation.

[17]  T. Walsh,et al.  Exact exchange and Wilson-Levy correlation: a pragmatic device for studying complex weakly-bonded systems. , 2005, Physical chemistry chemical physics : PCCP.

[18]  Pavel Hobza,et al.  Accurate interaction energies of hydrogen-bonded nucleic acid base pairs. , 2004, Journal of the American Chemical Society.

[19]  Dmitrii E. Makarov,et al.  van der Waals Energies in Density Functional Theory , 1998 .

[20]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[21]  Evert Jan Baerends,et al.  Density functional results for isotropic and anisotropic multipole polarizabilities and C6, C7, and C8 Van der Waals dispersion coefficients for molecules , 1997 .

[22]  J. Vondrášek,et al.  Structure and IR spectrum of phenylalanyl-glycyl-glycine tripetide in the gas-phase: IR/UV experiments, ab initio quantum chemical calculations, and molecular dynamic simulations. , 2005, Chemistry.

[23]  M Elstner,et al.  Quantum mechanics simulation of protein dynamics on long timescale , 2001, Proteins.

[24]  Xin Xu,et al.  From The Cover: The X3LYP extended density functional for accurate descriptions of nonbond interactions, spin states, and thermochemical properties. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

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

[26]  D. Svozil,et al.  Ab initio electronic structure of thymine anions. , 2005, Physical chemistry chemical physics : PCCP.

[27]  Jirí Cerný,et al.  Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs. , 2006, Physical chemistry chemical physics : PCCP.

[28]  Peter Pulay,et al.  CAN (SEMI) LOCAL DENSITY FUNCTIONAL THEORY ACCOUNT FOR THE LONDON DISPERSION FORCES , 1994 .

[29]  A. Becke,et al.  Exchange-hole dipole moment and the dispersion interaction. , 2005, The Journal of chemical physics.

[30]  D. Langreth,et al.  Van Der Waals Interactions In Density Functional Theory , 2007 .

[31]  Michal Hocek,et al.  Exceptional thermodynamic stability of DNA duplexes modified by nonpolar base analogues is due to increased stacking interactions and favorable solvation: Correlated ab initio calculations and molecular dynamics simulations. , 2006, Chemistry.

[32]  P. Hobza,et al.  Nonempirical calculations on all the 29 possible DNA base pairs , 1987 .

[33]  Yingkai Zhang,et al.  Describing van der Waals interaction in diatomic molecules with generalized gradient approximations: The role of the exchange functional , 1997 .

[34]  Marco Häser,et al.  Auxiliary basis sets to approximate Coulomb potentials (Chem. Phys. Letters 240 (1995) 283-290) , 1995 .

[35]  Kimihiko Hirao,et al.  A density-functional study on pi-aromatic interaction: benzene dimer and naphthalene dimer. , 2005, The Journal of chemical physics.

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

[37]  Michael J. Frisch,et al.  Self‐consistent molecular orbital methods 25. Supplementary functions for Gaussian basis sets , 1984 .

[38]  Saroj K. Nayak,et al.  Towards extending the applicability of density functional theory to weakly bound systems , 2001 .

[39]  A. Lesk,et al.  Determinants of a protein fold. Unique features of the globin amino acid sequences. , 1987, Journal of molecular biology.

[40]  Sandro Scandolo,et al.  Dynamical and thermal properties of polyethylene by ab initio simulation , 2000 .

[41]  Thomas A. Halgren,et al.  The representation of van der Waals (vdW) interactions in molecular mechanics force fields: potential form, combination rules, and vdW parameters , 1992 .

[42]  Ehud Gazit,et al.  A possible role for π‐stacking in the self‐assembly of amyloid fibrils , 2002, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[43]  Pavel Hobza,et al.  Potential Energy Surface of the Cytosine Dimer: MP2 Complete Basis Set Limit Interaction Energies, CCSD(T) Correction Term, and Comparison with the AMBER Force Field , 2004 .

[44]  Petros Koumoutsakos,et al.  Dispersion corrections to density functionals for water aromatic interactions. , 2004, The Journal of chemical physics.

[45]  P. Hobza,et al.  Quantum chemical and statistical thermodynamic investigations of anesthetic activity. 3. The interaction between CH4, CH3Cl, CH2Cl2, CHCl3, CCl4, and an O—H … O hydrogen bond , 1984 .

[46]  Ashcroft,et al.  Fluctuation attraction in condensed matter: A nonlocal functional approach. , 1991, Physical review. B, Condensed matter.

[47]  J. C. Slater.  A Simplification of the Hartree-Fock Method , 1951 .

[48]  Jan Kroon,et al.  TRANSFERABLE AB INITIO INTERMOLECULAR POTENTIALS. 2. VALIDATION AND APPLICATION TO CRYSTAL STRUCTURE PREDICTION , 1999 .

[49]  Donald G Truhlar,et al.  How well can new-generation density functional methods describe stacking interactions in biological systems? , 2005, Physical chemistry chemical physics : PCCP.

[50]  F. J. Luque,et al.  Nature of base stacking: reference quantum-chemical stacking energies in ten unique B-DNA base-pair steps. , 2006, Chemistry.

[51]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[52]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

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

[54]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[55]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[56]  I. Adamovic,et al.  Dynamic polarizability, dispersion coefficient C6 and dispersion energy in the effective fragment potential method , 2005 .

[57]  Michael J. Frisch,et al.  The performance of the Becke-Lee-Yang-Parr (B-LYP) density functional theory with various basis sets , 1992 .

[58]  M. Dion,et al.  van der Waals density functional for general geometries. , 2004, Physical review letters.

[59]  Pavel Hobza,et al.  On geometries of stacked and H-bonded nucleic acid base pairs determined at various DFT, MP2, and CCSD(T) levels up to the CCSD(T)/complete basis set limit level. , 2005, The Journal of chemical physics.

[60]  Dennis R. Salahub,et al.  Gaussian density functional calculations on hydrogen-bonded systems , 1992 .

[61]  J. Vondrášek,et al.  Sexual attraction in the silkworm moth. Nature of binding of bombykol in pheromone binding protein--an ab initio study. , 2003, Chemistry & biology.

[62]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[63]  A. Daniel Boese,et al.  New exchange-correlation density functionals: The role of the kinetic-energy density , 2002 .

[64]  Timothy Clark,et al.  Efficient diffuse function‐augmented basis sets for anion calculations. III. The 3‐21+G basis set for first‐row elements, Li–F , 1983 .

[65]  Trygve Helgaker,et al.  Basis-set convergence in correlated calculations on Ne, N2, and H2O , 1998 .

[66]  G. Folkers,et al.  Density functional studies on herpes simplex virus type 1 thymidine kinase–substrate interactions: The role of Tyr‐172 and Met‐128 in thymine fixation , 1998, Proteins.

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

[68]  Dennis R. Salahub,et al.  Hydrogen-bonding in glycine and malonaldehyde: Performance of the Lap1 correlation functional , 1997 .

[69]  Pavel Hobza,et al.  Unexpectedly strong energy stabilization inside the hydrophobic core of small protein rubredoxin mediated by aromatic residues: correlated ab initio quantum chemical calculations. , 2005 .

[70]  P. Hobza,et al.  Quantum chemical and statistical thermodynamic investigations of anesthetic activity. 2. The interaction between chloroform, fluoroform, and an NH...O:C hydrogen bond , 1982 .

[71]  T. Wesołowski,et al.  Gradient-free and gradient-dependent approximations in the total energy bifunctional for weakly overlapping electron densities , 2003 .

[72]  Michiel Sprik,et al.  A density‐functional study of the intermolecular interactions of benzene , 1996 .

[73]  A. Daniel Boese,et al.  The role of the basis set: Assessing density functional theory , 2003, physics/0305060.

[74]  Jirí Cerný,et al.  The X3LYP extended density functional accurately describes H-bonding but fails completely for stacking. , 2005, Physical chemistry chemical physics : PCCP.

[75]  J. Dobson.  Prospects for a van der Waals density functional , 1998 .

[76]  Giacinto Scoles,et al.  Intermolecular forces in simple systems , 1977 .

[77]  K. Borrelli,et al.  A binding mechanism in protein-nucleotide interactions: implication for U1A RNA binding. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[78]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[79]  A. Bondi.  van der Waals Volumes and Radii , 1964 .

[80]  F. Weigend,et al.  RI-MP2: first derivatives and global consistency , 1997 .

[81]  José M. Pérez-Jordá,et al.  A density-functional study of van der Waals forces: rare gas diatomics. , 1995 .

[82]  Ivano Tavernelli,et al.  Optimization of effective atom centered potentials for london dispersion forces in density functional theory. , 2004, Physical review letters.

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

[84]  Pavel Hobza,et al.  True stabilization energies for the optimal planar hydrogen-bonded and stacked structures of guanine...cytosine, adenine...thymine, and their 9- and 1-methyl derivatives: complete basis set calculations at the MP2 and CCSD(T) levels and comparison with experiment. , 2003, Journal of the American Chemical Society.

[85]  Georg Jansen,et al.  Intermolecular dispersion energies from time-dependent density functional theory , 2003 .

[86]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[87]  K. Jordan,et al.  Electronic Structure of Dipole-Bound Anions , 1998 .

[88]  John W. Hepburn,et al.  A simple but reliable method for the prediction of intermolecular potentials , 1975 .

[89]  T. Kubař,et al.  Binding of cationic and neutral phenanthridine intercalators to a DNA oligomer is controlled by dispersion energy: quantum chemical calculations and molecular mechanics simulations. , 2005, Chemistry.

[90]  Trygve Helgaker,et al.  Basis-set convergence of the energy in molecular Hartree–Fock calculations , 1999 .