第一性原理计算软件攻略-利用VESTA绘制差分电荷密度图解

1. 软件介绍

VESTA全称为Visualization for electronic and structural analysis. 是一款免费的可视化软件，可以帮助工作者利用第一性原理计算得到的数据制作出赏心悦目的图片。并且VESTA可以运行Materials Studio建模软件中的部分输出文件，通过VESTA和MS（Materials Studio）两款软件，可以完成第一性原理计算中的大部分模型处理和作图分析。VESTA可以读取MS生成的CIF文件，并将CIF文件直接转化为计算所用的POSCAR文件。一定程度下可以节省在MS里通过自己建设基矢导出原子位置的时间。另外，VESTA也可直接将POSCAR导出格式为*.xyz的文件（也就是MS识别的原子坐标文件），将*.xyz文件导入MS里可以调控POSCAR里的原子位置。添加吸附原子、分子等工作。图1为运行VESTA软件的界面显示。

图1

如图所示，左面栏里黑色箭头代表模型三维空间的立体展现，通过点击黑色箭头，在右方结构模型处移动鼠标可以实现模型的三维旋转。而第二个白色标号箭头，则可以实现原子的选取。其中十字标号可以使模型在可视化窗口中进行上下左右的移动。而图2这几个标号，分别代表两个原子之间的距离、三个原子之间呈现的夹角、平面角和界面之间的夹角。如图3所示，界面栏上方的a, b, c和a*, b*和c*则可以直接实现模型在x、y和z轴上面的视角。

图2

图3

2. 差分电荷密度图

差分电荷密度图是通过成键后的电荷密度与对应的点的原子电荷密度相减获得。通过差分电荷密度的计算和分析，可以清楚地得到在成键和成键电子耦合过程中的电荷移动以及成键极化方向等性质。差分电荷密度图在第一性原理计算中的应用多集中于分析结构优化后模型的原子间成键情况，尤其在界面计算中，通过差分电荷密度分析，可以初步直观的判断界面的结合方式。举两种材料组成的界面为例，差分电荷密度的计算公式为： ∆ρ = ρAB − ρA – ρB。ρAB代表了界面优化后的结构电荷密度，ρA代表组成界面的材料A的电荷密度，ρB代表组成界面的材料B的电荷密度。将它们作减，即可得到材料A和材料B组成界面AB后的电荷密度变化，进而可以分析界面中的电荷移动等性质。如图4所示，显示了石墨烯、掺杂石墨烯和氧化锂组成界面后的差分电荷密度图。其中黄色区域代表电子集聚的区域，蓝色区域表示电子损失的区域。通过差分电荷密度图，我们可以看出石墨烯和氧化锂组成界面后，石墨烯的电子发生重排，集聚在石墨烯和氧化锂的界面中。

图4

在计算研究体系的差分电荷密度时，界面AB，材料A和材料B的模型应保持同一，也就是说需要放在相同的大小格子中。当得到差分电荷密度的输出文件后，可以在各单独体系中的OUTCAR检验界面AB，材料A和材料B的NGX, NGY, NGZ是否一致，避免因为格点不同而导致差分电荷密度图无法正常绘制的情况发生。

如果在OUTCAR里发现参数并不相同时，需要在INCAR里通过设置NGX，NGY，NGZ保持一致使数据正常输出。使用VESTA绘制差分电荷密度图时，需要输出文件CHGCAR，例如要得到A、B两相组成的界面体系的差分电荷密度图，则需要得到界面AB、材料A和材料B各自的CHGCAR。然后利用上述公式∆ρ = ρAB − ρA – ρB 将对应的CHGCAR在VESTA里进行相减则可以达成目的。

3. 如何利用VESTA作图差分电荷密度图

在利用VESTA作图差分电荷密度图时，为了方便用图示解释过程，我们依旧采用石墨烯/氧化锂的界面模型作图。

（1）首先将总体系的CHGCAR导入VESTA里，在这里，直接将文件夹里的CHGCAR拖入VESTA就可以，CHGCAR文件一般为一百兆或者更大，所以，在VESTA导入CHGCAR时，如果响应比较慢的话，请耐心等待，如果强行运行，可能导致软件崩溃。

图5

（2）在导入界面总体系的CHGCAR后，导入材料A的CHGCAR；在导入材料A的CHGCAR时，要在上一个步骤的基础上进行导入，依次点击上栏中的Edit-Edit Data-Volumetric Data。

图6

此后，界面显示如下图所示，然后点击Import，导入材料A的CHGCAR

图7

