计算化学公社

标题: 请教高斯优化过渡态未收敛，振动方式正确的结构算正确吗？ [打印本页]

作者
Author: 阿秋Autumn 时间: 2016-12-6 16:33
标题: 请教高斯优化过渡态未收敛，振动方式正确的结构算正确吗？
请教一下，用高斯优化过渡态的时候，优化后出现l9999报错，但是振动方式和猜测的相同，而且在log文件中也给出了能量，请问这个结构可以使用吗？
也用过其他方法尝试继续收敛，可是收敛之后优化出的结构和初猜结构相差很大，而且也不是想要的振动方式~
希望前辈们能够帮帮忙，十分感谢！！！！

作者
Author: sobereva 时间: 2016-12-7 00:21
如果你是指优化过渡态没收敛，然后用最后一帧的结构直接做振动分析，发现虚频模式是想要的，那么能否用要看freq任务最后对收敛的判断情况。如果受力和位移比默认收敛限大不超过2倍，勉强可以接受，否则必须重新优化过渡态，calcfc/calcall结合maxstep减小步长上限一般都能解决。

作者
Author: 阿秋Autumn 时间: 2016-12-7 08:49

sobereva 发表于 2016-12-7 00:21
如果你是指优化过渡态没收敛，然后用最后一帧的结构直接做振动分析，发现虚频模式是想要的，那么能否用要看 ...

sob老师，我的是在l9999的时候虚频模式是想要的，而继续用最后优化的结构继续优化则出现了错误的振动模式，使用的是calcall，有两个结构都是这样，这是因为我的初猜错了吗？

作者
Author: sobereva 时间: 2016-12-7 12:51

阿秋Autumn 发表于 2016-12-7 08:49
sob老师，我的是在l9999的时候虚频模式是想要的，而继续用最后优化的结构继续优化则出现了错误的振动模式 ...

用opt(TS,maxstep=5,calcall)基于之前最后的结构继续优化过渡态，如果确实跑到其它地方去了，那即便之前那个结构虚频正确也白搭

有唯一的虚频且振动方向正确也不代表结构就是过渡态结构，必须结合受力为0的条件（做振动分析时末尾输出的受力接近0）

作者
Author: 阿秋Autumn 时间: 2016-12-8 08:48

sobereva 发表于 2016-12-7 12:51
用opt(TS,maxstep=5,calcall)基于之前最后的结构继续优化过渡态，如果确实跑到其它地方去了，那即便之前 ...

sob老师，我用您的方法对最后的结构继续进行了优化，可是最后还是出现了l9999的错误，信息如下：
Error termination request processed by link 9999.
Error termination via Lnk1e in /lustre/software/g09/l9999.exe at Thu Dec  8 00:32:58 2016.
Job cpu time:    0 days  7 hours 15 minutes 46.9 seconds.
File lengths (MBytes):  RWF= 285 Int=    0 D2E=    0 Chk=    8 Scr=    1

希望老师指点，感谢~

作者
Author: sobereva 时间: 2016-12-8 17:23

阿秋Autumn 发表于 2016-12-8 08:48
sob老师，我用您的方法对最后的结构继续进行了优化，可是最后还是出现了l9999的错误，信息如下：
Error ...

别光看l9999，要用gview观看优化过程结构和能量、受力变化趋势，判断怎么处理
凡是可能有用的做法都在这里提了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

作者
Author: 阿秋Autumn 时间: 2016-12-8 18:09

sobereva 发表于 2016-12-8 17:23
别光看l9999，要用gview观看优化过程结构和能量、受力变化趋势，判断怎么处理
凡是可能有用的做法都在这 ...

好的，谢谢sob老师~

作者
Author: laoman 时间: 2017-10-10 00:05

sobereva 发表于 2016-12-7 12:51
用opt(TS,maxstep=5,calcall)基于之前最后的结构继续优化过渡态，如果确实跑到其它地方去了，那即便之前 ...

Sobereva老师，请问什么情况下会出现TS优化出只有一个虚频且是反应方向上的虚频，而freq时的Forces远大于0的情况呢？

作者
Author: sobereva 时间: 2017-10-10 00:15

laoman 发表于 2017-10-10 00:05
Sobereva老师，请问什么情况下会出现TS优化出只有一个虚频且是反应方向上的虚频，而freq时的Forces远大于 ...

不会出现这种情况。force远大于0的情况opt任务根本不会正常结束。

作者
Author: laoman 时间: 2017-10-10 06:05

sobereva 发表于 2017-10-10 00:15
不会出现这种情况。force远大于0的情况opt任务根本不会正常结束。

明白了，感谢社长的回答。最近发现原来从freq的chk（来自PES扫描后的最高点，有反应方向的虚频，-1600到-1800的范围）读取geom和guess会比导出坐标算的快了一点（虽然也用了guess），读取chk的前提应该是freq和后面的ts都用的同样的基组吧。

作者
Author: sobereva 时间: 2017-10-10 14:01

laoman 发表于 2017-10-10 06:05
明白了，感谢社长的回答。最近发现原来从freq的chk（来自PES扫描后的最高点，有反应方向的虚频，-1600到- ...

快一点是因为从chk里读取了初猜，导致SCF部分基本不耗什么时间。
guess=read前后的基组可以不同，程序自动会做投影

作者
Author: laoman 时间: 2017-10-10 15:38

sobereva 发表于 2017-10-10 14:01
快一点是因为从chk里读取了初猜，导致SCF部分基本不耗什么时间。
guess=read前后的基组可以不同，程序 ...

好的，谢谢社长～

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