计算化学公社

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

[综合交流] 激发能算出来是负值怎么处理?

[复制链接 Copy URL]

49

帖子

0

威望

465

eV
积分
514

Level 4 (黑子)

各位老师好。最近在算一个有机分子,想计算得到他的绝热能级图,但是在计算T2时得到的发射能为负值。
我尝试了社长在博文里面提到的解决方法,但是结果还是优化着又变为负值了。

首先我是在基态优化好的结构上去做的S1,T1的发射能,数值都为正,很正常,用TDDFT计算T2时,出现如图一的负值结果。
然后我使用了#p stable=opt b3lyp/6-31g(d) scrf=(solvent=dichloromethane) 关键词找出真正的基态波函数,激发态计算的时候写#p guess=read b3lyp/6-31g(d) scrf=(solvent=dichloromethane) opt td(nstates=5,triplet,root=2)进行计算,得到图二结果。
是我对社长博文理解的不对吗?请给各位老师多多指教,万分谢谢。

XKB82)94@M)2B]$]~J{BW)Q.png (20.46 KB, 下载次数 Times of downloads: 36)

图一

图一

Y)FG`B%RFA)K~M8%Z%RDXBK.png (22.04 KB, 下载次数 Times of downloads: 35)

图二

图二

1236

帖子

1

威望

3495

eV
积分
4751

Level 6 (一方通行)

2#
发表于 Post on 2019-8-28 10:09:06 | 只看该作者 Only view this author
可能stable=opt并未找到稳定波函数,跑完这个你要看一下输出里是否找到了稳定波函数

49

帖子

0

威望

465

eV
积分
514

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2019-8-28 13:58:27 | 只看该作者 Only view this author
biogon 发表于 2019-8-28 10:09
可能stable=opt并未找到稳定波函数,跑完这个你要看一下输出里是否找到了稳定波函数

就是我用stable=opt关键词计算结果显示波函数已经稳定了,那您的意思是,我要去后面计算激发态的输出文件里再看是否波函数稳定吗?

_]ZMI}T6AD20[TOE$6P5W7Q.png (21.81 KB, 下载次数 Times of downloads: 41)

_]ZMI}T6AD20[TOE$6P5W7Q.png

190

帖子

0

威望

6100

eV
积分
6290

Level 6 (一方通行)

4#
发表于 Post on 2019-8-28 16:55:30 | 只看该作者 Only view this author
可能是你计算的分子具有明显的电荷转移性质,B3LYP泛函明显低估此类分子的发射能,建议尝试wB97XD或者M062X泛函。
https://www.x-mol.com/groups/weijie_chi

1236

帖子

1

威望

3495

eV
积分
4751

Level 6 (一方通行)

5#
发表于 Post on 2019-8-28 17:03:33 | 只看该作者 Only view this author
liyuanfu7 发表于 2019-8-28 13:58
就是我用stable=opt关键词计算结果显示波函数已经稳定了,那您的意思是,我要去后面计算激发态的输出文件 ...

我现在就好奇这啥分子

5万

帖子

99

威望

5万

eV
积分
112373

管理员

公社社长

6#
发表于 Post on 2019-8-28 18:45:01 | 只看该作者 Only view this author
如果你是指优化读取stable=opt的波函数刚开始优化T2的时候仍然所有激发能都是正值,但是后来才有的成为负值,用最后一步的结构再次做stable=opt,然后再guess=read进一步继续优化。
如果是第一步输出的激发能就有负值,说明stable=opt没真正奏效,尝试通过改变自旋多重度、调换轨道等方式得到真正基态波函数
北京科音自然科学研究中心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!

49

帖子

0

威望

465

eV
积分
514

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2019-8-28 19:04:36 | 只看该作者 Only view this author
sobereva 发表于 2019-8-28 18:45
如果你是指优化读取stable=opt的波函数刚开始优化T2的时候仍然所有激发能都是正值,但是后来才有的成为负值 ...

是的老师,我是在优化读取stable=opt的波函数刚开始优化T2的时候仍然所有激发能都是正值,但是后来才有的成为负值。它是到了第八次迭代出现的负值,我查看了一下,之前用TDDTF的方法直接计算T2也是到了第八次出现的负值,且得到的负值和用stable=opt之后得到的值一样。这该怎么解释呢?
老师那我可以用读取stable=opt的波函数用TDA进行计算T2来尝试一下吗?

49

帖子

0

威望

465

eV
积分
514

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2019-8-28 19:06:45 | 只看该作者 Only view this author
biogon 发表于 2019-8-28 17:03
我现在就好奇这啥分子

这是一个能发激光的TADF分子,和它同一系列的其它分子计算出来数值也都很正常,搞不懂为啥就它不一样

49

帖子

0

威望

465

eV
积分
514

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2019-8-28 19:09:28 | 只看该作者 Only view this author
chiweijie 发表于 2019-8-28 16:55
可能是你计算的分子具有明显的电荷转移性质,B3LYP泛函明显低估此类分子的发射能,建议尝试wB97XD或者M062X ...

也不是很明显的CT,就是我用该基组和泛函计算出来做得到的吸收发射值能比较好的与实验吻合,并且我计算该分子同一系列的其他几个计算结果都挺正常,所以我才没考虑去换泛函。

1236

帖子

1

威望

3495

eV
积分
4751

Level 6 (一方通行)

10#
发表于 Post on 2019-8-28 20:02:03 | 只看该作者 Only view this author
liyuanfu7 发表于 2019-8-28 19:09
也不是很明显的CT,就是我用该基组和泛函计算出来做得到的吸收发射值能比较好的与实验吻合,并且我计算该 ...

激发能为负和用什么泛函是没什么关系,我在想是不是有什么BUG存在,你把这个分子的结构微微改动一下再算看看是是不是还会出现这个问题

151

帖子

0

威望

3277

eV
积分
3428

Level 5 (御坂)

11#
发表于 Post on 2019-8-28 21:25:54 | 只看该作者 Only view this author
可能是三重态不稳定,用CAM-B3LYP加TDA。

49

帖子

0

威望

465

eV
积分
514

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2019-8-28 21:45:00 | 只看该作者 Only view this author
superrice 发表于 2019-8-28 21:25
可能是三重态不稳定,用CAM-B3LYP加TDA。

嗯嗯,我也想过用TDA,那这样的话,是不是我的T1也应该用TDA呢?

49

帖子

0

威望

465

eV
积分
514

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2019-8-28 21:46:20 | 只看该作者 Only view this author
biogon 发表于 2019-8-28 20:02
激发能为负和用什么泛函是没什么关系,我在想是不是有什么BUG存在,你把这个分子的结构微微改动一下再算 ...

嗯嗯,那我试一试改变一点儿二面角,谢谢老师

13

帖子

0

威望

203

eV
积分
216

Level 3 能力者

14#
发表于 Post on 2022-4-28 17:17:53 | 只看该作者 Only view this author
同学 想问下你这个问题解决了吗 我在用TDDFT优化S1的时候也遇到了激发能为负值的情况,想请教一下

1万

帖子

0

威望

7396

eV
积分
18149

Level 6 (一方通行)

15#
发表于 Post on 2022-4-28 18:02:06 | 只看该作者 Only view this author
happiness 发表于 2022-4-28 10:17
同学 想问下你这个问题解决了吗 我在用TDDFT优化S1的时候也遇到了激发能为负值的情况,想请教一下

2楼已经有人回答了
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 21:01 , Processed in 0.225702 second(s), 23 queries , Gzip On.

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