自旋轨道耦合（SOC）时，能带计算疑惑

yeahhanpei · 发表于 Post on 2019-8-27 22:16:45

最近看了篇拓扑绝缘体的模拟文章，其中提到了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设置有需要注意的吗，还是和常规的一样即可？

恳请老师指教，谢谢

Light · 发表于 Post on 2019-9-21 19:25:55

麻烦楼主留个联系方式，一起学习

Anni · 发表于 Post on 2019-10-28 19:03:36

您好，我想请教一下您这个考虑SOC的图是用什么软件怎么画的呢？

Anni · 发表于 Post on 2019-10-29 18:38:06

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

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

yeahhanpei · 发表于 Post on 2019-11-4 11:47:04

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

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

wcwlk22 · 发表于 Post on 2020-10-4 10:18:02

第一个问题：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β，至于取向我就真没办法了；本质上应该要做一个超交换的计算算出来，但那太麻烦了。我自己是根据文章等效出来的一个磁矩方向取的。

qqq21 · 发表于 Post on 2022-4-20 21:34:19

楼主问题解决了吗

weilu · 发表于 Post on 2024-4-8 20:11:49

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

汪杰 · 发表于 Post on 2024-12-2 20:04:21

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

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

saxis的理解好像有点不对

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] 自旋轨道耦合（SOC）时，能带计算疑惑

评分 Rate