计算化学公社

标题: 求助l9999错误 [打印本页]

作者
Author: 小菜鸟求关照 时间: 2014-12-24 19:09
标题: 求助l9999错误
我在优化过渡态的时候出现Error termination request processed by link 9999.错误，我知道是优化圈数不够，#p b3lyp genecp opt=(qst2 Noeigentest) freq out=wfx iop(1/8=1)，这个是输入文件题目，这个错误怎么解决呢？谢谢各位老师啦，我在导出log文件把它转化为gjf的时候出现这个错误

作者
Author: diaolanxinyu 时间: 2014-12-24 21:02
L9999不一定就是迭代圈数不够，有时候增加到足够大也还是会出错，有可能是构型或是方法，都要尝试修改一下，过渡态本来就不好找，如果只局限在一个错误，会卡住的，呵呵

还有，要查看高斯输出文件最后有ERROR行，并多查看至少10行，详细查找原因。

对于初学过渡态计算的，我建议用QST3，我之前帮别人简单算过一个，用QST2算，总是出错，各种错误，然后各种修改后还是出错，用QST3就解决了，呵呵

一点拙见，希望对你有用！

作者
Author: diaolanxinyu 时间: 2014-12-24 21:09
“我在导出log文件把它转化为gjf的时候出现这个错误”
这是什么意思呢？为什么要把log文件转化成gjf呢，过渡态计算完成后不就是要查看和分析log文件吗？
转化做什么呢？
如果说需要保存l输出文件中其中某一步结构的话，直接点gaussview上的保存按钮就可以了，也不会出错的。

这个图中的错误应该是你打开log文件时出现的吧，我觉得还是要修改输入文件重新计算。

作者
Author: 小菜鸟求关照 时间: 2014-12-24 21:21

diaolanxinyu 发表于 2014-12-24 21:09
“我在导出log文件把它转化为gjf的时候出现这个错误”
这是什么意思呢？为什么要把log文件转化成gjf呢，过 ...

我就是像你那样保存的，这个时候出现这个错误的

作者
Author: 小菜鸟求关照 时间: 2014-12-24 21:25

diaolanxinyu 发表于 2014-12-24 21:02
L9999不一定就是迭代圈数不够，有时候增加到足够大也还是会出错，有可能是构型或是方法，都要尝试修改一下 ...

学到了我刚开始学过渡态，谢谢您大神

作者
Author: zhuhongkunky 时间: 2014-12-24 23:09
本帖最后由 zhuhongkunky 于 2014-12-24 23:22 编辑

我最近也初学过渡态，这个很有可能是优化到一半了，你可以看看构型是不是介于反应物与生成物之间（主要看看反应位置的键长键角的变化），如果是的话，把这个构型变为反应物接着算，或者再从开始增加圈数后在算。
还有可能构型之间的键长太长了，可适当缩短（一般比正常的多30%左右比较好）。
再不济就是某些原子进攻的位置不对，可以多试几个不同的位置。
建议加上calcfc，找到TS概率比较高。
本人拙见，要有不对还多多交流哈:P

作者
Author: 小菜鸟求关照 时间: 2014-12-25 09:23

zhuhongkunky 发表于 2014-12-24 23:09
我最近也初学过渡态，这个很有可能是优化到一半了，你可以看看构型是不是介于反应物与生成物之间（主要看看 ...

好的我试一试哈

作者
Author: 小菜鸟求关照 时间: 2014-12-25 09:24

diaolanxinyu 发表于 2014-12-24 21:02
L9999不一定就是迭代圈数不够，有时候增加到足够大也还是会出错，有可能是构型或是方法，都要尝试修改一下 ...

用qst3的时候需要加入过渡态初始构型这个好难构造的，很难猜测出它的初始构型吧

作者
Author: diaolanxinyu 时间: 2014-12-25 09:43
可以把你用QST2计算的结果作为QST3的过渡态结构或者根据你的反应物和产物猜测一个过渡态结构

我当时用QST2时老是出错，虽然出错但还是会得到一个过渡态结构，然后将这个结构作为QST3的过渡态结构，然后就成功了，呵呵

