用户工具

站点工具

本页面的其他翻译:
  • zh

adf:2d-peda

这是本文档旧的修订版!


如何进行二维周期性体系的结合能、能量分解pEDA

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窗口即可。

第二步:设置计算参数

严格的说,需要先对吸附模型进行结构优化,但本例仅仅为了演示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能量各项的含义

adf/2d-peda.1575871701.txt.gz · 最后更改: 2019/12/09 14:08 由 liu.jun

© 2014-2022 费米科技(京ICP备14023855号