计算化学公社

标题: TS计算正常收敛，虚频正常；但是IRC曲线错误的原因 [打印本页]

作者
Author: chenxin199261 时间: 2014-11-7 20:11
标题: TS计算正常收敛，虚频正常；但是IRC曲线错误的原因
我用TS对一个过渡态进行优化，结果收敛，检查虚频，虚频方向正确。
我用这个过渡态结果计算IRC。但结果显然是错误的。通过看动画发现，原子向反应物方向走几步，又反返回了。
IRC输入文件输入如下

################################################
%nprocshared=24
%chk=FeONH36-TS-5-IRC-p.chk
# irc=(maxpoints=10,calcall) ub3lyp/genecp geom=connectivity scf=xqc

5 2
....

Fe 0
lanl2dz
****
C H O N 0
6-31g(d,p)
****

Fe 0
lanl2dz
###############################################

由于多重度是5，计算SCF收敛有点困难，故加上XQC

麻烦各位了，如果遇到同样问题，麻烦指导一下。

作者
Author: sobereva 时间: 2014-11-7 20:54
你的多重度设的是2，电荷是+5，和你描述不符。

先确认过渡态结构定位得很准确，用int=ultrafine结合tight优化收敛限，并且最好calcall，然后再做IRC试试。

如果问题依旧，把IRC的步长调小点，或者尝试别的IRC生成算法。

作者
Author: chenxin199261 时间: 2014-11-7 20:59

sobereva 发表于 2014-11-7 20:54
你的多重度设的是2，电荷是+5，和你描述不符。

先确认过渡态结构定位得很准确，用int=ultrafine结合tigh ...

好的，我在论坛上打错了，gaussian计算时多重度跟电荷没反。

我TS优化calcall已经用上了。
接下来。我打算按您说的int=ultrafine 并结合tight重新计算。
谢谢您

作者
Author: chenxin199261 时间: 2014-11-11 12:22

sobereva 发表于 2014-11-7 20:54
你的多重度设的是2，电荷是+5，和你描述不符。

先确认过渡态结构定位得很准确，用int=ultrafine结合tigh ...

按您说的，我的问题解决了，谢谢你

作者
Author: wyy015020 时间: 2019-1-14 19:43

chenxin199261 发表于 2014-11-11 12:22
按您说的，我的问题解决了，谢谢你

您好，我现在遇到了和您一样的问题，就是用TS找过渡态，正常收敛，有一个虚频，虚频模式正确，但是跑irc的时候只能往产物方向跑，往反应物跑几步就会折返或者终止。所以想请问您最后是怎么解决的呢？重新找到的过渡态和之前的结构差别大吗？

作者
Author: sobereva 时间: 2019-1-15 02:36

wyy015020 发表于 2019-1-14 19:43
您好，我现在遇到了和您一样的问题，就是用TS找过渡态，正常收敛，有一个虚频，虚频模式正确，但是跑irc ...

仔细看
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

作者
Author: wyy015020 时间: 2019-1-17 20:13

sobereva 发表于 2019-1-15 02:36
仔细看
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/threa ...

谢谢sob老师

作者
Author: wyy015020 时间: 2019-1-17 20:37

sobereva 发表于 2019-1-15 02:36
仔细看
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/threa ...

sob老师，我现在加了tight和int=ultrafine优化过渡态，然后进行irc分析两边曲线走出来还是一样的，调大步长也不行。现在用calcall对过渡态重新优化，发现过渡态的结构和之前的结构没有区别，键长键角以及能量都差不多，请问我还有别的尝试方法吗？还是这个过渡态结构压根就是不对的？

作者
Author: sobereva 时间: 2019-1-18 00:18

wyy015020 发表于 2019-1-17 20:37
sob老师，我现在加了tight和int=ultrafine优化过渡态，然后进行irc分析两边曲线走出来还是一样的，调大步 ...

看虚频模式，足够80%确定过渡态是不是你想要的
即便过渡态不是你想要的，只要过渡态是准确的势能面的一阶鞍点，原理上IRC总能正常跑下去
走IRC相关问题该说的我在博文里都说了，没有可以补充的

作者
Author: wyy015020 时间: 2019-1-20 11:27

sobereva 发表于 2019-1-18 00:18
看虚频模式，足够80%确定过渡态是不是你想要的
即便过渡态不是你想要的，只要过渡态是准确的势能面的一 ...

好的，谢谢您，我再试试

作者
Author: wyy015020 时间: 2019-1-20 11:50
本帖最后由 wyy015020 于 2019-1-21 20:35 编辑

sob老师，不好意思再打扰您一下，我刚刚看了优化过渡态的能量是E=-1221.36780221,然后将这个过渡态结构跑IRC，对应坐标为0的能量却是-1221.37933361，请问这是什么原因呢？我是将优化出的TS结构直接写在IRC输入文件里的，其他的过渡态我都是这样操作的，都能正常的跑IRC，就这个出现了这种问题

作者
Author: sobereva 时间: 2019-1-20 15:19

wyy015020 发表于 2019-1-20 11:50
sob老师，不好意思再打扰您一下，我刚刚看了优化过渡态的能量是E=-1227.36780221,然后将这个过渡态结构跑IR ...

8成是做优化和跑IRC用的影响能量的关键词不同

作者
Author: wyy015020 时间: 2019-1-21 09:24
老师，我把关键词发给您看一下，优化过渡态：# M062X/6-311G(d,p) opt(calcfc,ts,noeigen,maxstep=3,notrust) freq scrf=smd int=ultrafine
跑IRC: # M062X/6-311G(d,p) irc(calcfc,maxpoint=20,stepsize=15,LQA) int=ultrafine scrf=smd
还请您指导一下

作者
Author: sobereva 时间: 2019-1-21 16:33

wyy015020 发表于 2019-1-21 09:24
老师，我把关键词发给您看一下，优化过渡态：# M062X/6-311G(d,p) opt(calcfc,ts,noeigen,maxstep=3,notrus ...

光从这上看不出什么来。走IRC时候加上guess=read读取优化任务的chk再试

作者
Author: wyy015020 时间: 2019-1-22 17:18
感谢sob老师，按您说的方法，我的问题解决了，谢谢您

作者
Author: 哈哈哈 时间: 2019-4-28 15:17

wyy015020 发表于 2019-1-17 20:37
sob老师，我现在加了tight和int=ultrafine优化过渡态，然后进行irc分析两边曲线走出来还是一样的，调大步 ...

请问就你说的这种情况你最后怎么解决的安?

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