Shana 发表于 2024-4-25 09:09
老师您好。我找到一个TBSOC软件，能够计算SOC矩阵元。这个是基于瓦尼尔轨道的。老师，我想请问一下，如果 ...

这个其实是基于紧束缚模型算的，要求基函数必须是px/py/pz这样的原子轨道。Wannier90给出的最大居域化函数并不是原子轨道。

万里云 发表于 2024-4-22 12:52
1. VASP中的自旋轨道耦合应该是lambda*L*S形式的。自旋轨道耦合是电子由于运动感受到的磁场，与自身自旋磁 ...

老师您好。我找到一个TBSOC软件，能够计算SOC矩阵元。这个是基于瓦尼尔轨道的。老师，我想请问一下，如果我想计算两个orbital之间的SOC，使用TBSOC可行吗？

万里云 发表于 2024-4-22 12:52
1. VASP中的自旋轨道耦合应该是lambda*L*S形式的。自旋轨道耦合是电子由于运动感受到的磁场，与自身自旋磁 ...

谢谢您的回复。我在大师兄群里看到跟您同名的，是同一个人吗？方便加一下qq吗？

1. VASP中的自旋轨道耦合应该是lambda*L*S形式的。自旋轨道耦合是电子由于运动感受到的磁场，与自身自旋磁矩间耦合导致的能量附加项，是一个类似势能项的单体算符，应该可以求出矩阵元后直接加到总能里。
2. 要计算lambda*L*S形式的自旋轨道耦合项，必须先把KS轨道投影到原子波函数上。所以OUTCAR中是按原子给出的。
3. 展开系数一般情况下是复数，PROCAR中给出的应该是其模方，不太好算。材料不复杂的话，可以用紧束缚模型来算。

PS：L/S型的自旋轨道耦合中的sx和sy项会使自旋期望值在x和y方向也有分量，形成自旋纹理。但这个纹理远比Rashba或者Dresselhaus型的复杂。而且后两者对体系对称性有要求。很多文章中用VASP计算上下表面对称材料的自旋纹理，还声称有Rashba效应，就不知道是怎么算出来的。