计算化学公社

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

[Gaussian/gview] Gaussian 16 用TDDFT优化激发态结构收敛和断点续跑问题

[复制链接 Copy URL]

4

帖子

0

威望

25

eV
积分
29

Level 2 能力者

本帖最后由 莱万少夫斯基 于 2023-3-22 07:40 编辑

大家好。我正在使用Gaussian 16通过TDDFT优化激发态结构,模拟UV-Vis。收敛问题。我发现在out输出文件中,似乎已经显示在某一个interation全部收敛,然而从下一个interation开始又变成全部不收敛了。我比较困惑这种情况是否正常。另外,我发现使用现在的out文件,Multiwfn已经可以模拟出UV-Vis图谱。我不知道是否可以停止优化?#P opt=restart B3LYP/gen pseudo=read TD(nstates=203,restart) guess=read。希望可以获得大佬的指点!不胜感激!

202303210339279119..png (86.71 KB, 下载次数 Times of downloads: 7)

从输出文件中搜索“converged”关键词显示

从输出文件中搜索“converged”关键词显示

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2023-3-21 04:19:26 | 只看该作者 Only view this author
你犯了一个关键性错误,就是你当前优化的是第一激发态,你计算最低3~5个激发态就足够了,算203个态完全是无意义地严重白浪费时间。而且模拟UV-Vis并不需要优化激发态。算发射光谱才需要优化激发态,而此时明显又用不着算200个态。
仔细把此文看了,该说的基本知识都说了:
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

4

帖子

0

威望

25

eV
积分
29

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2023-3-21 05:28:20 | 只看该作者 Only view this author
本帖最后由 莱万少夫斯基 于 2023-3-21 06:00 编辑
sobereva 发表于 2023-3-21 04:19
你犯了一个关键性错误,就是你当前优化的是第一激发态,你计算最低3~5个激发态就足够了,算203个态完全是无 ...

谢谢社长回复。可是我是按照这个体系的相关文献来的。他们模拟了S0到S200的TDDFT去获得UV-Vis光谱。而且,在此之前我优化过类似最低的5个态,那个时候获得的光谱区间离UV-Vis相当远。当我增加到200个态的时候,就已经和文献报道的光谱区间比较吻合了。我看了那个贴文,不知道我理解的对不对。我们需要计算垂直吸收能,要保持能量基态最低结构。如果要计算到S200的结构不优化激发态的话,正确的设置是否是#P B3LYP/gen pseudo=read TD(nstates=203, root=203) guess=read? 请您明示,谢谢。

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2023-3-21 06:19:05 | 只看该作者 Only view this author
莱万少夫斯基 发表于 2023-3-21 05:28
谢谢社长回复。可是我是按照这个体系的相关文献来的。他们模拟了S0到S200的TDDFT去获得UV-Vis光谱。而且 ...

完全错误,基本概念和逻辑还没搞明白。没好好看我的文章。
北京科音自然科学研究中心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

4

帖子

0

威望

25

eV
积分
29

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2023-3-21 07:07:01 | 只看该作者 Only view this author
sobereva 发表于 2023-3-21 06:19
完全错误,基本概念和逻辑还没搞明白。没好好看我的文章。

感谢您的回复。按照您4.1的例子,我使用基态结构
#P B3LYP/gen pseudo=read TD(nstates=203,root=200),这样的设置是否可以正确的用于考察S0到S200的垂直吸收能,并用于UV-Vis的模拟?

245

帖子

0

威望

780

eV
积分
1025

Level 4 (黑子)

6#
发表于 Post on 2023-3-21 13:06:21 | 只看该作者 Only view this author
请问楼主方便提供下相关的文章吗?这吸收光谱对应到S0→S200的激发我还真是挺好奇是什么样的体系… 平常如果我算的低激发态这么红的不可接受我就换泛函了,从来没想过(也没条件)一次算203个激发能的……

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

7#
发表于 Post on 2023-3-21 16:56:41 | 只看该作者 Only view this author
莱万少夫斯基 发表于 2023-3-21 00:07
感谢您的回复。按照您4.1的例子,我使用基态结构
#P B3LYP/gen pseudo=read TD(nstates=203,root=200), ...

可以,这样会算203个激发态,其中除了最高的几个激发态可能例外以外,第n个激发态就是能量第n低的激发态。虽然不能保证前200个激发态都严格满足”第n个激发态就是能量第n低的激发态“,但是这一点引入的不确定性对于画UV-Vis图而言是可以接受的。
另外注意此时无需写root,但写了也不影响结果的可用性
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
本团队长期招收研究生,有意者可私信联系

4

帖子

0

威望

25

eV
积分
29

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2023-3-21 19:46:18 | 只看该作者 Only view this author
wzkchem5 发表于 2023-3-21 16:56
可以,这样会算203个激发态,其中除了最高的几个激发态可能例外以外,第n个激发态就是能量第n低的激发态 ...

感谢您的回复!这样是不是还得考虑三重激发态的可能,然后把三重激发态和单重态一起对比才能说谁是能量最低的n个激发态呢?

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

9#
发表于 Post on 2023-3-21 20:16:39 | 只看该作者 Only view this author
莱万少夫斯基 发表于 2023-3-21 12:46
感谢您的回复!这样是不是还得考虑三重激发态的可能,然后把三重激发态和单重态一起对比才能说谁是能量最 ...

我说得不准确,我的意思是”第n个激发态就是能量第n低的单重态激发态“,也就是你的计算即使连这个也不能保证,但是这一点对计算UV-Vis的影响可以忽略

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
莱万少夫斯基 + 1 谢谢

查看全部评分 View all ratings

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
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 22:34 , Processed in 0.193324 second(s), 30 queries , Gzip On.

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