计算化学公社

 找回密码 Forget password
 注册 Register
Views: 489|回复 Reply: 10
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 怎么才能得到transition quadrupole moments呢

[复制链接 Copy URL]

28

帖子

0

威望

453

eV
积分
481

Level 3 能力者

高斯输出文件中只有transition velocity quadrupole moments,为什么没有transition quadrupole moments。怎么才能得到transition quadrupole moments呢,或者哪些量化软件可以计算。

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

2#
发表于 Post on 2025-4-3 00:13:45 | 只看该作者 Only view this author
如果是要诸如<psi_A|xy|psi_B>、<psi_A|zz|psi_B>这种积分,其中psi_A和psi_B是TDDFT激发态或其一是基态,改Multiwfn源代码可以容易地实现
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

28

帖子

0

威望

453

eV
积分
481

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-4-3 15:54:56 | 只看该作者 Only view this author
本帖最后由 cdw32156 于 2025-4-3 15:59 编辑
sobereva 发表于 2025-4-3 00:13
如果是要诸如、这种积分,其中psi_A和psi_B是TDDFT激发态或其一是基态,改Multiwfn源代码可以容易地实现

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

1万

帖子

0

威望

8969

eV
积分
20713

Level 6 (一方通行)

4#
发表于 Post on 2025-4-3 16:20:16 | 只看该作者 Only view this author
长度表象的transition quadrupole moments不满足gauge invariance,也就是你把分子平移一下,结果就变了,所以很少有人需要用。这个和transition dipole moment的情况是不一样的,后者不管是长度表象还是速度表象,都满足gauge invariance,但速度表象的结果对基组更敏感,所以大家普遍用长度表象。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

6万

帖子

99

威望

5万

eV
积分
120081

管理员

公社社长

5#
发表于 Post on 2025-4-4 03:30:58 | 只看该作者 Only view this author
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。

先按楼上说的想清楚算这个有没有意义再去改。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

28

帖子

0

威望

453

eV
积分
481

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-4-6 13:13:00 | 只看该作者 Only view this author
wzkchem5 发表于 2025-4-3 16:20
长度表象的transition quadrupole moments不满足gauge invariance,也就是你把分子平移一下,结果就变了, ...

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

28

帖子

0

威望

453

eV
积分
481

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2025-4-6 13:13:51 | 只看该作者 Only view this author
sobereva 发表于 2025-4-4 03:30
改excittrans.f90里的subroutine exctransdip,对应Multiwfn主功能18的子功能5,原本是用来计算态之间的 ...

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

1万

帖子

0

威望

8969

eV
积分
20713

Level 6 (一方通行)

8#
发表于 Post on 2025-4-6 13:30:25 | 只看该作者 Only view this author
cdw32156 发表于 2025-4-6 13:13
感谢老师,我的理解太浅显了,我是想计算一下由单分子到二聚体,属于四极矩跃迁的激发态分裂为两个激发态 ...

XY面指的是和屏幕平行的面吗?如果是的话,可以认为是主要由电四极-电四极相互作用贡献的相互作用能
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

28

帖子

0

威望

453

eV
积分
481

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2025-4-6 15:01:09 | 只看该作者 Only view this author
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。请问这是合理的吗?

1万

帖子

0

威望

8969

eV
积分
20713

Level 6 (一方通行)

10#
发表于 Post on 2025-4-6 18:32:37 | 只看该作者 Only view this author
cdw32156 发表于 2025-4-6 15:01
是的,与屏幕平行。好的,感谢。我现在卡在电四极矩跃迁的计算上了。计算只能得到的transition velocity q ...

当基组足够大时就可以,至于多大算大取决于你需要的精度,也可参考文献确定,transition velocity dipole moments的结论可以大体适用于transition velocity quadrupole moments的情形
你的体系Qzz=0是因为分子对称性
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

28

帖子

0

威望

453

eV
积分
481

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2025-4-7 10:42:59 | 只看该作者 Only view this author
wzkchem5 发表于 2025-4-6 18:32
当基组足够大时就可以,至于多大算大取决于你需要的精度,也可参考文献确定,transition velocity dipole ...

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-13 05:27 , Processed in 0.172965 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list