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

Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

A new density functional (DF) of the generalized gradient approximation (GGA) type for general chemistry applications termed B97-D is proposed. It is based on Becke's power-series ansatz from 1997 and is explicitly parameterized by including damped atom-pairwise dispersion corrections of the form C(6) x R(-6). A general computational scheme for the parameters used in this correction has been established and parameters for elements up to xenon and a scaling factor for the dispersion part for several common density functionals (BLYP, PBE, TPSS, B3LYP) are reported. The new functional is tested in comparison with other GGAs and the B3LYP hybrid functional on standard thermochemical benchmark sets, for 40 noncovalently bound complexes, including large stacked aromatic molecules and group II element clusters, and for the computation of molecular geometries. Further cross-validation tests were performed for organometallic reactions and other difficult problems for standard functionals. In summary, it is found that B97-D belongs to one of the most accurate general purpose GGAs, reaching, for example for the G97/2 set of heat of formations, a mean absolute deviation of only 3.8 kcal mol(-1). The performance for noncovalently bound systems including many pure van der Waals complexes is exceptionally good, reaching on the average CCSD(T) accuracy. The basic strategy in the development to restrict the density functional description to shorter electron correlation lengths scales and to describe situations with medium to large interatomic distances by damped C(6) x R(-6) terms seems to be very successful, as demonstrated for some notoriously difficult reactions. As an example, for the isomerization of larger branched to linear alkanes, B97-D is the only DF available that yields the right sign for the energy difference. From a practical point of view, the new functional seems to be quite robust and it is thus suggested as an efficient and accurate quantum chemical method for large systems where dispersion forces are of general importance.

Semiempirical GGA-type density functional constructed with a long-range dispersion correction.

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)...

Density functionals with broad applicability in chemistry.

The hydrogen bond is the most important of all directional intermolecular interactions. It is operative in determining molecular conformation, molecular aggregation, and the function of a vast number of chemical systems ranging from inorganic to biological. Research into hydrogen bonds experienced a stagnant period in the 1980s, but re-opened around 1990, and has been in rapid development since then. In terms of modern concepts, the hydrogen bond is understood as a very broad phenomenon, and it is accepted that there are open borders to other effects. There are dozens of different types of X-H.A hydrogen bonds that occur commonly in the condensed phases, and in addition there are innumerable less common ones. Dissociation energies span more than two orders of magnitude (about 0.2-40 kcal mol(-1)). Within this range, the nature of the interaction is not constant, but its electrostatic, covalent, and dispersion contributions vary in their relative weights. The hydrogen bond has broad transition regions that merge continuously with the covalent bond, the van der Waals interaction, the ionic interaction, and also the cation-pi interaction. All hydrogen bonds can be considered as incipient proton transfer reactions, and for strong hydrogen bonds, this reaction can be in a very advanced state. In this review, a coherent survey is given on all these matters.

The hydrogen bond in the solid state.

Recent gas phase measurements and high level ab initio calculations of the CH/π interactions show that the nature of the CH/π interactions is significantly different from conventional hydrogen bonds, although the CH/π interactions were often regarded as the weakest class of hydrogen bonds. A clear correlation between the binding energies of the benzene clusters with alkanes in the gas phase and the polarizabilities of alkanes was reported. The correlation suggests that the dispersion interaction is the major source of the attraction of the CH/π interactions. Ab initio calculations also show that the major source of attraction in the CH/π interactions is the dispersion interaction and the electrostatic contribution is small. The nature of the “typical CH/π interactions is close to that of van der Waals interactions, if some exceptional “activated” CH/π interactions of highly acidic C–H bonds (C–H bonds of acetylene and chloroform) are excluded. Conventional hydrogen bonds are sufficiently strong and directional due to the large contribution of the highly orientation dependent electrostatic interactions to the attraction. The hydrogen bond is important in controlling structures of molecular assemblies, since the hydrogen bond is sufficiently strong and directional. On the other hand the “typical” CH/π interactions are weak and their directionality is very weak. Although the “typical” CH/π interactions are often regarded as important interactions in controlling the structures of molecular assemblies as in the cases of conventional hydrogen bonds, the importance of the “typical” CH/π interactions is questionable.

CH/π-Interactions

Physicochemical Properties and Application of Ionic Liquids with N-P Bonds as Lithium Secondary Battery Electrolytes

The halogen bond in solution: general discussion

Preference of the monodentate contact in the CH/π interaction between an alkyl group and a single phenyl ring: Stable structures of benzene–ethane clusters

Gas-phase energy differences between the Z and E rotamers and the rotational barrier heights of methyl formate and trifluromethyl formate: an ab initio study

Seiji Tsuzuki

Papers

CH/π-Interactions

Physicochemical Properties and Application of Ionic Liquids with N-P Bonds as Lithium Secondary Battery Electrolytes

The halogen bond in solution: general discussion

Preference of the monodentate contact in the CH/π interaction between an alkyl group and a single phenyl ring: Stable structures of benzene–ethane clusters

Gas-phase energy differences between the Z and E rotamers and the rotational barrier heights of methyl formate and trifluromethyl formate: an ab initio study