在导入材料A的CHGCAR后，界面显示如图8所示，点击Subtract from current data和Raw Data，最后点击确定ok。备注：点击OK后反应时间较长，请耐心等待。

图8

在出现如下图所示的界面后，则导入材料A的CHGCAR已经成功，按照上述程序，再重复导入材料B的CHGCAR。

图9

导入材料B的CHGCAR后界面如图10所示，此时点击OK得到该体系的差分电荷密度图的初步显示。

图10

图11

点击Show Section对绘制的图进行简单处理，可以看到图11的小蓝色区域消失，如图12所示。

图12

此时，调整等值面(isosurface level)的数值，美化图形就可以得到完整的差分电荷密度图了。不同材料的等值面设置不同，根据自己的体系进行调整。

图13

图14

调整完后，就可以导出图形了，此时点击File-Export-Raster Image就可以导出到自己数据图形放置的文件夹了。之后采用PS对图形进行进一步美化，就可以得到漂亮的差分电荷密度密度图了，如图15所示。

图15

配图中所用的参考文献
1. Wang, T.; Zhao, N.; Shi, C.; Ma, L.; He, F.; He, C.; Li, J.; Liu, E., Interface and Doping Effects on Li Ion Storage Behavior of Graphene/Li2O. The Journal of Physical Chemistry C 2017.

也可参照

http://blog.sciencenet.cn/blog-567091-736154.html

第一性原理计算软件攻略-利用VESTA绘制差分电荷密度图解相关推荐

第一性原理计算软件FLEUR: The Jülich FLAPW code family
第一性原理计算软件FLEUR: The Jülich FLAPW code family 链接为:http://www.flapw.de/pm/index.php QE和Abinit都是采用平面波加赝 ...
JAVA程序员装机必备软件攻略
因为自己的电脑年久失修,今天用tomcat跑了同时跑了两个项目电脑卡死不动了,于是用了单位的电脑,装了一下午软件,晚上终于成功的将项目跑了起来,我把我自己的软件攻略分享给大家,希望可以帮助你们更加有效 ...
3岁宝宝记忆力测试软件,攻略：2-3岁宝宝记忆力大作战！
原标题:攻略:2-3岁宝宝记忆力大作战! 小宝贝们,好久不见. 受了些挫折,病也还没好,就要跟大家一起奔赴战场啦!一边受挫折,一边迎难而上:一边疲惫不堪,一边咬牙坚持,这不就是妈妈们的常态吗?怎么听上 ...
怎样画流程图攻略：流程图绘制软件使用方法
我们在日常的办公中,经常需要绘制一些流程图,借助这些流程图,除了可以让自己更加清晰地认识到事情的每个细节之外,也可以让其他人知道你的想法.那么,对于流程图的绘制,借助软件会更加方便.接下来就要来和大家 ...
matlab生成计算软件,利用Matlab从第一性原理计算软件Siesta读取和生成能带图
Siesta关于E-k关系计算结果保存于systemlabel.bands文件内,故Matlab程序只需要读取此文档即可.经笔者测试,在siesta-2.0.2与siesta-3.2版本内均可完美运行 ...
针对计算机网络或者计算,北京自考“计算机网络与通信”计算题攻略（2）
35.海明码中的信息位为7位,在接收端收到的报文为11110111011,求: (1)问传输是否有错? (2)求发送端发送的信息位. CRC码 1.设信息码为1011,采用的生成多项式为G(X)=X4 ...
iOS应用开发攻略 pdf电子书
尊敬的用户您好,由于该iOS应用开发攻略pdf书受百度网盘影响无法做公共分享,只能私密分享,有不到之处请多多谅解! 百度链接: http://pan.baidu.com/s/1xPovn 密码: ...
Wordhero AI智能写作工具：购买、攻略、资源 -揭秘如何根据关键词和标题生成段落
Wordhero AI是我经过仔细对比,实际应用之后获得相应效果之后,坚持使用的AI写作工具.为了持续获得这方面的提升和反馈,我建立了两个渠道去让自己不断获得国内外同样工具使用者的反馈和建议: AI ...
国内外量子计算软件和量子计算语言
1.量子计算软件国外量子计算软件 1.pyquil pyquil 是美国的量子计算公司 Rigetti 开发的,可以实现量子线路的计算模拟.含噪声的量子逻辑门计算模拟,量子芯片的云端运行等功能.P ...

第一性原理计算软件攻略-利用VESTA绘制差分电荷密度图解

第一性原理计算软件攻略-利用VESTA绘制差分电荷密度图解相关推荐

最新文章

热门文章