用高斯，在找到过渡态之后，怎么快速地找到相应的反应物和产物？

michaelm

请教大家一下，用高斯，在找到过渡态之后，怎么快速地找到相应的反应物和产物。
我的作法是这样的：
1）用过渡态做IRC，同时往前往后各做一些。一般我会把步长设大些，步数设少些。
2）然后用得到的在IRC里的最前和最后的结构，分别再优化，得到相应的反应物和产物结构。

这个好像略有些麻烦，两步不能自动进行，必须在运行完第一步之后，手工提取结构再做第二步的两个计算。

我想请教一下大家，有没有更简单的方法；或者高斯里有没有简单的能直接读取chk文件里的irc的某个结构？你们一般都是怎么做的呢？

谢谢。

sobereva

优化时候用geom=check，则从IRC的chk里会读取IRC任务最后一个点的结构。对于双方向跑，这样读取的是最接近反应物的那个IRC点的结构，但是最接近产物那端的结构没法读取。做法只能是自己写个程序，从输出文件中或者从fch文件中提取出来
代码可参考
将Gaussian的IRC任务输出转换为.xyz轨迹文件的工具
http://sobereva.com/285

michaelm

谢谢！
从irc里读取的是irc最后一个点，那就足够了。把双方向分成单方跑两次，配合oldchk，就可以解决问题了。