请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 76|回复: 6

[新手求助] irc两端结构优化收敛问题

[复制链接]

13

帖子

0

威望

65

eV
积分
78

Level 2 能力者

发表于 7 天前 | 显示全部楼层 |阅读模式
各位老师好:
  附件是我进行IRC计算得到附件中的曲线,其中的两端结构以及过渡态结构都符合预期,计算所用命令为:#p IRC=(CALCFC,STEPSIZE=5,MAXPOINT=25)  M062X/6-31g(d,p),采用相同方法和基组,取右端反应物结构进行几何优化时,出现SCF不收敛的情况,我的问题如下:  1:对irc曲线端点结构进行几何优化时,采用的方法与基组是否必须与IRC计算时相同?
  2:我按照帖子中的方法,尝试过加上int=ultrafine,或者更改为B3LYP泛函进行优化,能够解决收敛问题,请问这时的结构是否有意义,势能面是否不同了呢?



IRC.jpg

1万

帖子

25

威望

2万

eV
积分
40571

管理员

公社社长

发表于 6 天前 | 显示全部楼层
1 原理上不需要。但如果用的不同,道理上说不通
2 此时相当于和IRC不在一个势能面上,如果要把IRC的数据和极小点的数据拿到一起说事,说不通
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群1号:18616395,2号:466017436。超过4000人,用于交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

13

帖子

0

威望

65

eV
积分
78

Level 2 能力者

 楼主| 发表于 6 天前 | 显示全部楼层
sobereva 发表于 2019-4-16 00:15
1 原理上不需要。但如果用的不同,道理上说不通
2 此时相当于和IRC不在一个势能面上,如果要把IRC的数据和 ...

那请问老师,我可以不用IRC两端结构来优化吗,比如我重新再gview里构建反应物初始结构,再选用与IRC计算相同的理论方法进行优化,这样可以吗?但我感觉这样得到的势垒就不在一个反应路径上了。

29

帖子

0

威望

418

eV
积分
447

Level 3 能力者

发表于 6 天前 | 显示全部楼层
sobereva 发表于 2019-4-16 00:15
1 原理上不需要。但如果用的不同,道理上说不通
2 此时相当于和IRC不在一个势能面上,如果要把IRC的数据和 ...

那遇到2这种情况如何处理呢?(类似于楼主的情况,优化的右端反应物有虚频,加tight int=ultrafine 优化后解决。但和IRC就不在一个势能面上了)
a 先int=ultrafine 优化rc 再拿该结构去int=ultrafine优化,再和IRC数据做比较。(这种有时去int=ultrafine 能得到正常的结构,但有时仍有虚频)
b ts 更精确 所有的都在int=ultrafine 下做优化
我想了这两种,不确定能不能这样做。(但除了当下优化的这部分还有其他正常优化的数据,一起比较感觉还是有些不合理)

1万

帖子

25

威望

2万

eV
积分
40571

管理员

公社社长

发表于 6 天前 | 显示全部楼层
jing-quan63 发表于 2019-4-16 08:36
那请问老师,我可以不用IRC两端结构来优化吗,比如我重新再gview里构建反应物初始结构,再选用与IRC计算 ...

IRC两端的结构肯定比你自己随便摆的更接近于极小点,这么做没实际意义
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群1号:18616395,2号:466017436。超过4000人,用于交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

1万

帖子

25

威望

2万

eV
积分
40571

管理员

公社社长

发表于 6 天前 | 显示全部楼层
哈哈哈 发表于 2019-4-16 08:39
那遇到2这种情况如何处理呢?(类似于楼主的情况,优化的右端反应物有虚频,加tight int=ultrafine 优化 ...

找过渡态、走IRC也用ultrafine
如果需要和IRC上的点对比能量,至少优化完了重新用ultrafine算个单点
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群1号:18616395,2号:466017436。超过4000人,用于交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

13

帖子

0

威望

65

eV
积分
78

Level 2 能力者

 楼主| 发表于 6 天前 | 显示全部楼层
sobereva 发表于 2019-4-16 11:13
找过渡态、走IRC也用ultrafine
如果需要和IRC上的点对比能量,至少优化完了重新用ultrafine算个单点

好的,谢谢sober老师,非常感谢您解惑!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-4-22 21:17 , Processed in 0.154007 second(s), 26 queries .

快速回复 返回顶部 返回列表