Unlike earlier DFT-D methods that were strictly (atomic) pairwise-additive,
DFT-D3 includes three-body (triatomic) corrections. These terms are
significant for non-covalent complexes assembled from large
monomers,
^{
617
}
J. Chem. Theory Comput.
(2018),
14,
pp. 2955.
Link
especially those that contain a large number of
polarizable centers.
^{
279
}
Int. J. Quantum Chem.
(2014),
114,
pp. 1157.
Link
The many-body dispersion (MBD)
method of Tkatchenko et al.
^{
1091
}
Phys. Rev. Lett.
(2012),
108,
pp. 236402.
Link
^{,}
^{
38
}
J. Chem. Phys.
(2014),
140,
pp. 18A508.
Link
represents a
more general and less empirical approach that goes beyond the pairwise-additive
treatment of dispersion. This is accomplished by including $n$-body
contributions to the dispersion energy up to the number of atoms, and
polarization screening contributions to infinite order. Even in small systems
such as benzene dimer, the MBD approach consistently outperforms other popular
vdW methods.
^{
110
}
Chem. Sci.
(2016),
7,
pp. 1712.
Link
The essential idea behind MBD is to approximate the dynamic response of a
system by that of dipole-coupled quantum harmonic oscillators (QHOs), each of
which represents a fragment of the system of interest. The correlation energy
of such a system can then be evaluated exactly by diagonalizing the
corresponding Hamiltonian:
^{
449
}
Chem. Rev.
(2017),
117,
pp. 4714.
Link
$${\widehat{H}}_{\text{MBD}}=\frac{1}{2}\sum _{A}^{\text{atoms}}{\widehat{\nabla}}_{{\xi}_{A}}^{2}+\frac{1}{2}\sum _{A}^{\text{atoms}}{\omega}_{A}^{2}{\xi}_{A}^{2}+\frac{1}{2}\sum _{A,B}^{\text{atoms}}{\omega}_{A}{\omega}_{B}{({\alpha}_{A}^{0}{\alpha}_{B}^{0})}^{1/2}{\xi}_{A}{T}_{AB}{\xi}_{B}.$$ | (5.50) |
Here, ${\xi}_{A}={m}_{A}^{1/2}|{\mathbf{r}}_{A}-{\mathbf{R}}_{A}|$ is the mass-weight displacement of oscillator $A$ from its center ${\mathbf{R}}_{A}$, ${\omega}_{A}$ is the characteristic frequency, and ${\alpha}_{A}^{0}$ is the static polarizability. ${T}_{AB}$ is the dipole potential between the oscillators $A$ and $B$. The MBD Hamiltonian is obtained through coarse-graining of the long-range correlation (through the long-range dipole tensor ${\mathbf{T}}_{\mathrm{lr}}$) and approximating the short-range polarizability via the adiabatic connection fluctuation-dissipation formula:
$${E}_{c,\text{lr}}^{\mathrm{MBD}}=-\sum _{n=2}^{\mathrm{\infty}}\frac{{(-1)}^{n}}{n}{\int}_{0}^{\mathrm{\infty}}\frac{\mathrm{d}u}{2\pi}\sum _{AB}\text{tr}[\u27e8{({\alpha}_{\mathrm{eff}}{\mathbf{T}}_{\mathrm{lr}})}^{n}{\u27e9}_{AB}(iu)].$$ | (5.51) |
This approximation expresses the dynamic polarizability ${\alpha}_{\mathrm{eff}}$ (of a given fragment) in terms of the polarizability of the corresponding QHO,
$${\alpha}_{A}^{\mathrm{QHO}}(u)=\frac{{q}_{A}^{2}}{{m}_{A}({\omega}_{A}^{2}-{u}^{2}-i\delta u)}$$ | (5.52) |
in which ${q}_{A}$ is the charge, ${m}_{A}$ the mass, and ${\omega}_{A}$ the characteristic frequency of the oscillator. The integration in the frequency domain in Eq. (5.51) can be done analytically, leading to the so-called plasmon pole formula for the correlation energy,
$${E}_{c}=\frac{1}{2}\sum _{p=1}^{3N}({\overline{\omega}}_{p}-{\omega}_{p})$$ | (5.53) |
in which $N$ is the number of fragments, ${\overline{\omega}}_{p}$ are the frequencies of the interacting (dipole-coupled) system, and ${\omega}_{p}$ are the frequencies of the non-interacting system (i.e., the collection of independent QHOs). The sum runs over all $3N$ characteristic frequencies of the system.
A particular method within the MBD framework is defined by the models for the
static polarizability (${\alpha}_{\text{eff},A}^{0}$), the non-interacting
characteristic frequencies (${\omega}_{A}$), and the damping function [$f(R)$] used
to define ${\mathbf{T}}_{\text{lr}}$. In Q-Chem, the MBD method is implemented
following the “MBD@rsSCS” approach, where “rsSCS” stands for
range-separated self-consistent screening.
^{
38
}
J. Chem. Phys.
(2014),
140,
pp. 18A508.
Link
In this
approach, ${\alpha}_{\text{eff},A}^{0}$ is obtained in a two-step process:
The free-volume scaling approach is applied to the free-atom polarizabilities, using the Hirshfeld-partitioned molecular electron density. This is the same procedure used in the TS-vdW method described in Section 5.7.4.
The short-range atomic polarizabilities ${\alpha}_{\mathrm{sr},AB}(iu)$ are obtained by applying a Dyson-like screening on only the short range part of the polarizabilities. The same range-separation will later be used to define ${\mathbf{T}}_{\mathrm{lr}}$.
The short-range atomic polarizabilities are summed up along one fragment coordinate to obtain the local effective dynamic polarizability, i.e., ${\alpha}_{\mathrm{eff},A}^{0}={\sum}_{B}{\alpha}_{\mathrm{sr},AB}$, and are then spherically averaged. The range-separation (damping) function $f(R)$ used to construct ${\alpha}_{\mathrm{sr},AB}(iu)$ and ${\mathbf{T}}_{\mathrm{lr}}$ is the same as that in Eq. (5.45), except with $d=6$ instead of $d=20$, and again ${s}_{R}$ for a given functional obtained by fitting to interaction energies for non-bonded complexes. The MBD energy is than calculated by diagonalizing the Hamiltonian Eq. (5.50) and using the plasmon-pole formula, Eq. (5.53).
BOP | B97-3 | B97 | BPW91 | revPBE | TPSS | BP86 | mPW91 | PBE | PBE0 | PW91 | SVWN | |
${s}_{R}$ | 0.42 | 0.45 | 0.50 | 0.53 | 0.54 | 0.56 | 0.64 | 0.68 | 0.83 | 0.85 | 0.91 | 1.19 |
The MBD-vdW approach greatly improves the accuracy of the interaction energies
for S66
^{
929
}
J. Chem. Theory Comput.
(2011),
7,
pp. 2427.
Link
test set, even if a simple functional like PBE is
used, with a mean absolute error of 0.3 kcal/mol and a maximum error of
1.3 kcal/mol, as compared to 2.3 kcal/mol (mean) and 7.2 kcal/mol (max) for
plain PBE. In general, the MBD-vdW method is superior to pairwise a
posteriori dispersion corrections.
^{
449
}
Chem. Rev.
(2017),
117,
pp. 4714.
Link
As mentioned above in the context of the TS-vdW method (Section 5.7.4), the FHI-aims or Quantum Espresso codes cannot perform exact unrestricted SCF calculations for the atoms and this leads to inconsistent free-atom volumes as compared to the spherical ones computed in Q-Chem, and thus inconsistent values for the vdW correction. Since the parameters of the TS-vdW and MBD-vdW models were fitted for use with FHI-aims and Quantum Espresso, results obtained using these codes are slightly closer to S66 benchmarks and thus scalar modifiers are available for the internally-computed Hirshfeld volume ratios in Q-Chem. For S66, the use of these modifiers leads to negligible differences between results obtained from all three codes.
The MBD-vdW correction is requested by setting MBDVDM = TRUE in the $rem section. Other job control variables, including the aforementioned modifiers for the free-atom volume ratios, are the same as those for the TS-vdW method and are described in Section 5.7.4.
MBDVDW
Flag to switch on the MBD-vdW method
TYPE:
INTEGER
DEFAULT:
0
OPTIONS:
0
Do not calculate MBD.
1
Calculate the MBD-vdW contribution to the energy.
2
Calculate the MBD-vdW contribution to the energy and the gradient.
RECOMMENDATION:
NONE
MBDVDW_BETA
Set custom value of the ${s}_{R}$ ($\beta $) damping parameter
TYPE:
INTEGER
DEFAULT:
no default value defined
OPTIONS:
$n$
Corresponding to $n\cdot {10}^{-4}$
RECOMMENDATION:
Use predefined values for supported functionals, otherwise consult
Ref. 38 and other relevant literature.
$molecule 0 1 O H 1 0.95 H 1 0.95 2 104.5 $end $rem BASIS 6-31G* METHOD PBE MBDVDW TRUE !vdw settings HIRSHFELD_CONV 6 ! sets SCF_CONVERGENCE for single atom calculations HIRSHMOD 4 ! Apply modifiers to the free-atom volumes for H, C, N, and O $end