Van der Waals Density Functional for General Geometries
TL;DR: The proposed functional, which includes van der Waals forces in a seamless fashion, is applied to rare gas and benzene dimers, where it is shown to give a realistic description.
Abstract: A scheme within density functional theory is proposed that provides a practical way to generalize to unrestricted geometries the method applied with some success to layered geometries [Phys. Rev. Lett. 91, 126402 (2003)]]. It includes van der Waals forces in a seamless fashion. By expansion to second order in a carefully chosen quantity contained in the long-range part of the correlation functional, the nonlocal correlations are expressed in terms of a density-density interaction formula. It contains a relatively simple parametrized kernel, with parameters determined by the local density and its gradient. The proposed functional is applied to rare gas and benzene dimers, where it is shown to give a realistic description.
Citations
More filters
TL;DR: It is shown by an extensive benchmark on molecular energy data that the mathematical form of the damping function in DFT‐D methods has only a minor impact on the quality of the results and BJ‐damping seems to provide a physically correct short‐range behavior of correlation/dispersion even with unmodified standard functionals.
Abstract: It is shown by an extensive benchmark on molecular energy data that the mathematical form of the damping function in DFT-D methods has only a minor impact on the quality of the results. For 12 different functionals, a standard "zero-damping" formula and rational damping to finite values for small interatomic distances according to Becke and Johnson (BJ-damping) has been tested. The same (DFT-D3) scheme for the computation of the dispersion coefficients is used. The BJ-damping requires one fit parameter more for each functional (three instead of two) but has the advantage of avoiding repulsive interatomic forces at shorter distances. With BJ-damping better results for nonbonded distances and more clear effects of intramolecular dispersion in four representative molecular structures are found. For the noncovalently-bonded structures in the S22 set, both schemes lead to very similar intermolecular distances. For noncovalent interaction energies BJ-damping performs slightly better but both variants can be recommended in general. The exception to this is Hartree-Fock that can be recommended only in the BJ-variant and which is then close to the accuracy of corrected GGAs for non-covalent interactions. According to the thermodynamic benchmarks BJ-damping is more accurate especially for medium-range electron correlation problems and only small and practically insignificant double-counting effects are observed. It seems to provide a physically correct short-range behavior of correlation/dispersion even with unmodified standard functionals. In any case, the differences between the two methods are much smaller than the overall dispersion effect and often also smaller than the influence of the underlying density functional.
14,151 citations
TL;DR: The CASTEP program as mentioned in this paper is a computer program for first principles electro-Nic structure calculations, and some of its features and capabilities are described and near-future development plans outlined.
Abstract: CASTEP Computer program / Density functional theory / Pseudopotentials / ab initio study / Plane-wave method / Computational crystallography Abstract. The CASTEP code for first principles electro- nic structure calculations will be described. A brief, non- technical overview will be given and some of the features and capabilities highlighted. Some features which are un- ique to CASTEP will be described and near-future devel- opment plans outlined.
9,884 citations
TL;DR: The re-optimization of a recently proposed long-range corrected hybrid density functional, omegaB97X-D, to include empirical atom-atom dispersion corrections yields satisfactory accuracy for thermochemistry, kinetics, and non-covalent interactions.
Abstract: We report re-optimization of a recently proposed long-range corrected (LC) hybrid density functional [J.-D. Chai and M. Head-Gordon, J. Chem. Phys., 2008, 128, 084106] to include empirical atom–atom dispersion corrections. The resulting functional, ωB97X-D yields satisfactory accuracy for thermochemistry, kinetics, and non-covalent interactions. Tests show that for non-covalent systems, ωB97X-D shows slight improvement over other empirical dispersion-corrected density functionals, while for covalent systems and kinetics it performs noticeably better. Relative to our previous functionals, such as ωB97X, the new functional is significantly superior for non-bonded interactions, and very similar in performance for bonded interactions.
9,184 citations
Journal Article•
TL;DR: Chai and Head-Gordon as discussed by the authors proposed a long-range corrected (LC) hybrid density functional with Damped Atom-Atom Dispersion corrections, which is called ωB97X-D.
Abstract: Long-Range Corrected Hybrid Density Functionals with Damped Atom-Atom Dispersion Corrections Jeng-Da Chai ∗ and Martin Head-Gordon † Department of Chemistry, University of California and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA (Dated: June 14, 2008) We report re-optimization of a recently proposed long-range corrected (LC) hybrid density func- tionals [J.-D. Chai and M. Head-Gordon, J. Chem. Phys. 128, 084106 (2008)] to include empirical atom-atom dispersion corrections. The resulting functional, ωB97X-D yields satisfactory accuracy for thermochemistry, kinetics, and non-covalent interactions. Tests show that for non-covalent sys- tems, ωB97X-D shows slight improvement over other empirical dispersion-corrected density func- tionals, while for covalent systems and kinetics, it performs noticeably better. Relative to our previous functionals, such as ωB97X, the new functional is significantly superior for non-bonded interactions, and very similar in performance for bonded interactions. I. INTRODUCTION Due to its favorable cost-to-performance ratio, Kohn- Sham density-functional theory (KS-DFT) [1, 2] has be- come the most popular electronic structure theory for large-scale ground-state systems [3–5]. Its extension for treating excited-state systems [6, 7], time-dependent den- sity functional theory (TDDFT), has also been developed to the stage where it is now very widely used. The essential ingredient of KS-DFT, the exchange- correlation energy functional E xc , remains unknown and needs to be approximated. Semi-local gradient-corrected density functionals, though successful in many applica- tions, lead to qualitative failures in some circumstances, where the accurate treatment of non-locality of exchange- correlation hole becomes crucial. These situations occur mostly in the asymptotic regions of molecular systems, such as spurious self-interaction effects upon dissociation [8, 9] and dramatic failures for long-range charge-transfer excitations [10–12]. Widely used hybrid density function- als, like B3LYP [13, 14], do not qualitatively resolve these problems. These self-interaction errors can be qualitatively re- solved using the long-range corrected (LC) hybrid density functionals [15, 16, 18], which employ 100% Hartree-Fock (HF) exchange for long-range electron-electron interac- tions. This is accomplished by a partition of unity, using erf(ωr)/r for long-range (treated by HF exchange) and erfc(ωr)/r for short-range (treated by an exchange func- tional), with the parameter ω controlling the partition- ing. Over the past five years, the LC hybrid scheme has been attracting increasing attention [15] since its compu- tational cost is comparable with standard hybrid func- tionals [13]. However, LC functionals have tended to be inferior to the best hybrids for properties such as ther- mochemistry. ∗ Electronic † Author address: jdchai@berkeley.edu to whom correspondence should be addressed. Electronic address: mhg@cchem.berkeley.edu Recently we have improved the overall accuracy at- tainable with the LC functionals by using a systematic optimization procedure [18]. One important conclusion is that optimizing LC and hybrid functionals with identical numbers of parameters in their GGA exchange and cor- relation terms leads to noticeably better results for all properties using the LC form. The resulting LC func- tional is called ωB97. Further statistically significant improvement results from re-optimizing the entire func- tional with one extra parameter corresponding to an ad- justable fraction of short-range exact exchange, defining the ωB97X functional. Independent test sets covering thermochemistry and non-covalent interactions support these conclusions. However, problems associated with the lack of non-locality of the correlation hole, such as the lack of dispersion interactions (London forces), still remain, as the semi-local correlation functionals cannot capture long-range correlation effects [19, 20]. There have been significant efforts to develop a frame- work that can account for long-range dispersion effects within DFT. Zaremba and Kohn (ZK) [21] derived an exact expression for the second-order dispersion energy in terms of the exact density-density response functions of the two separate systems. To obtain a tractable non- local dispersion functional, Dobson and Dinite (DD) [22] made local density approximations to the ZK response functions. DD’s non-local correlation functional was ob- tained independently [23] by modifying the effective den- sity defined in the earlier work of Rapcewicz and Ashcroft Starting from the formally exact expression of KS- DFT, the adiabatic connection fluctuation-dissipation theorem (ACFDT), for the ground-state exchange- correlation energy, Langreth and co-workers [25] devel- oped a so-called van der Waals density functional (vdW- DF) by making a series of reasonable approximations to yield a computationally tractable scheme. Recently, Becke and Johnson (BJ) developed a series of post-HF correlation models with a novel treatment for dispersion interactions based on the exchange-hole dipole moment [26]. The origin of dispersion claimed in the BJ models was recently questioned by Alonso, and A.
6,345 citations
TL;DR: This Account compared the performance of the M06-class functionals and one M05-class functional (M05-2X) to that of some popular functionals for diverse databases and their performance on several difficult cases.
Abstract: Although density functional theory is widely used in the computational chemistry community, the most popular density functional, B3LYP, has some serious shortcomings: (i) it is better for main-group chemistry than for transition metals; (ii) it systematically underestimates reaction barrier heights; (iii) it is inaccurate for interactions dominated by medium-range correlation energy, such as van der Waals attraction, aromatic−aromatic stacking, and alkane isomerization energies. We have developed a variety of databases for testing and designing new density functionals. We used these data to design new density functionals, called M06-class (and, earlier, M05-class) functionals, for which we enforced some fundamental exact constraints such as the uniform-electron-gas limit and the absence of self-correlation energy. Our M06-class functionals depend on spin-up and spin-down electron densities (i.e., spin densities), spin density gradients, spin kinetic energy densities, and, for nonlocal (also called hybrid)...
5,876 citations