计算化学公社

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

[Gaussian/gview] 关于IRC同时用lqa后的结果是否正确求助各位老师

[复制链接 Copy URL]

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

我找到了一个过渡态,有一个虚频,振动强度为-142,振动模式是我想要的结果,但是在做IRC分析时总是出现一个突变点,结果报错无法收敛(Delta-x Convergence NOT Met),除了此突变点reverse和forward趋势都是正确的。(forward、reverse两个方向分别计算)
后来我在论坛上看到可以加上lqa的命令,重新计算正常结束,找到了反应物和产物,但是路径分析的图中出现很多拐点,不平滑。
这样可以认为正确吗?
我又用同样的方法计算了另外一个我得到的过渡态,有的拐点甚至高于初始点,是不是我找的这个过渡态存在问题?谢谢大家!

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2017-8-24 15:49:08 | 只看该作者 Only view this author
用stepsize减小步长(步长越小越平滑),或者用recalc=N设定每N步重新精确计算一次Hessian(N越小越平滑但也越耗时)
北京科音自然科学研究中心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

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2017-8-24 15:53:59 | 只看该作者 Only view this author
sobereva 发表于 2017-8-24 15:49
用stepsize减小步长(步长越小越平滑),或者用recalc=N设定每N步重新精确计算一次Hessian(N越小越平滑但 ...

谢谢sob老师,我尝试一下减小步长。老师,我得到的第二个图那样的结果算不算是过渡态存在问题呀?我需要重新找吗?

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

4#
发表于 Post on 2017-8-24 16:19:35 | 只看该作者 Only view this author
chrishuo 发表于 2017-8-24 15:53
谢谢sob老师,我尝试一下减小步长。老师,我得到的第二个图那样的结果算不算是过渡态存在问题呀?我需要 ...

这曲线看着像数值问题,先做平滑了再说。

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2017-8-24 16:24:36 | 只看该作者 Only view this author
liyuanhe211 发表于 2017-8-24 16:19
这曲线看着像数值问题,先做平滑了再说。

好的,谢谢您!

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2017-8-24 19:50:36 | 只看该作者 Only view this author
liyuanhe211 发表于 2017-8-24 16:19
这曲线看着像数值问题,先做平滑了再说。

我针对第二个图减小了步长(之前的步长是30),#p irc=(forward,calcfc,LQA,stepsize=5) b3lyp/6-31g(d)任务正常结束,输出文件中出现了This type of calculation cannot be archived.(我看您的之前的帖子说这不是报错)正常结束。但是我在查看路径时发现只有两个点,这样可以吗?我应该怎么办?谢谢您!

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

7#
发表于 Post on 2017-8-24 20:22:43 | 只看该作者 Only view this author
chrishuo 发表于 2017-8-24 19:50
我针对第二个图减小了步长(之前的步长是30),#p irc=(forward,calcfc,LQA,stepsize=5) b3lyp/6-31g(d) ...

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

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2017-8-25 20:35:16 | 只看该作者 Only view this author

好的,我再试试看,谢谢老师!

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2017-8-26 20:04:17 | 只看该作者 Only view this author

sob老师,我用这样的命令#p opt=(ts,calcall,noeigentest,tight) int=ultrafine freq b3lyp/6-31g(d)重新计算了过渡态,对得到的过渡态进行IRC分析时,forward(给的命令是#p irc=(forward,calcfc,recalc=3,LQA,stepsize=10) int=ultrafine b3lyp/6-31g(d))方向仍旧只有两个点,reverse方向得到了平滑的曲线。我接下来应该怎么办?是不是说明我的过渡态存在问题?谢谢老师

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

10#
发表于 Post on 2017-8-26 20:14:41 | 只看该作者 Only view this author
chrishuo 发表于 2017-8-26 20:04
sob老师,我用这样的命令#p opt=(ts,calcall,noeigentest,tight) int=ultrafine freq b3lyp/6-31g(d)重新 ...


stepsize设到15或20,如果能走下去但是还不平滑,把LQA去掉试试(由于用了recalc=3了,重校正步不收敛的错误几率已经降低不少)看是否已经足够平滑,如果不行,再把recalc设2乃至1(等价于calcall)。
或者尝试用GS2算法产生IRC
北京科音自然科学研究中心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

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

11#
 楼主 Author| 发表于 Post on 2017-8-26 20:38:59 | 只看该作者 Only view this author
sobereva 发表于 2017-8-26 20:14
stepsize设到15或20,如果能走下去但是还不平滑,把LQA去掉试试(由于用了recalc=3了,重校正步不收敛 ...

好的,谢谢老师

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

12#
 楼主 Author| 发表于 Post on 2017-8-29 09:20:51 | 只看该作者 Only view this author
sobereva 发表于 2017-8-26 20:14
stepsize设到15或20,如果能走下去但是还不平滑,把LQA去掉试试(由于用了recalc=3了,重校正步不收敛 ...

sob老师,我在用#p opt=(ts,calcall,noeigentest,tight) int=ultrafine freq b3lyp/6-31g(d)对已经找到的过渡态重新计算时,计算了将近4天还没有结束,我用chemcraft看能量已经维持在一个固定值优化多次没有变化,遇到这样的问题我应该怎么处理?谢谢老师!

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

13#
发表于 Post on 2017-8-29 09:48:22 | 只看该作者 Only view this author
chrishuo 发表于 2017-8-29 09:20
sob老师,我在用#p opt=(ts,calcall,noeigentest,tight) int=ultrafine freq b3lyp/6-31g(d)对已经找到的 ...

用gview看,甭用chemcraft,直接就能给你能量和受力变化,看有无收敛的迹象。
calcall显然慢。用了calcall就没有太大必要用tight,这导致难收敛,而且也没必要用ultrafine
北京科音自然科学研究中心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

32

帖子

0

威望

83

eV
积分
115

Level 2 能力者

14#
 楼主 Author| 发表于 Post on 2017-8-29 09:58:02 | 只看该作者 Only view this author
本帖最后由 chrishuo 于 2017-8-29 10:10 编辑
sobereva 发表于 2017-8-29 09:48
用gview看,甭用chemcraft,直接就能给你能量和受力变化,看有无收敛的迹象。
calcall显然慢。用了calca ...

老师,我用gaussview看了一下,发现能量也不再变化了

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 06:48 , Processed in 0.212496 second(s), 23 queries , Gzip On.

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