Hi,
I have spent most of the afternoon improving the debug output for the
GTO calculation code, and checking into an issue Michael spotted (thanks
for adding Molpro support too). We have three files that should
essentially give the same output now - Gaussian, Q-Chem and Molpro. The
only problem is that they do not agree!
http://avogadro.openmolecules.net/wiki/Basis_sets
As you can see, all of the exponents match (a). The contraction
coefficients © for Gaussian and Molpro are identical, but the Q-Chem
contraction coefficients differ quite significantly. The shape and size
of the isosurfaces produced using Gaussian and Molpro both look good,
but I still think that there is a problem with the scaling of the
Gaussian values that I need to look into more.
This is more a status update on the GTO code. I would like to produce
some reference cube files for each supported file type ideally, and then
add some unit tests. It seems like methane HF/6-31G is a reasonable test
molecule to use for this purpose. A cube rms deviation type value might
be the most appropriate method of comparison. We should expect some
variation between different packages due to convergence criteria etc,
but not that much with such a simple molecule/basis set as far as I can
tell.
Thoughts?
Thanks,
Marcus