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
J. Chem. Theory Comput.
(2018), 14, pp. 2955. especially those that contain a large number of polarizable centers. 284 Int. J. Quantum Chem.
(2014), 114, pp. 1157. The many-body dispersion (MBD) method of Tkatchenko et al. 1123 Phys. Rev. Lett.
(2012), 108, pp. 236402. , 38 J. Chem. Phys.
(2014), 140, pp. 18A508. represents a more general and less empirical approach that goes beyond the pairwise-additive treatment of dispersion. This is accomplished by including -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. 112 Chem. Sci.
(2016), 7, pp. 1712.
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
(2017), 117, pp. 4714.
Here, is the mass-weight displacement of oscillator from its center , is the characteristic frequency, and is the static polarizability. is the dipole potential between the oscillators and . The MBD Hamiltonian is obtained through coarse-graining of the long-range correlation (through the long-range dipole tensor ) and approximating the short-range polarizability via the adiabatic connection fluctuation-dissipation formula:
This approximation expresses the dynamic polarizability (of a given fragment) in terms of the polarizability of the corresponding QHO,
in which is the charge, the mass, and 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,
in which is the number of fragments, are the frequencies of the interacting (dipole-coupled) system, and are the frequencies of the non-interacting system (i.e., the collection of independent QHOs). The sum runs over all characteristic frequencies of the system.
A particular method within the MBD framework is defined by the models for the
static polarizability (), the non-interacting
characteristic frequencies (), and the damping function  used
to define . In Q-Chem, the MBD method is implemented
following the “MBD@rsSCS” approach, where “rsSCS” stands for
range-separated self-consistent screening.
J. Chem. Phys.
(2014), 140, pp. 18A508. In this approach, 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.5.
The short-range atomic polarizabilities 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 .
The short-range atomic polarizabilities are summed up along one fragment coordinate to obtain the local effective dynamic polarizability, i.e., , and are then spherically averaged. The range-separation (damping) function used to construct and is the same as that in Eq. (5.45), except with instead of , and again 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).
The MBD-vdW approach greatly improves the accuracy of the interaction energies
J. Chem. Theory Comput.
(2011), 7, pp. 2427. 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. 463 Chem. Rev.
(2017), 117, pp. 4714.
As mentioned above in the context of the TS-vdW method (Section 5.7.5), 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.5.
$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