跳到主要内容
版本:2024.05-cpu

NEB 计算过渡态

整个计算流程有以下几步:

Note

step1: 得到稳定的初末态结构;

step2: 自洽计算得到初态能量与末态能量;

step3: 过渡态计算。

NOTE: 初末态计算的参数设置要完全一致

以 Li 在 LiF 中的间隙运输为例子。建模准备 LiF 原胞,分别进行晶格优化+原子位置优化,如何进行优化参考结构优化部分。 优化完成后点击结构优化组件查看结果,并导出结构到结构库中。

从结构库中找到刚才保存的结构,并点击编辑进行结构操作:因为之前优化时使用的是原胞,需要寻找对称性变成惯用晶胞再进行扩胞。

然后为 LiF 构建超胞及 Li 间隙的初末态结构。

初末态结构准备后,再次进入 Q-Flow 设置工作流分别进行结构优化及自洽计算,此时只需要进行原子位置优化。即点击组件系统组件→ 拖拽结构优化自洽计算组件到工作流区域 → 添加从结构导入结构优化再到自洽计算的连线。同样在优化完成后点击结构优化组件查看结果,并导出结构到结构库中。

最后是过渡态的计算:点击组件系统组件→ 拖拽过渡态计算组件到工作流区域 → 添加从结构导入过渡态计算的连线。其中结构导入组件中导入的是初态优化后的结构。

过渡态计算可以设置的基础参数除了有常规的 GPU 卡数、是否 k 点并行、赝势类型、泛函类型、截断能、k 点网格、原子优化算法、控制弛豫的步数及收敛条件等参数。还有一些对 NEB 计算的收敛标准、方法设置等。并在该组件中导入优化好的末态结构。其中初末态总能从前面自洽计算组件的结果中获取。

计算进度栏内包含一些图表数据,表格内会给出体系的基本参数。图里给出了平均总能变化、原子受力变化、各 image 总能、各 image 间距及自洽步数等。

计算完成后点击过渡态计算组件查看结果,可以查看过渡态势垒。