请教计算小分子RESP电荷时，为什么高斯计算已收敛，但是结果显示为星号

784539690

各位老师好，我在使用高斯与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

sobereva

我不方便从这个网盘上下载

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

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


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

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

牧生

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

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