Abstract:A recently developed empirical dispersion correction (Grimme et al., J. Chem. Phys. 2010, 132, 154104) to standard density functional theory (DFT-D3) is implemented in the plane-wave program package VASP. The DFT-D3 implementation is compared with an implementation of the earlier DFT-D2 version (Grimme, J. Comput. Chem. 2004, 25, 1463; Grimme, J. Comput. Chem. 2006, 27, 1787). Summation of empirical pair potential terms is performed over all atom pairs in the reference cell and over atoms in shells of neighboring cells until convergence of the dispersion energy is obtained. For DFT-D3, the definition of coordination numbers has to be modified with respect to the molecular version to ensure convergence. The effect of three-center terms as implemented in the original molecular DFT-D3 version is investigated. The empirical parameters are taken from the original DFT-D3 version where they had been optimized for a reference set of small molecules. As the coordination numbers of atoms in bulk and surfaces are much larger than in the reference compounds, this effect has to be discussed. The results of test calculations for bulk properties of metals, metal oxides, benzene, and graphite indicate that the original parameters are also suitable for solid-state systems. In particular, the interlayer distance in bulk graphite and lattice constants of molecular crystals is considerably improved over standard functionals. With the molecular standard parameters (Grimme et al., J. Chem. Phys. 2010, 132, 154104; Grimme, J. Comput. Chem. 2006, 27, 1787) a slight overbinding is observed for ionic oxides where dispersion should not contribute to the bond. For simple adsorbate systems, such as Xe atoms and benzene on Ag(111), the DFT-D implementations reproduce experimental results with a similar accuracy as more sophisticated approaches based on perturbation theory (Rohlfing and Bredow, Phys. Rev. Lett. 2008, 101, 266106).

Gaussian and Plane-Wave Mixed Density Fitting for Periodic Systems.

Gaussian Basis Density Functional Theory for Systems Periodic in Two or Three Dimensions.

From plane waves to local Gaussians for the simulation of correlated periodic systems

Efficient Implementation of Density Functional Theory Based Embedding for Molecular and Periodic Systems Using Gaussian Basis Functions

Range-separated double-hybrid density-functional theory applied to periodic systems

Time-Dependent Density Functional Theory Excited-State Forces with the Gaussian and Augmented Plane Wave Method

Realizing Effective Cubic-Scaling Coulomb Hole Plus Screened Exchange Approximation in Periodic Systems Via Interpolative Separable Density Fitting with a Plane-Wave Basis Set.

Excited-State Forces with the Gaussian and Augmented Plane Wave Method for the Tamm-Dancoff Approximation of Time-Dependent Density Functional Theory

Gaussian Plane-Wave Neural Operator for Electron Density Estimation

Efficient periodic resolution-of-the-identity Hartree–Fock exchange method with k -point sampling and Gaussian basis sets

Density fitting in periodic systems: application to TDHF in diamond and oxides

Efficient Treatment of Relativistic Effects with Periodic Density Functional Methods: Energies, Gradients, and Stress Tensors

Efficient evaluation of two-center Gaussian integrals in periodic systems

Plane-Wave-Based Stochastic-Deterministic Density Functional Theory for Extended Systems

Automated optimization and uncertainty quantification of convergence parameters in plane wave density functional theory calculations

Resolution-of-identity approach to Hartree–Fock, hybrid density functionals, RPA, MP2 and GW with numeric atom-centered orbital basis functions

Exact Exchange with Range-separated Algorithm for Thermodynamic Limit of Periodic Hartree-Fock Theory

Hybrid functionals for periodic systems in the density functional tight-binding method

Implementation of empirical dispersion corrections to density functional theory for periodic systems

First Principles Excitons in Periodic Systems with Gaussian Density Fitting and Ewald Potential Functions

A Mixed Basis Density Functional Approach for One-Dimensional Systems with B-splines