计算化学公社

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

[量化理论] 求教TD输出文件里跃迁电偶极矩、跃迁速度偶极矩、跃迁磁偶极矩三者的区别

[复制链接 Copy URL]

8

帖子

0

威望

105

eV
积分
113

Level 2 能力者

求教TD输出文件里跃迁电偶极矩、跃迁速度偶极矩、跃迁磁偶极矩三者的区别以及两个振子强度的区别

Image_1713454402999_edit_290975888779562.png (77.13 KB, 下载次数 Times of downloads: 28)

Image_1713454402999_edit_290975888779562.png

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

2#
发表于 Post on 2024-4-19 03:02:03 | 只看该作者 Only view this author
相当于<G|t|E>积分,G和E分别是基态和激发态波函数,三个跃迁矩对应的算符t不一样。不说具体问这个的目的没法进一步回答你
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

285

帖子

2

威望

2852

eV
积分
3177

Level 5 (御坂)

3#
发表于 Post on 2024-5-5 18:07:02 | 只看该作者 Only view this author
如何得到S0到T1跃迁的跃迁电偶极矩和跃迁磁偶极矩?用高斯可以吗?

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

4#
发表于 Post on 2024-5-5 18:28:26 | 只看该作者 Only view this author
虎王 发表于 2024-5-5 11:07
如何得到S0到T1跃迁的跃迁电偶极矩和跃迁磁偶极矩?用高斯可以吗?

需要用支持在TDDFT里考虑SOC贡献的软件。高斯本身不行,高斯搭配其他能算SOC的程序可能行,但是数据后处理可能比较麻烦(需要用高斯算的标量态跃迁偶极矩和其他程序算的SOC来计算旋量态跃迁偶极矩,也就是程序算完以后很可能还需要一些手工计算)。
建议用ORCA、BDF等软件算
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

285

帖子

2

威望

2852

eV
积分
3177

Level 5 (御坂)

5#
发表于 Post on 2024-5-5 20:22:49 | 只看该作者 Only view this author
wzkchem5 发表于 2024-5-5 18:28
需要用支持在TDDFT里考虑SOC贡献的软件。高斯本身不行,高斯搭配其他能算SOC的程序可能行,但是数据后处 ...

ORCA似乎不提供跃迁磁偶极矩。

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

6#
发表于 Post on 2024-5-5 20:26:42 | 只看该作者 Only view this author
虎王 发表于 2024-5-5 13:22
ORCA似乎不提供跃迁磁偶极矩。

即使程序不打印跃迁磁偶极矩,由CD强度和跃迁电偶极矩,应该能反推出跃迁磁偶极矩
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

5

帖子

0

威望

139

eV
积分
144

Level 2 能力者

7#
发表于 Post on 2024-5-6 15:14:48 | 只看该作者 Only view this author
sobereva 发表于 2024-4-19 03:02
相当于积分,G和E分别是基态和激发态波函数,三个跃迁矩对应的算符t不一样。不说具体问这个的目的没法进一 ...

老师 我想通过振子强度去相对说明从基态跃迁到某个激发态的概率的大小 那我应该用跃迁偶极矩还是跃迁速度偶极矩上的振子强度呢

195

帖子

0

威望

4029

eV
积分
4224

Level 6 (一方通行)

8#
发表于 Post on 2024-5-7 11:05:44 | 只看该作者 Only view this author
wzkchem5 发表于 2024-5-5 20:26
即使程序不打印跃迁磁偶极矩,由CD强度和跃迁电偶极矩,应该能反推出跃迁磁偶极矩

老师,请问下通过强度算出来的磁偶极矩是不是只有数值,不是向量;请问如何获得跃迁磁偶极矩向量,dalton或者orca可以吗,谢谢老师

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

9#
发表于 Post on 2024-5-7 15:29:17 | 只看该作者 Only view this author
tiandikuoyuan 发表于 2024-5-7 04:05
老师,请问下通过强度算出来的磁偶极矩是不是只有数值,不是向量;请问如何获得跃迁磁偶极矩向量,dalton ...

ORCA给出的CD强度是分x,y,z分量的,所以应该可以得到跃迁磁偶极矩的向量

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
tiandikuoyuan + 5 谢谢

查看全部评分 View all ratings

BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

10#
发表于 Post on 2024-5-7 15:31:15 | 只看该作者 Only view this author
ppzhang 发表于 2024-5-6 08:14
老师 我想通过振子强度去相对说明从基态跃迁到某个激发态的概率的大小 那我应该用跃迁偶极矩还是跃迁速度 ...

