转载自  http://blog.sciencenet.cn/blog-567091-736154.html  仅供学习交流  侵权必删

处理步骤：【需要商议的步骤】---另一种定义电荷密度差

http://muchong.com/html/201512/9799282.html  【是这样做的吗？ ask Ruizhi】

用VESTA 做VASP差分电荷（C-A-B）图方法：
1、导入 C 的CHGCAR（拖拽到vesta）面板框中
2、Edit——>Edit data——>Volumetric Data————>Import (选择A 的CHGCAR)——>Subtract
from current data——>OK

3、Edit——>Edit data——>Volumetric Data————>Import (选择B 的CHGCAR)——>Subtract
from current data——>OK
注意：A,B,C三个文件的格点要一一对应。

文献分析：

(1)Cr合金化Mg_2Ni氢化物能量与电子结构的第一性原理研究

      为了更直观更形象地理解Mg2NiH4和Mg(I)CrNiH4中原子间的成键情况，选取2 个晶胞(002)面
的差分电荷密度图(是指原子---组成体系(团簇)之后电荷的重新分布， 可以很清楚地看出体系中原子的成
键情况)如图4 所示。

     发现在Mg2NiH4晶胞中，H-Ni间的电荷密度高，而H-Mg 及Ni-Mg 间的电荷密度低，故H-Ni 之间存在较强的相互作用，而H-Mg 及Ni-Mg 间的作用不明显；Cr 合金化Mg2Ni 氢化物(见图4b)后，发现H-Ni 间的电荷密度降低，这说明H-Ni 键的作用减弱, 故有利于释氢， 体系的解氢能力增强。

http://image.sciencenet.cn/album/201608/18/110838ddxwii5dsq5a8sld.png

（２）　放热型金属合金化对钒基贮氢材料性能影响的理论研究

　　在氢化物 VH₂ 中加入放热型金属 Mg、 Ca、 Ti、 Zr 以后,V₅₁M₁₂H₆₄(M=Mg、 Ca、 Ti、 Zr)中最 中心的 V和 H 之间电荷密度重叠性增强, 表明 V-H 之间的相互作用增加; 同时 M(M=Mg、 Ca、 Ti、 Zr)和 H 之间的电荷密度重叠性较大, 表明 M-H(M=Mg、 Ca、 Ti、 Zr)之间也有较强的相互作用。所以, 放热型 金属合金化后, 由于放热型金属本身 Mg、 Ca、 Ti、 Zr 可以吸氢, 同时又使 V-H 之间的相互作用增强, 氢 化物可以结合更多的氢, 从而增加了贮氢量。但氢 化物不容易放出氢, 即吸放氢性能减弱。这与实验 结果及前面的结论一致。

（３）金属氢化物力学性能的第一性原理研究

　    金属氢化物 的稳定性与Fermi能级处电子态密度的分布密切相关。 不稳定的金属氢化物如TiH₂和HfH₂在Ferm能级处一般有电子态密度尖峰出现。

       如果电子被束缚在H原子的局域势场附近，则体系的稳定性变差(可由差分电荷密度图判断)。

       当金属氢化物体系总的电子态密度在Fermi能级处出现奇异性的尖峰，或者电子被束缚在H原子的局域势场附近时，体系的稳定性差， 有可能发生相变。

     在低温下MH2(M=La,Nd,Gd,Tb,Y,Dy,Ho,Er,Lu,Sc)具有稳定的氟化钙结构，而MH2(M=Ti,Zr,Hf)的氟化钙结构不稳定。对于稳定的金属氢化物而言， 抵抗体积形变、切应变和拉伸(或压缩)形变的能力由弱递增的顺序为LaH₂，NdH₂,GdH₂,TbH₂,YH₂,DyH₂,HoH₂,ErH₂,LuH₂,ScH₂。

http://image.sciencenet.cn/album/201608/18/114300bqzkf9zfkq4iqk9k.png

http://image.sciencenet.cn/album/201608/18/115701yodetz06osydjetc.png

（4） Hydride-Based Electride Material, LnH2 (Ln = La, Ce, or Y)

   

