mutron 发表于 2018-5-13 17:07
留意到把坐标提取出来单独做freq居然是基态。无论怎么变动设置：opt=vtight、int=superfine、提取内坐标、g ...

谢谢您

mutron 发表于 2018-5-13 17:11
奇怪了，怎么最近传的gif都被重新编辑过导致动不了

因为直接上传的图会自动被加水印
压缩成压缩包后上传即可解决，但gif略大的时候还是建议用网盘以节约论坛空间

mutron 发表于 2018-5-13 17:07
留意到把坐标提取出来单独做freq居然是基态。无论怎么变动设置：opt=vtight、int=superfine、提取内坐标、g ...

看来涉及到O的时候一定要注意有没有势能面交叉

奇怪了，怎么最近传的gif都被重新编辑过导致动不了

留意到把坐标提取出来单独做freq居然是基态。无论怎么变动设置：opt=vtight、int=superfine、提取内坐标、geom=allcheck，ts优化结果有一个虚频，但是单独做freq还是基态


而且留意到freq（ultrafinegrid）的E（UB3LYP）=-2386.18773517，ts结果的E（UB3LYP）=-2386.15877802，很有可能freq（包括IRC）任务计算收敛到能量更低的态上面，导致题主IRC做的其实是基态


用multiwfn查看自旋密度：
     ts计算结果E（UB3LYP）=-2386.15877802
     提取坐标后的计算结果E（UB3LYP）=-2386.18773517，基本对应1楼IRC最高点的能量



解决办法就是加上guess=tcheck，结果就正常了：
     

liyuanhe211 发表于 2018-5-9 22:03
百度盘传一下opt和几个IRC的输入输出文件吧

链接：https://pan.baidu.com/s/1RtIr74YBcjxkCDQOQIuikQ 密码：saom

zc000 发表于 2018-5-9 17:24
李老师，我做过渡态的命令是#p opt=(calcall,tight,ts) int=ultrafine freq b3lyp/6-311g(d) geom=connec ...

百度盘传一下opt和几个IRC的输入输出文件吧

我本是个娃娃 发表于 2018-5-9 20:25
很有可能你的IRC曲线两侧走的是同一个路径，也就是说，不是对应的TS的两端，只是其中一端而已

老师，像我这种情况应该怎么处理？

很有可能你的IRC曲线两侧走的是同一个路径，也就是说，不是对应的TS的两端，只是其中一端而已

liyuanhe211 发表于 2018-5-9 16:52
“不精确”不是“不对”，和试了多少种没关系。比如优化时加opt=tight（DFT可能需要int=ultrafine）

李老师，我做过渡态的命令是#p opt=(calcall,tight,ts) int=ultrafine freq b3lyp/6-311g(d) geom=connectivity

你的IRC明显是错的，正常的IRC不可能有个尖峰
仔细阅读此文，确定该注意的都注意了。maxpoints文中有详细说明
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400

zc000 发表于 2018-5-9 16:48
我的这个体系比较简单，就三个原子，过渡态我试了很多种，然后对每种过渡态都进行了分析，发现出来的结果 ...

“不精确”不是“不对”，和试了多少种没关系。比如优化时加opt=tight（DFT可能需要int=ultrafine）

liyuanhe211 发表于 2018-5-9 15:46
考虑过渡态不精确地可能

我的这个体系比较简单，就三个原子，过渡态我试了很多种，然后对每种过渡态都进行了分析，发现出来的结果左右两边的结构基本一样。我就想是不是我的这个反应不成立？

考虑过渡态不精确地可能