要不你试试，或者网上搜一下会有很多讨论关于过渡态计算的

作者
Author: sobereva 时间: 2014-12-25 09:47
l9999这种提示没有什么意义，N多种错误都会导致同样的报错，要把输出文件从后往前看找到和报错直接对应的提示语句来判断原因。
通常没必要用QST2/3，直接用opt=TS最省事。
如果从输出文件看就是没收敛导致的，而且用gview观看TS优化过程发现体系发生了震荡，那么尝试调整初始结构、加上int=ultrafine、减小步长上限、切换优化算法、加上calcall、换泛函或基组等方式去解决。没有万能的解决方法，可行的做法在此帖有讨论（虽然主要谈几何优化，但是处理方法对于搜索过渡态大多也同样适用）：《量子化学计算中帮助几何优化收敛的常用方法》（http://sobereva.com/164）。

顺带一提，iop(1/8=x)这个调整步长的关键词不宜总是加上，有必要的时候才需要加，通常都是优化到后期发生了震荡的时候基于后期的结构继续优化时才用。尤其是你让x=1，会令优化过程变得很缓慢，白费很长时间。如果确实要限制步长上限，通常也应该结合notrust关键词，否则可能起不到效果。

作者
Author: diaolanxinyu 时间: 2014-12-25 09:48
http://hi.baidu.com/sobereva/item/c5a8de0b88de1c10acdc70d7

http://hi.baidu.com/sobereva/item/f6adaa62841af809a1cf0f9a

http://emuch.net/bbs/viewthread.php?tid=1527579
这是sobereva老师和小卒老师的帖子，原理和经验，希望对你有用

作者
Author: diaolanxinyu 时间: 2014-12-25 09:50
http://emuch.net/bbs/viewthread.php?tid=1620532
还有一个，呵呵，多多使用你手中的搜索引擎，会有惊喜等你奥，呵呵

作者
Author: diaolanxinyu 时间: 2014-12-25 09:55
http://hi.baidu.com/sobereva/item/c23f2edeba7400312b35c7dd
这是sobereva老师的blog，现在已经成了我每天必登的网站。

成长可以从这里开始！:)

作者
Author: 小菜鸟求关照 时间: 2014-12-25 10:08

diaolanxinyu 发表于 2014-12-25 09:55
http://hi.baidu.com/sobereva/item/c23f2edeba7400312b35c7dd
这是sobereva老师的blog，现在已经成了我每 ...

谢谢

作者
Author: 小菜鸟求关照 时间: 2014-12-25 10:14

diaolanxinyu 发表于 2014-12-25 09:55
http://hi.baidu.com/sobereva/item/c23f2edeba7400312b35c7dd
这是sobereva老师的blog，现在已经成了我每 ...

谢谢您，您推荐的那几个帖子对我很有用，您也是一位大神！

作者
Author: diaolanxinyu 时间: 2014-12-25 10:23
不管谁回复，只要对你有用，有收获就行，sobereva老师才是我们的大神，不要盲目崇拜！

作者
Author: 小菜鸟求关照 时间: 2014-12-25 11:43
我用的qst3,然后出现l9999错误以后接着提取分子坐标，这样连续几次提取分子坐标作为过渡态，如此进行，最后收敛正常，存在一个虚频。成功！

作者
Author: 小菜鸟求关照 时间: 2014-12-25 11:45

diaolanxinyu 发表于 2014-12-25 10:23
不管谁回复，只要对你有用，有收获就行，sobereva老师才是我们的大神，不要盲目崇拜！

sobereva老师是我们的大神，这是谁也无法替代的，从老师那里我也学到了好多，您帮我解决了问题，您也是我的老师，真的非常感谢sobereva老师和您。

作者
Author: zhuhongkunky 时间: 2014-12-25 14:54

小菜鸟求关照发表于 2014-12-25 09:24
用qst3的时候需要加入过渡态初始构型这个好难构造的，很难猜测出它的初始构型吧

那用qst2也行啊

作者
Author: 小菜鸟求关照 时间: 2014-12-25 15:30

zhuhongkunky 发表于 2014-12-25 14:54
那用qst2也行啊

嗯呢已经成功解决了

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