计算化学公社

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

[新手求助] IRC计算始终不完全

[复制链接 Copy URL]

355

帖子

0

威望

1917

eV
积分
2272

Level 5 (御坂)

在M062X/6-311+g(2d,p)下算IRC,结果只有两个点
把基组换成cc-pVTZ,def2-TZVP重新优化再算IRC结果也一样,加int=superfine或者irc=tight也没用
但是看虚频的振动模式感觉像是对的
之后也试过换理论级别,PM6算下来64个点(stepsize=5),有一个反应过程的雏形,但是也不完整,感觉还没到反应物或者产物就停了

IRC.rar

170.31 KB, 下载次数 Times of downloads: 8

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

2#
发表于 Post on 2021-11-22 05:51:46 | 只看该作者 Only view this author
仔细看
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400http://bbs.keinsci.com/thread-7736-1-1.html

加“int=superfine或者irc=tight”本来就没用
最简单粗暴做法是加大stepsize

IRC再怎么跑也不可能精确跑到极小点去,搞清楚原理

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

355

帖子

0

威望

1917

eV
积分
2272

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2021-11-23 08:14:52 | 只看该作者 Only view this author
sobereva 发表于 2021-11-22 05:51
仔细看
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400(http://bbs.keinsci.com/thread- ...

S老师,我看了教程,把stepsize调成20,但是跑HPC迭代不收敛,GS2虽然满足收敛标准还是报告不收敛,LQA总算算出来了,但是还是不完整,尤其是右边,明显不平坦,不知为何能够达到收敛限我看教程,这个是因为鞍点不精确导致的问题吗?我重新优化了结构,opt=verytight,IRC计算下来还是只有3个点


Snipaste_2021-11-23_08-01-50.png (20.16 KB, 下载次数 Times of downloads: 37)

Snipaste_2021-11-23_08-01-50.png

Snipaste_2021-11-23_08-12-37.png (24.54 KB, 下载次数 Times of downloads: 39)

Snipaste_2021-11-23_08-12-37.png

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

4#
发表于 Post on 2021-11-23 08:19:12 | 只看该作者 Only view this author
mfdsrax2 发表于 2021-11-23 08:14
S老师,我看了教程,把stepsize调成20,但是跑HPC迭代不收敛,GS2虽然满足收敛标准还是报告不收敛,LQA总 ...

没必要跑完整就不需要跑完整。搞清楚跑IRC的目的是什么
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

355

帖子

0

威望

1917

eV
积分
2272

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2021-11-23 08:44:30 | 只看该作者 Only view this author
sobereva 发表于 2021-11-23 08:19
没必要跑完整就不需要跑完整。搞清楚跑IRC的目的是什么

IRC是确定过渡态是否连接对应的反应物和产物吧,没有跑完怎么确定呢

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

6#
发表于 Post on 2021-11-23 09:07:13 | 只看该作者 Only view this author
mfdsrax2 发表于 2021-11-23 08:44
IRC是确定过渡态是否连接对应的反应物和产物吧,没有跑完怎么确定呢

对两端优化啊

何况只观看过渡态开始的IRC轨迹走的趋势,从肉眼和直觉上也自然能判断出TS连通的是什么结构

好好把此文看了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400http://bbs.keinsci.com/thread-7736-1-1.html

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

355

帖子

0

威望

1917

eV
积分
2272

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2021-11-23 10:21:36 | 只看该作者 Only view this author
sobereva 发表于 2021-11-23 09:07
对两端优化啊

何况只观看过渡态开始的IRC轨迹走的趋势,从肉眼和直觉上也自然能判断出TS连通的是什么 ...

S老师,我认真看了您的文章,知道IRC只是一个定性分析的工具,但是现在项目需要高质量的IRC动画达到演示目的,所以无论如何都想把这个IRC做好,现在这个LQA+stepsize=20的动画很不连贯,长度也不够,无法达到要求,所以来请教您问题出在哪里。文章里面写的可能原因是过渡态优化精度不够,IRC步长太大,没有精确计算Hessian等原因,我已经尽力去排除了,最后都用上m062x/cc-pvtz stepsize=3 calcall int=superfine,仍然得不到好的IRC图,麻烦您再指点指点

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

8#
发表于 Post on 2021-11-23 10:26:04 | 只看该作者 Only view this author
mfdsrax2 发表于 2021-11-23 10:21
S老师,我认真看了您的文章,知道IRC只是一个定性分析的工具,但是现在项目需要高质量的IRC动画达到演示 ...

LQA+stepsize=20的时候加上calcall
IRC往往无法跑得特别完整是Gaussian自身判断结束标准的缺陷,不改源代码没法解决
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

355

帖子

0

威望

1917

eV
积分
2272

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2021-11-23 14:49:36 | 只看该作者 Only view this author
sobereva 发表于 2021-11-23 10:26
LQA+stepsize=20的时候加上calcall
IRC往往无法跑得特别完整是Gaussian自身判断结束标准的缺陷,不改源 ...

跑这段IRC的时候已经用了calcall

再请问一下,还有其他软件能够算IRC吗

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

10#
发表于 Post on 2021-11-24 07:35:43 | 只看该作者 Only view this author
mfdsrax2 发表于 2021-11-23 14:49
跑这段IRC的时候已经用了calcall

再请问一下,还有其他软件能够算IRC吗

Gaussian已经是做的最好的了

没尝试其它程序的必要
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

355

帖子

0

威望

1917

eV
积分
2272

Level 5 (御坂)

11#
 楼主 Author| 发表于 Post on 2021-11-28 17:20:11 | 只看该作者 Only view this author
sobereva 发表于 2021-11-24 07:35
Gaussian已经是做的最好的了

没尝试其它程序的必要

S老师,正好借着这个帖子,问一下论坛附件的问题。我在编辑之前帖子的时候,根本没有编辑图片和删除附件的选项,只能上传新图片,可能你的权限比较高可以编辑,但是我就不行,换了好几个浏览器都是一样的,这应该是论坛的BUG。





5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

12#
发表于 Post on 2021-12-2 05:58:42 | 只看该作者 Only view this author
mfdsrax2 发表于 2021-11-28 17:20
S老师,正好借着这个帖子,问一下论坛附件的问题。我在编辑之前帖子的时候,根本没有编辑图片和删除附件 ...

跟权限没关系。不是bug。用搜狗浏览器,即便是普通用户也能修改之前上传的图片。
提问时候别借帖子问,论坛功能的问题发到公社大厅板块问
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 06:21 , Processed in 0.200974 second(s), 29 queries , Gzip On.

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