24小时热门版块排行榜    

查看: 497  |  回复: 3
当前主题已经存档。

dralexandero

木虫 (小有名气)

[交流] 【求助】过渡态优化的问题

对CH3OCl分解途径的一个过渡态优化,计算结果显示l9999出错。
Item                              Value          Threshold  Converged
Maximum Force               0.042266     0.000450     NO
RMS     Force                  0.019617     0.000300     NO
Maximum Displacement     0.021614     0.001800     NO
RMS     Displacement        0.008988     0.001200     NO
Predicted change in Energy=-1.631015D-03
Optimization stopped.
    -- Number of steps exceeded,  NStep= 100
    -- Flag reset to prevent archiving.
四个收敛条件没一个符合的,我的输入结构是照文献上来的,应该比较接近预想的优化后构型,事实上在计算过程中用Gview看了下算着算着就开始向不是预想的方向偏离了。请教下这可能是什么原因?
输入文件如下:
%mem=200mb
#p ump2(full)/6-31g(d,p) opt(ts,calcfc,noeigen,maxcyc=400) scf(maxcyc=400) iop(1/8=2,2/16=1,5/13=1) freq

Title Card Required

0 1
C
H                    1    B1
H                    1    B2       2    A1
H                    1    B3       2    A2       3    D1       0
O                    1    B4       4    A3       2    D2       0
Cl                   5    B5       1    A4       4    D3       0
      
  B1                    1.389                    
  B2                    1.07                     
  B3                    1.07                     
  B4                    1.322                    
  B5                    2.233                    
  A1                  120.80489                  
  A2                  120.80487                  
  A3                  120.80488                  
  A4                   84.5                     
  D1                 -143.83572                  
  D2                   65.04006                  
  D3                   71.91775

[ Last edited by dralexandero on 2009-3-24 at 15:40 ]
回复此楼
来 是 空 言 去 绝 踪
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)

★ ★
xuefei06(金币+2,VIP+0):3Q! 3-24 23:58
你把中间优化方向没有错误之前的结构取出来,按同样的关键字再接着算吧。
如果还不行,观察一下离收敛还有多远时开始偏离你想要的方向,如果离收敛不远,同样把结构再取出来,加calcall参数优化,应该就行了。
另外我对noeigen参数没多少把握,你确定你这个优化加上它合适吗?
2楼2009-03-24 22:18:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dralexandero

木虫 (小有名气)

谢谢,我明天试试看吧,noeigen应该是要加的,我算它的另一个途径的过渡态时不加就会出错,应该不是这个问题。
来 是 空 言 去 绝 踪
3楼2009-03-24 23:20:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sshadow

铁虫 (初入文坛)

★ ★
yuhuobuku(金币+2,VIP+0):欢迎参加讨论 4-8 09:12
noeigen这个按照说明书上的说法是说不要程序检查负本征值的数量,如果数目不正确就终止,如果要用就最好把maxcycle设置成最小5,我对这个关键字一直也很疑惑,有高手指点下就好了。
4楼2009-04-04 22:28:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dralexandero 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见