Use the PM6-DH2 procedure of: M. Korth, M. Pitonák, J. Rezác, and P. Hobza, "A Transferable H-bonding Correction For Semiempirical Quantum-Chemical Methods", J. Chem. Theory Comp. 2010, 6, 344–352 and J. Rezác, J. Fanfrlik, D. Salahub and P. Hobza, " Semiempirical Quantum Chemical PM6 Method Augmented by Dispersion and H-Bonding Correction Terms Reliably Describes Various Types of Noncovalent Complexes " J. Chem. Theory Comp. 2009, 5, 1749-1760. See: Abstract.
The PM6-DH2 method was parameterized to reproduce interaction energies for geometries obtained from high-level quantum mechanical calculations, see accuracy. While it is possible to use it for geometry optimization, please be aware of the following limitation and that the method should be used with extra care. Users are recommended to optimize the geometry with the PM6 or other PM6-DXX methods and calculate the final energy or interaction energy using PM6-DH2.
KNOWN LIMITATION OF PM6-DH2: The gradients obtained by the current implementation of the -H correction do not include the term containing the derivative of atomic charges with change of the coordinates. The structure obtained by minimization using this gradient is not the exact minimum of PM6-DH2 energy. This error is negligible in the case of weaker H-bonds, but in case of a strongly bound formic acid dimer, the error is 0.30 kcal/mol. This error can be avoided by using NOANCI, but then the calculations will then take much longer.
The PM6-DH2 procedure corrects binding errors in the PM6 method. It can be used with geometry optimization or with a single point (1SCF) calculation. Normally, two or three calculations would be needed to get the binding energy.
To print the individual components of DH2, add DISP