计算化学公社

标题: 怎么才能得到transition quadrupole moments呢 [打印本页]

作者
Author: cdw32156 时间: 2025-4-2 23:54
标题: 怎么才能得到transition quadrupole moments呢
高斯输出文件中只有transition velocity quadrupole moments，为什么没有transition quadrupole moments。怎么才能得到transition quadrupole moments呢，或者哪些量化软件可以计算。

作者
Author: sobereva 时间: 2025-4-3 00:13
如果是要诸如<psi_A|xy|psi_B>、<psi_A|zz|psi_B>这种积分，其中psi_A和psi_B是TDDFT激发态或其一是基态，改Multiwfn源代码可以容易地实现

作者
Author: cdw32156 时间: 2025-4-3 15:54
本帖最后由 cdw32156 于 2025-4-3 15:59 编辑

sobereva 发表于 2025-4-3 00:13
如果是要诸如、这种积分，其中psi_A和psi_B是TDDFT激发态或其一是基态，改Multiwfn源代码可以容易地实现

感谢！就是想得到基态到激发态的跃迁四极矩。XX、YY、ZZ、、、方向的数值。但是不知道怎么去实现

作者
Author: wzkchem5 时间: 2025-4-3 16:20
长度表象的transition quadrupole moments不满足gauge invariance，也就是你把分子平移一下，结果就变了，所以很少有人需要用。这个和transition dipole moment的情况是不一样的，后者不管是长度表象还是速度表象，都满足gauge invariance，但速度表象的结果对基组更敏感，所以大家普遍用长度表象。

作者
Author: sobereva 时间: 2025-4-4 03:30

cdw32156 发表于 2025-4-3 15:54
感谢！就是想得到基态到激发态的跃迁四极矩。XX、YY、ZZ、、、方向的数值。但是不知道怎么去实现

改excittrans.f90里的subroutine exctransdip，对应Multiwfn主功能18的子功能5，原本是用来计算态之间的跃迁电/磁偶极矩的。里面有call genDbas_curr和call genDorb，这两个子程序分别用来产生基函数和轨道之间的电偶极矩算符的积分矩阵，可以参照integral.f90里其源代码新写个比如genQbas_curr、genQorb子程序用来产生四极矩积分矩阵，用它们分别替换掉exctransdip里对genDbas_curr和genDorb的调用。genDbas_curr和genDorb中都利用到genDbas子程序，里面的产生GTF之间的偶极矩积分的子程序dodipoleint可直接改为现成的产生GTF之间多极矩积分的子程序domultipoleint。

先按楼上说的想清楚算这个有没有意义再去改。

作者
Author: cdw32156 时间: 2025-4-6 13:13

wzkchem5 发表于 2025-4-3 16:20
长度表象的transition quadrupole moments不满足gauge invariance，也就是你把分子平移一下，结果就变了， ...

感谢老师，我的理解太浅显了，我是想计算一下由单分子到二聚体，属于四极矩跃迁的激发态分裂为两个激发态的情况。
(, 下载次数 Times of downloads: 5) 比如这个D4h对称性的分子，这是第四激发态的跃迁密度。具有两个与X轴呈45°的相位相反的跃迁偶极，总的跃迁偶极为0。当两个同样的分子放于XY面，这个激发态会耦合分裂为两个激发态。其相互作用能是由什么决定的呀。对于跃迁偶极来说，利用激子耦合模型，知道跃迁偶极矩的大小就能计算。但是对于这种类型，也可以通过跃迁偶极矩计算呢，还是要利用跃迁四极矩呀，我先用公式拟合一下

作者
Author: cdw32156 时间: 2025-4-6 13:13

sobereva 发表于 2025-4-4 03:30
改excittrans.f90里的subroutine exctransdip，对应Multiwfn主功能18的子功能5，原本是用来计算态之间的 ...

感谢sob老师，我会尝试学习一下的，非常感谢

作者
Author: wzkchem5 时间: 2025-4-6 13:30

cdw32156 发表于 2025-4-6 13:13
感谢老师，我的理解太浅显了，我是想计算一下由单分子到二聚体，属于四极矩跃迁的激发态分裂为两个激发态 ...

XY面指的是和屏幕平行的面吗？如果是的话，可以认为是主要由电四极-电四极相互作用贡献的相互作用能

作者
Author: cdw32156 时间: 2025-4-6 15:01

wzkchem5 发表于 2025-4-6 13:30
XY面指的是和屏幕平行的面吗？如果是的话，可以认为是主要由电四极-电四极相互作用贡献的相互作用能

是的，与屏幕平行。好的,感谢。我现在卡在电四极矩跃迁的计算上了。计算只能得到的transition velocity quadrupole moments。可是相互作用可以用transition velocity quadrupole moments计算吗，我不太懂这个。还有四极矩满足Qxx+Qyy+Qzz=0。我这个体系的这个激发态的transition velocity quadrupole moments的Qxx+Qyy=0，Qzz=0。请问这是合理的吗?

作者
Author: wzkchem5 时间: 2025-4-6 18:32

cdw32156 发表于 2025-4-6 15:01
是的，与屏幕平行。好的,感谢。我现在卡在电四极矩跃迁的计算上了。计算只能得到的transition velocity q ...

当基组足够大时就可以，至于多大算大取决于你需要的精度，也可参考文献确定，transition velocity dipole moments的结论可以大体适用于transition velocity quadrupole moments的情形
你的体系Qzz=0是因为分子对称性

作者
Author: cdw32156 时间: 2025-4-7 10:42

wzkchem5 发表于 2025-4-6 18:32
当基组足够大时就可以，至于多大算大取决于你需要的精度，也可参考文献确定，transition velocity dipole ...

好的，非常感谢。您说的基组足够大，相互作用可以用transition velocity quadrupole moments计算吗。我再去找找相关的文献

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