计算化学公社

标题: 关于qst2算过渡态时遇到的问题求助 [打印本页]

作者
Author: Kellyyyy 时间: 2019-9-21 15:16
标题: 关于qst2算过渡态时遇到的问题求助
在进行qst2计算的时候（已经保证了输入文件的两帧分子的标号完全一致且经过了一定的优化），运算了很久最后error termination。
看了一下log文件，最后一次结果显示如下

Maximum Force          0.000128    0.000450    YES
RMS    Force          0.000027    0.000300    YES
Maximum Displacement    0.019984    0.001800    NO
RMS    Displacement    0.004499    0.001200    NO
Predicted change in Energy=-1.528013D-06
Optimization stopped.
-- Number of steps exceeded,  NStep= 132
-- Flag reset to prevent archiving.

看起来四个标准只满足了两个，应该还需要继续跑。但是不知道是不是因为系统默认最多只能走132次，达到了最大步数之后系统自动终止还是什么原因就停止了。看了看生成的fchk文件，有一个与反应物和生成物比较相关的分子出现，而且有且只有一个虚频，不知道这个是不是过渡态呢？

就想请教一下大家如果是跑到了系统默认的最大步数，那么可以如何继续下去呢？是添加某个指令增大步数还是在这次计算的结果上继续ts或者restart呢？

作者
Author: sobereva 时间: 2019-9-21 17:17
原因已经提示得很清楚了Number of steps exceeded, NStep= 132

如果最后一步的结构看起来像你找的过渡态的样子，做振动分析可以看出对应反应路径的方向有显著的虚频，就用这个结构用opt=TS找过渡态
如果届时出现震荡之类导致难收敛，用此文的方法
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

作者
Author: Kellyyyy 时间: 2019-9-21 17:40

sobereva 发表于 2019-9-21 17:17
原因已经提示得很清楚了Number of steps exceeded, NStep= 132

如果最后一步的结构看起来像你找的过渡 ...

确实是有一个虚频，不过要怎么才能看得出来这个虚频与反应路径方向是否相关呢？

作者
Author: Kellyyyy 时间: 2019-9-21 17:46

Kellyyyy 发表于 2019-9-21 17:40
确实是有一个虚频，不过要怎么才能看得出来这个虚频与反应路径方向是否相关呢？

哦哦哦，点开frequency里面看那个虚频对应的动画确实能看出来。好像是对的

作者
Author: Kellyyyy 时间: 2019-9-21 19:21

sobereva 发表于 2019-9-21 17:17
原因已经提示得很清楚了Number of steps exceeded, NStep= 132

如果最后一步的结构看起来像你找的过渡 ...

以其作为过渡态进行设置步数为50的irc计算，算了大半天之后，error termination。。。
出现的错误如下：

Cartesian Forces:  Max    0.002770011 RMS    0.000642311
IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC
Error in corrector energy =       -0.0000058440
Magnitude of corrector gradient =    0.0049232741
Magnitude of analytic gradient =    0.0052181590
Magnitude of difference =          0.0022517538
Angle between gradients (degrees)= 25.4434
Pt 16 Step number  20 out of a maximum of  20
CORRECTOR INTEGRATION CONVERGENCE:
Recorrection delta-x convergence threshold: 0.010000
Delta-x Convergence NOT Met
Maximum number of corrector steps exceded.
Error termination

一个是校正能出错，一个是校正步数超过了最大值。。。
您看这样的话还有必要restart吗？或者说需要添加指令增大校正步数吗？

作者
Author: Kellyyyy 时间: 2019-9-21 19:53

Kellyyyy 发表于 2019-9-21 19:21
以其作为过渡态进行设置步数为50的irc计算，算了大半天之后，error termination。。。
出现的错误如下： ...

我在您之前的教程中已经看到相关的调试方法，谢谢！

在另一个计算文件中，同样出现了超过了最大步数。。。
Maximum Force          0.010487    0.000450    NO
RMS    Force          0.001550    0.000300    NO
Maximum Displacement    0.090045    0.001800    NO
RMS    Displacement    0.020281    0.001200    NO
Predicted change in Energy=-9.002281D-04
Optimization stopped.
-- Number of steps exceeded,  NStep= 111
-- Flag reset to prevent archiving.

这个是qst3计算，不知道有没有什么合适的指令可以把设置的步数设置得更大啊？或者说这个是系统默认的还是受电脑的计算能力及内存的限制呢？

作者
Author: sobereva 时间: 2019-9-22 07:15

Kellyyyy 发表于 2019-9-21 19:53
我在您之前的教程中已经看到相关的调试方法，谢谢！

在另一个计算文件中，同样出现了超过了最大步数。 ...

maxcyc设优化步数上限
别盲目设这个，先看收敛趋势再说

作者
Author: Kellyyyy 时间: 2019-9-22 12:21
好的，谢谢您！

我在跑另外一个qst2的时候，遇到了一个奇怪的问题，刚一开始跑就报错了：
New curvilinear step not converged.
***** Convergence failure in GTrans *****
Try using 3 structures as input for
QST transition state search
Error termination

可是我把它与我之前跑的qst2的文件对比了一下，gjf的前几行几乎都是差不多的，那为什么这个曲线不收敛？还建议跑qst3呢？

附gjf：
# opt=(qst2,calcfc) freq ub3lyp/6-311++g(d,p) scrf=(cpcm,solvent=water)
geom=connectivity

请解答一下哈

作者
Author: Kellyyyy 时间: 2019-9-22 12:23
而我不在集群上面在我自己电脑安装的g09上面跑就不会出现这个错误不会终止，这是为什么呢？

作者
Author: sobereva 时间: 2019-9-23 07:54

Kellyyyy 发表于 2019-9-22 12:21
好的，谢谢您！

我在跑另外一个qst2的时候，遇到了一个奇怪的问题，刚一开始跑就报错了：

QST2本来就很不好用，我的观点在这里写明了
简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460（http://bbs.keinsci.com/thread-12093-1-1.html）

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