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.