计算化学公社

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

[新手求助] 关于TDDFT激发态计算的一些问题

[复制链接 Copy URL]

17

帖子

0

威望

119

eV
积分
136

Level 2 能力者

各位老师好,最近一直在学习TDDFT激发态的一些知识,其中也遇到了一些问题,肯请各位老师指点解惑,谢谢!
1.我不懂TDDFT计算激发态的信息 高斯软件是怎么计算的,理论依据是什么?是如何得到激发态的?
2.做TDDFT计算的时候不同的泛函运用得到的激发态信息不同,如何确认哪种泛函靠谱?(最好就是和实验去比对?)
3.TDDFT计算的时候 td(nstates=5),假如在PBE0泛函计算下得到S1是价层局域激发,S2是CT态,S3是里德堡态,那么此时S2,S3的真实性怎么样?那么需要换成M06-2X泛函对S2,S3重新计算?
4.在做激发态结构优化时,如果此激发态的势能面是排斥态,那么势能面上没有极小点,是不是就一直优化不出来该激发态的结构?那么该怎么解决?
5.势能面的扫描,分刚性和柔性的,那么该选择哪种扫描方式,得到的激发态势能面才是可靠的?
6.关于opt结构优化,其过程就是从输入的input 文件中的结构信息进行初猜,找到临近的极小点?还是对该结构进行全局的势能面扫描得到最小值?
以上问题有无知的地方,请老师们给予海涵,谢谢各位老师,谢谢!

418

帖子

4

威望

2910

eV
积分
3408

Level 5 (御坂)

2#
发表于 Post on 2018-12-17 20:48:01 | 只看该作者 Only view this author
TDDFT根据Mark-Casida方程求解激发能,依据好像是激发能是极化函数的奇点。

17

帖子

0

威望

119

eV
积分
136

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2018-12-17 20:57:28 | 只看该作者 Only view this author
万里云 发表于 2018-12-17 20:48
TDDFT根据Mark-Casida方程求解激发能,依据好像是激发能是极化函数的奇点。

谢谢您 您那有相关的学习资料? 我先去搜一波。谢谢!

418

帖子

4

威望

2910

eV
积分
3408

Level 5 (御坂)

4#
发表于 Post on 2018-12-17 21:10:28 | 只看该作者 Only view this author
ZMINER 发表于 2018-12-17 20:57
谢谢您 您那有相关的学习资料? 我先去搜一波。谢谢!

Springer出版过一本Fundamentals of Time-Dependent Density Functional Theory.

6万

帖子

99

威望

5万

eV
积分
120101

管理员

公社社长

5#
发表于 Post on 2018-12-18 03:44:27 | 只看该作者 Only view this author
1 你若真想了解深入,去看
含时密度泛函理论(TDDFT)综述与专著合集 http://bbs.keinsci.com/thread-6365-1-1.html
如果只是搞应用的,别去花过多时间去钻,没必要

2 看
乱谈激发态的计算方法
http://sobereva.com/265
如果计算模型和实验有可比性,一般去和实验对比来说事。没有实验就和更高精度计算级别结果对比

3 S2可能不理想,S3一定很糟。不代表M06-2X就一定合适

4 是。看具体情况

5 一般都用柔性扫描,实际化学意义比刚性扫描强得多。但绝对不能说刚性扫描的势能面是“不可靠”,这是逻辑错误。势能面完全是由计算级别定义的

6 opt只会收敛到当前势能面上与初猜结构最近的极小点
北京科音自然科学研究中心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

17

帖子

0

威望

119

eV
积分
136

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2018-12-18 09:15:31 | 只看该作者 Only view this author
sobereva 发表于 2018-12-18 03:44
1 你若真想了解深入,去看
含时密度泛函理论(TDDFT)综述与专著合集 http://bbs.keinsci.com/thread-6365-1 ...

感谢sob老师的指导。
1.谢谢老师分享理论知识连接;
2.更高精度指的是更好更高的计算水平?
3.那么怎么确保算出来的每个激发态信息才是可靠的?与问题2一样,和更高的精度计算进行比较?
4.那么激发态结构优化不出来,怎么扫势能面(思路是,从基态出发优化出激发态结构,再从优化的激发态结构出发扫描其势能面,不知道我的做法对不对?)
5,6.谢谢sob老师解答疑惑。
谢谢老师

6万

帖子

99

威望

5万

eV
积分
120101

管理员

公社社长

7#
发表于 Post on 2018-12-18 15:25:27 | 只看该作者 Only view this author
ZMINER 发表于 2018-12-18 09:15
感谢sob老师的指导。
1.谢谢老师分享理论知识连接;
2.更高精度指的是更好更高的计算水平?

2 y
3 可以考虑TD-双杂化泛函、EOM-CCSD、CC2等,普适性和可靠性都比单一普通泛函要强
4 什么情况都可以做扫描,怎么扫描取决于具体目的,和激发态极小点能否优化出来是两码事
北京科音自然科学研究中心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

