计算化学公社

标题: 固定二面角优化，优化一定圈数后，之后的收敛趋势基本不变了，求大神帮忙看看！ [打印本页]

作者
Author: 风影月瞳 时间: 2016-1-14 20:15
标题: 固定二面角优化，优化一定圈数后，之后的收敛趋势基本不变了，求大神帮忙看看！
输入命令： #p opt=ModRedundant b3lyp/6-31G(d) iop(1/8=10) td(nstates=5,root=1) sc f=(tight,maxcycle=500)

第26262行
  D54    -0.02709  -0.00067 0.00000 0.00000 0.00000  -0.02709
               Item             Value    Threshold  Converged?
Maximum Force                 0.062960    0.000450    NO
RMS    Force                   0.014983    0.000300    NO
Maximum Displacement    0.000000    0.001800    YES
RMS    Displacement          0.000000    0.001200    YES
Predicted change in Energy=-5.602965D-03
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

之后的收敛趋势基本都是这样，如

26983行
      Item             Value    Threshold  Converged?
Maximum Force                0.062960    0.000450    NO
RMS    Force                   0.014983    0.000300    NO
Maximum Displacement    0.000000    0.001800    YES
RMS    Displacement       0.000000    0.001200    YES
Predicted change in Energy=-5.602860D-03

27639行
      Item             Value    Threshold  Converged?
Maximum Force                0.062960    0.000450    NO
RMS    Force                   0.014983    0.000300    NO
Maximum Displacement    0.000000    0.001800    YES
RMS    Displacement       0.000000    0.001200    YES
Predicted change in Energy=-5.602863D-03

一直到9万多行都是差不多这个趋势，是不是死循环了

最后得出 ! Non-Optimized Parameters !
                        ! (Angstroms and Degrees)  !
--------------------------                         --------------------------
! Name  Definition             Value       Derivative Info.             !
--------------------------------------------------------------------------------
! R1 R(1,2)                1.3928       -DE/DX = 0.0131             !
! R2 R(1,5)                1.4252       -DE/DX = -0.0045             !
! R3 R(1,9)                1.0857       -DE/DX = 0.0014             !
! R4 R(2,3)                1.3836       -DE/DX = 0.0076             !
! R5 R(2,10)                1.0846       -DE/DX = 0.0025             !
! R6 R(3,12)                1.4224       -DE/DX = -0.007             !
! R7 R(3,16)                1.0835       -DE/DX = 0.0016             !
! R8 R(4,5)                1.4524       -DE/DX = -0.0307             !
! R9 R(4,6)                1.5272       -DE/DX = -0.0322             !
! R10 R(4,12)                1.3696       -DE/DX = 0.0136             !
! R11 R(5,21)                1.3443       -DE/DX = 0.012             !
! R12 R(6,8)                1.3879       -DE/DX = -0.0491             !
! R13 R(6,14)                1.3511       -DE/DX = 0.0239             !
! R14 R(7,11)                1.0803       -DE/DX = 0.0032             !
! R15 R(7,14)                1.4466       -DE/DX = -0.0069             !
谁有经验，怎样才能解决呢？祝大家工作顺利，生活愉快！求帮助！

作者
Author: sobereva 时间: 2016-1-14 23:06

位移都已经是0了，如果观看优化轨迹，应该结构都已经不动了，所以继续跑下去也是徒劳。
尝试GDIIS之类优化算法看看能不能继续走，或者基组、泛函换换试试。

BTW：

tight是多余的，本身就是默认的
不要写maxcyc=500，毫无意义

1/8=10这种设定步长上限的设定，用的时候必须和notrust一起用，否则起不到任何效果

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