计算化学公社

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

[新手求助] IRC为何一步还没有跑完就正常结束了?

[复制链接 Copy URL]

146

帖子

0

威望

1319

eV
积分
1465

Level 4 (黑子)

一个活泼质子氢的转移反应,为何跑IRC的时候作业正常结束了,但一步都没有跑完?附件是输出文件,还请各位前辈给点建议,是势垒太小的原因么?还是哪个参数设置的不合适?

ts1_RO_Ar_irc_p.log

1.21 MB, 下载次数 Times of downloads: 29

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

2#
发表于 Post on 2016-10-10 20:28:09 | 只看该作者 Only view this author
胡说,分明跑了30步:


146

帖子

0

威望

1319

eV
积分
1465

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2016-10-10 21:09:45 | 只看该作者 Only view this author
liyuanhe211 发表于 2016-10-10 20:28
胡说,分明跑了30步:

不好意思啊李老师,是我言不达意了,

我是想问这个react coord为何等于问号?是哪里出问题了?这个是用chemcraft打开Log文件后显示的结果

3097

帖子

29

威望

1万

eV
积分
17221

Level 6 (一方通行)

4#
发表于 Post on 2016-10-10 22:00:44 | 只看该作者 Only view this author
珊珊来迟 发表于 2016-10-10 21:09
不好意思啊李老师,是我言不达意了,

我是想问这个react coord为何等于问号?是哪里出问题了?这个是 ...

用其他软件尝试。

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

5#
发表于 Post on 2016-10-11 01:27:12 | 只看该作者 Only view this author
明明都已经正常跑完了啊,还想怎么跑,gview看到的IRC曲线也很正常
Maximum number of steps reached.
Calculation of FORWARD path complete.
Reaction path calculation complete.

Energies reported relative to the TS energy of       -1076.212813
--------------------------------------------------------------------------
    Summary of reaction path following
--------------------------------------------------------------------------
                        Energy    RxCoord
   1                    0.00000   0.00000
   2                   -0.00017   0.14107
   3                   -0.00154   0.28600
   4                   -0.00543   0.43270
   5                   -0.01148   0.57951
   6                   -0.01783   0.72626
   7                   -0.02225   0.87196
   8                   -0.02433   1.01289
   9                   -0.02568   1.15820
  10                   -0.02665   1.30405
  11                   -0.02745   1.45035
  12                   -0.02812   1.59692
  13                   -0.02870   1.74348
  14                   -0.02920   1.89011
  15                   -0.02964   2.03677
  16                   -0.03004   2.18340
  17                   -0.03040   2.33015
  18                   -0.03073   2.47687
  19                   -0.03103   2.62353
  20                   -0.03130   2.77026
  21                   -0.03156   2.91699
  22                   -0.03179   3.06379
  23                   -0.03200   3.21053
  24                   -0.03219   3.35720
  25                   -0.03237   3.50391
  26                   -0.03253   3.65062
  27                   -0.03267   3.79740
  28                   -0.03280   3.94413
  29                   -0.03291   4.09079
  30                   -0.03302   4.23752
  31                   -0.03311   4.38425
北京科音自然科学研究中心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

146

帖子

0

威望

1319

eV
积分
1465

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2016-10-11 09:24:20 | 只看该作者 Only view this author
sobereva 发表于 2016-10-11 01:27
明明都已经正常跑完了啊,还想怎么跑,gview看到的IRC曲线也很正常
Maximum number of steps reached.
Ca ...

恩恩,谢谢sob老师,做完IRC之后不知道怎么正确的分析结果,一直纠结于上个用chemcraft的截图显示的信息,这下用GVIEW看到IRC曲线了,谢谢老师啦

146

帖子

0

威望

1319

eV
积分
1465

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2016-10-11 09:30:37 | 只看该作者 Only view this author
sobereva 发表于 2016-10-11 01:27
明明都已经正常跑完了啊,还想怎么跑,gview看到的IRC曲线也很正常
Maximum number of steps reached.
Ca ...


另外sob老师,还有一个小问题请教您,请问文献中这种反应路径图是用什么软件做的呢?

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

8#
发表于 Post on 2016-10-11 13:12:28 | 只看该作者 Only view this author
珊珊来迟 发表于 2016-10-11 09:30
另外sob老师,还有一个小问题请教您,请问文献中这种反应路径图是用什么软件做的呢?

在Origin中绘制能量折线图的方法
http://sobereva.com/320
北京科音自然科学研究中心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

146

帖子

0

威望

1319

eV
积分
1465

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2016-10-11 18:02:55 | 只看该作者 Only view this author
sobereva 发表于 2016-10-11 13:12
在Origin中绘制能量折线图的方法
http://sobereva.com/320

太感谢老师啦!!!

146

帖子

0

威望

1319

eV
积分
1465

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2016-10-12 22:47:15 | 只看该作者 Only view this author
sobereva 发表于 2016-10-11 01:27
明明都已经正常跑完了啊,还想怎么跑,gview看到的IRC曲线也很正常
Maximum number of steps reached.
Ca ...

sob老师,产物方向IRC正常结束了,可为什么产物方向IRC报错,然后结果是这样的?第二个的能量突然降到这么低?附件是我的log文件,还请老师给点建议,麻烦您啦

ts1_RO_Ar_irc_r.log

1.19 MB, 下载次数 Times of downloads: 6

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

11#
发表于 Post on 2016-10-12 22:57:53 | 只看该作者 Only view this author
珊珊来迟 发表于 2016-10-12 22:47
sob老师,产物方向IRC正常结束了,可为什么产物方向IRC报错,然后结果是这样的?第二个的能量突然降到这 ...

IRC里写LQA,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

146

帖子

0

威望

1319

eV
积分
1465

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2016-10-12 23:05:29 | 只看该作者 Only view this author
sobereva 发表于 2016-10-12 22:57
IRC里写LQA,IRC正常跑完了看到的就正常了。

恩恩,谢谢老师!

146

帖子

0

威望

1319

eV
积分
1465

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2016-10-13 18:39:21 | 只看该作者 Only view this author
sobereva 发表于 2016-10-12 22:57
IRC里写LQA,IRC正常跑完了看到的就正常了。

sob老师您好,加了LQA关键词后IRC正常结束了。太谢谢您啦!
老师,我是要比较A ,B两条路径哪个势垒比较低,现在IRC跑完了,证明两个过渡态找的是正确的,但其中PathB的过渡态△U和△H都为负值,比反应物还要低,但△G是比反应物高的,请问老师这种情况怎么处理呢?

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

14#
发表于 Post on 2016-10-13 18:42:55 | 只看该作者 Only view this author
珊珊来迟 发表于 2016-10-13 18:39
sob老师您好,加了LQA关键词后IRC正常结束了。太谢谢您啦!
老师,我是要比较A ,B两条路径哪个势垒比较 ...

只要你的计算过程、数据读取都确保无误,自由能垒为正,就说明是严格意义的过渡态。H和U就甭管它了
北京科音自然科学研究中心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

146

帖子

0

威望

1319

eV
积分
1465

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2016-10-13 18:52:01 | 只看该作者 Only view this author
sobereva 发表于 2016-10-13 18:42
只要你的计算过程、数据读取都确保无误,自由能垒为正,就说明是严格意义的过渡态。H和U就甭管它了

有没有可能是因为找到的过渡态结构不够准确?

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 05:11 , Processed in 0.188624 second(s), 23 queries , Gzip On.

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