17

帖子

0

威望

119

eV
积分
136

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2018-12-18 16:30:04 | 只看该作者 Only view this author
sobereva 发表于 2018-12-18 15:25
2 y
3 可以考虑TD-双杂化泛函、EOM-CCSD、CC2等,普适性和可靠性都比单一普通泛函要强
4 什么情况都可 ...

谢谢sob老师耐心指导解惑。谢谢您

4.那我做S2态势能面柔性扫描 比如 #p opt=modredundant td(nstates=6,singlets,root=2) b3lyp/6-311+g(d,p) geom=connectivity 设置反应坐标(扫描步长等),那么结构信息可以直接用基态优化好的结构?这么做可行?

谢谢老师!

6万

帖子

99

威望

5万

eV
积分
120101

管理员

公社社长

9#
发表于 Post on 2018-12-19 14:46:45 | 只看该作者 Only view this author
ZMINER 发表于 2018-12-18 16:30
谢谢sob老师耐心指导解惑。谢谢您

4.那我做S2态势能面柔性扫描 比如 #p opt=modredundant td(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

17

帖子

0

威望

119

eV
积分
136

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2018-12-19 16:36:10 | 只看该作者 Only view this author

谢谢老师解答疑惑,受益匪浅。谢谢您!

82

帖子

0

威望

512

eV
积分
594

Level 4 (黑子)

11#
发表于 Post on 2019-4-1 16:55:18 | 只看该作者 Only view this author
本帖最后由 huangyoyo 于 2019-4-1 17:22 编辑
sobereva 发表于 2018-12-18 03:44
1 你若真想了解深入,去看
含时密度泛函理论(TDDFT)综述与专著合集 http://bbs.keinsci.com/thread-6365-1 ...

请问我用ORCA计算了一个Ne原子的激发态,第7个激发态跟第12个激发态轨道跃迁的成分差不多为什么激发能会差那么多呢?哪个数值才表示基态到2P轨道的跃迁能?

9

帖子

0

威望

73

eV
积分
82

Level 2 能力者

12#
发表于 Post on 2020-10-19 09:14:02 | 只看该作者 Only view this author
老师们好,我做了一个锕系内嵌富勒烯结构的TDDFT计算,计算的泛函用的是m062x,基组为混合基组。内嵌元素用的基组是genecp赝势基组,C原子使用的是6-31G*基组。审稿人提出了质疑基组使用的问题:“Very large basis set has been used for the center atom ((14s13p10d8f6g)/ [10s9p5d4f3g]), however, the basis set used for the C atom may not be adequate to investigate weak superatom-superatom interaction, and particularly for its excited states. So large size basis set must be used for the C atom.”是否我的基组使用确实有误,我该如何有利的回复?谢谢大家~

1万

帖子

0

威望

8983

eV
积分
20731

Level 6 (一方通行)

13#
发表于 Post on 2020-10-19 10:14:15 | 只看该作者 Only view this author
Believe 发表于 2020-10-19 09:14
老师们好,我做了一个锕系内嵌富勒烯结构的TDDFT计算,计算的泛函用的是m062x,基组为混合基组。内嵌元素用 ...

不存在“genecp赝势基组”这种说法,genecp只是读入赝势的一个关键字,不是赝势基组的名字。真正的赝势基组是你在坐标部分后面写的那些东西。
从审稿人的回复看,你应该是在算两个内嵌富勒烯之间的色散相互作用?那至少应该用triple zeta或加弥散,不然BSSE很明显。
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?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

1657

帖子

5

威望

4556

eV
积分
6313

Level 6 (一方通行)

喵星人

14#
发表于 Post on 2020-10-19 10:21:14 | 只看该作者 Only view this author
锕系你为啥用M062X啊。。。。这种金属+C团簇的用PBE0多稳
然后赝势基组用的啥?
另外审稿人的意思很明确,C用的基组太小了,至少上TZVP吧

16

帖子

0

威望

53

eV
积分
69

Level 2 能力者

15#
发表于 Post on 2025-5-16 09:48:26 | 只看该作者 Only view this author
wzkchem5 发表于 2020-10-19 10:14
不存在“genecp赝势基组”这种说法,genecp只是读入赝势的一个关键字,不是赝势基组的名字。真正的赝势基 ...

王老师,我想请问我使用ORCA 计算异核双原子分子激发(含过渡金属),用TD-DFT cc-PVTZ基组,但是我使用某些泛函的时候本该双重简并的态,这两个能量数值出现小数点后两位的偏差(单位是eV),而使用另外一些可能层次高一些的泛函,却也成功算出了能量相同,我以为是泛函水平低的问题,但是另外一个老师告诉我水平低的泛函也不会把简并的态算的差别如此之大,我也使用对称性,那么请问这是数值积分的问题吗,我在其他软件上用波函数方法其实并没有出现此问题

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 22:28 , Processed in 0.233060 second(s), 21 queries , Gzip On.

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