|
发布时间: 2025-4-2 23:54
正文摘要:高斯输出文件中只有transition velocity quadrupole moments，为什么没有transition quadrupole moments。怎么才能得到transition quadrupole moments呢，或者哪些量化软件可以计算。 |
wzkchem5 发表于 2025-4-6 18:32 好的，非常感谢。您说的基组足够大，相互作用可以用transition velocity quadrupole moments计算吗。我再去找找相关的文献 |
cdw32156 发表于 2025-4-6 15:01 当基组足够大时就可以，至于多大算大取决于你需要的精度，也可参考文献确定，transition velocity dipole moments的结论可以大体适用于transition velocity quadrupole moments的情形 你的体系Qzz=0是因为分子对称性 |
wzkchem5 发表于 2025-4-6 13:30 是的，与屏幕平行。好的,感谢。我现在卡在电四极矩跃迁的计算上了。计算只能得到的transition velocity quadrupole moments。可是相互作用可以用transition velocity quadrupole moments计算吗，我不太懂这个。还有四极矩满足Qxx+Qyy+Qzz=0。我这个体系的这个激发态的transition velocity quadrupole moments的Qxx+Qyy=0，Qzz=0。请问这是合理的吗? |
cdw32156 发表于 2025-4-6 13:13 XY面指的是和屏幕平行的面吗？如果是的话，可以认为是主要由电四极-电四极相互作用贡献的相互作用能 |
sobereva 发表于 2025-4-4 03:30 感谢sob老师，我会尝试学习一下的，非常感谢 |
cdw32156 发表于 2025-4-3 15:54 改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。 先按楼上说的想清楚算这个有没有意义再去改。 |
| 长度表象的transition quadrupole moments不满足gauge invariance，也就是你把分子平移一下，结果就变了，所以很少有人需要用。这个和transition dipole moment的情况是不一样的，后者不管是长度表象还是速度表象，都满足gauge invariance，但速度表象的结果对基组更敏感，所以大家普遍用长度表象。 |
本帖最后由 cdw32156 于 2025-4-3 15:59 编辑 sobereva 发表于 2025-4-3 00:13 感谢！就是想得到基态到激发态的跃迁四极矩。XX、YY、ZZ、、、方向的数值。但是不知道怎么去实现 |
| 如果是要诸如<psi_A|xy|psi_B>、<psi_A|zz|psi_B>这种积分，其中psi_A和psi_B是TDDFT激发态或其一是基态，改Multiwfn源代码可以容易地实现 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-13 22:29 , Processed in 0.985778 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites