atk:利用鬼原子计算表面功函数
差别
这里会显示出您选择的修订版和当前版本之间的差别。
| 两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
| atk:利用鬼原子计算表面功函数 [2019/02/25 21:40] – [计算和分析] xie.congwei | atk:利用鬼原子计算表面功函数 [2019/02/28 10:02] (当前版本) – [注解] xie.congwei | ||
|---|---|---|---|
| 行 147: | 行 147: | ||
| | Tl (100) | 3.84 | 3.74 | | | Tl (100) | 3.84 | 3.74 | | ||
| | W (100) | 4.32–5.22 | | W (100) | 4.32–5.22 | ||
| + | |||
| + | |||
| + | |||
| + | [1] 所有计算均使用 DoubleZetaPolarized 基组和 9x9x1 的k点网格。其余参数如上例所示设置。请注意，这意味着所有值不一定会在这些参数下收敛。表面也不需要弛豫。 | ||
| + | |||
| + | [2] 已知 Ag（110） 表面经历了重建；而此处显示的计算值不考虑这一点。 | ||
| + | |||
| + | [3] 意外地，用 SingleZetaPolarized 计算得到的金的功函数为 5.28 eV。 | ||
| + | |||
| + | |||
| + | |||
| + | |||
| ===== 注解 ===== | ===== 注解 ===== | ||
| + | * 我们在上文的一些地方说明了仔细检查某些数值的收敛和几何参数的重要性。功函数对所有这些参数都非常敏感，因此，为方便起见，我们总结一下最重要的几个： | ||
| + | * **Number of slab layers**：对于银，测试表明，如上所述的 13 层是足够的，但任何较小的数字都将会产生不同的结果； | ||
| + | * **Extension of vacuum on both sides**：在本例中使用的右晶胞边界上的 20 Å 是安全的；至少对于 Ag（100） 来说，即使只有 10 Å 也足够了； | ||
| + | * **K-point sampling**：9x9x1 可以给出合理的结果，但一如既往地，k 点取样这个参数的收敛不是简单的单调函数； | ||
| + | * **Basis set**：对于过渡金属，需要极化轨道来获得正确合适的原子函数，因此我们为银选择 “SingleZetaPolarized”。但对其他类别的元素，可能适用其他规则。与元素无关的话，“DoubleZeta” 基组是一个很好的选择。“DoubleZetaPolarized” 总会是一个更安全的选择，虽然它需要更多的计算时间； | ||
| + | * **Exchange-correlation functional**：LDA 通常能够很好地描述“简单”的金属，但显然并不适用所有材料。 | ||
| + | * 由于 k 点（使 QuantumATK 可以很好并行化的参数之一）的数量大，脚本可以使用 MPI 并行运行以获得相当大的执行优势。 | ||
| + | * 通过比较上表中的 Ag 的（100），（110）和（111）晶面的值，显然功函数取决于所考虑的晶面。 | ||
| ===== 参考 ===== | ===== 参考 ===== | ||
| + | * 英文原文：[[https:// | ||
atk/利用鬼原子计算表面功函数.1551102018.txt.gz · 最后更改: 2019/02/25 21:40 由 xie.congwei