Figure 3a shows the electron density map of LaH₂ in the energy region from −1.5 eV to EF.

   Surprisingly, conduction electrons populate the cage site and show strong anisotropic shape into the six neighboring La ions.

        As seen in Figure 3b, the (110) plane obviously contains a large contribution from La 5dz2 (eg, red part). Thus, we may regard LaH₂V as an inorganic electride in which the electrons reside in the center of the Ht8 cube. However, the electrides reported to date are usually metal-rich compounds having cages surrounded by cations.2,3 

   Thus, we need to check the calculated electron distribution carefully by estimation of the electron localization function (ELF) map in order to get information about unshared electrons.

  Figure 3c shows the (110) plane in the ELF map. It shows clearly the ionic bonding in the compound and that there are no lobes【叶片】 in the cage originated from lone pair (LP) electrons. Td site symmetry of the anion【阴离子，H- ions】 site prevents the appearance of the LP state. 【为什么】

     http://image.sciencenet.cn/album/201608/20/1132449eey9r9ad0a8z909.png

（5） 差分电荷密度图--透明钠：Transparent dense sodium

   http://image.sciencenet.cn/album/201608/20/120056pdm1muqkqjkxmki1.png

 

Difference charge density (e Å^-3）; crystal density minus superposition of isolated atomic densities) plotted in the (110) plane at 320GPa (6.50Å^-3 per atom)

Thus, the insulating state in Na-hP4 arises from the strong localization of valence electrons in the interstices of Na-hP4 (Fig. 4c and Supplementary Fig. 6) enabled by p–d hybridization.

（5） Pressure-stabilized superconductive yttrium hydrides

http://image.sciencenet.cn/album/201608/20/124155ivrr98arrv0vapks.png

 

Valence electrons localization was found between the two neighbouring H2 atoms while absent between H1 and H2 atoms (Fig. 2e). This indicates the presence of both molecular “H2” and monoatomic H in tI10-YH4.

Despite the longer distance, covalent interaction between H atoms is clearly visible from the localized valence electrons between the H atoms (Fig. 2f).

http://image.sciencenet.cn/album/201608/20/125606klwz945swetbeb9w.png

  The empirical consideration is support from quantitative calculations of the difference of the electron density of tI10-YH4 to that of a hypothetical structure consisting only H sublattice.

  It is clearly shown in Fig. 3(a,b), that there is no electron density (covalent bond) between the two H₂ atoms in the pure H structure. However, when the Y atoms were present, localized electrons are found between two H2 atoms. Therefore, the charge transfer  from Y to H₂ is responsible to the formation of molecular “H₂” in YH₄.

    Similarly, the foramtion of “H₆” hexagons in YH₆ results from the accommodation by H of excess electrons from Y atom (Fig. 3c,d).

   For comparison, a survey of molecular “H₂” in most hydrogen-rich compounds only show a slightly elongated H-H bond length than pure solid H₂. Examples are 0.87Å in GeH₄ , 0.79Å  in SnH₄ , 0.84Å  in SiH₂(H₂)₂ , 0.76 Å  in LiH₂ and 0.8 Å  in NaH₉ .

http://image.sciencenet.cn/album/201608/20/13354186tkkzcmkmmxxr11.png

网络问答：

(1) AA:: 如何在vaspview中画电荷密度等高线，图是我画出的电荷密度，想在上面添加等高线，求大能帮助
 QQ:::

vesta:
Utilities->2D Data Display->Coutour

(2) AA::

怎么调节切片方向呢？ 我切出来是这样的，调节ranges of fractional coordinates中参数都不管用，是不是调节其他参数啊

QQ:::

Utilities->2D Data Display->Coutours->Slice->hkl
然后还有distance from origin可以调 [距离原点的距离？]

QQ:::

用lev00来分析也很好。。。比较方便。

对的，你安装好lev00后，里面会有很多选项，包括可以任意指定切面。。
然后再用origin画。

(3) vaspview和vesta比较

比较了一下，感觉还是vesta出来的图美观和操作方便些。vaspview里面的功能太少了。vesta切slice也很方便，不管z轴还是x轴还是任意方面的面都可以切。vesta中你的线消失是因为你切的面没落在原子上当然等高线就出不来了。

       切时先用shift选中三个原子，切面就容易些。貌似切出的面只能存成图片格式，export 2D data不知为什么是灰色，懂vesta的讲下

QQ:::可切任意面，注意切面到原点的距离，就那么几项，仔细看看，你一定能搞定

