计算化学公社

标题: 一堆的L9999怎么办啊怎么办 [打印本页]

作者
Author: alystone 时间: 2015-6-22 11:47
标题: 一堆的L9999怎么办啊怎么办
请教：我也知道L9999是优化步数不够，可是我都read上一个很多次了，还是L9999，并且我的输入结构本来就是在gas phase里优化过的坐标，然后加上关键词scrf就总是L9999。
      这是输入：
      #p opt freq rb3lyp/gen pseudo=read int=noxctest pop=(full,nbo,savenbos) scrf=pcm
      这是输出：
      Reordering of NBOs for storage: 177  363  293  187  158  409  454  446  401  233
Reordering of NBOs for storage: 319  364  341  357  312  332  399  444  377  467
Reordering of NBOs for storage: 422  386  354  309  170  232  301  321  366  433
Reordering of NBOs for storage: 247  478  298  244  385  214  340  456  411  378
Reordering of NBOs for storage: 215  333  423  387  468  222  457  412  231  220
Reordering of NBOs for storage: 269  367  322  342  403  448  213  358  313  447
Reordering of NBOs for storage: 402  173  287  223  271  289  145  311  356  224
Reordering of NBOs for storage: 376  331  421  466  455  410  263  300  246  274
Reordering of NBOs for storage: 430  475  292  212  320  365  160  257  284  229
Reordering of NBOs for storage: 193  203  264  302  248  291  273  159  198  189
Reordering of NBOs for storage: 166  290  272  165  167  169  183  178  157  168
Reordering of NBOs for storage: 179  172  176  188  182  186  171
Labels of output orbitals:  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR
Labels of output orbitals:  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR  CR
Labels of output orbitals:  CR  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD
Labels of output orbitals:  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD
Labels of output orbitals:  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD
Labels of output orbitals:  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  BD  LP  LP  LP  BD  BD  BD
Labels of output orbitals:  BD  BD  BD  BD  BD  LP  LP  BD  BD  BD* BD* LP* BD* BD* BD* BD* BD* BD* BD* BD*
Labels of output orbitals:  BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD*
Labels of output orbitals:  BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD*
Labels of output orbitals:  BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD*
Labels of output orbitals:  BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* BD* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY* RY*
Labels of output orbitals:  RY* RY* RY* RY* RY* RY* RY*
Leave Link  607 at Mon Apr 13 12:53:12 2015, MaxMem= 2621440000 cpu:    162.0
(Enter /public/soft/software/g09-c01/g09/l9999.exe)

LOVE IS BLIND, THAT'S WHY ALL THE WORLD LOVES A LOUVER.
Error termination request processed by link 9999.
Error termination via Lnk1e in /public/soft/software/g09-c01/g09/l9999.exe at Mon Apr 13 12:53:12 2015.
Job cpu time: 15 days 23 hours 43 minutes  8.8 seconds.
File lengths (MBytes):  RWF= 1041 Int=    0 D2E=    0 Chk=    79 Scr=    1

作者
Author: alystone 时间: 2015-6-22 11:52
本帖最后由 alystone 于 2015-6-22 11:59 编辑

这是优化步数C:\Users\liuzy\Desktop\QQ图片20150622114836.png

作者
Author: sobereva 时间: 2015-6-22 13:24
几何优化出现L9999绝对不是优化的步数上限不够（这是初学者犯的最严重误区之一，并因此总是盲目加大步数上限到好几百，导致不仅完全没解决问题，还白浪费大量计算时间），而是默认的优化方法和相关设定对当前体系优化不合适，需要调整。

合理的解决方法此文有详述。最常用的几种处理见文末。
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

另外，int=noxctest这个关键词绝对不要写，没有任何好处，只会增加结果错误的风险。

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