AM1 parameter bug
Letif Mones
letif... at gmail.com
Mon Apr 22 13:27:02 UTC 2013
Dear CP2K users,
I found a bug related to an AM1 parameter of the phosphorus atom. It is in file semi_empirical_parameters.F line 1258:
sep%fn3(2) = 0.03329000_dp
However, the corresponding correct parameter should be:
sep%fn3(2) = 1.77935400_dp
With this correction the energy and forces of a test system (some phosphate) are consistent with those obtained from other programs (sander/gaussian/mopac).
Kind regards,
Letif
More information about the CP2K-user
mailing list