终于搞定了，原因出在CHG与CHGCAR格式的差异上，原来做差时用的是CHGCAR，差值也按CHGCAR的格式保存的，但是VESTA认的是CHG的格式，用CHGCAR的格式就出现了上述情况，现在按CHG格式保存电荷密度差，结果就正常了。
(5)  与MS的比较

   其实没有 ms 方便 ，但是我经常在 linux 下，所以我感觉  VESTA+isotropy 是目前最好的选择了

     建议调整模型就不要用VESTA了， 周期性结构建模也不要用gaussview，用MS。
第一个问题就可以在MS中把晶胞调整参数后rebuilt一下就可以了；第二个问题建模用MS的话就更方便了，建晶胞--切表面--做真空层--建超胞

QQ:::   boundary建立模型，然后再POSCAR中将相应原子删除

（6） 能不能告诉我具体第几章第几节谈到了为什么不能导出2D电荷密度数据的问题？
我看了Manual上面13章仅仅一笔带过说这个选项是导出数据。

QQQ:::"15.7 Exporting 2D data
To export 2D data of the specied slice as a text le, choose File" menu ) Export 2D Data. . . "
in the 2D Data Display window. No 2D data les can be output on selection of style (hkl) plane
in the bounding box" in the Slice Properties dialog box."

（7） 我现在想用脚本来处理态密度，之前用过论坛上的dos-procar.f，但是都是处理s,p,和 d的，这个脚本没有针对s,py,pz,px,dxy,dyz,dz2,dxz,dx2的，所以请大家帮帮忙，我知道送金币很俗，但也只有金币能聊表谢意了

(8)   VESTA中怎么能调出等高线图的标尺呢，就是右边的标尺怎么获得呢？谢谢

【如何把Slice图放到结构当中，像下面这样？】
http://blog.sciencenet.cn/static/js/grey.gif

QQQ:::可以呀 在2D视图里面，可以设置等高线和色彩的

QQQQ3::::

http://blog.sciencenet.cn/static/js/grey.gif
看到那个saturation levels 吗？最大的是2.98023（100%），最小的是0.5871（0%），而右边的标尺是根据那个把最大和最小分成十份，你自己算一下标上去就行了

（9） 立体图中，就是Isosurface，比如at a value of 0.75；画出来之后没找到标尺的选项啊，请问能给出立体图的标尺吗？谢谢比如说下图

【如何把Slice图放到结构当中，像下面这样？】
http://blog.sciencenet.cn/static/js/grey.gif

QQQQ::

这个图可能不是vesta做的
        我只能在截取一个晶面时显示相应的isosurface【isosurface的物理含义？】，但此时section中小于该值的就不会显示。要做可能也是两张合成的。对于标尺是针对section的，这个是不变的，你在2D里面看好数据，把标尺弄出来，就可以了。
对于这个图可以这么做：
1，先从2D里面吧标尺弄出来放在图片一中，,
2，然后在boundary里面设置晶体的晶面，返回带properties里面选择isosurfaces，在这里修改isosurface level为0（里面的颜色不会随你的这个值而改变，设置为零（你的电荷的最低值，也就是你标尺最下端得值）一般能把所有的数据显示出来），这样可以得到图二，
3，返回主界面，把Volumetric data 里面的 show section前面的钩去掉，再到properties里面选择isosurfaces，在这里修改isosurface level为设置为你想要的值，但此时这个颜色会和那个标尺不一样（默认的是黄色）。你可以现在标尺中找出你想要值的颜色，使用其他的程序读取RGB值，然后再在color中填上就可以。如果你都不了那个RGB值的话，就靠你的眼睛去判断，点后面那个按钮。可以选择你认为是一样的值。得到图三
4，在画图版中，使用那个白色透明模式（工具栏中最下面有两个大的黄色的下面的一个），把这几张图合并起来就行。
在整个过程中你不要去调图像看看图像怎么样，或者你看好之后用orientation，回到以前的视图角度，这样可以保证两张图是可以完全重合的

（10）空位、间隙原子 密度图的处理

   刚刚接触VESTA，想用VESTA处理VASP的运算结果，有两个问题求教虫友：
（a）vasp算出CHGCAR后可以导入VESTA，给出整个超胞的电荷密度图，但我只关心其中某几个原子附近的电荷密度分布，求教高手，能否只显示我关心的这几个原子上的电荷密度，如果能够，如何实现？
（b）超胞计算中引入了空位，能否在图示中用虚线球给出该空位的理想位置，并且用虚线表示出其本来的成键情况，如果能够，如何操作？
非常感谢！

