NPA电荷不同版本高斯运行报错

lonelycvavc

各位老师和前辈们好。计算NPA电荷时我遇到一个问题，在A.02版和D.01版（均为windows系统下）提交附件中的gjf文件会报错：Bad arguments to LAPack or BLAS routine. log文件是同一个输入文件提交到linux系统下的E.01版高斯中的运行结果，虽然计算正常终止但开始NPA计算之前有一行 Computed number of electrons  139.99822915 does not match expected value of   140. 请问这样计算出来的结果正确吗？以及报错的原因是高斯版本或系统版本的原因吗？有没有解决办法呢？谢谢！

我本是个娃娃

用不着这么大的基组，把基组降到6-31G(d)就可以了

helpme

又是win又是linux，又是A又是E，这样一下就把问题搞复杂了。
可以固定一个参数试试结果，比如linux下试试A和E的，看到底什么问题？

我觉得linux的结果应该没有什么问题。

另外 Linux E01版本的高斯都已经用上了。

lonelycvavc

谢谢娃娃，刚才试了用您给的基组计算没有问题，也没有出现Computed number of electrons  139.99822915 does not match expected value of   140.的提示。但是我分子是用大基组优化的，不知NPA电荷改用小基组计算是否严谨。
也谢谢helpme老师，Linux系统下的输出文件是我拜托一位老师帮我算的，E.01版的高斯也是装在那位老师自己的服务器里，我手里只有D.01版（A.02版因之前怀疑是版本问题已卸），所以您说的试验有点难以实现。
我更想知道出现这个问题的原因和有没有能让我在自己电脑（windows系统Gaussian D.01版）中正常计算得出结果而不用更改基组的方法。再次谢谢各位。

sobereva

bad argument那个报错是程序bug。
139.99822915偏离140属于可忽略的数值误差，可以无视。
既然去掉弥散函数可以算，你就在6-31G*下面做就行了。弥散函数对于NBO分析影响甚微，完全可忽略。

lonelycvavc

好的，谢谢Sob老师！