计算化学公社

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

[新手求助] 寻找过渡态/柔性扫描中SCF不收敛求助

[复制链接 Copy URL]

175

帖子

0

威望

2136

eV
积分
2311

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
计算一个气体小分子在催化剂表面吸附后断键过程的过渡态。用QST2法生成的初猜结果不是很好,第一个结构的单点能都收敛失败;而用TS法根据末态手动摆的结构可能也不是太好,无法计算到过渡态。无奈下想用柔性扫描找个初猜。用PM6法计算第一个点就SCF不收敛;提高到HF/3-21G后,跑了10步左右再次出现SCF不收敛(共18步)。之前一般在结构很差时才会碰到SCF不收敛的情况。请问大家,在过渡态/柔性扫描过程中,碰到这样的情况,有什么好的办法吗?谢谢大家。
It's my life and it's now or never

6万

帖子

99

威望

5万

eV
积分
124706

管理员

公社社长

9#
发表于 Post on 2020-1-27 02:21:07 | 只看该作者 Only view this author
wrtgcn 发表于 2020-1-27 01:49
感谢您的回复。如果确实想要计算O3在材料表面吸附后断键过程的过渡态,您在方法或相应关键词上有什么建议 ...

用稍微像样的级别,诸如至少B3LYP/6-31G*。HF完全不能接受
解决SCF不收敛的做法我博文里都写了,没有可补充的

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
wrtgcn + 3 谢谢

查看全部评分 View all ratings

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

175

帖子

0

威望

2136

eV
积分
2311

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2020-1-27 01:49:17 | 只看该作者 Only view this author
sobereva 发表于 2020-1-26 18:54
O3属于很难算的体系,静态相关强,HF的结果是离谱的,要是再牵扯其断键,那事情就更复杂了。
SCF难收敛 ...

感谢您的回复。如果确实想要计算O3在材料表面吸附后断键过程的过渡态,您在方法或相应关键词上有什么建议吗?我在这一步上卡了好长时间了,一直没能找到。
It's my life and it's now or never

6万

帖子

99

威望

5万

eV
积分
124706

管理员

公社社长

7#
发表于 Post on 2020-1-26 18:54:33 | 只看该作者 Only view this author
wrtgcn 发表于 2020-1-26 06:39
谢谢社长。体系中没有过渡金属,但是气体分子是臭氧,有没有可能是其断键过程中的电子分布比较特殊,导致 ...

O3属于很难算的体系,静态相关强,HF的结果是离谱的,要是再牵扯其断键,那事情就更复杂了。
SCF难收敛跟这个不无关系
北京科音自然科学研究中心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

175

帖子

0

威望

2136

eV
积分
2311

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2020-1-26 06:39:24 | 只看该作者 Only view this author
sobereva 发表于 2020-1-25 23:16
所有可能的原因和解决办法在http://sobereva.com/61那篇文章里都提了。
没有具体文件没法给出确切的建议。 ...

谢谢社长。体系中没有过渡金属,但是气体分子是臭氧,有没有可能是其断键过程中的电子分布比较特殊,导致SCF难以收敛?
It's my life and it's now or never

175

帖子

0

威望

2136

eV
积分
2311

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2020-1-26 06:21:12 | 只看该作者 Only view this author
Machaophy 发表于 2020-1-25 15:09
HF/3-21G 计算的方法不可信吧,没有用过HF。。。。。。。。。。。。
套用Sob老师的话,最穷的级别B3LYP/6- ...

想用HF/3-21G扫描个靠谱的初始结构,再提高级别做过渡态计算。之前用B3LYP/6-31g*做柔性扫描,28个小时才跑了10个点不到,每个点的SCF收敛的很慢。
It's my life and it's now or never

6万

帖子

99

威望

5万

eV
积分
124706

管理员

公社社长

4#
发表于 Post on 2020-1-25 23:16:42 | 只看该作者 Only view this author
所有可能的原因和解决办法在http://sobereva.com/61那篇文章里都提了。
没有具体文件没法给出确切的建议。
HF/3-21G算了也白算
PM6都SCF不收敛,有较大可能输入文件就存在显著硬伤。如果催化剂里有过渡金属,改用xtb靠谱得多,Gaussian结合xtb找过渡态又快又容易,见《将Gaussian与Grimme的xtb程序联用搜索过渡态、产生IRC、做振动分析》(http://sobereva.com/421
北京科音自然科学研究中心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

15

帖子

0

威望

70

eV
积分
85

Level 2 能力者

3#
发表于 Post on 2020-1-25 15:09:31 | 只看该作者 Only view this author
HF/3-21G 计算的方法不可信吧,没有用过HF。。。。。。。。。。。。
套用Sob老师的话,最穷的级别B3LYP/6-31g*,这个都用不起就别做量化计算了
推荐M06-2X,这个找过渡态性能好些
你不贴文件出来,不知道你怎么做的呀

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
wrtgcn + 1 谢谢

查看全部评分 View all ratings

205

帖子

0

威望

1730

eV
积分
1935

Level 5 (御坂)

2#
发表于 Post on 2020-1-25 12:27:44 | 只看该作者 Only view this author
解决SCF不收敛问题的方法  http://sobereva.com/61  看完基本搞定

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
wrtgcn + 1 谢谢

查看全部评分 View all ratings

Shanghai JiaoTong University
欢迎关注我的推送号:AI计算生化

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

GMT+8, 2026-1-25 19:51 , Processed in 0.436385 second(s), 22 queries , Gzip On.

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