计算化学公社

标题: 请教计算小分子RESP电荷时，为什么高斯计算已收敛，但是结果显示为星号 [打印本页]

作者
Author: 784539690 时间: 2026-3-31 19:19
标题: 请教计算小分子RESP电荷时，为什么高斯计算已收敛，但是结果显示为星号
本帖最后由 784539690 于 2026-3-31 19:37 编辑

各位老师好，我在使用高斯与Multiwfn计算丙二酰辅酶A的RESP电荷时，
先使用# B3LYP/6-311G** em=GD3BJ opt freq，在linux系统（g16）上计算，显示不收敛报错。
后使用# B3LYP/6-311G** em=GD3BJ opt=GDIIS freq int=ultrafine SCF=XQC，计算12小时后，结束，但是freq不收敛。
Maximum Force          0.000166    0.000450    YES
RMS    Force          0.000014    0.000300    YES
Maximum Displacement    0.025960    0.001800    NO
RMS    Displacement    0.003697    0.001200    NO

提取chk文件的最后一帧结构，使用# B3LYP/6-311G** em=GD3BJ opt=(GDIIS,tight) freq int=ultrafine SCF=QC计算，结果基本收敛
Maximum Force          0.000001    0.000015    YES
RMS    Force          0.000000    0.000010    YES
Maximum Displacement    0.000125    0.000060    NO
RMS    Displacement    0.000027    0.000040    YES

，随后在linux系统上使用formchk result.chk命令将chk转换成fch，chgcons.txt 设置为1-87 -5
在multiwfn程序中载入.fch，然后输入
7  //布居分析与原子电荷计算
18  //RESP模块
6  //设置做拟合静电势计算时用的原子电荷约束
1  //从外部文件中读取电荷约束
chgcons.txt  //电荷约束文件的实际路径
2  //做一步方式的拟合静电势计算
计算约900秒后，结果都带星号。请问原因是什么。

本人这方面很不专业，上面有什么信息写错了，还望老师批评指正。
文件上传到网盘了，[color=rgba(6, 10, 38, 0.9)][backcolor=rgba(255, 255, 255, 0.5)]链接：https://pan.quark.cn/s/a3a517119736?pwd=6Ac5
提取码：6Ac5

作者
Author: sobereva 时间: 2026-4-1 05:27
我不方便从这个网盘上下载

如果不做约束能正常用Multiwfn得到RESP电荷，只可能是约束设得不当
只有87个原子，干嘛约束为-5？Gaussian计算时候净电荷设成-5，Multiwfn直接算RESP电荷就完了，根本不需要任何约束

当前是去质子化态，必须带溶剂模型

没事甭用危险的SCF=QC，下文明确说了
解决SCF不收敛问题的方法
http://sobereva.com/61

没必要用tight收敛限，算RESP电荷对几何优化收敛精度没什么要求，loose收敛限都可以

作者
Author: 牧生 时间: 2026-4-1 07:55
本帖最后由牧生于 2026-4-1 08:05 编辑

你这个结构式本身就画的不对，原子间的距离有显著的不合理之处，后面的优化那一系列结果全都不能用
你打开result.gjf，应该一眼就能看到不合理所在

体系的净电荷也不应该是-5

欢迎光临计算化学公社 (http://bbs.keinsci.com/)