计算化学公社

标题: 大基组转化fchk报Warning [打印本页]

作者
Author:
Dave666    时间: 3 hour ago
标题: 大基组转化fchk报Warning
本帖最后由 Dave666 于 2026-2-26 14:36 编辑

有个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,比较一下是否结果有所不同。

Mulliken电荷本来就对基组比较敏感,变化明显也算是正常的,无法用来说明fch文件转化问题
作者
Author:
rugals    时间: halfhour ago
感觉像是chk本身有问题,
条件允许的话最好再算一遍单点对比下。
另外基组都到qz了,泛函精度是不是也得跟上。

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

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3