计算化学公社

 找回密码 Forget password
 注册 Register
Views: 26816|回复 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: 76)

微信图片_20190506151514.png

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

2#
发表于 Post on 2019-5-7 08:02:48 | 只看该作者 Only view this author
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出“提问”、“求助”要素(仔细看http://bbs.keinsci.com/thread-9348-1-1.html),我已把你的标题“高斯算激发态偶极矩和跃迁偶极矩”改了,以后务必注意
北京科音自然科学研究中心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
积分
124671

管理员

公社社长

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

怎么获得激发态偶极矩此文3.4节说了,必须写density
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314http://bbs.keinsci.com/thread-2413-1-1.html

是基组不是机组
北京科音自然科学研究中心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

163

帖子

0

威望

1025

eV
积分
1188

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2019-5-7 08:56:24 | 只看该作者 Only view this author
sobereva 发表于 2019-5-7 08:02
如置顶的新社员必读贴和论坛首页的公告栏所示,求助帖必须在帖子标题明确体现出“提问”、“求助”要素(仔 ...

好的,谢谢社长提醒

163

帖子

0

威望

1025

eV
积分
1188

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2019-5-7 08:57:45 | 只看该作者 Only view this author
sobereva 发表于 2019-5-7 08:07
TD=(Read  到底read了什么没交代清楚

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

谢谢sob老师,请问如果是需要在上一步的基础上继续算的话,并且读取20个激发态中的其中一个,那么read中的关键词应该怎么写?

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

6#
发表于 Post on 2019-5-7 09:08:18 | 只看该作者 Only view this author
Abandon-fmt 发表于 2019-5-7 08:57
谢谢sob老师,请问如果是需要在上一步的基础上继续算的话,并且读取20个激发态中的其中一个,那么read中 ...

什么叫继续算?如果之前的态数算少了,算更多态的时候用read虽然名义上能节约时间,但往往实际起到反效果
北京科音自然科学研究中心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

163

帖子

0

威望

1025

eV
积分
1188

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2019-5-7 09:59:34 | 只看该作者 Only view this author
sobereva 发表于 2019-5-7 09:08
什么叫继续算?如果之前的态数算少了,算更多态的时候用read虽然名义上能节约时间,但往往实际起到反效果

好的,谢谢sob老师!

163

帖子

0

威望

1025

eV
积分
1188

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2019-5-7 10:12:01 | 只看该作者 Only view this author
sobereva 发表于 2019-5-7 08:07
TD=(Read  到底read了什么没交代清楚

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

加了density之后,激发态偶极矩我算出来了,再请问一下sob老师,跃迁偶极矩是找到log文件中最后一个Ground to excited state transition electric dipole moments (Au),根据第二步的f最大值对应的激发态,找到该激发态对应的数值么?

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

9#
发表于 Post on 2019-5-8 16:40:41 | 只看该作者 Only view this author
Abandon-fmt 发表于 2019-5-7 10:12
加了density之后,激发态偶极矩我算出来了,再请问一下sob老师,跃迁偶极矩是找到log文件中最后一个Groun ...

考察哪个态完全取决于自己,没有非得看f最大的这一条
北京科音自然科学研究中心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

163

帖子

0

威望

1025

eV
积分
1188

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2019-5-8 20:55:17 | 只看该作者 Only view this author
sobereva 发表于 2019-5-8 16:40
考察哪个态完全取决于自己,没有非得看f最大的这一条

好的,谢谢sob老师

330

帖子

0

威望

3006

eV
积分
3336

Level 5 (御坂)

计算化学路人甲

11#
发表于 Post on 2019-8-23 22:16:43 | 只看该作者 Only view this author
sobereva 发表于 2019-5-8 16:40
考察哪个态完全取决于自己,没有非得看f最大的这一条

社长,我有个疑惑,写上density后,对应的电偶极矩,是仅仅分子构型为激发态构型,但是仍是未激发的分子(即分子的电子排布为HOMO及以下轨道全占据,LUMO及以上轨道全为空)时的偶极矩吗?
如果是这种情况,那么我怎么获得分子处于激发状态下(LUMO及以上电子轨道有电子占据)的偶极矩呢?

6万

帖子

99

威望

5万

eV
积分
124671

管理员

公社社长

12#
发表于 Post on 2019-8-23 22:26:00 | 只看该作者 Only view this author
Novice 发表于 2019-8-23 22:16
社长,我有个疑惑,写上density后,对应的电偶极矩,是仅仅分子构型为激发态构型,但是仍是未激发的分子 ...

root设成几,density关键词产生的波函数及相关信息(包括偶极矩)就对应第几激发态

另见
使用Multiwfn计算激发态间的跃迁偶极矩和各个激发态的偶极矩
http://sobereva.com/227
北京科音自然科学研究中心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

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这个关键词后是不是激发态偶极矩和跃迁偶极矩都可以从其输出文件中读取,上面第三步计算就没有必要了?辛苦老师

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 (黑子)

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

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

本版积分规则 Credits rule

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

GMT+8, 2026-1-24 02:42 , Processed in 0.180247 second(s), 23 queries , Gzip On.

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