计算化学公社

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

[Gaussian/gview] TDDFT优化激发能变成了负值

[复制链接 Copy URL]

227

帖子

0

威望

2800

eV
积分
3027

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
g09,tddft优化第一三重态,从第三步几何开始,激发能变成了负值。还有继续优化下去的意义么?

谢谢

2494

帖子

11

威望

7054

eV
积分
9768

Level 6 (一方通行)

2#
发表于 Post on 2017-6-30 18:19:39 | 只看该作者 Only view this author
请问你前两步做啥了?

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

3#
发表于 Post on 2017-6-30 19:57:06 | 只看该作者 Only view this author

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
小范范1989 + 5 冰老师得到sob老师的真传了,资料太丰富了

查看全部评分 View all ratings

Stand on the shoulders of giants

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

4#
发表于 Post on 2017-7-1 01:00:45 | 只看该作者 Only view this author
你优化T1直接用0 3在UDFT下优化就完了,别用TDDFT,这里提了
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314

为负说明当前参考态不是真正基态波函数,很可能S0和T1发生了交叉,在第三步的结构开始T1成了真正的基态。
北京科音自然科学研究中心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

227

帖子

0

威望

2800

eV
积分
3027

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2017-7-2 13:51:33 | 只看该作者 Only view this author
sobereva 发表于 2017-7-1 01:00
你优化T1直接用0 3在UDFT下优化就完了,别用TDDFT,这里提了
Gaussian中用TDDFT计算激发态和吸收、荧光、 ...

需要和T2比较

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

6#
发表于 Post on 2017-7-2 22:40:19 | 只看该作者 Only view this author

T2再说T2的事,先把T1弄完了再说
北京科音自然科学研究中心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

10

帖子

0

威望

49

eV
积分
59

Level 2 能力者

7#
发表于 Post on 2019-3-7 19:10:19 | 只看该作者 Only view this author
本帖最后由 藍天許 于 2019-3-7 19:15 编辑

請問如果優化T1,激發能為負值;優化T2時出現l9999.exe 的錯誤;讀進T1結構優化T3,激發能為正,這樣也代表當前的參考態不是真正的S0嗎?

---------------
#P opt td=(triplet, root=3) wb97xd/6-31G(d) int=acc2e=11 scf=(novaracc,noincfock)
---------------

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

8#
发表于 Post on 2019-3-8 06:01:36 | 只看该作者 Only view this author
藍天許 发表于 2019-3-7 19:10
請問如果優化T1,激發能為負值;優化T2時出現l9999.exe 的錯誤;讀進T1結構優化T3,激發能為正,這樣也代表 ...

不同结构下的情况不一样,可能有的结构下默认设置就能收敛到稳定的波函数,但有的结构下可能不能。
只能说在你优化T3最后的结构下,你的参考态波函数已经是这个结构下的真实的S0了。
北京科音自然科学研究中心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

10

帖子

0

威望

49

eV
积分
59

Level 2 能力者

9#
发表于 Post on 2019-3-9 14:41:57 | 只看该作者 Only view this author
sobereva 发表于 2019-3-8 06:01
不同结构下的情况不一样,可能有的结构下默认设置就能收敛到稳定的波函数,但有的结构下可能不能。
只能 ...

好的 謝謝您

230

帖子

0

威望

790

eV
积分
1020

Level 4 (黑子)

10#
发表于 Post on 2019-6-16 17:01:39 | 只看该作者 Only view this author
sobereva 发表于 2017-7-1 01:00
你优化T1直接用0 3在UDFT下优化就完了,别用TDDFT,这里提了
Gaussian中用TDDFT计算激发态和吸收、荧光、 ...

