24小时热门版块排行榜    

Znn3bq.jpeg
查看: 581  |  回复: 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的回帖

sshadow

铁虫 (初入文坛)

★ ★
yuhuobuku(金币+2,VIP+0):欢迎参加讨论 4-8 09:12
noeigen这个按照说明书上的说法是说不要程序检查负本征值的数量,如果数目不正确就终止,如果要用就最好把maxcycle设置成最小5,我对这个关键字一直也很疑惑,有高手指点下就好了。
4楼2009-04-04 22:28:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +6 1624447980 2026-04-08 7/350 2026-04-10 08:06 by 探123
[考研] 314求调剂 +16 xhhdjdjsjks 2026-04-09 17/850 2026-04-09 18:26 by 喝醉酒N
[材料工程] 材料调剂推荐 +3 蛋糕x2 2026-04-07 3/150 2026-04-09 17:01 by Lilly_Li
[考研] 085404 293求调剂 +7 勇远库爱314 2026-04-08 7/350 2026-04-09 16:02 by 猪会飞
[考研] 机械调剂 +3 zzzbcb 2026-04-07 3/150 2026-04-07 22:19 by hemengdong
[考博] 博士申请 +3 IQwQl 2026-04-05 3/150 2026-04-07 20:31 by greychen00
[考研] 专硕085403,291分,有两篇专利,一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 材料求调剂 +18 一样YWY 2026-04-05 18/900 2026-04-07 15:49 by dxlg
[考研] 316求调剂 +7 yyx想调剂 2026-04-05 7/350 2026-04-07 14:31 by shdgaomin
[考研] 293求调剂 +3 勇远库爱314 2026-04-06 3/150 2026-04-07 11:15 by hugr
[考研] 考研调剂 +7 15615482637 2026-04-04 7/350 2026-04-06 22:56 by chenzhimin
[考研] 308求调剂 +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by 蓝云思雨
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666!
[考研] 296求调剂 +3 汪!?! 2026-04-05 5/250 2026-04-05 17:38 by 蓝云思雨
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 22408求调剂 354分 可跨专业 +3 hannnnnnn 2026-04-04 3/150 2026-04-04 14:35 by 土木硕士招生
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-03 16:46 by wxiongid
信息提示
请填处理意见