计算化学公社

标题: 走IRC的时候forward方向只走几十个点求助 [打印本页]

作者
Author: 平娜英 时间: 2020-11-30 21:41
标题: 走IRC的时候forward方向只走几十个点求助
老师，您好，我想问一下，就是我在走IRC的时候，一个方向可以走几百个点，但是另一个方向只走几十个点，这是为什么吖？我尝试过换步长，尝试过读取chk文件，也尝试过分开走reverse和forward都是这样，并且有3个过渡态走IRC时都出现同样的问题。现在还有一个IRC正在计算。我是分着走的，reverse方向可以走几百个点还没有结束，forward还是走几十个点就结束了，实在不知道是为什么。3个图分别如下，第4个forward方向只走了28个点。谢谢您啦！

作者
Author: hebrewsnabla 时间: 2020-11-30 23:08
收敛了不就结束了，你还想怎样。还是你觉得没收敛

作者
Author: liyuanhe211 时间: 2020-12-1 00:24
IRC的终点是极小点（附近），超过它IRC没法知道下面再往哪个方向走。所以这个选项叫“maxpoints”不是“number of points”

作者
Author: 平娜英 时间: 2020-12-1 08:39

hebrewsnabla 发表于 2020-11-30 23:08
收敛了不就结束了，你还想怎样。还是你觉得没收敛

收敛了呀，但是计算能垒的时候和过渡态的能量差不多

作者
Author: 平娜英 时间: 2020-12-1 08:40

liyuanhe211 发表于 2020-12-1 00:24
IRC的终点是极小点（附近），超过它IRC没法知道下面再往哪个方向走。所以这个选项叫“maxpoints”不是“num ...

那怎么设置会让它继续往下走呢

作者
Author: biogon 时间: 2020-12-1 09:23
那个位置就在一个极小点附近了，再跑也跑不动了

作者
Author: FETEYA 时间: 2020-12-1 09:47
取出forward方向终点的结构做opt,看看优化得到的稳定结构和终点的结构是不是同一个

作者
Author: 平娜英 时间: 2020-12-1 10:57

FETEYA 发表于 2020-12-1 09:47
取出forward方向终点的结构做opt,看看优化得到的稳定结构和终点的结构是不是同一个

两者的构型差不多，在键长、键角上会有0.1或0.2的区别

作者
Author: 平娜英 时间: 2020-12-1 11:00

biogon 发表于 2020-12-1 09:23
那个位置就在一个极小点附近了，再跑也跑不动了

那计算出来的能垒和TS的能垒相差0.1或0.2，这数据感觉有问题啊

作者
Author: zjxitcc 时间: 2020-12-1 11:04

平娜英发表于 2020-12-1 11:00
那计算出来的能垒和TS的能垒相差0.1或0.2，这数据感觉有问题啊

IRC正常停掉的末端点，只是很接近local minimum，不一定就是。还需要提取结构再优化。

“计算出来的能垒和TS的能垒相差0.1或0.2”描述大概率有误，你想说TS的能量？

作者
Author: 平娜英 时间: 2020-12-1 11:14
我已拿出IRC末点优化了。我说的是下图中TS和PC1的能量差值

作者
Author: FETEYA 时间: 2020-12-1 20:49

平娜英发表于 2020-12-1 11:14
我已拿出IRC末点优化了。我说的是下图中TS和PC1的能量差值

这种情况也是可能的,算出来能垒是多少就是多少

作者
Author: naoki 时间: 2020-12-1 21:45
我之前TS找错，走IRC就像你图2一样，后来发现我TS虚频只有30多，于是重新摆初猜重新来过就正常了

作者
Author: 平娜英 时间: 2020-12-2 09:52

naoki 发表于 2020-12-1 21:45
我之前TS找错，走IRC就像你图2一样，后来发现我TS虚频只有30多，于是重新摆初猜重新来过就正常了

我的TS虚频是80多，并且震动趋势也是对的

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