计算化学公社

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

[综合交流] 垂直与绝热激发能差异大,基态激发态结构差异也大,选基态还是激发态结构做TDDFT

[复制链接 Copy URL]

18

帖子

0

威望

689

eV
积分
707

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 zjx-leslie 于 2022-11-16 20:22 编辑

算一个激发FAD转移电子到聚合单体的体系,对基态和感兴趣的激发态先后做了优化。结果FAD的体系基态与激发态结构差异很大,计算得到的垂直激发能和绝热激发能也差将近0.7eV,同时绝热激发能与实验的到的数据更加吻合。

想请教一下各位,遇到这种情况,如果要去分析空穴-电子,我应该选取基态结构去做TDDFT,还是用感兴趣的激发态结构去做TDDFT?


此外也想问一下大家有没有讨论绝热激发能与垂直激发能的文章,我去拜读一下。


附件里我也附上了包含部分结构的gjf文件,都是已经优化好的结构。



Excited_state.gjf

2.15 KB, 下载次数 Times of downloads: 10

激发态优化后结构

Ground_state.gjf

2.11 KB, 下载次数 Times of downloads: 8

基态优化后结构

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2022-11-16 21:34:25 | 只看该作者 Only view this author
都可以分析。基态结构的空穴-电子分析用来讨论吸收,激发态结构的空穴-电子分析用来讨论从相应态的发射。如果讨论垂直激发后弛豫过程的电子结构变化,可以对激发态的FC点结构和极小点结构分别做布居分析、键级分析等来比较差异。
北京科音自然科学研究中心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

18

帖子

0

威望

689

eV
积分
707

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-11-17 14:35:04 | 只看该作者 Only view this author
本帖最后由 zjx-leslie 于 2022-11-17 14:36 编辑
sobereva 发表于 2022-11-16 21:34
都可以分析。基态结构的空穴-电子分析用来讨论吸收,激发态结构的空穴-电子分析用来讨论从相应态的发射。如 ...

好的,感谢sob老师的回复,目前我对基态结构做了TDDFT,但垂直激发能比实验测得的数据大0.7eV,但用激发态结构算得的绝热激发能和实验值很接近,对实验现象进行解释时,能否认为是在激发过程中,结构发生较大变化导致垂直激发能与实验不吻合

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

4#
发表于 Post on 2022-11-17 16:12:34 | 只看该作者 Only view this author
zjx-leslie 发表于 2022-11-17 14:35
好的,感谢sob老师的回复,目前我对基态结构做了TDDFT,但垂直激发能比实验测得的数据大0.7eV,但用激发 ...

这种事要全面考虑各种可能的因素



至于你提到的结构问题,实验上没法直接测垂直或是绝热激发能。精确跟实验对照需要做振动耦合电子光谱计算。不考虑振动耦合的情况,看




北京科音自然科学研究中心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

18

帖子

0

威望

689

eV
积分
707

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-11-17 17:05:12 | 只看该作者 Only view this author
sobereva 发表于 2022-11-17 16:12
这种事要全面考虑各种可能的因素

好的,谢谢sob老师
目前我用的激发态优化关键词是#p opt freq td=(nstates=6,root=2,restart) cam-b3lyp/6-31g(d) scrf=(smd,solvent=water) nosymm em=gd3bj
TDDFT计算的关键词是#P CAM-B3LYP/6-311g(d) scrf=(smd,solvent=water) nosymm em=gd3bj TD(nstates=30) IOp(9/40=4)
基态优化关键词是#p opt freq b3lyp/6-31g(d) scrf=(smd,solvent=water) nosymm em=gd3bj
体系是从酶活性口袋的簇模型,大约200个原子,其中结构变化较大的FAD部分也放在了附件中,请您帮忙看一下关键词是否合理,由于计算量较大,目前暂时没有做构象平均

我计算得到的垂直激发能是3.8eV,绝热激发能是2.8eV,实际测得反应在光照2.7-3.0eV能够发生,垂直激发能与绝热激发能差异很大,与实验值差异也很大,您看一下是否是我的计算方法不对或者是数据应该做一些校正或者进一步处理

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

6#
发表于 Post on 2022-11-17 17:12:13 | 只看该作者 Only view this author
zjx-leslie 发表于 2022-11-17 10:05
好的,谢谢sob老师
目前我用的激发态优化关键词是#p opt freq td=(nstates=6,root=2,restart) cam-b3lyp ...

CAM-B3LYP有高估局域激发能的倾向。此外酶口袋里的介电常数(一般在4~8左右)远小于水的介电常数(80左右),用水作溶剂未必是最合适的选择
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
本团队长期招收研究生,有意者可私信联系

18

帖子

0

威望

689

eV
积分
707

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-11-17 18:17:37 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-17 17:12
CAM-B3LYP有高估局域激发能的倾向。此外酶口袋里的介电常数(一般在4~8左右)远小于水的介电常数(80左右 ...

感谢您的回复,请问您是否有相关的文章推荐,我去拜读一下,您提到的介电常数的问题确实很关键,我也用TDDFT算了不加溶剂化模型的结构,激发能反而与实验结果比较吻合,因此我也猜想出现这种情况的原因,是否是由于酶活性口袋本身是一个疏水环境,所以加了溶剂模型反而会导致激发能被高估?

18

帖子

0

威望

689

eV
积分
707

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2022-11-17 18:22:16 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-17 17:12
CAM-B3LYP有高估局域激发能的倾向。此外酶口袋里的介电常数(一般在4~8左右)远小于水的介电常数(80左右 ...

像酶活性口袋这种体系,酶的外部是水环境,但口袋内部又是疏水的,拿来做簇模型的话,溶剂模型怎么加也是个问题,请问您是否有好的建议或者相关的参考文献,非常感谢您~

1万

帖子

0

威望

9859

eV
积分
22096

Level 6 (一方通行)

9#
发表于 Post on 2022-11-17 18:32:25 | 只看该作者 Only view this author
zjx-leslie 发表于 2022-11-17 11:17
感谢您的回复,请问您是否有相关的文章推荐,我去拜读一下,您提到的介电常数的问题确实很关键,我也用TD ...

很多用团簇模型算酶的人会用氯苯模拟蛋白环境,我一时找不到相关文献,你可以搜一下
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 01:09 , Processed in 0.169204 second(s), 22 queries , Gzip On.

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