求助，DFT计算自旋轨道耦合情形的跃迁偶极矩遇到困难

Shana · 发表于 Post on 2024-4-17 10:24:42

各位老师好，

小弟最近在学习一篇文关于计算自旋轨道耦合下跃迁偶极矩的文献，参考这篇npj的文献，在一阶微扰下，可以将跃迁偶极矩写成，
其中1A是单态，3A是三重态。小弟对于这种计算方法有几处疑问，恳请各位老师解答。

1、Δμ应该怎么计算？<3A|er|3A>应该怎么计算？这是否可以理解为r的平均值乘上电荷？那么又应该怎么计算这个平均值呢？

2、自旋轨道耦合部分应该怎么计算？在文献中有这么一段描述，这里的矩阵元在vasp中是怎么计算的呢？

3、在文献附录中，将自旋轨道耦合拆分成如下的形式，这里面的推导过程我能够理解，但是最后应该计算各个量呢？例如，<a''↑|lz·sz|a''↑>的值怎么计算？小弟设想的是LS=J^2-L^2-S^2，然后假定微扰前的态是J^2,L^2和S^2的本征态。但是作者将其分成了单电子态的处理，那么这里的值应该怎么算呢？

恳请各位老师和同学指导！

Shana · 发表于 Post on 2024-4-17 15:13:14

@卡开发发卡老师有时间指导一下吗？

卡开发发 · 发表于 Post on 2024-4-17 17:38:34

Shana 发表于 2024-4-17 15:13
@卡开发发卡老师有时间指导一下吗？

旋轨耦合这方面我了解有限，建议看看别人的说法。对于问题1，郑奇靖老师写了一个vasptdm的程序，你可以参考他写的原理和程序。

Shana · 发表于 Post on 2024-4-17 19:38:05

卡开发发发表于 2024-4-17 17:38
旋轨耦合这方面我了解有限，建议看看别人的说法。对于问题1，郑奇靖老师写了一个vasptdm的程序，你可以参 ...

感谢卡老师！打扰了

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

[综合讨论] 求助，DFT计算自旋轨道耦合情形的跃迁偶极矩遇到困难