使用cp2k计算TDDFT的Transition dipole为0

wangxc · 发表于 Post on 2025-9-15 11:17:14

老师们，我使用cp2k对一个结构进行了TDDFT计算，nstate=5，计算三重态，开启SOC。在输出文件中发现multiplicity 3处的Transition dipole为0。我想请问我的设置是否有问题？以下是我提供的input和log文件，我删除了&coord 的部分。

pwzhou · 发表于 Post on 2025-9-15 11:25:24

不考虑悬轨耦合的情况下，单重态到三重态的跃迁是禁阻的，跃迁偶极矩为0是完全正常的结果。

sobereva · 发表于 Post on 2025-9-16 04:11:32

你读的是考虑SOC之前的，考虑SOC后的结果在这后面，SOC CORRECTED SPECTRUM

wangxc · 发表于 Post on 2025-9-18 10:22:34

sobereva 发表于 2025-9-16 04:11
你读的是考虑SOC之前的，考虑SOC后的结果在这后面，SOC CORRECTED SPECTRUM

感谢老师指导，我想获得数据是T1的Transition dipole，在SOC的部分没有提供，我应该如何获得呢？

sobereva · 发表于 Post on 2025-9-18 23:38:35

wangxc 发表于 2025-9-18 10:22
感谢老师指导，我想获得数据是T1的Transition dipole，在SOC的部分没有提供，我应该如何获得呢？

不输出跃迁偶极矩
如果你的目的只是绘制电子光谱，有振子强度就够了，Multiwfn直接支持基于CP2K的SOC-TDDFT输出文件绘制考虑旋轨耦合的光谱。
如果非要跃迁偶极矩，需要自己改CP2K源代码，从输出振子强度的地方开始顺藤摸瓜找到给出跃迁偶极矩的地方

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

[CP2K] 使用cp2k计算TDDFT的Transition dipole为0