计算化学公社

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

[综合交流] 关于优化激发态发生势能面交叉,如何得到所希望态的优化结果

[复制链接 Copy URL]

338

帖子

0

威望

1107

eV
积分
1445

Level 4 (黑子)

本帖最后由 cokie 于 2023-7-28 11:51 编辑

想请教一下各位老师,我在做分子优化时,出现了势能面交叉,我希望优化得到T3的极小点(即基态下排序Excited state2这条势能面的极小点)。
首先我分别优化了T2和T3,方法如下:
1. 基于优化后的S0,做tddft triplet nstates=5, root=3来优化
2. 基于优化后的S0,做tddft triplet nstates=5, root=2来优化
结果是,优化T2正常,跟踪过程没有发生势能面交叉。但优化T3过程中发生了势能面交叉,ROOT=3得到的结果T3的能量小于T2且得到的结构以及能量均与ROOT=2得到的结果几乎一致


在拜读了sob老师的谈谈势能面交叉对激发态优化的影响 - 思想家公社的门口:量子化学·分子模拟·二次元 (sobereva.com)一文后,我做出了如下尝试:


1. 基于方法一优化后的结构,做tddft triplet nstates=5, root=3来优化
2. 跟踪优化路径,找到发生势能面交叉的那一步,选择那一步的下一步(已交叉)作为初态,做tddft triplet nstates=5, root=3来优化
3. 减小步长,追踪ROOT=2

但这三种尝试结果均得到的依旧是T2,也就是Excitd state=3(黄线)这条线的极小点。

请问老师们我该如何优化得到Excited state =2 (白线)这条势能面的极小点呢?

微信截图_20230728113300.png (29.84 KB, 下载次数 Times of downloads: 12)

势能面示意图

势能面示意图

这是选择root=2.png (54.29 KB, 下载次数 Times of downloads: 10)

这是选择root=2.png

这是选择root=3.png (51.03 KB, 下载次数 Times of downloads: 13)

这是选择root=3.png

19

帖子

0

威望

489

eV
积分
508

Level 4 (黑子)

5#
发表于 Post on 2024-8-31 10:42:49 | 只看该作者 Only view this author
cokie 发表于 2024-8-22 16:03
我把优化步长改到2解决的

谢谢回复 祝科研顺利

338

帖子

0

威望

1107

eV
积分
1445

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2024-8-22 16:03:16 | 只看该作者 Only view this author
本帖最后由 cokie 于 2024-8-22 16:04 编辑
jjm 发表于 2024-8-22 15:44
请问你后来这个问题是怎么解决的呢

我把优化步长改到2解决的

19

帖子

0

威望

489

eV
积分
508

Level 4 (黑子)

3#
发表于 Post on 2024-8-22 15:44:13 | 只看该作者 Only view this author
请问你后来这个问题是怎么解决的呢

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

2#
发表于 Post on 2023-7-28 16:01:50 | 只看该作者 Only view this author
首先,你确定excited state 2一定有极小值点吗?
激发态势能面本来就不一定存在极小值点
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

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

GMT+8, 2024-11-23 22:14 , Processed in 0.282473 second(s), 25 queries , Gzip On.

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