这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
atk:块体金的相对论效应 [2016/08/08 14:33] – [旋轨耦合相对论GGA能带计算] liu.jun | atk:块体金的相对论效应 [2022/11/18 11:32] (当前版本) – [旋轨耦合相对论GGA能带计算] liu.jun | ||
---|---|---|---|
行 3: | 行 3: | ||
脚本文件下载: | 脚本文件下载: | ||
- | * [[docs.quantumwise.com/ | + | * [[http://docs.quantumwise.com/ |
- | * [[docs.quantumwise.com/ | + | * [[http://docs.quantumwise.com/ |
- | * [[docs.quantumwise.com/ | + | * [[http://docs.quantumwise.com/ |
金是较重的元素之一。体相的金,其电子结构受到相对论效应,包括自旋轨道耦合(SOC)的显著影响。 | 金是较重的元素之一。体相的金,其电子结构受到相对论效应,包括自旋轨道耦合(SOC)的显著影响。 | ||
行 19: | 行 19: | ||
首先,双击添加一个New Calculator{{: | 首先,双击添加一个New Calculator{{: | ||
- | {{: | + | {{ : |
**编辑New | **编辑New | ||
行 25: | 行 25: | ||
在Basis set/ | 在Basis set/ | ||
- | {{: | + | {{ : |
在Basis菜单选择Density mesh-cutoff到300Hartree,k-point采样增加到9x9x9: | 在Basis菜单选择Density mesh-cutoff到300Hartree,k-point采样增加到9x9x9: | ||
- | {{: | + | {{ : |
**编辑Bandstructure设置:** | **编辑Bandstructure设置:** | ||
行 35: | 行 35: | ||
设置Points per segment为101,这样高对称点之间的能带很光滑;Brillouin zone route设置为G, | 设置Points per segment为101,这样高对称点之间的能带很光滑;Brillouin zone route设置为G, | ||
- | {{: | + | {{ : |
保存脚本的名为gold_gga.py,通过Job Manager{{: | 保存脚本的名为gold_gga.py,通过Job Manager{{: | ||
行 43: | 行 43: | ||
在VNL LabFloor中找到输出文件的Bandstructure项,并使用Bandstructure Analyzer作图。结果如下图所示: | 在VNL LabFloor中找到输出文件的Bandstructure项,并使用Bandstructure Analyzer作图。结果如下图所示: | ||
- | {{: | + | {{ : |
与文献< | 与文献< | ||
=====旋轨耦合相对论GGA能带计算===== | =====旋轨耦合相对论GGA能带计算===== | ||
- | 下面演示使用GGA泛函计算包含自旋轨道耦合的能带结构。电子的自旋必须设置为非共线。这样比标准的GGA计算耗时更大,并且不容易收敛。可以使用两个方法来解决:a) 从自旋激化的计算结果(SGGA)开始计算;b)使用电子密度混合的方法。 | + | 下面演示使用GGA泛函计算包含自旋轨道耦合的能带结构。电子的自旋必须设置为非共线。这样比标准的GGA计算耗时更大,并且不容易收敛。可以使用两个方法来解决:a) 从自旋极化的计算结果(SGGA)开始计算;b)使用电子密度混合的方法。 |
<WRAP center round box 100%> | <WRAP center round box 100%> | ||
注意: | 注意: | ||
- | 我们知道金是没有磁性的,因此对金进行自旋激化的计算实际上看起来有些傻瓜。这样得到的结果和共线、自旋补偿计算的到的结果是一致的。 | + | 我们知道金是没有磁性的,因此对金进行自旋极化的计算实际上看起来有些傻瓜。这样得到的结果和共线、自旋补偿计算的到的结果是一致的。 |
- | 实际上,和自旋补偿态的结果相比,收敛的自旋激化的计算结果,能够给非共线旋轨耦合态提供一个好得多的初始猜测。如果从自旋激化电子态出发,非共线的自洽迭代只需要少量几步就能收敛,从而能够大大节省计算时间。 | + | 实际上,和自旋补偿态的结果相比,收敛的自旋极化的计算结果,能够给非共线旋轨耦合态提供一个好得多的初始猜测。如果从自旋极化电子态出发,非共线的自洽迭代只需要少量几步就能收敛,从而能够大大节省计算时间。 |
</ | </ | ||
**SGGA初始态:** | **SGGA初始态:** | ||
行 178: | 行 178: | ||
</ | </ | ||
- | 保存该脚本名为gold_sogga.py,并运行计算。使用4核进行并行化计算,大约需要6分钟。注意ATK并行计算需要使用MPI库,具体参考[[atk: | + | 保存该脚本名为gold_sogga.py,并运行计算。使用4核进行并行化计算,大约需要6分钟。注意QuantumATK并行计算需要使用MPI库,具体参考[[atk: |
结果: | 结果: | ||
行 184: | 行 184: | ||
下图为SOGGA计算得到的能带结构。与文献< | 下图为SOGGA计算得到的能带结构。与文献< | ||
- | {{: | + | {{ : |
你可以使用Compare Data插件更直接的对比这两个能带(如下图所示)。同时高亮显示GGA能带和SOGGA能带,可以看到相对论效应对金的价带有显著影响。下图与文献< | 你可以使用Compare Data插件更直接的对比这两个能带(如下图所示)。同时高亮显示GGA能带和SOGGA能带,可以看到相对论效应对金的价带有显著影响。下图与文献< | ||
- | {{: | + | {{ : |