目录

沉积超薄氧化层调控银表面功函数

版本：2015.1

沉积在金属基底上的绝缘材料超薄膜构成了一类特殊的材料，在不同的领域具有可调特性和增长潜力 [GP11] [Pac12] [PH13]。在本教程中，您将学习如何通过沉积另一种材料的覆盖层来更改基材的功函数。

注意

功函数的调整

在金属基板上沉积薄绝缘膜的一个重要结果就是金属载体功函数的诱导变化，会根据界面的性质而降低或增加。在文献中已经报道过这种功函数变化的例子。利用开尔文探针力显微镜或扫描隧道显微镜研究 Au(111) 和 Ag(100) 上的碱金属氯化物薄膜，显示功函数有 0.5-1.2 eV [LZE04] [PPS + 05] 的降低。另一项基于场发射共振的研究发现，高达3 ML的NaCl岛状物的功函数降低了 1.3 eV [PBP + 07]。理论计算预测了各种金属上的 NaCl，MgO 和其他氧化物功函数的降低 [BTI04] [GCP06] [PMP08]。

您将计算由于沉积 1 至 3 层绝缘 MgO 形成金属的 Ag(100) 表面功函数变化。计算功函数的过程遵循教程 Computing the work function of a metal surface using ghost atoms 中给出的指示。，

特别地，您将：

  1. 创建并优化 Ag 和 MgO 的块体结构；
  2. 构造 Ag(100) 和 MgO(100) 表面；
  3. 构建 MgO(100)/Ag(100) 界面；
  4. 设置功函数的计算并运行；
  5. 分析结果，与文献做比较。

警告

计算设置

在本教程中，您将使用一组特定的计算设置（网格截断，k 点取样，基组，交换关联，金属层数等），这些设置的选择使得 QuantumATK 的结果可以与文献比较。但请务必记住，您应该检查您的结果在这些设置下是否适当地收敛。

Ag(100) 和 MgO(100) 表面

按照以下步骤创建银和氧化镁块体：

在本教程中，您将使用 PW91 交换关联函数进行 DFT 模拟，以便与文献 [PMP08] 中 Prada 等人的结果进行比较。PW91 中 Ag 的晶格常数为 4.16Å，因此将此值用于银的块体构型：

提示

您还可以使用 PW9 1和 11x11x11 的 k 点简单地执行银块体构型的 DFT 几何优化，从而弛豫晶胞。对于晶格常数，这种更一般的方法的结果约为 4.16 Å。

接下来，构建 Ag(100) 表面：

按照相同的步骤创建一个 4 层的 MgO(100) 表面，然后将两个新的 Stash 项目分别重命名为 “Ag(100)” 和 “MgO(100)”。

Ag/MgO 界面

您现在应该使用 Builders Interface 插件构建 MgO(100)/Ag(100) 界面：

提示

您可以在教程 The Interface Builder in QuantumATK 的技术说明中了解更过有关创建界面的更多信息。

DFT 计算

如教程 Computing the work function of a metal surface using ghost atoms 中所述，您需要在表面或界面上方添加“鬼原子”。功函数的计算需要对延伸到真空中的电荷密度进行非常好地描述，鬼原子正好可以保障这点。

添加鬼原子

ATK-DFT 计算

当表面之上晶胞边界的有效势为零时，功函数就被认为是化学势。Dirichlet 边界条件 （BC） 则用于强制执行这点。您还将计算有效势，以便直观地检验。按照以下所列方法设置所需的 DFT 计算：

保存脚本为 MgO3LAg.py，并用 Job Manager 执行。如有需要，您也可以在此处 ↓ MgO3LAg.py 下载最终的脚本。如果是在一个现代的笔记本电脑上连续执行计算，应该仅需要 5 分钟就可以完成。如果是用 4 个 MPI 进程并行执行，则可减少到 2.5 分钟。

分析结果

OptimizeGeometry，ChemicalPotential 和 EffectivePotential 数据块现在应该已出现在 QuantumATK 的 LabFloor 上。

