计算化学公社

标题: 关于resp计算小分子电荷的误差大小的问题 [打印本页]

作者
Author: yyfdemajia 时间: 2019-9-4 22:32
标题: 关于resp计算小分子电荷的误差大小的问题
各位大神好，我最近在重复一篇论文用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

作者
Author: sobereva 时间: 2019-9-5 01:09
各种数值细节设定都可能影响结果，包括计算级别、构象、拟合格点分布和密度等等

反正你就按照下文的做法算，结果肯定可靠
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）

作者
Author: yyfdemajia 时间: 2019-9-5 02:21
好的多谢

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