这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:三维周期性体系的能量分解peda [2016/03/24 13:33] – [第四步:结果查看] liu.jun | adf:三维周期性体系的能量分解peda [2016/05/18 15:57] (当前版本) – 移除 liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ======三维周期性体系的能量分解pEDA====== | ||
- | ADF2016以后的版本,支持周期性体系的能量分解。以下以体心立方的铁晶体孔状结构吸附H< | ||
- | |||
- | =====第一步:建模===== | ||
- | |||
- | 在[[|http:// | ||
- | |||
- | 创建孔洞结构的过程,参考[[adf: | ||
- | |||
- | 点击如下图黄色圆圈所示的H符号,点击孔洞内的空间两次,则创建好了H< | ||
- | |||
- | {{: | ||
- | |||
- | 当然,真正进行计算的时候,需要进行结构优化。但本例中仅演示能量分解的功能,因此省去了结构优化的步骤。 | ||
- | |||
- | ======第二步:设置参数====== | ||
- | |||
- | 设置参数如下: | ||
- | |||
- | 泛函、基组的选择,可以参考[[adf: | ||
- | |||
- | {{: | ||
- | |||
- | **为能量分解设置片段区域,设置方法,如果没有做过的话,务必参考[[adf: | ||
- | |||
- | {{: | ||
- | |||
- | {{: | ||
- | |||
- | **对于pEDA而言,只需要关心Gamma点就可以了**,所以k点设置如下(或者K-Space选择Gamma Only也可以。实际上默认情况下,pEDA也会自动取gamma点,而不读取非1 1 1的k点设置): | ||
- | |||
- | {{: | ||
- | |||
- | {{: | ||
- | |||
- | ======第三步:运行任务====== | ||
- | |||
- | 提交任务的方式,参考[[adf: | ||
- | |||
- | ======第四步:结果查看====== | ||
- | |||
- | 查看结果文件,参考[[adf: | ||
- | |||
- | 打开*.out文件,点击Properties > PEDA即可看到能量分解的情况。PEDA的分解方式和分子中的分解方式很类似。分子的能量分解物理意义,可以参考[[adf: | ||
- | |||
- | △E< | ||
- | |||
- | |||
- | 其中: | ||
- | |||
- | * △E< | ||
- | * △E< | ||
- | * △E< | ||
- | * △E< | ||
- | *** △E< | ||
- | |||
- | 因此用BAND的*.out的书写方式,重新写一遍上面的公式(本例省去E_disp部分): | ||
- | |||
- | E_int=E_elstat+E_Pauli+E_orb | ||
- | |||
- | 上面E_orb,细分为T_orb、Elst_orb、XC_orb,也就是轨道相互作用,对动能T、静电作用能Elst、交换相关能的贡献; | ||
- | |||
- | 与轨道无关的部分有两种成分:E_elstat、E_Pauli是与轨道相互作用无关的部分。注意静电作用能既包括与轨道作用有关的部分,也包括与轨道作用无关的部分; | ||
- | |||
- | 与轨道无关的部分也可以按另一种方式拆分:T^0、Elst^0、XC^0,即动能、静电作用能、交换相关能。注意交换相关能中也是既包括与轨道作用有关的部分,也包括与轨道作用无关的部分; | ||
- | |||
- | 因此上面的公式,可以用其他形式展开: | ||
- | |||
- | E_int | ||
- | |||
- | =(E_elstat+E_Pauli)+E_orb | ||
- | |||
- | =(E_elstat+E_Pauli)+(T_orb+Elst_orb+XC_orb) | ||
- | |||
- | =(T^0+Elst^0+XC^0)+(T_orb+Elst_orb+XC_orb) |