QQQ:::vesta 似乎是没有办法只显示几个原子的电荷密度的，只关心几个原子，是不是可以考虑差分电荷分析？就是把你要关注的几个原子作为一个体系（去掉其他原子,仍用原来晶格）保持坐标不动，单独算一下，得到CHGCAR1， 然后把其他部分的原子（去掉关心的原子）保持坐标不变再单独算一下，得到CHGCAR2， 然后用总体系的CHGCAR-CHGCAR1-CHCAR2 得到CHGCAR——difference。 然后再用vesta 做CHGCAR-difference。但和你的要求可能有点区别。

引入空位的话，vesta 似乎不能做到你的要求，但是你可以在vesta里面作结构图，然后导出图片，在ppt或者Photoshop 中修改图片。可以得到满足你要求的图片

QQQ：：：：

只显示几个原子而不显示其他原子的电荷密度很奇怪的，没有太多意义。我建议你采用截面的方法来显示某几个原子附近的电荷密度分布

（11）再用VESTA画部分电荷密度（partial charge density）的时候，isosurface level 一般取值为多少。有什么原则吗,比如最大值的多少分之一。我刚刚学不清楚。请高手指点下，提前谢

（12）A：：：选择适当的平面使要观察的原子的ELF落在这个平面上

在用VESTA做ELFCAR的2D图的时候要选择适当的平面使要观察的原子的ELF落在这个平面上，我的问题是用VESTA作2D图时显示不出原子，那么怎么自己加原子呢？比如下面两个图：图1是3D图，其中紫色的是Zn，灰色是O，黄色是cu，图2选取的001面，那么在图2的那些区域，怎样确定那个是In产生的哪个是粗cu产生的呢？
又或者：怎样让001面上只显示出cu及其周围的Zn和O？
Q：：：记得slice里面不仅可以选哪个面，还可以定义面的位置，比如在最底下（z=0），或者中间（z=0.5）……
应该是可以定义的。

你图片上看到“slice”这个键了吗？点击就出现那个对话框里，h、k、l 三个指数定义面，下面有个“Range of projection”，计算定义原点到面的距离，一个的以angstrom为单位，另一个是分数坐标吧。
不知道怎么回事，上传不了图，不好意思啊。

（13）电荷密度图单位问题：CHGCAR的单位应该是e, 要获得真正的电荷密度需要再除以晶包体积【怎样除以晶胞体积】

  VESTA读取的CHGCAR单位是e/bohr^3 如果想要改成e/A^3 就要在Edit--> Edit data--> volumetric  data --> import...  ---> convert the unit 里面选择bohr^3 to  angstrom^3  
我想问的是有没有简便的设置不用每次都这样改？太麻烦了！

QQ：：

问一下楼主，你从哪里知道CHGCAR的单位就是e/Bohr^3呢？
我最近看了一些相关资料知道，CHGCAR的单位应该是e, 要获得真正的电荷密度需要再除以晶包体积，如果你晶胞体积采用的是Bohr为单位，那么得到的电荷密度就是以e/Bohr^3为单位i，反之，如果晶胞采用A为单位，那么电荷密度就应该是以e/A^3为单位。

VASP说明书上有这样一句话：
“CHGCAR file
This file contains the lattice vectors, atomic coordinates, the total charge density multiplied by the volume  on the fine FFT-grid (NG(X,Y,Z)F)”
“ Please do not forget to divide by the volume before visualizing the file!”

(14) 差分电荷密度的另一种算法：Δρ=ρ(AB)-ρ(A)-ρ(B)

Q：：：

会不会是因为INCAR里面某些参数是使用默认值，而这些值是根据POTCAR读来的，结果因为POTCAR不同而导致INCAR也发生了变化？比如ENCUT是可以不写的，如果3个INCAR里面都没有这一行，而是采用PREC=Accurate的设置的话，那么AB那一次计算的ENCUT会跟ENMAX比较高的那个（比如B）的一致，而ENMAX比较低的那个（比如A）算出来的就是另一回事了。不妨把那些可以手动设置的参数自行设置一下，比如ENCUT、NGX、NGY、 NGZ和NGXF、 NGYF、NGZF之类的…当然我理解也不是很深，不确定是不是这种情况。

