这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版上一修订版两侧同时换到之后的修订记录 | ||
atk:si_100_表面重构_用atk进行几何结构优化 [2016/08/08 18:42] – [结果] liu.jun | atk:si_100_表面重构_用atk进行几何结构优化 [2016/08/13 10:07] – [结果] dong.dong | ||
---|---|---|---|
行 3: | 行 3: | ||
=====简介===== | =====简介===== | ||
- | 虽然LDA和GGA计算,对于Si(以及大部分其他半导体材料)得不到合适的带隙,但这些“简单”的泛函却能够很好的预测几何结构性质。在本文中,将演示如何使用ATK研究所谓的Si(100)表面不对称二聚重构。 | + | 虽然 LDA 和 GGA 计算,对于 Si(以及大部分其他半导体材料)得不到合适的带隙,但这些“简单”的泛函却能够很好的预测几何结构性质。在本文中,将演示如何使用 ATK 研究所谓的Si(100)表面不对称二聚重构。 |
- | 本文着重关注物理以及如何正确的设置模型,而非如何操作VNL。本文假定读者已经有了建模队经验(例如切开表面),以及设置计算、查看结果,因此提到具体步骤,不再非常详尽的解释。 | + | 本文着重关注物理以及如何正确的设置模型,而非如何操作 VNL。本文假定读者已经有了建模的经验(例如如何切开表面),以及设置计算、查看结果,因此提到具体步骤,不再详尽的解释。 |
=====创建结构===== | =====创建结构===== | ||
- | 1,启动VNL并打开Builder | + | 1,启动 VNL 并打开 Builder |
- | 2,从数据库添加alpha硅晶体结构(Add ‣ From Database) | + | 2,从数据库添加 alpha-硅晶体结构(Add ‣ From Database) |
3,打开Builders ‣ Surface (Cleave) | 3,打开Builders ‣ Surface (Cleave) | ||
- | | + | |
- | * 因为二聚重构只发生在2x1超胞,所以需要创建较大的表面单胞,而不是用最小表面单胞。因此在第二页窗口中设置V< | + | * 因为二聚重构只发生在 2x1 超胞,所以需要创建较大的表面单胞,而不是用最小表面单胞。因此在第二页窗口中设置$V_2=2u_2$ |
{{ : | {{ : | ||
- | 在表面建模工具的最后一页选择Slab厚度为4.5层(默认真空厚度是合适的) | + | 在表面建模工具的最后一页选择 Slab 厚度为 4.5 层(默认真空厚度是合适的) |
{{ : | {{ : | ||
行 26: | 行 26: | ||
4,将结构沿所有方向居中(Coordinate Tools ‣ Center) | 4,将结构沿所有方向居中(Coordinate Tools ‣ Center) | ||
- | 5,选择z坐标值最小的两个Si原子,点击几次Rattle工具{{: | + | 5,选择 z 坐标值最小的两个 Si 原子,点击几次 Rattle 工具 {{: |
- | 6,选择Z坐标值最大的两个Si原子,使用{{: | + | 6,选择 Z 坐标值最大的两个 Si 原子,使用 {{: |
生成的结构如下图所示: | 生成的结构如下图所示: | ||
行 36: | 行 36: | ||
=====设置计算===== | =====设置计算===== | ||
- | 1,依次插入New Calculator和OptimizeGeometry。 | + | * 依次插入 New Calculator 和 OptimizeGeometry。 |
- | + | | |
- | 2,Calculator中除了k-point sampling之外所有参数使用默认值。设置k-point sampling为9x9(当然c方向一个k点就够了,因为是层状结构)——精度非常重要,因为对称、不对称二聚状态的能量差很小((A. Ramstad, G. Brocks, and P.J. Kelly: Theoretical study of the Si(100) surface reconstruction. Phys. Rev. B 51, 14504, 1995))。 | + | |
- | + | | |
- | 3,优化,限制住底部的两个Si原子以及4个H原子。如果可能的话,限制弛豫的自由度总是很重要的,这样能够投影为纯的平移和转动。 | + | |
- | + | ||
- | 4,确保勾选了Save trajectory,并设置轨迹文件的名字 | + | |
{{ : | {{ : | ||
+ | * 设置输出文件的名字,现在已经准备好运行了。保存脚本! | ||
- | 5,设置输出文件的名字,现在已经准备好运行了。保存脚本! | + | 将脚本送入 Job Manager,或者将其上传到集群之后运行。该计算并行化非常好,根据 MPI 并行节点数差别,计算耗时约为 1 个小时或几个小时。 |
- | + | ||
- | 将脚本送入Job Manager,或者将其上传到集群之后运行。该计算并行化非常好,根据MPI并行节点数差别,计算耗时约为1个小时或几个小时。 | + | |
=====结果===== | =====结果===== | ||
行 57: | 行 52: | ||
{{ : | {{ : | ||
- | 分析结构弛豫的轨迹也很有意义。选择LabFloor中,你之前保存的轨迹文件,点击Viewer大卡,可以看到优化的动画过程: | + | 分析结构弛豫的轨迹也很有意义。选择LabFloor中,你之前保存的轨迹文件,点击Viewer大卡,可以看到优化的动画过程(点开大图): |
- | [[http:// | + | {{ :atk: |
如果我们对弛豫过程对每一步中的总能量和力作图,可以看到总能量或多或少是在持续下降的,这符合我们采用的BFGS能量最小化方法的期望。BFGS的初衷是为了能量的最小化,而不是仅仅跟随力的方向变化。下面的小脚本能够读取轨迹数据,并为每一步的总能量和力最大值作图: | 如果我们对弛豫过程对每一步中的总能量和力作图,可以看到总能量或多或少是在持续下降的,这符合我们采用的BFGS能量最小化方法的期望。BFGS的初衷是为了能量的最小化,而不是仅仅跟随力的方向变化。下面的小脚本能够读取轨迹数据,并为每一步的总能量和力最大值作图: | ||
行 131: | 行 126: | ||
有一点非常重要:不对称的出现,是由长程作用导致的,这也是短程方法,比如tight-binding或者经典力场只能弛豫到对称二聚体到原因。 | 有一点非常重要:不对称的出现,是由长程作用导致的,这也是短程方法,比如tight-binding或者经典力场只能弛豫到对称二聚体到原因。 | ||
=====总结===== | =====总结===== | ||
+ | 本例演示了ATK能够通过DFT结构优化,成功地预测相当复杂的表面结构重构。Si(100)面首先形成一个对称二聚体,它反过来诱导出电子层面的不对称,从而我们在再次优化结构之后,找到基态不对称结构。这种不对称是一种长程效应,来自表面以下的少量几层原子。二聚体Si晶格在二聚体下面被向下压缩。 |