计算化学公社

标题: 使用cp2k计算TDDFT的Transition dipole为0 [打印本页]

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

作者
Author: pwzhou 时间: 2025-9-15 11:25
不考虑悬轨耦合的情况下，单重态到三重态的跃迁是禁阻的，跃迁偶极矩为0是完全正常的结果。

作者
Author: sobereva 时间: 2025-9-16 04:11
你读的是考虑SOC之前的，考虑SOC后的结果在这后面，SOC CORRECTED SPECTRUM

作者
Author: wangxc 时间: 2025-9-18 10:22

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

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

作者
Author: sobereva 时间: 2025-9-18 23:38

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

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

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