这是本文档旧的修订版!
BAND中的二维体系是真正的二维体系,并不需要加真空层。在Main菜单中,Periodicity一栏显示的是Slab。因此BAND对二维、一维体系的计算效率高于VASP等平面波程序。
由于pEDA分析支持正常k空间布点,而NOCV则只支持Gamma点,因此建议二者分别计算。本文使用AMS2019.301完成。
Mg 0.59538358 0.59538226 0.20750000 Mg 2.08384335 -0.89307752 -1.89750000 Mg -2.38153597 0.59538226 0.20750000 Mg -0.89307620 -0.89307752 -1.89750000 O 2.08384335 -0.89307752 0.20750000 O -0.89307620 -0.89307752 0.20750000 O 0.59538358 0.59538226 -1.89750000 O -2.38153597 0.59538226 -1.89750000 O 0.59538553 0.59539200 3.95250000 C 0.59538493 0.59538903 2.80750000 VEC1 2.97691955 -2.97691955 0.00000000 VEC2 2.97691955 2.97691955 0.00000000
把这段数据直接复制到Input窗口即可。注意,此时Main窗口的Periodicity自动修改为Slab了,也就是一个真正的二维材料,上下表面为半无限大的真空。
严格的说,需要先对吸附模型进行结构优化,但本例仅仅为了演示pEDA的功能,因此省去了结构优化,而假定上述得到的结构是优化好的结构。设置参数如下:
注意,泛函、基组的选择,请参考不同泛函的特点、ADF:一般情况下如何选择基组Basis Set、对同一种元素的不同原子指定不同基组
如果有原子不在Cell范围内,通过Edit > Crystal > Map Atoms to Cell。因为进行pEDA计算的时候,如果有原子不在Cell范围内可能会报错。
将H2和铁表面分成2个区,分区方法参考如何创建分区,如下图:
设置其他参数:
关于k点的设置:2019版可以用默认k点设置,也可以在Details → K-Space Integration → Number of Points输入具体k点个数,例如7 7;之前的版本需要选择Gamma Only(也就是1 1,只有一个1*1=1个k点)
设置完成后,File > Save As保存任务。
提交任务的方式,参考正式版的安装、维护与升级
查看结果文件,参考如何查看结果。
打开*.out文件,点击Properties > PEDA Energy Terms即可看到能量分解的情况。其中△Eint表示两个分区之间的结合能。其他各项能量的含义,请参考pEDA能量各项的含义。