计算化学公社

标题: 求助，DFT计算自旋轨道耦合情形的跃迁偶极矩遇到困难 [打印本页]

作者
Author: Shana 时间: 2024-4-17 10:24
标题: 求助，DFT计算自旋轨道耦合情形的跃迁偶极矩遇到困难
各位老师好，

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

1、Δμ应该怎么计算？<3A|er|3A>应该怎么计算？这是否可以理解为r的平均值乘上电荷？那么又应该怎么计算这个平均值呢？
(, 下载次数 Times of downloads: 17)

2、自旋轨道耦合部分应该怎么计算？在文献中有这么一段描述，这里的矩阵元在vasp中是怎么计算的呢？
(, 下载次数 Times of downloads: 14)

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

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

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

作者
Author: 卡开发发 时间: 2024-4-17 17:38

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

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

作者
Author: Shana 时间: 2024-4-17 19:38

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

感谢卡老师！打扰了

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