大基组转化fchk报Warning

Dave666

有个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转化有问题呢？

zjxitcc

g16A/B/C也不能混用formchk和unfchk，最好再检查一下版本，比如在提交任务脚本里写which formchk，与当前终端运行which formchk，比较一下是否结果有所不同。

Mulliken电荷本来就对基组比较敏感，变化明显也算是正常的，无法用来说明fch文件转化问题

rugals

感觉像是chk本身有问题，
条件允许的话最好再算一遍单点对比下。
另外基组都到qz了，泛函精度是不是也得跟上。

sobereva

下文明确说过Mulliken电荷对基组敏感度极高，不同基组下没有可比性
一篇深入浅出、完整全面介绍原子电荷的综述文章已发表！
http://sobereva.com/714（http://bbs.keinsci.com/thread-46067-1-1.html）
《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）

不同基组下轨道顺序往往有所不同，导致HOMO、LUMO显著变化。并且QZ基组比TZ弥散程度明显更高，这一点上对空轨道的特征影响也往往很明显