计算化学公社

标题: 自旋轨道耦合（SOC）时，能带计算疑惑 [打印本页]

作者
Author: yeahhanpei 时间: 2019-8-27 22:16
标题: 自旋轨道耦合（SOC）时，能带计算疑惑
最近看了篇拓扑绝缘体的模拟文章，其中提到了GeTe/Sb2Te3的自旋轨道耦合的能带结构，打算自己做一遍。查询资料后得出计算步骤如下：

第一步：进行常规的结构优化
第二步：提取上步的结构文件，INCAR中加上SOC，进行静态计算以得到CHGCAR和WAVCAR。添加的参数为
         LOSRBIT=.T.
         GGA_COMPAT=.F.
         NBANDS= ***             !上一步OUTCAR中的两倍
         LMAXMIX=4
         ISYM=0
         SAXIS= ? ? ?
         MAGMOM= ? ? ?
第三步：设置ICHARG=11读取上一步的CHGCAR，添加常规能带计算时的参数设置，上一步列出的参数在这一步同样添加，再进行能带计算即可。

但是还有些疑惑：
1）以上步骤中又不合理的地方吗
2）第一步中INCAR需要添加 ISYM=0；ISPIN=2吗？
3）SAXIS参数网上说通常设为 0 0 1，这对所有体系都合适吗？
4）MAGMOM是体系中所有元素的磁矩，元素的磁矩有什么标准吗，在 x y z三轴如何分？
5）第二、第三步的KPOINTS设置有需要注意的吗，还是和常规的一样即可？

恳请老师指教，谢谢

作者
Author: Light 时间: 2019-9-21 19:25
麻烦楼主留个联系方式，一起学习

作者
Author: Anni 时间: 2019-10-28 19:03
您好，我想请教一下您这个考虑SOC的图是用什么软件怎么画的呢？

作者
Author: Anni 时间: 2019-10-29 18:38

Anni 发表于 2019-10-28 19:03
您好，我想请教一下您这个考虑SOC的图是用什么软件怎么画的呢？

您好，我刚刚看到您写错了，是LSORBIT，ISPIN=1，ISYM=0要添加，

作者
Author: yeahhanpei 时间: 2019-11-4 11:47

Anni 发表于 2019-10-29 18:38
您好，我刚刚看到您写错了，是LSORBIT，ISPIN=1，ISYM=0要添加，

嗯嗯，多谢多谢。帖子中的图是文献里的

作者
Author: wcwlk22 时间: 2020-10-4 10:18
第一个问题：LMAXMIX=4 应该是d电子体系的，但我看了你这个主要是5p..

2 saxis = 0 0 1 意思是磁矩平行z轴, 我的理解是这是一个假设平均场近似的设置. 比如你有一个反铁磁, 你把所有磁矩等效到z上，然后看成z+和z-的叠加，这个应该是可以的。
但是在SOC计算里我认为这是不合理的。
因为你p电子原子轨道磁矩至少分px py pz，原本你的体系里自旋可能跟px和py耦合的，你全部align到z方向就变成所有电子和pz耦合了，耦合强度仅取决于各自自旋在z方向上的分量。这是不合理的。

最后，magmon的单位是玻尔磁子，你的文章里应该有这个体系的ueff=uβ，至于取向我就真没办法了；本质上应该要做一个超交换的计算算出来，但那太麻烦了。我自己是根据文章等效出来的一个磁矩方向取的。

作者
Author: qqq21 时间: 2022-4-20 21:34
楼主问题解决了吗

作者
Author: weilu 时间: 2024-4-8 20:11
新手刚入，像是CrI3求能带的话，需要考虑磁性吗，考虑的话结构优化需要考虑soc磁性这些吗，那SAXIS= MAGMOM= （扩胞24个I 8个Cr）如何设置呢，静态计算需要接着考虑这些吗，然后第三步再进行常规的能带计算吗，谢谢各位指教

作者
Author: 汪杰 时间: 2024-12-2 20:04

wcwlk22 发表于 2020-10-4 10:18
第一个问题：LMAXMIX=4 应该是d电子体系的，但我看了你这个主要是5p..

2 saxis = 0 0 1 意思是磁矩平行z ...

saxis的理解好像有点不对

欢迎光临计算化学公社 (http://bbs.keinsci.com/)