计算化学公社

标题: 过渡态做完IRC后，优化得到的中间体比过渡态能量高 [打印本页]

作者
Author: radical 时间: 2022-1-7 09:30
标题: 过渡态做完IRC后，优化得到的中间体比过渡态能量高
得到过渡态后，读取chk文件做IRC,这上面ts的能量还是最高的，但读取irc最后一步做为初始构型去优化，最终得到的能量比过渡态还高，而且还发现明明irc最后一步的的能量是-724.44498996a.u.,而优化的第一步能量变成-724.42202488a.u.，明明是同样的构型不知道为什么能量差别这么大

过渡态关键词：opt=(calcfc,ts,noeigen) nosymm freq ub3lyp/6-311+g(d) optcyc=1000 guess=mix
IRC关键词：irc=(calcfc,maxpoint=50,LQA,forward) ub3lyp/6-311+g(d) geom=allcheck guess=read
优化中间体的关键词：opt freq ub3lyp/6-311+g(d) optcyc=1000 guess=mix nosymm

体系是双自由基反应

作者
Author: 喵星大佬 时间: 2022-1-7 09:42
stable=opt

作者
Author: radical 时间: 2022-1-7 09:56

喵星大佬发表于 2022-1-7 09:42
stable=opt

请问我可以读取优化得到的过渡态chk文件，写上stable=opt；还是要从头算

作者
Author: sobereva 时间: 2022-1-7 10:54

radical 发表于 2022-1-7 09:56
请问我可以读取优化得到的过渡态chk文件，写上stable=opt；还是要从头算

如果IRC最后一个点正好是你要优化的初始结构的那个点，优化时用guess=read读取IRC的chk保证能收敛到相同波函数
如果不是，用stable=opt得到当前结构最稳定波函数，之后优化时guess=read读取这个任务的chk

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