The method of dispersion correction as an add-on to standard Kohn-Sham density functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coefficients and cutoff radii that are both computed from first principles. The coefficients for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination numbers (CN). They are used to interpolate between dispersion coefficients of atoms in different chemical environments. The method only requires adjustment of two global parameters for each density functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of atomic forces. Three-body nonadditivity terms are considered. The method has been assessed on standard benchmark sets for inter- and intramolecular noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean absolute deviations for the S22 benchmark set of noncovalent interactions for 11 standard density functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C(6) coefficients also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in molecules and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems.

http://dns2.asia.edu.tw/~ysho/YSHO-English/1000%20WC/PDF/J%20Che%20Phy132,%20154104.pdf

A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu

抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1（PfPMP1）与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用，在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员，通过启动转录不同的var基因变异体为抗原变异提供了分子基础。

疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A

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.

Effect of the damping function in dispersion corrected density functional theory

We present an analysis of the performances of a parameter free density functional model (PBE0) obtained combining the so called PBE generalized gradient functional with a predefined amount of exact exchange. The results obtained for structural, thermodynamic, kinetic and spectroscopic (magnetic, infrared and electronic) properties are satisfactory and not far from those delivered by the most reliable functionals including heavy parameterization. The way in which the functional is derived and the lack of empirical parameters fitted to specific properties make the PBE0 model a widely applicable method for both quantum chemistry and condensed matter physics.

Toward reliable density functional methods without adjustable parameters: The PBE0 model

Recent extensions of the DMol3 local orbital density functional method for band structure calculations of insulating and metallic solids are described. Furthermore the method for calculating semilocal pseudopotential matrix elements and basis functions are detailed together with other unpublished parts of the methodology pertaining to gradient functionals and local orbital basis sets. The method is applied to calculations of the enthalpy of formation of a set of molecules and solids. We find that the present numerical localized basis sets yield improved results as compared to previous results for the same functionals. Enthalpies for the formation of H, N, O, F, Cl, and C, Si, S atoms from the thermodynamic reference states are calculated at the same level of theory. It is found that the performance in predicting molecular enthalpies of formation is markedly improved for the Perdew–Burke–Ernzerhof [Phys. Rev. Lett. 77, 3865 (1996)] functional.

From molecules to solids with the DMol3 approach

A Comment on the Letter by John P. Perdew, Kieron Burke, and Matthias Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). The authors of the Letter offer a Reply.

Comment on “Generalized Gradient Approximation Made Simple”

The difficulty of widely used density functionals in describing the dissociation behavior of some homonuclear and heteronuclear diatomic radicals is analyzed. It is shown that the self-interaction error of these functionals accounts for the problem—it is much larger for a system with a noninteger number of electrons than a system with an integer number of electrons. We find the condition for the erroneous dissociation behavior described by approximate density functionals: when the ionization energy of one dissociation partner differs from the electron affinity of the other partner by a small amount, the self-interaction error will lead to wrong dissociation limit. Systems with a noninteger number of electrons and hence the large amount of self-interaction error in approximate density functionals arise also in the transition states of some chemical reactions and in some charge-transfer complexes. In the course of analysis, we derive a scaling relation necessary for an exchange-correlation functional to be self-interaction free.

A challenge for density functionals: Self-interaction error increases for systems with a noninteger number of electrons

A major challenge for combined quantum mechanical and molecular mechanical methods (QM/MM) to study large molecules is how to treat the QM/MM boundary that bisects some covalent bonds. Here a pseudobond approach has been developed to solve this problem for ab initio QM/MM calculations: a one-free-valence atom with an effective core potential is constructed to replace the boundary atom of the environment part and to form a pseudobond with the boundary atom of the active part. This pseudobond, which is described only by the QM method, is designed to mimic the original bond with similar bond length and strength, and similar effects on the rest of the active part. With this pseudobond approach, some well-known deficiencies of the link atom approach have been circumvented and a well-defined potential energy surface of the whole QM/MM system has been provided. The construction of the effective core potential for the pseudobond is independent of the molecular mechanical force field and the same effective core potential is applicable to both Hartree–Fock and density functional methods. Tests on a series of molecules yield very good structural, electronic, and energetic results in comparison with the corresponding full ab initio quantum mechanical calculations.

A pseudobond approach to combining quantum mechanical and molecular mechanical methods

A new practical approach to studying enzyme reactions by combining ab initio QM/MM calculations with free energy perturbation is presented An efficient iterative optimization procedure has been developed to determine optimized structures and minimum energy paths for a system with thousands of atoms on the ab initio QM/MM potential: the small QM sub-system is optimized using a quasi-Newton minimizer in redundant internal coordinates with ab initio QM/MM calculations, while the large MM sub-system is minimized by the truncated Newton method in Cartesian coordinates with only molecular mechanical calculations The above two optimization procedures are performed iteratively until they converge With the determined minimum energy paths, free energy perturbation calculations are carried out to determine the change in free energy along the reaction coordinate Critical to the success of the iterative optimization procedure and the free energy calculations is the smooth connection between the QM and MM regions p

Free energy calculation on enzyme reactions with an efficient iterative procedure to determine minimum energy paths on a combined ab initio QM/MM potential energy surface

For a linear combination of electron densities of degenerate ground states, it is shown that the value of any energy functional is the ground state energy, if the energy functional is exact for ground state densities, size consistent, and translational invariant. The corresponding functional of kinetic and interaction energy is the linear combination of the functionals of the degenerate densities. Without invoking ensembles, it is shown that the energy functional of fractional number electrons is a series of straight lines interpolating its values at integers. These results underscore the importance of grand canonical ensemble formulation in density functional theory.

Yingkai Zhang

Papers

Comment on “Generalized Gradient Approximation Made Simple”

A challenge for density functionals: Self-interaction error increases for systems with a noninteger number of electrons

A pseudobond approach to combining quantum mechanical and molecular mechanical methods

Free energy calculation on enzyme reactions with an efficient iterative procedure to determine minimum energy paths on a combined ab initio QM/MM potential energy surface

Degenerate Ground States and a Fractional Number of Electrons in Density and Reduced Density Matrix Functional Theory