尝试选择 OptimizeGeometry 数据块，并通过单击 LabFloor 右侧的 Viewer 插件可视化弛豫轨迹。单击 按钮启动视频。确认受约束的 Ag 和鬼原子在弛豫期间确实是固定的。

功函数

选择 ChemicalPotential 数据块，然后单击 Show Text Representation 插件以读取出计算得到的化学势为 -2.99 eV。

因此，这种 3 层 MgO 在 Ag(100) 上的功函数为 2.99 eV，这与 Prada 等人在文献中报道的计算值 2.96 eV 非常一致[PMP08]。

您还可以按照上述步骤构建 2L-MgO/Ag(100)，1L-MgO/Ag(100) 和 Ag(100) 系统。另一种方法是使用弛豫过的 3L-MgO/Ag(100) 构型作为其他系统的起点。这将减少弛豫这些系统所需的 BFGS 步数，从而节省计算时间。

例如，创建 2L-MgO/Ag 系统，计算功函数：

表 4 计算得到的 PW91 功函数（eV）。对应的功函数差异为 wrt 是为了清除括号中的 Ag(100)。

QuantumATK Pada et al. [PMP08]
Ag(100) 4.22 4.23
1L-MgO/Ag 3.31 (-0.91) 3.29 (-0.94)
2L-MgO/Ag 2.97 (-1.25) 2.95 (-1.28)
3L-MgO/Ag 2.99 (-1.23) 2.96 (-1.27)

可以在此处下载计算表中所有标记为 QuantumATK 的功函数所需要的脚本：↓ Ag100.py，↓ MgO1LAg.py，↓ MgO2LAg.py，↓ MgO3LAg.py。

提示

上表显示了 QuantumATK 与文献 （VASP） 中功函数计算之间具有良好的一致性。但是，如果计算设置发生变化，结果可能会改变。例如，SZP 基组是用于 Ag 原子的—— DZP 基组可能会给出略微不同的结果。所使用的赝势的类型也可能影响结果，并且在某些情况下可能需要更多的鬼原子。

有效势

您已使用一组特定的边界条件进行功函数的计算——左 C 面上的 Neumann 和右 C 面上的 Dirichlet。您现在可以使用 1D Projector 插件显示计算中的平均有效势：

Ag(100) 和 MgO(100) 区域的有效势明显不同。此外，两个不同 BC 的影响从超胞两端的势值和斜率就可以非常清楚地看出：

1D 投影插件

1D 投影可用于将各种 3D 网格数据投影到 1D 表示。这非常有助于实现可视化的目的，适用于 QuantumATK 网格对象（参见方框）的较广范围。

常用的 QuantumATK 网格对象： BlochState，EffectivePotential，Eigenstate，ElectronDensity, ElectronDifferenceDensity，ElectrostaticDifferencePotential，ExchangeCorrelationPotential，ExternalPotential，LocalDeviceDensityOfStates，TransmissionEigenstate，ElectronLocalizationFunction。

插件小程序中提供了几个选项：

Grid

您可以打开投影工具，在 LabFloor 上选择多个对象，将它们彼此相邻绘制。在这里，您可以选择要绘制哪一个。

Axis

选择您想要要投影 3D 数据网格的方向。

Projection type

对垂直于所选方向的平面中的所有数据求和或求平均值。您还可以沿着穿过特定投影点的直线绘制单个值。

Projection point

在分数坐标中指定选择投影类型为 Through point 时使用的投影点。

Spin projection

在自旋极化计算的情况下，共线或非共线，您可以选择特定的旋转投影。

Add line

上述选项指定后，单击此按钮可在窗口右侧绘制投影。您可以在同一个图中添加更多投影。

Remove line

在 Projection Plot 窗口中选择一行，然后单击以从图中删除这条线。

Clear plot

从绘图中删除所有线条。

Line Info

显示了与当前所选绘图线/点相关的一些有用信息。请注意，该图是交互式的。单击绘图的任意点可以输出相应的信息。

Projection Plot

右键单击可缩放，自定义或将数据导出到文件。

参考