有个150原子左右的有机分子,之前用def2tzvp基组转化fchk,都是正常结束。
审稿人说基组有些小,于是对def2tzvp的优化得到的结构,将中心的几个关键原子的基组提高到def2qzvpp再计算了次单点能。
看了log也是正常结束,但是输出的Mulliken电荷已经很大不同。
在同一服务器上转化fchk的时候,报了个Warning
Warning: fchk array ILSW element 12 has value 1649267441664 set to -1.
机子里就装了一个gaussian,版本是ES64L-G16RevC.01,应该不是其他帖子提到的混用不同版本gaussian导致的报错。
gview或者multiwfn还是可以正常打开和查看分子轨道的,但关键结论和之前已经不同了。比如HOMO,LUMO的形状已经和def2tzvp算出来很大差异了。
请问这种情况,是正常的吗?还是fchk转化有问题呢?
作者Author: zjxitcc 时间: 1 hour ago
g16A/B/C也不能混用formchk和unfchk,最好再检查一下版本,比如在提交任务脚本里写which formchk,与当前终端运行which formchk,比较一下是否结果有所不同。