计算化学公社

标题: 求助：BSSE计算输出的电荷不正确 [打印本页]

作者
Author: SpongeBob 时间: 2021-3-30 09:46
标题: 求助：BSSE计算输出的电荷不正确
我再计算结合能过程中发现电荷不正确，请问是我的输入文件有错误吗？
这是输入文件
%chk=test-n1-mg.chk
%nprocshared=32
%mem=14GB
# opt b3lyp/6-311++g(d,p)  counterpoise=2

Title Card Required

3, 1, 1, 1, 2, 1
C(Fragment=1)    1.68038900 -1.51864900 -0.23093500
H(Fragment=1)    1.39502200 -1.66519900 -1.27534300
H(Fragment=1)    1.00658100 -2.09242200 0.40950600
C(Fragment=1)    3.11257000 -1.97638700 0.00444500
H(Fragment=1)    3.44552900 -1.83305100 1.03318700
H(Fragment=1)    3.13394100 -3.05292200 -0.18570700
H(Fragment=1)    3.82818000 -1.52320000 -0.68052800
C(Fragment=1)    2.32290700 0.85729400 -0.79433000
H(Fragment=1)    3.32419200 0.57960900 -0.47326800
H(Fragment=1)    2.18424700 0.50764500 -1.81919500
C(Fragment=1)    2.12763500 2.35969100 -0.67959700
H(Fragment=1)    2.90549100 2.82209400 -1.29379400
H(Fragment=1)    1.17562500 2.71593600 -1.07873600
H(Fragment=1)    2.27124700 2.73900400 0.33333600
C(Fragment=1)    1.61790400 0.26335700 1.53998900
H(Fragment=1)    1.24582800 1.25709800 1.77999200
H(Fragment=1)    1.12878600 -0.49187800 2.15375200
H(Fragment=1)    2.68743100 0.22926400 1.73212300
C(Fragment=1)    -0.05166800 0.29476900 -0.32481100
H(Fragment=1)    -0.20274800 1.34428500 -0.07804000
H(Fragment=1)    -0.10919200 0.18630900 -1.40926400
C(Fragment=1)    -1.12999900 -0.56875400 0.35344100
H(Fragment=1)    -1.15805000 -1.58607000 -0.03372300
H(Fragment=1)    -1.05204100 -0.59523200 1.44026800
O(Fragment=1)    -2.51388600 -0.01467900 0.07013400
N(Fragment=1)    1.37896700 -0.02909700 0.07403800
C(Fragment=1)    -2.81788300 1.28548800 0.77127300
H(Fragment=1)    -3.89611000 1.46364600 0.68109900
H(Fragment=1)    -2.55457800 1.17949300 1.82292900
H(Fragment=1)    -2.27970800 2.09117100 0.27864700
Mg(Fragment=2) -4.09937300 -0.82794400 -0.68778500

最终log文件中电荷却不正确
ESP charges:
            1
   1  C 0.000000
   2  H 0.000000
   3  H 0.000000
   4  C 0.000000
   5  H 0.000000
   6  H 0.000000
   7  H 0.000000
   8  C 0.000000
   9  H 0.000000
10  H 0.000000
11  C 0.000000
12  H 0.000000
13  H 0.000000
14  H 0.000000
15  C 0.000000
16  H 0.000000
17  H 0.000000
18  H 0.000000
19  C 0.000000
20  H 0.000000
21  H 0.000000
22  C 0.000000
23  H 0.000000
24  H 0.000000
25  O 0.000000
26  N 0.000000
27  C 0.000000
28  H 0.000000
29  H 0.000000
30  H 0.000000
31  Mg 2.000000
Sum of ESP charges = 2.00000
ESP charges with hydrogens summed into heavy atoms:
            1
   1  C 0.000000
   4  C 0.000000
   8  C 0.000000
11  C 0.000000
15  C 0.000000
19  C 0.000000
22  C 0.000000
25  O 0.000000
26  N 0.000000
27  C 0.000000
31  Mg 2.000000
      ESP charges condensed to fragments:
   1    0.000000
   2    2.000000

作者
Author: snljty 时间: 2021-3-30 10:26
本帖最后由 snljty 于 2021-3-30 10:27 编辑

不要带着CounterPoise做优化和频率计算，没有解析导数，要做有限差分，计算量大的多的多还可能有数值问题。电荷这个问题应该是由于缺乏相应的解析导数，Gaussian不给算。看卢老师博文http://sobereva.com/46 《谈谈BSSE校正与Gaussian对它的处理》

作者
Author: SpongeBob 时间: 2021-3-30 12:11

snljty 发表于 2021-3-30 10:26
不要带着CounterPoise做优化和频率计算，没有解析导数，要做有限差分，计算量大的多的多还可能有数值问题。 ...

好的，谢谢

作者
Author: sobereva 时间: 2021-3-31 13:12
没事别给氢加弥散函数，仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

B3LYP都不带色散校正还指望算结合能？
谈谈“计算时是否需要加DFT-D3色散校正？”
http://sobereva.com/413（http://bbs.keinsci.com/thread-9772-1-1.html）

在counterpoise任务中看原子电荷也是完全莫名其妙。counterpoise任务仅适合用于算结合能、BSSE校正量，不要同时干别的事、读取任何其它信息

作者
Author: SpongeBob 时间: 2021-3-31 14:42

sobereva 发表于 2021-3-31 13:12
没事别给氢加弥散函数，仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

我看到老师的博文已经进行了修改，谢谢老师

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