PM6-DH2 - Accuracy

PM6-DH2 gives much improved interaction energies compared to PM6.  Whereas in PM6, the average unsigned error for interaction energies for the complexes in the S22 set is 3.27 Kcal/mol, for PM6-DH2, the equivalent average unsigned error is 0.37 Kcal/mol. For details of PM6-DH2, see: M. Korth, M. Pitonak, J. Rezac, and P. Hobza, "A Transferable H-bonding Correction For Semiempirical Quantum-Chemical Methods",  J. Chem. Theory Comp. 2010, 6, 344-352

Some representative complexes are shown in the Table.

Table of Representative Complexes

Complex  CCSD(T)/CBS  PM6  Error PM6-DH2 Error 
A-T WC  -16.37 -9.06 7.31 -16.45 -0.08
G-C WC  -31.40 -21.35 10.05 -30.75 +0.65
A-T Stack  -12.23 -4.94 7.39 -11.67 +0.54
G-C Stack  -19.02 -13.07 5.93 -19.66 -0.64
Water dimer  -5.02 -3.94 1.08 -4.90 +0.12
Formic acid dimer  -18.61 -11.14 7.47 -18.65 -0.04
Formamide dimer  -15.96 -12.55 3.41 -15.87 +0.10
Lys+ - Glu-  -108.40 -98.83 9.57 -104.63 +3.77
Arg+ - Asp-  -110.80 -100.83 9.97 -109.80 +1.01

Calculations were done using the high-quality geometries from the S22 database, see: Jurecka, P.; Sponer, J.; Cerny, J.; Hobza, P. Phys. Chem. Chem. Phys. 2006, 8 (17), 1985-1993
Reference interaction energies (CCSD(T) from the  [CCSD(T)/complete basis set (CBS) limit] of the dispersion-bonded complexes.

Worked example:

A. PM6-DH2 Formic acid dimer ΔHf = -188.158
B. PM6-DH2 two separated formic acid molecules
ΔHf = -169.504

Interaction energy = B - A = -18.654

Note: Heats of formation predicted by PM6-DH2 should not be compared to reference heats of formation.  This is because PM6-DH2 is designed for predicting energies of interaction.  If heats of formation are wanted, use PM6, for which the predicted  ΔHcan be compared to experimental values.

See also PM6-DH+ accuracy.