5.7 DFT Methods for van der Waals Interactions

5.7.3 Exchange-Dipole Model (XDM)

Becke and Johnson have proposed an exchange dipole model (XDM) of dispersion.47, 442 The attractive dispersion energy arises in this model via the interaction between the instantaneous dipole moment of the exchange hole in one molecule, and the induced dipole moment in another. This is a conceptually simple yet powerful approach that has been shown to yield very accurate dispersion coefficients without fitting parameters. This allows the calculation of both intermolecular and intramolecular dispersion interactions within a single DFT framework. The implementation and validation of this method in the Q-Chem code is described in Ref. 496.

The dipole moment of the exchange hole function hσ(𝐫,𝐫) is given at point 𝐫 by

dσ(𝐫)=-𝐫-hσ(𝐫,𝐫)𝐫𝑑𝐫, (5.31)

where σ=α,β. This depends on a model of the exchange hole, and the implementation in Q-Chem uses the Becke-Roussel (BR) model.49 In most implementations the BR model, hσ is not available in analytic form and its value must be numerically at each grid point. Q-Chem developed for the first time an analytical expression for this function,496 based on non-linear interpolation and spline techniques, which greatly improves efficiency as well as the numerical stability.

Two different damping functions have been used with XDM. One of them relies only the intermolecular C6 coefficient, and its implementation in Q-Chem is denoted as “XDM6”. In this version the dispersion energy is

EvdW=AatomsB<AatomsEvdW,AB=-AatomsB<AatomsC6,ABRAB6+kC6,AB/EABcorr (5.32)

where k is a universal parameter, and EABcorr is the sum of the absolute values of the correlation energies of the free atoms A and B. The dispersion coefficients C6,AB is computed according to

C6,ij=dX2AdX2BαAαBdX2AαB+dX2BαA (5.33)

where dX2A is the square of the exchange-hole dipole moment of atom A, whose effective polarizability (in the molecule) is αA.

The XDM6 scheme can be further generalized to include higher-order dispersion coefficients, which leads to the “XDM10” model in Q-Chem:

EvdW=-AatomsB<Aatoms(C6,ABRvdW,AB6+RAB6+C8,ABRvdW,AB8+RAB8+C10,ABRvdW,AB10+RAB10). (5.34)

The higher-order dispersion coefficients are computed using higher-order multipole moments of the exchange hole.443 The quantity RvdW,AB is the sum of the effective van der Waals radii of atoms A and B,

RvdW,AB=a1Rcrit,AB+a2 (5.35)

with a critical distance

Rcrit,AB=13[(C8,ABC6,AB)1/2+(C10,ABC6,AB)1/4+(C10,ABC8,AB)1/2]. (5.36)

XDM10 contains two universal parameters, a1 and a2, whose default values of 0.83 and 1.35, respectively, were fit to reproduce intermolecular interaction energies.442 Becke later suggested several other XC functional combinations with XDM, which employ different values of a1 and a2. The user is advised to consult the recent literature for details.48, 458

As in DFT-D, the van der Waals energy is added as a post-SCF correction. Analytic gradients and Hessians are available for both XDM6 and XDM10. Additional job control and customization options are listed below.

DFTVDW_JOBNUMBER
       Basic vdW job control
TYPE:
       INTEGER
DEFAULT:
       0
OPTIONS:
       0 Do not apply the XDM scheme. 1 Add vdW as energy/gradient correction to SCF. 2 Add vDW as a DFT functional and do full SCF (this option only works with XDM6).
RECOMMENDATION:
       None

DFTVDW_METHOD
       Choose the damping function used in XDM
TYPE:
       INTEGER
DEFAULT:
       1
OPTIONS:
       1 Use Becke’s damping function including C6 term only. 2 Use Becke’s damping function with higher-order (C8 and C10) terms.
RECOMMENDATION:
       None

DFTVDW_MOL1NATOMS
       The number of atoms in the first monomer in dimer calculation
TYPE:
       INTEGER
DEFAULT:
       0
OPTIONS:
       0–Natoms
RECOMMENDATION:
       None

DFTVDW_KAI
       Damping factor k for C6-only damping function
TYPE:
       INTEGER
DEFAULT:
       800
OPTIONS:
       10–1000
RECOMMENDATION:
       None

DFTVDW_ALPHA1
       Parameter in XDM calculation with higher-order terms
TYPE:
       INTEGER
DEFAULT:
       83
OPTIONS:
       10-1000
RECOMMENDATION:
       None

DFTVDW_ALPHA2
       Parameter in XDM calculation with higher-order terms.
TYPE:
       INTEGER
DEFAULT:
       155
OPTIONS:
       10-1000
RECOMMENDATION:
       None

DFTVDW_USE_ELE_DRV
       Specify whether to add the gradient correction to the XDM energy. only valid with Becke’s C6 damping function using the interpolated BR89 model.
TYPE:
       LOGICAL
DEFAULT:
       1
OPTIONS:
       1 Use density correction when applicable. 0 Do not use this correction (for debugging purposes).
RECOMMENDATION:
       None

DFTVDW_PRINT
       Printing control for VDW code
TYPE:
       INTEGER
DEFAULT:
       1
OPTIONS:
       0 No printing. 1 Minimum printing (default) 2 Debug printing
RECOMMENDATION:
       None

Example 5.12  Sample input illustrating a frequency calculation of a vdW complex consisted of He atom and N2 molecule.

$molecule
   0 1
   He 0.000000  0.00000   3.800000
   N  0.000000  0.000000  0.546986
   N  0.000000  0.000000 -0.546986
$end

$rem
   JOBTYPE             FREQ
   IDERIV              2
   EXCHANGE            B3LYP
   INCDFT              0
   SCF_CONVERGENCE     8
   BASIS               6-31G*
   !vdw parameters settings
   DFTVDW_JOBNUMBER    1
   DFTVDW_METHOD       1
   DFTVDW_PRINT        0
   DFTVDW_KAI          800
   DFTVDW_USE_ELE_DRV  0
$end

The original XDM implementation by Becke and Johnson used Hartree-Fock exchange but XDM can be used in conjunction with GGA, meta-GGA, or hybrid functionals, or with a specific meta-GGA exchange and correlation (the BR89 exchange and BR94 correlation functionals, for example). Encouraging results have been obtained using XDM with B3LYP.496 Becke has found more recently that this model can be efficiently combined with the P86 exchange functional, with the hyper-GGA functional B05. Using XDM together with PBE exchange plus LYP correlation, or PBE exchange plus BR94 correlation, has been also found fruitful. See Refs. 458 and 716 for some recent choices in this regard.