计算化学公社

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

[新手求助] Gaussian算激发态偶极矩和跃迁偶极矩和文献对比问题

[复制链接 Copy URL]

163

帖子

0

威望

1025

eV
积分
1188

Level 4 (黑子)

最近在重复一篇文献(Design, Synthesis, and Reaction of π-Extended Coumarin-based New Caged Compounds with Two-photon Absorption Character in the Near-IR Region)的双光子计算部分,用了相同的机组和泛函,基态的偶极矩是差不多的,但是激发态和跃迁偶极矩差距极大。不明白问题出在哪儿,求大神指点。
以下是我的高斯输入文件:
第一步:基态结构和频率优化
#n B3LYP/6-31G(d) Opt freq
第二步:计算了20个激发态
# B3LYP/6-31G(d) TD=NStates=20
Geom=Check Guess=Read
第三步:从第二步中找到f最大值对应的激发态,做结构和频率优化
# B3LYP/6-31G(d) TD=(Read,NStates=20,Root=1)
Geom=Modify Guess=Read Opt freq


从第一步的log文件中找最后一个 Dipole moment (field-independent basis, Debye),对应的数值作为基态偶极矩
从第三步的log文件中找最后一个 Dipole moment (field-independent basis, Debye),对应的数值作为激发态偶极矩; 找到最后一个Ground to excited state transition electric dipole moments (Au),根据第二步的f最大值对应的激发态,找到该激发态对应的数值作为跃迁偶极矩。




微信图片_20190506151622.png (149.79 KB, 下载次数 Times of downloads: 88)

微信图片_20190506151622.png

微信图片_20190506151514.png (66.77 KB, 下载次数 Times of downloads: 74)

微信图片_20190506151514.png

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

26#
发表于 Post on 2022-5-14 11:25:50 | 只看该作者 Only view this author

看置顶的新社员必读贴了解怎么正确贴图
北京科音自然科学研究中心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

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

25#
发表于 Post on 2022-5-12 15:27:24 | 只看该作者 Only view this author

图片没有上传对
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

2

帖子

0

威望

83

eV
积分
85

Level 2 能力者