当基组完备时,两种做法等价;当基组并非完备时,跃迁偶极矩得到的振子强度往往更准
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

195

帖子

0

威望

4029

eV
积分
4224

Level 6 (一方通行)

11#
发表于 Post on 2024-5-7 17:22:15 | 只看该作者 Only view this author
wzkchem5 发表于 2024-5-7 15:29
ORCA给出的CD强度是分x,y,z分量的,所以应该可以得到跃迁磁偶极矩的向量

老师,我查了下ECD强度的计算公式(J. Chem. Phys. 150, 174121 (2019)),我想问下orca计算结果CD SPECTRUM中输出的“MX        MY        MZ”是不是就是公式中m跃迁磁偶极矩?另外输出的SOC CORRECTED CD SPECTRUM* 中给出的x,y,z分量是复数形式,请问应该怎么转换为实数向量。非常感谢老师!


1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

12#
发表于 Post on 2024-5-7 19:10:03 | 只看该作者 Only view this author
本帖最后由 wzkchem5 于 2024-5-7 12:26 编辑
tiandikuoyuan 发表于 2024-5-7 10:22
老师,我查了下ECD强度的计算公式(J. Chem. Phys. 150, 174121 (2019)),我想问下orca计算结果CD SPECT ...

我看了一下代码,确实是跃迁磁偶极矩。这一点手册没有写清楚,我们会补上。抱歉这一点我之前的帖子也说错了,我想当然以为是转子强度的x,y,z分量了,因为很早的ORCA版本里这里写的是RX, RY, RZ(但我问了一下负责这块代码的人,即使在当时,这里实际打印的仍然是跃迁磁偶极矩)。
SOC校正的CD谱的跃迁磁偶极矩本来就应该是复数,谈不上转化为实数。如果直接报道,报道复数就行;如果要做计算,用复数乘法的规则套实数的公式就行,但是该取复共轭的地方不要忘了取复共轭;如果要作图表示跃迁磁偶极矩的方向,我不太清楚怎么画,但估计肯定有文献画过,可以查查文献是怎么画的

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
tiandikuoyuan + 5 谢谢

查看全部评分 View all ratings

BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

195

帖子

0

威望

4029

eV
积分
4224

Level 6 (一方通行)

13#
发表于 Post on 2024-9-3 15:02:10 | 只看该作者 Only view this author
wzkchem5 发表于 2024-5-7 19:10
我看了一下代码,确实是跃迁磁偶极矩。这一点手册没有写清楚,我们会补上。抱歉这一点我之前的帖子也说错 ...

老师,我尝试使用orca6计算了SOC校正的CD强度,发现给出的形式和orca5有很大变化,目前给出了跃迁电偶极矩和跃迁速度偶极矩形式的CD强度,请问下这边的MXMYMZ是否是跃迁磁偶极矩?
“SOC CORRECTED ABSORPTION SPECTRUM VIA TRANSITION ELECTRIC DIPOLE MOMENTS”和“SOC CORRECTED ABSORPTION SPECTRUM VIA TRANSITION VELOCITY DIPOLE MOMENTS ”中的|DX|和|PX|能否可以认为分别是相应激发态SOC校正后的跃迁电偶极矩和跃迁速度偶极矩?
请问下目前给出的偶极矩数据都是复数的‌模,应该如何知道偶极矩方向?两个版本计算的CD强度有很大区别,请问是软件内换了方法吗,计算文件中只有“DoQuad true”不同,因为在orca6中该关键词被替换了。非常感谢!


orca5计算结果

orca6计算结果



1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

14#
发表于 Post on 2024-9-7 15:13:46 | 只看该作者 Only view this author
tiandikuoyuan 发表于 2024-9-3 08:02
老师,我尝试使用orca6计算了SOC校正的CD强度,发现给出的形式和orca5有很大变化,目前给出了跃迁电偶极 ...



当偶极矩是复数的时候,原则上没法谈论偶极矩的方向,这是考虑了SOC的计算和没考虑SOC的计算的一个很大区别。
去掉DoQuad true再试试,如果还是不一样,建议在orca论坛上发帖问

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
tiandikuoyuan + 5 谢谢

查看全部评分 View all ratings

BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 06:00 , Processed in 0.214492 second(s), 35 queries , Gzip On.

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