本帖最后由 sige 于 2019-9-26 10:07 编辑
管理员注：此文有大量严重错误，勿轻信。斧正见2L
注：文章系转载http://blog.sciencenet.cn/blog-567091-836986.html 原文有大量严重错误，勿轻信。斧正见2楼sob老师回复。
Mulliken电荷是一种计算原子partial charge（局部电荷）的方法。本来，电子是在空间分布的一团电子云，但是广大玩计算的同志希望能有一种更简洁的方法来描述这种混沌不清的分布，看看谁多谁少，谁剥削了谁，因此计算原子电荷就非常流行，称为布局分析(population analysis)，通俗说，如何在原子间分割大家共同的家产。大致说先后提出了三种分割家产的思路。
第一类是以原子轨道对分子轨道的贡献划分。其中最著名的是Mulliken童鞋1955年提出的Mulliken population analuysis. Mullken认为，如果两个原子共用了一些电子对，那么电子可以分成3部分，属于A的(净原子布居)， 属于B的， 和共同的(重叠布居)。最公平的方法是，把大家共同拥有的一分为二，各拿一半，加上自己原有的，就是你的局部电荷。怎么计算净原子布居和重叠布居？根据分子轨道中的原子轨道系数。 因为分子轨道是原子轨道的线性组合，假设两原子核外各有一条电子轨道A, B组成了一个电子对分子轨道C=xA+yB，总电荷等于波函数的平方：|C|^2=2个电子；
仔细看这个公式，谁是谁的不用说了吧。。。类似依轨道划分的还有所谓NBO电荷，Born电荷。
第二种分割思路是按波函数节点划分，这就是Bader charge， 依据所谓的AIM理论，atom in molecule。Bader童鞋非常愤青，经常用很犀利的语言驳斥各种电荷划分。他说依轨道划分类似于在两原子间砍了一刀，各拿一半，这是没有任何物理意义的。那么什么叫有物理意义？他说你这一刀砍下去，要顺着波函数的节点砍下去。我们知道波函数的节点是波函数变号的地方，对应的电荷密度为0. 他说这是电子云空间分割最理想最天然的分割点，比你那轨道划分好5倍，呵呵。不可否认，分子轨道的空间节点可以显示成键的中心位置等信息，的确是一种比较好的划分方法。
第三种思路，都是搞力场的童鞋感兴趣的，叫ESP电荷，包括原装ESP, 改进版的ChelpG, Hershifiled等方法。搞力场的童鞋说，你们搞切割的，虽然思路清晰，但是我拿到我的力场中一用，计算静电作用，发现不行，连分子偶极距都是错的。肿么办呢？俺们也不搞切割了，俺们拿一个探测电荷放到分子外头探，探探受到了多大的静电作用；把空间均匀划分成几千个格点，探几千次；然后呢，我就在原子核位点尝试放置一些电荷，来拟合探测电荷受到的静电作用。不断尝试，最后最能够重现分子静电场的那一套原子电荷就是俺要的电荷。这就是原始的ESP电荷分析。后来发展的ChelpG什么，就是在拟合中加各种约束，比如增加偶极距的权重。这种电荷分析方法基本上是经验拟合，但是能够正确重现偶极距，分子间作用等性质，搞力场的最喜欢。对他们来说，只要结果正确即可，何须在意经由何种途径？唯一不爽的是屏蔽效应很严重，即包裹在分子里面的原子的电荷不容易探到它们的静电作用，拟合时你电荷给大点给小点，拟合结果都差不多，郁闷。。。 Bader分析之前需要完成的步骤和注意事项：第一步，Bader分析之前的静态计算参数：LAECHG；在静态计算的INCAR中，设置LAECHG=.TRUE.（同时注意：ICHARG=2，LCHARG =.TRUE.）.计算结束后，会产生AECCAR0，AECCAR1, AECCAR2这3个文件，以便后续分析使用。第二步，Bader分析需要的两个后处理工具：chgsum.pl和bader；（1）进入“VASP TSTtools”页面下载chgsum.pl文件，链接： http://theory.cm.utexas.edu/vasp/scripts.html （2）进入“BaderCharge Analysis”页面下载bader程序，链接 http://theory.cm.utexas.edu/henkelman/code/bader/推荐使用Linux版本，下载后用如下命令解压：tar-zxvf bader_lnx_64.tar.gz第三步，后处理工具的使用。（1）chgsum.pl的使用方法如下：运行：chgsum.pl AECCAR0 AECCAR2，得到CHGCAR_sum文件2）bader的使用方法如下：运行：bader -ref CHGCAR_sum CHGCAR，得到ACF.dat（包含价电子电荷信息）、BCF.dat、AVF.dat这3个文件。一些帖子中使用：bader CHGCAR -ref CHGCAR_sum，这样也是可以的。 注意事项： （1）Bader 只支持PAW 的方法，需要vasp4.6.31 以及以上版本。 （2）计算电荷布局分析时，需要比较密的FFT 网格，可以通过NG(X,Y,Z) F来调整，增大这些量，直到总的电荷数目正确。一般情况下，使用默认的FFT网格设置不会出现问题。 （3）Bader 电荷布居分析基于电荷密度。和Mulliken ,Hirhfeld 分析比较，原子电荷：Hirshfeld<Mulliken<Bader.（这个比较是定性的）。 （4）现在还不能分析某一个轨道的电荷布居情况，不过很快就可以了，新版本将具有这个功能（暂时没有验证新版本是否具有此功能，待后续更新补充）。 ACF.dat含义： 我一般只看ACF.dat
# X Y Z CHARGE MIN DIST ATOMIC VOL
--------------------------------------------------------------------------------
1 4.2634 1.4211 1.4211 8.1310 1.1369 9.7368
，1是第一个原子， x,y,z是原子的笛卡尔坐标。CHARGE是原子的bader电荷，这个是最主要的。。。然后是最小距离，原子体积。。。
NUMBER OF ELECTRONS: 64.0000 是你算的这个体系总共64个电子。
AVF.dat 我也不懂啥意思，不过，从名字上看，就是atom volume file
BCF.dat 就是bondcharge file。。你就知道意思了。。。 问：你好 我想问一下 我算的O是7.2左右 Ti是1.4左右，Ag 是10.3左右，这个是什么情况呀
POTCAR中 Os2p4 Ti d3s1 Ag s1d10 答：O 得到1.2e 约为-2价，Ti失去2.6e 约为+3/+4价，Ag失去0.7e 约为+1价，这只能作为趋势判断，定量效果并不好。
---------------------------------------------
VASP小白，潜心学习ing
| 
发表于 Post on 2019-9-24 09:08:38
收藏 Add to favorites
发表于 Post on 2019-9-25 02:52:21
楼主 Author