计算化学公社

标题: 求助CO在类似石墨结构上脱附IRC过程的产物优化到了反应物 [打印本页]

作者
Author: wch990122 时间: 2021-7-13 10:56
标题: 求助CO在类似石墨结构上脱附IRC过程的产物优化到了反应物
各位老师好，本人在分析一个脱附过程时遇到了问题，具体是CO在类似石墨结构上脱附，如图1。
首先根据感觉进行过渡态的判断，但始终得不到正确的过渡态。于是对键长做了柔性扫描，发现存在最大值，如图2，于是将最大值结构设置成初猜的过渡态，得到的过渡态有且至于一个虚频，并且振动方向也是正确的，进一步进行IRC分析，因为默认步长始终是走了几步就停止，按照相关帖子将步长放大到30，得到IRC曲线如图3，反应物和产物都是我想要的，但是对产物进行进一步的优化时，却得到了反应物，不知为何所以特来求救。
扫描，过渡态，IRC的输入文件如下：

%nprocshared=4
%mem=8GB
%chk=scan2.chk
#p opt=modredundant ub3lyp/6-31g(d) nosymm

.......
12 15 s 20 0.1

%nprocshared=4
%mem=8GB
%chk=TS3.chk
#p opt=(TS,calcfc,noeigen) freq ub3lyp/6-31g(d) em=gd3

%nprocshared=4
%mem=8GB
%chk=IRC3.chk
#p irc=(calcfc,maxpoint=150,stepsize=30,lqa) ub3lyp/6-31g(d)

图片：
(, 下载次数 Times of downloads: 9) (, 下载次数 Times of downloads: 20) (, 下载次数 Times of downloads: 18)

作者
Author: wch990122 时间: 2021-7-13 11:07
(, 下载次数 Times of downloads: 7)

作者
Author: sobereva 时间: 2021-7-13 17:47
IRC加上recalc=3，然后再对末端优化。而且优化时opt加上maxstep=5,notrust

还有这问题的话，改用M06-2X重新算一遍

欢迎光临计算化学公社 (http://bbs.keinsci.com/)