双自由基反应IRC关键词怎么写

radical

请问各位大神，一个加成了氢氧自由基的偶氮苯与另一个氢氧自由基发生脱氢反应，生成物都是单重态闭壳层，在对这个反应的过渡态做irc的时候，还要用ub3lyp和guess=mix吗，关键词应该怎么写呢

zjxitcc

不需要考虑这些问题，无论对于任何反应，跑IRC总是建议使用guess=read读取过渡态算好的波函数，否则过渡态波函数岂不算好了放那里浪费。

冰释之川

除了楼上提到的从TS读取初猜的建议，还得注意：算双自由基体系的IRC要将 forward 和reverse 方向分开跑（就是创建两个任务），不然就算跑IRC读取了TS波函数，当跑完一边的时候，另一边的初始波函数不会从*.chk文件中读取的。

radical

冰释之川 发表于 2021-12-21 16:13
除了楼上提到的从TS读取初猜的建议，还得注意：算双自由基体系的IRC要将 forward 和reverse 方向分开跑（就 ...

谢谢两位老师，没有像老师说的创两个任务跑果然曲线有问题，我现在重新跑一下

radical

冰释之川 发表于 2021-12-21 16:13
除了楼上提到的从TS读取初猜的建议，还得注意：算双自由基体系的IRC要将 forward 和reverse 方向分开跑（就 ...

老师我还有一个问题想问一下，请问guess=read与geom=check有什么区别，我想读取ts的chk文件做irc应该用哪个呢

zjxitcc

radical 发表于 2021-12-22 09:37
老师我还有一个问题想问一下，请问guess=read与geom=check有什么区别，我想读取ts的chk文件做irc应该用哪 ...

不矛盾，这些是不同的功能，guess=read读取轨道，geom=check读取几何结构，geom=allcheck读取电荷+自旋多重度+几何结构。做IRC计算最简单的输入文件，应该写guess=read geom=allcheck（当然，方法、基组也都要写），底下不需要写Title、电荷、自旋多重度、坐标等。

radical

zjxitcc 发表于 2021-12-22 09:46
不矛盾，这些是不同的功能，guess=read读取轨道，geom=check读取几何结构，geom=allcheck读取电荷+自旋多 ...

谢谢老师，那我的关键词这么写您看可以吗，irc=(calcfc,maxpoint=50,LQA,reverse) ub3lyp/6-311+g(d) guess=mix geom=allcheck guess=read

zjxitcc

radical 发表于 2021-12-22 15:15
谢谢老师，那我的关键词这么写您看可以吗，irc=(calcfc,maxpoint=50,LQA,reverse) ub3lyp/6-311+g(d) gue ...

不可以。guess=mix和guess=read混用，这是读进来以后再混合轨道，辛苦算好的过渡态波函数瞬间被搅乱。

没有必要用calcfc，因为绝大多数时候过渡态肯定做过频率分析计算（因为要看虚频），过渡态的chk文件里有力常数，直接读取即可，因此应该用rcfc，省时间。

radical

zjxitcc 发表于 2021-12-22 15:18
不可以。guess=mix和guess=read混用，这是读进来以后再混合轨道，辛苦算好的过渡态波函数瞬间被搅乱。

...

好的老师，我马上改掉，结合老师您之前在答疑专帖上的介绍，总算搞清楚这几个关键词的关系了

cheng1234

你好！麻烦问一下你在找过渡态的时候都加了ub3lyp和guess=mix吗

NisangPro

那这个例子的关键词就是这样了？
irc=(rcfc,maxpoint=50,LQA,reverse) ub3lyp/6-311+g(d) geom=allcheck guess=read