对AIM electron density的 basin做Laplacian积分，结果不为零

Charles_Gao

请教Sob老师：
我想求BCP点对各个Basin的Source function积分，得到各原子对该BCP的电子密度的贡献。文献中说Laplacian积分应当接近零，以证明我的积分数值是正确的。然而我的Laplacian积分的值比较大，可以和该BCP的电子密度相比了。提高grid的密度并没有改善。请问这是什么原因呢？
PS，我看到积分的时候似乎忽略电子密度小于0.001的空间了，请问这会有影响吗？

Charles_Gao

我的输出结果：
Total result:
     Atom       Basin       Integral(a.u.)   Vol(Bohr^3)   Vol(rho>0.001)
      1 (C )       6          0.00269775        39.507        35.776
      2 (F )       8         -0.00117154       611.648       102.938
      3 (F )       3         -0.00063455       670.528       102.864
      4 (F )      17         -0.00100000       562.031       105.629
      5 (C )       7          0.00186629        39.199        35.675
      6 (F )       9         -0.00121250       615.782       102.639
      7 (F )       1          0.00003214       674.040       102.849
      8 (F )      15         -0.00104788       583.927       105.443
      9 (C )       4          0.00093024        39.192        35.748
     10 (F )      14         -0.00005318       645.543       102.693
     11 (F )      16         -0.00097808       613.460       102.921
     12 (F )       2         -0.00040709       600.274       105.566
     13 (C )      10         -0.00308807        39.367        35.669
     14 (F )      12         -0.00088592       659.451       102.918
     15 (F )      11         -0.00075808       569.770       105.556
     16 (F )      13         -0.00123067       630.475       102.720
     17 (Cu)       5          0.07449841       191.016        99.024
Sum of above integrals:             0.06755726
Sum of basin volumes (rho>0.001):    1486.629 Bohr^3

这个BCP的电子密度为0.1079

sobereva

2L的数据是直接用混合格点对电子密度拉普拉斯函数进行的积分？这种积分精度是很正常的，毕竟拉普拉斯函数是分布特征波动非常剧烈的函数，积分误差会大于诸如电子密度这样分布较平缓的函数。

你看看利用源函数算出的BCP的密度和直接算出来的BCP的密度相差多少，如果相差不大，结果完全可以接受。

并没有忽略电子密度小于0.001的空间，只不过输出体积的时候把密度>0.001的区域记入原子体积

Charles_Gao

sobereva 发表于 2018-8-17 22:13
2L的数据是直接用混合格点对电子密度拉普拉斯函数进行的积分？这种积分精度是很正常的，毕竟拉普拉斯函数是 ...

谢谢！source function积分结果约为0.1056，似乎是已经比较接近了