计算化学公社

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

[新手求助] 高斯寻找过渡态几何不收敛(已经尝试过社长帖子的大部分方法)该怎么办?

[复制链接 Copy URL]

22

帖子

0

威望

919

eV
积分
941

Level 4 (黑子)

学生在寻找五价铁和小分子污染的过渡态中,遇到了几何不收敛的问题,基于社长的帖子http://sobereva.com/164在之前的过渡态寻找过程中得到了极大的帮助,但是在这计算里面,尝试了许多方法(包括gdiis,maxstep=1都用过了,笛卡尔坐标等....)都无法收敛到我想要的的过渡态(铁氧中的氧被小分子污染物拔走),因此来发第一次贴求助,恳求老师的帮忙。先谢过大家~

5H3TSd-gdiis.rar

654.59 KB, 下载次数 Times of downloads: 17

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

2#
发表于 Post on 2020-12-5 11:48:04 | 只看该作者 Only view this author
先用opt=loose 试试看吧

5H3TSd-new.gjf

1.67 KB, 下载次数 Times of downloads: 26

Stand on the shoulders of giants

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

3#
发表于 Post on 2020-12-5 12:36:32 | 只看该作者 Only view this author
本帖最后由 liyuanhe211 于 2020-12-5 12:40 编辑

含铁的注意检查一下每一步的电子结构对不对,我遇到过优化过程中它震荡的两个结构其实是在两个不同的电子态之间不断切换。

既然都能接受calcall的耗时了,可以看看每一步的虚频是不是对的,你的虚频数值比较小。

可以做做柔性扫描,不只是键长,还包括从什么方向接近之类的。

1万

帖子

0

威望

9876

eV
积分
22118

Level 6 (一方通行)

4#
发表于 Post on 2020-12-5 12:37:20 | 只看该作者 Only view this author
没准这个反应就是没有过渡态的,活化能为0。
可以做个柔性扫描,看看是不是真的没有势垒
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
本团队长期招收研究生,有意者可私信联系

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

5#
发表于 Post on 2020-12-5 13:05:40 | 只看该作者 Only view this author
liyuanhe211 发表于 2020-12-5 12:36
含铁的注意检查一下每一步的电子结构对不对,我遇到过优化过程中它震荡的两个结构其实是在两个不同的电子态 ...

我在想是不是应该对波函数稳定性检测一下……
Stand on the shoulders of giants

22

帖子

0

威望

919

eV
积分
941

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2020-12-5 13:24:01 | 只看该作者 Only view this author
liyuanhe211 发表于 2020-12-5 12:36
含铁的注意检查一下每一步的电子结构对不对,我遇到过优化过程中它震荡的两个结构其实是在两个不同的电子态 ...

学生在气相中能找出来,虚频也合适,一加隐性水就跑不出来了(虽然水的影响肯定会很大,但总是觉得气相中能找出来在水中应该也能找出来)。谢谢老师,我用柔性扫描试试

22

帖子

0

威望

919

eV
积分
941

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2020-12-5 13:24:33 | 只看该作者 Only view this author
wzkchem5 发表于 2020-12-5 12:37
没准这个反应就是没有过渡态的,活化能为0。
可以做个柔性扫描,看看是不是真的没有势垒

谢谢您的回复,学生待会试试

22

帖子

0

威望

919

eV
积分
941

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2020-12-5 13:30:12 | 只看该作者 Only view this author
冰释之川 发表于 2020-12-5 11:48
先用opt=loose 试试看吧

直接上loose吗,学生先试一下,谢谢您~

6万

帖子

99

威望

6万

eV
积分
125150

管理员

公社社长

9#
发表于 Post on 2020-12-6 14:56:57 | 只看该作者 Only view this author
Vincentears 发表于 2020-12-5 13:24
学生在气相中能找出来,虚频也合适,一加隐性水就跑不出来了(虽然水的影响肯定会很大,但总是觉得气相中 ...

对于反应物或产物带显著电荷的情况,极性溶剂环境可能极为显著地影响势能面,甚至决定过渡态能否出现。例如

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

1万

帖子

0

威望

9876

eV
积分
22118

Level 6 (一方通行)

10#
发表于 Post on 2020-12-6 17:10:24 | 只看该作者 Only view this author
冰释之川 发表于 2020-12-5 13:05
我在想是不是应该对波函数稳定性检测一下……

同意,不过波函数稳定性检测也不够,含铁体系经常有两个以上的稳定波函数。。。
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
本团队长期招收研究生,有意者可私信联系

22

帖子

0

威望

919

eV
积分
941

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2020-12-10 10:52:26 | 只看该作者 Only view this author
liyuanhe211 发表于 2020-12-5 12:36
含铁的注意检查一下每一步的电子结构对不对,我遇到过优化过程中它震荡的两个结构其实是在两个不同的电子态 ...

老师,再请教一下您,如何确定优化有没有收敛到正确的电子态?

1万

帖子

0

威望

9876

eV
积分
22118

Level 6 (一方通行)

12#
发表于 Post on 2020-12-10 12:09:49 | 只看该作者 Only view this author
Vincentears 发表于 2020-12-10 10:52
老师,再请教一下您,如何确定优化有没有收敛到正确的电子态?

你做计算之前必然是本来就知道正确的电子态是哪个的,比如铁是高自旋还是低自旋,其他的一些带自旋的分子和铁是铁磁耦合还是反铁磁耦合。所以一般看一下自旋布居是不是和你的预期相符就行了,有时还需要作出前线轨道的图。
如果不知道哪个电子态是正确的,那就穷举出所有可能的电子态,分别找到这些电子态,取能量最低的那个。
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-23 05:41 , Processed in 0.172806 second(s), 23 queries , Gzip On.

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