计算化学公社

标题: 对AIM electron density的 basin做Laplacian积分，结果不为零 [打印本页]

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

作者
Author: Charles_Gao 时间: 2018-8-17 21:41
我的输出结果：
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

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

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

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

作者
Author: Charles_Gao 时间: 2018-8-20 16:13

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

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

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