关于resp计算小分子电荷的误差大小的问题

yyfdemajia

各位大神好， 我最近在重复一篇论文用resp计算Acetophenone小分子partial charge的结果。 我用了和论文一样的泛函和基组，发现结果和论文发表的数值大致相同，但是最大误差却差了0.08。

比如，同样的碳原子我算了 -0.058，而论文的结果确实是-0.0085，这个原子就差了0.05。我想问一下这么大的误差合理吗？ 还是我自己哪里算错了？  


我用的高斯输入文件是：

%chk=./g09.chk
%mem=2GB
%NProcShared=20
# opt=tight b3lyp/6-311++g** nosymm scf=qc pop=mk iop(6/33=2,6/42=6,6/50=1)

Acetophenone

0 1
C                 -0.70411985   -0.21325049    0.22444779
C                  0.69485405   -0.23835167    0.30295828
C                  1.42777349    0.94327388    0.12828111
C                  0.76171882    2.15000129   -0.12490283
C                 -0.63725504    2.17510232   -0.20341414
C                 -1.37017433    0.99347626   -0.02873978
H                 -1.26372012   -1.11544780    0.35781781
H                  1.20340142   -1.15971431    0.49627272
H                  1.32131899    3.05219893   -0.25827098
H                 -1.14580240    3.09646494   -0.39672866
H                 -2.43832202    1.01264178   -0.08868311
C                  3.69703557   -0.41038360    0.49278220
H                  3.91897492   -0.89657754   -0.43417970
H                  4.60766733   -0.21170281    1.01830645
H                  3.07301187   -1.04423322    1.08753347
C                  2.96510752    0.91568994    0.21455513
O                  3.62323926    1.97674144    0.05770207

AP_ini.gesp

AP.gesp

sobereva

各种数值细节设定都可能影响结果，包括计算级别、构象、拟合格点分布和密度等等

反正你就按照下文的做法算，结果肯定可靠
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）
计算RESP原子电荷的超级懒人脚本（一行命令就算出结果）
http://sobereva.com/476（http://bbs.keinsci.com/thread-12858-1-1.html）

yyfdemajia

好的 多谢