Sob老师好,我最近在进行激发态计算的时候也遇到了激发能变为负值的问题。想比较的是T1和S1的能量差,所以直接用了50-50的选项。关键词如下:
#p B3LYP/TZVP opt td(50-50,nstates=5) em=gd3bj
结果计算过程中,T1的激发能就变成了负的。
想问您一下,比较T1和S1的能量值我这样写关键词合理吗,是否是由于我写了opt,所以对三重态结构进行优化了,导致的三重态能量低于基态呢?

是不是先优化一个基态结构,然后用优化好的基态结构,在计算激发态能量时把opt去掉更合理一点呢?

或者是先优化一个基态结构,用stable=opt找到稳定的波函数,再把关键词改为:#p B3LYP/TZVP opt td(50-50,nstates=5) em=gd3bj guess=read 这样进行操作可以吗?

还有想问您的就是这个作业最后其实并没能正常结束,报了一个没有遇到过的l914错误:

想问您一下,这个错误是什么原因导致的,怎么解决比较好?

第一次遇到比较S1和T1能量的问题,不是很会操作,辛苦老师解答一下了,谢谢Sob老师

最后我把我之前的输入文件附在后面,或许您能够用的到。
sun_01a_b3d3_50_50_td_TZVP (1.62 KB, 下载次数 Times of downloads: 17)

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

11#
发表于 Post on 2019-6-17 01:22:46 | 只看该作者 Only view this author
Jameslzy 发表于 2019-6-16 17:01
Sob老师好,我最近在进行激发态计算的时候也遇到了激发能变为负值的问题。想比较的是T1和S1的能量差,所 ...

优化激发态的时候甭用50-50,否则你都不知道被优化的那个态到底是S1还是T1
北京科音自然科学研究中心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

230

帖子

0

威望

790

eV
积分
1020

Level 4 (黑子)

12#
发表于 Post on 2019-6-17 09:21:53 | 只看该作者 Only view this author
sobereva 发表于 2019-6-17 01:22
优化激发态的时候甭用50-50,否则你都不知道被优化的那个态到底是S1还是T1

谢谢老师回复。那么我先优化基态结构,然后用优化好的基态结构,在计算激发态能量时把opt去掉,这时用50-50比较T1和S1的能量,这样操作可以吗?还有就是想问下老师图里的l914错误是什么原因产生的,怎么解决比较好呢?谢谢老师

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

13#
发表于 Post on 2019-6-17 21:56:11 | 只看该作者 Only view this author
Jameslzy 发表于 2019-6-17 09:21
谢谢老师回复。那么我先优化基态结构,然后用优化好的基态结构,在计算激发态能量时把opt去掉,这时用50- ...

如果你想比基态结构下S1和T1能量,这样可以

没完整的输入输出文件没法说
北京科音自然科学研究中心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

230

帖子

0

威望

790

eV
积分
1020

Level 4 (黑子)

14#
发表于 Post on 2019-6-17 22:52:47 | 只看该作者 Only view this author
sobereva 发表于 2019-6-17 21:56
如果你想比基态结构下S1和T1能量,这样可以

没完整的输入输出文件没法说

谢谢老师回复。如果想比较S1和T1最稳定结构下S1和T1能量就要分别用singlet和triplet关键词以及opt关键词分别交两个作业来比较电子能了是吧?之前说报错的输入和输出文件我放在下面附件里了,麻烦老师帮忙看一下,谢谢老师
sun_03_out.zip (280.06 KB, 下载次数 Times of downloads: 16)
sun_03_b3d3_50_50_td_TZVP (2.77 KB, 下载次数 Times of downloads: 11)




6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

15#
发表于 Post on 2019-6-18 02:51:01 | 只看该作者 Only view this author
Jameslzy 发表于 2019-6-17 22:52
谢谢老师回复。如果想比较S1和T1最稳定结构下S1和T1能量就要分别用singlet和triplet关键词以及opt关键词 ...


对于T1结构尽量用UKS而非TDDFT方式优化

把50-50去了再说,还有问题把nstates增加再试
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 05:03 , Processed in 1.130931 second(s), 24 queries , Gzip On.

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