计算化学公社

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

[ORCA] ORCA优化激发态能量出现极大跳变的结果是否可信?

[复制链接 Copy URL]

245

帖子

0

威望

780

eV
积分
1025

Level 4 (黑子)

各位老师好,我在ORCA上进行了一个双核金属有机配合物的S1激发态优化(PBE0-D3BJ,TDDFT,def2-SV(P)/def2-TZVP(-f)基组和对应的辅助基组),从优化过的S0构型开始优化(优化到第11次时机器断电退出,之后重新读取最后一帧的坐标和波函数继续优化),到某一帧突然出现了能量的巨大跳变,结构上没有明显不合理之处。请问各位老师出现这样的情况(还开始有一点震荡)时是否还应该继续跑下去?
out文件如下。 1a-S1.zip (2.41 MB, 下载次数 Times of downloads: 6)

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

2#
发表于 Post on 2023-3-29 22:16:00 | 只看该作者 Only view this author
可能是基态SCF突然收敛到了另一个波函数。
检查能量跳变前后的波函数是否定性符合,如果有定性区别,检查哪个更符合化学常识
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
本团队长期招收研究生,有意者可私信联系

245

帖子

0

威望

780

eV
积分
1025

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2023-3-29 22:53:00 | 只看该作者 Only view this author
wzkchem5 发表于 2023-3-29 22:16
可能是基态SCF突然收敛到了另一个波函数。
检查能量跳变前后的波函数是否定性符合,如果有定性区别,检查 ...

从原子坐标来看没有特别明显的差别。但我注意到在跳变前1步出现warning“ SERIOUS PROBLEM WITH INTERNALS - ANGLE IS APPROACHING 180 OR 0 DEGREES!! REBUILDING A NEW SET OF INTERNALS ” 是不是因为这个导致重新载入优化循环的结构已经不对了?

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

4#
发表于 Post on 2023-3-30 00:03:43 | 只看该作者 Only view this author
KSeGaSn 发表于 2023-3-29 15:53
从原子坐标来看没有特别明显的差别。但我注意到在跳变前1步出现warning“ SERIOUS PROBLEM WITH INTERNAL ...

那检查一下这一步前后结构是不是发生了跳变。可能虽然前后结构都合理,但是彼此差别比较大
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
本团队长期招收研究生,有意者可私信联系

245

帖子

0

威望

780

eV
积分
1025

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2023-3-30 09:02:41 | 只看该作者 Only view this author
wzkchem5 发表于 2023-3-30 00:03
那检查一下这一步前后结构是不是发生了跳变。可能虽然前后结构都合理,但是彼此差别比较大

我选择了跳变前最后一步和跳变后两步的结构进行比较,RMSD~0.04,从图像上来看也几乎完全重叠了。 详细问题是“dihedral-internal angle (14, 5, 3, 0) reaching linearity!  179.4317”,看起来也是一个很正常的情况。
如果我不放心的话,是不是最好读取跳变前的结构重新载入优化呢?毕竟在这个warning之后的经历了原子坐标的重建?

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

6#
发表于 Post on 2023-3-30 14:56:20 | 只看该作者 Only view this author
KSeGaSn 发表于 2023-3-30 02:02
我选择了跳变前最后一步和跳变后两步的结构进行比较,RMSD~0.04,从图像上来看也几乎完全重叠了。详细问 ...

可以试一下
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
本团队长期招收研究生,有意者可私信联系

245

帖子

0

威望

780

eV
积分
1025

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2023-3-30 22:18:18 | 只看该作者 Only view this author

好的老师,我去试试看吧。另外老师我在对优化好的基态单独进行频率计算时,CP-SCF计算中分成了14个IBatch,前面都成功收敛了,结果进行到IBatch 11的第七个Iteration突然无报错不再更新out文件,CPU占用也回到0了,请问这是什么问题呢?也是这个分子,所以计算量特别大,我还不敢随便kill掉。

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

8#
发表于 Post on 2023-3-30 23:37:38 | 只看该作者 Only view this author
KSeGaSn 发表于 2023-3-30 15:18
好的老师,我去试试看吧。另外老师我在对优化好的基态单独进行频率计算时,CP-SCF计算中分成了14个IBatch ...

可能是硬盘写满了,或者刚好遇到了偶然发生的硬件问题
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-22 06:28 , Processed in 0.171121 second(s), 23 queries , Gzip On.

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