VESTA软件的用法讨论：摘自大家identation：

http://emuch.net/bbs/viewthread.php?tid=6422652&fpage=1

请大家交流一下VESTA软件的用法

这个取值是你自己取得，就看画出的图能不能描述出问题的所在，比如你中心原子电荷是0.5，外部空间是0.2，那你就不能将0-1取为一个颜色，因为这样你无法分辨外部中心原子和外部空间。VESTA作图的单位我还真不清楚，我做密度图的时候都是用lev00然后导出数据用origin画的，不过好像VESTA的单位和VASP的单位不太一样，这个你自己查一下

a）首先，VESTA的作用可以是--画晶体结构，当你知道了晶体的空间群，原子的Wyckoff坐标之后，就可以用VESTA画结构了：

(1)File->New Structure->Unit cell->Space group处输入空间群的编号，Setting处可以对该空间群编号下的不同表述方式进行调整，例如No. 62号空间群属于正交晶系，输入62在Space group处之后，如果setting是1，则显示为Pnma，如果setting为3，则显示为Pbnm，实际上这两个空间群是等价的。

（2）接下来，structure parameters处，点击New，然后在Symbol处选择你想输入的元素，Label处输入你想显示的标签，同时把x，y和z后面的空格处填上该元素的坐标，则该原子输入完毕，再点击New，输入其他原子坐标，直至完成。

（3）完成之后，你将看到你所画的晶体结构图。

（4）对于该图，你可以有多种显示方式，例如键线式，八面体或四面体等等。这需要如下步骤来做到--点击Edit->Bonds->点击New，然后将A1和A2处填上成键的原子，例如考虑体系BiFeO3，如果想显示Bi-O八面体，则可把A1处的原子写成Bi，A2处写成O，然后注意调整Min. length和Max. length处的数值，这两个值的意思是，如果A2原子到A1原子的距离在Min. length和Max. length数值之间，则认为这两个原子成键，两原子直接连一线。然后点击OK即可。

（5）如果要八面体或者四面体显示该晶体，在完成前4步之后，点击Style按钮，选上Polyhedral，则可以以多面体显示。

（6）Objects->Properties处可以改变显示的风格--例如可以改变多面体的颜色，可以使用超晶报显示等等。

      完成上述步骤后，点击File->Export data，你将把你画的晶体以不同格式导出，包括cif格式以及vasp格式等等，vasp格式导出的结果，第一行是你造的晶体的名称，由你来取，默认为New structure；第二行一般为1.0，是基矢的缩放系数，接下来的三列是基矢，也就是格子的abc轴的描述矢量；接下来是元素种类，接下来是不同种类元素的原子个数，接下来是原子坐标--导出的这部分数据可以作为VASP的POSCAR，也可以经过拷贝黏贴之后修改格式作为PWSCF,ABINIT,ELK等软件的晶体坐标输入。

此外，也可以将所画晶体用图片格式导出，可以导出矢量图也可以导出普通的图片，请点击File->Export Raster Image或File->Export Vector Image

（b）VESTA可以显示电荷密度，

        VASP的计算结果，直接用VESTA打开CHGCAR或CHG文件，你将获得显示其原子以及电荷密度。这个也可以调整，（1）首先是3D显示，可以如下修改或者调整：Objects->Properties->Isosurfaces->Isosurface level处的数值做适当调整，你将显示出漂亮的电荷密度等高线图。

（2）Utilities->2D data display  将使你可以显示出电荷密度等高线的2维投影，输入你想显示的晶面指数，当然你可以实现选定一个原子，输入晶面指数之后，再点击Calculate the best plane for the selected atoms，你将获得通过这个原子，且是你选定晶面上的2D电荷密度。

（3）你同样可以试着在Edit->Lattice planes处做类似操作，获得在原晶体结构上显示的某一晶面的2D电荷密度投影。

（c）改变基矢：

     对于一个已经用VESTA打开的晶体结构，你可以点击Edit->Edit data->unit cell，然后点击Option按钮，在弹出对话框的Rotation matrix P处输入基矢的变换矩阵，你将对晶体的基矢作一定的变换。

最最简单的玩法是，当你把Rotation matrix P第三行第三列的那个数值由1改成2的时候，你将获得z方向的超晶报；同理如果把第二行第二列的数值由1改为2，你将或者y方向的超晶胞。

