这是本文档旧的修订版!
版本:2017
在本教程中,您将学习如何使用 ATK 计算金属表面功函数。您将创建 Ag(100)表面的平板模型,利用鬼原子改进其电子属性的描述,并使用合适的边界条件正确描述其上方真空的存在。假设您已有一些使用 ATK 的经验。
功函数的计算主要依赖于精确地模拟表面电荷密度进入到真空中的衰减,但短程的局部基组有时不足以满足要求。在这种情况下,在表面上方放置额外的 LCAO 基函数 - 这些称为鬼原子。
在本节中,您将设置含有鬼原子的 Ag(100) 表面平板模型,并计算其表面功函数。
1.打开 Builder(按钮 )。点击 Add Add From Database,搜索 “silver” 并将其添加到 Stash。
2.打开 Builders Surface (Cleave)….,利用 Surface (Cleave) 插件切割垂直于 [100] 方向上的块体银定义曲面。
3.点击 “Next” 两次。在 Surface (Cleave) 插件的最后一页,设置平面外晶胞矢量 V3 为 “Non-periodic and normal”。按照如下设置其他的参数:
这些参数设置将确保功函数的收敛。
在垂直于表面的方向上,系统的非周期性需要适当的边界条件。
4.点击 “Finish”。该表面结构将出现在 Stash。
5.为完成几何设置,请在 Builder 中选择 Z 坐标最大的原子,然后通过单击图标 将其转换为鬼原子。
6.结果,平板由 13 层组成。在表面上的鬼原子的作用是描述真空中电子密度的平滑衰减。
7.将结构发送到 Script Generator (按钮 )
将原子变为鬼原子意味着将它的势(如赝势或核电荷)设置为零。然而,部署在那里的基函数仍用于描述电子密度。另一种方法是扩大在每个原子位置处基函数的范围,但如果系统不像本例中那么简单,则会使计算成本高昂。
1.在 Script Generator,双击 Blocks 面板的图标添加一个 New Calculator。并点击 Analysis ChemicalPotential 计算平板的化学势。
如下所示在平板的顶部和底部施加了适当的边界条件使得化学势等于功函数。
2.为 NetCDF 文件设置一个合适的文件名,如 “Ag_100.nc”。
3.双击 Script 面板上的 New Calculator,设置如下参数:
谨记:一般情况下,要检验与 k 点网格精细度有关性质的收敛性。