求助：BSSE计算输出的电荷不正确

SpongeBob

我再计算结合能过程中发现电荷不正确，请问是我的输入文件有错误吗？
这是输入文件
%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

snljty

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

SpongeBob

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

好的，谢谢

sobereva

没事别给氢加弥散函数，仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

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

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

SpongeBob

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

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