Bader电荷分析时的格子划分

patience2014

在做第一原理计算后，需要对电子密度进行Bader电荷分析，但问题是如何划分格子能让分析后的电子总数与实际电子总数相等或很接近（0.1e以内）？比如说体系有1361个电子，经过Bader分析后的电子数变成了1361.78388个电子，这个时候电子总数就比实际电子数多了0.78388个电子了，导致原子所带电荷的误差比较大。很无奈，一直都是试验很多个格子，比如136*136*136，如果误差大，就换140*140*140，直至误差最小。请问有没有其他的办法可以准确快速地进行格子划分呢？

sobereva

这要看你的体系，以及用的什么程序做得量化计算。
如果你是指bader那个程序，那个程序用的是均匀立方网格，无论你把格点间距设得多小，对于全电子计算，由于内核附近电子密度非常大，积分误差也会很大。
只有Multiwfn这种原子中心+均匀格点混合方式做AIM分析来得到AIM电荷（也叫bader电荷）的程序才能准确积分全电子计算得到的电子密度。
如果你是用赝势产生的波函数，由于没有内核电子密度很大的区域，单纯用均匀立方网格积分还算说得过去，格点间距越小（且盒子距离边缘原子足够大）的时候可以得到越准的积分精度。

patience2014

感谢您的回复！
是在Cu表面放上有机分子这样的体系，使用SIESTA程序计算的。是使用的赝势波函数。是用bader那个程序做电荷分析的。Multiwfn能分析SIESTA产生的TOCH文件吗？SIESTA可以产生cube文件，产生的时候就需要输入格点数目了。
您说格点间距越小（且盒子距离边缘原子足够大）积分越准，整体趋势是这样的。但具体操作的时候，会有波动，摸不出规律，还是得一个个尝试。
就想知道有没有别的途径可以更好地解决这个问题。

sobereva

patience2014 发表于 2016-9-14 14:08
感谢您的回复！
是在Cu表面放上有机分子这样的体系，使用SIESTA程序计算的。是使用的赝势波函数。是用bade ...

Multiwfn分析不了siesta的结果。
Multiwfn也可以基于电子密度的cube文件做AIM分析，不过此时没法用原子中心格点，所以就和bader就没什么区别了

没什么更好途径，只能增加格点密度。
盒子距离边缘原子不用太远，只要不太近就行。