24#
发表于 Post on 2022-5-12 15:03:49 | 只看该作者 Only view this author
sobereva 发表于 2019-5-7 08:07
TD=(Read  到底read了什么没交代清楚

怎么获得激发态偶极矩此文3.4节说了,必须写density

C:\Users\97355\Desktop\图片1

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

23#
发表于 Post on 2020-11-13 08:06:10 | 只看该作者 Only view this author
wzkchem5 发表于 2020-11-12 21:47
不会直接影响吸收光谱,吸收光谱的强度只和跃迁偶极矩有关,吸收峰的位置倒是可能有一些间接关系,偶极矩 ...

嗯呢,知道啦,非常感谢,辛苦啦

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

22#
发表于 Post on 2020-11-12 21:47:15 | 只看该作者 Only view this author
xxzj 发表于 2020-11-12 21:22
对偶极矩的理解不太透彻,想问一下这个矢量差会对吸收光谱有影响吗? 我是研究有机太阳能电池方向的,然 ...

不会直接影响吸收光谱,吸收光谱的强度只和跃迁偶极矩有关,吸收峰的位置倒是可能有一些间接关系,偶极矩变化越大,溶剂对吸收峰位置的影响就越大。
至于激子解离,不是说这个值大有利于激子解离,而是说这个值大说明激子解离的程度大,换句话说激子解离是偶极矩变化大的原因,而不是偶极矩变化大是激子解离的原因。直观地想,激子解离程度明显的话,正负电荷分离得比较开,自然会对偶极矩有比较大的影响。
如果想衡量是否容易极化,应该算极化率,这和偶极矩没有任何直接关系。
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

21#
发表于 Post on 2020-11-12 21:22:20 | 只看该作者 Only view this author
wzkchem5 发表于 2020-11-12 10:29
哦,这个是矢量差的模。
你是打算算什么,算吸收光谱?荧光光谱?还是要考察溶剂对吸收和发射波长的影响 ...

对偶极矩的理解不太透彻,想问一下这个矢量差会对吸收光谱有影响吗? 我是研究有机太阳能电池方向的,然后看到有的文章说这个值大,可能有利于激子解离,从而提升器件的能量转换效率。所以它大,是代表容易极化,然后电子容易转移,从而利于电子-空穴分离吗? 辛苦您啦

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

20#
发表于 Post on 2020-11-12 10:29:40 | 只看该作者 Only view this author
xxzj 发表于 2020-11-12 08:26
谢谢您,就是按照上面步骤进行计算偶极矩,然后根据下面的公式求解的

哦,这个是矢量差的模。
你是打算算什么,算吸收光谱?荧光光谱?还是要考察溶剂对吸收和发射波长的影响?
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

19#
发表于 Post on 2020-11-12 08:26:36 | 只看该作者 Only view this author
本帖最后由 xxzj 于 2020-11-12 08:30 编辑
wzkchem5 发表于 2020-11-11 22:22
你这个差值指的是标量差还是矢量差?

谢谢您,就是按照上面步骤进行计算偶极矩,然后根据下面的公式求解的

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

18#
发表于 Post on 2020-11-11 22:51:41 | 只看该作者 Only view this author
xxzj 发表于 2020-11-11 20:18
老师,关于这篇帖子,我还有一个问题需要请教老师,我也想计算基态偶极矩,激发态偶极矩还有跃迁偶极矩, ...

关于第二步,你得说清楚你到底想要哪个激发态的偶极矩。你这样写的话默认root=1,给出的是第一激发态的偶极矩
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

17#
发表于 Post on 2020-11-11 22:50:40 | 只看该作者 Only view this author
xxzj 发表于 2020-11-11 20:57
嗯呢,谢谢您,非常感谢,还有一个问题就是计算出了基态和激发态偶极矩,然后得到二者差值uge,这个值大小 ...

显然就是电子激发过程中体系偶极矩的变化
北京科音自然科学研究中心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

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

16#
发表于 Post on 2020-11-11 22:22:03 | 只看该作者 Only view this author
xxzj 发表于 2020-11-11 20:57
嗯呢,谢谢您,非常感谢,还有一个问题就是计算出了基态和激发态偶极矩,然后得到二者差值uge,这个值大小 ...

你这个差值指的是标量差还是矢量差?
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

15#
发表于 Post on 2020-11-11 20:57:52 | 只看该作者 Only view this author
wzkchem5 发表于 2020-11-11 20:32
第二步算出来的是吸收的跃迁偶极矩,第三步算出来的是发射的跃迁偶极矩
如果你只关心吸收的跃迁偶极矩, ...

嗯呢,谢谢您,非常感谢,还有一个问题就是计算出了基态和激发态偶极矩,然后得到二者差值uge,这个值大小能说明什么问题呢?

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

14#
发表于 Post on 2020-11-11 20:32:13 | 只看该作者 Only view this author
xxzj 发表于 2020-11-11 20:18
老师,关于这篇帖子,我还有一个问题需要请教老师,我也想计算基态偶极矩,激发态偶极矩还有跃迁偶极矩, ...

第二步算出来的是吸收的跃迁偶极矩,第三步算出来的是发射的跃迁偶极矩
如果你只关心吸收的跃迁偶极矩,那就不需要做第三步了
另外还有几个理解上的问题,第三步root=1优化的是最低的激发态,而不是f最大的激发态;第二步不应该简单地找最后的Dipole moment (field-independent basis, Debye)和Ground to excited state transition electric dipole moments (Au),而应该取决于你想算哪个激发态
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

13#
发表于 Post on 2020-11-11 20:18:33 | 只看该作者 Only view this author
sobereva 发表于 2019-5-7 08:07
TD=(Read  到底read了什么没交代清楚

怎么获得激发态偶极矩此文3.4节说了,必须写density

老师,关于这篇帖子,我还有一个问题需要请教老师,我也想计算基态偶极矩,激发态偶极矩还有跃迁偶极矩,和上面帖子相似:
第一步:基态结构和频率优化
#n B3LYP/6-31G(d) Opt freq
第二步:计算了20个激发态
# B3LYP/6-31G(d) TD=NStates=20 density
第三步:从第二步中找到f最大值对应的激发态,做结构和频率优化
# B3LYP/6-31G(d) TD=(Read,NStates=20,Root=1)  Opt freq

从第一步的log文件中找最后一个 Dipole moment (field-independent basis, Debye),对应的数值作为基态偶极矩
从第二步步的log文件中找最后一个 Dipole moment (field-independent basis, Debye),对应的数值作为激发态偶极矩; 找到最后一个Ground to excited state transition electric dipole moments (Au),该激发态对应的数值作为跃迁偶极矩。
就是当第二步加了density这个关键词后是不是激发态偶极矩和跃迁偶极矩都可以从其输出文件中读取,上面第三步计算就没有必要了?辛苦老师

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

GMT+8, 2026-1-23 23:20 , Processed in 0.198008 second(s), 30 queries , Gzip On.

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