VESTA是一个功能很强大的晶体建模显示工具，我个人认为其基本上可以囊括MS的可视化模块，而且显示起来的图更好看。也希望大家多多交流，说说你所知道的VESTA的用法或作用。有金币悬赏！

(d) 講一些關於 isosurfaces 分析: 這是因為 VESTA 讀入 CHGCAR 後 會將所有的值除以體積，且單位是用 Bohr?
1) VESTA 可以直接讀取 vasp 的 CHGCAR，意思當然就是說只要檔案格式內容跟 CHGCAR 一樣
    VESTA都可以直接讀，比方 
              ELFCAR (electron localization function)
              LOCPOT (local potential)
              PARCHG (partial charge density)
    再延伸出去，對任何空間相關資料只要把格式轉成像 CHGCAR 般~即可給 VESTA 讀取
    (若發生讀取錯誤，請把副檔名改成 .vasp 即可)
2) 在 isosurfaces 分析過程有許多參數可以調整，東調西調~調完一大堆參數後
    建議將它存檔~存成VESTA 自己的 .vesta 檔 (File --> Save As)
    下次打開只要開這個 .vesta 檔就可以出現調整後的結果
    這個 .vesta 檔是可以讀取，裡面有一大堆資料
    第九行 IMPORT_DENSITY 1 的下方~第十行
    就是你所給予的 CHGCAR ~
   也就是說~如果你今天要分析 PARCHG，你可能有很多檔PARCHG.0001, ... PARCHG.xxxx
    每個檔案的結構是一樣，但電荷分部不一樣
    這時你只要先處理好一個檔案~接著把這個檔案複製成許多檔案
    然後修改每個檔案的第十行~來改變 import file position
    就可以非常快速地分析 partial charge
3) 在 VESTA , Properties --> Isosurfaces 介面
    會顯示 isosurfaces 的最大與最小值
    這個值跟你所給予的 CHGCAR 內容是不一致的
    這是因為 VESTA 讀入 CHGCAR 後
    會將所有的值除以體積，且單位是用 Bohr
    (VASP 所有長度單位都是用 Angstrom)
    因此 VESTA 顯示的最大最小值，是除以 體積(unit in Bohr^3) 後的最大最小值
4) 把 CHGCAR 丟入 VESTA 之後
    VESTA 很貼心的會立刻幫你選了某個 isovalue
    這個值不是隨便選的
    值的公式在手冊的 16.7 Default isosurface level
    在這公式內有個參數 n ，手冊上說明這個 n is a parameter to adjust d(iso)
    n 可以在 Edit --> Preferences 內的 Default isosurface level 調整
5) 做出 3D isosurface 可以配合 Edit --> Lattice Planes 做剖面分析
    這個切面也是可以調整設定
    在 Properties --> Sections

 

 

 

 

VESTA画局域电荷密度图 ELF

http://muchong.com/html/201508/9267921.html

VESTA可以切片
Utilities---2D DATE   选择面指数、厚度等，调整等值面，就可以了

就在VESTAT中调整啊  
1、体材料上，调整Isosurface level
2、切面：只能调整显示    调整sauration levels，一般按照后边的百分比调整（先0-100看看再微调）

电子局域函数源表征:电子的局域化分布特征:用于确定成键类型和找出孤对电子分布的理由。

http://blog.sina.com.cn/s/blog_6813cb430100tvka.html

     电子局域函数源于 J. Chem. Phys. 92(1990)这篇文章，作者定义了一个ELF（Electron Localization Function，电子局域密度函数），以此来表征电子的局域化分布特征。优点是无需计算局域分子轨道，计算量较小。
       具体的推导过程请大家自行查看文献，我们主要关心的是它的用途，也就是为什么它能被用于确定成键类型和找出孤对电子分布的理由。
        按照定义，ELF的值在０到1之间，取上限值1表示电子完全局域化，而0值可能表示电子完全离域化（当然也可能是表示该处没有电子http://blog.sciencenet.cn/static/js/grey.gif ）而ELF=1/2这个中间值则表示该处电子形成了类似于电子气的电子对分布(原文是electron-gas-like pair probability)。
        VASP可以计算ELF的分布情况，数据可以在ELFCAR中找到，可以用VASPVIEW或VESTA等软件作图查看。