24小时热门版块排行榜    

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

dralexandero

木虫 (小有名气)

谢谢,我明天试试看吧,noeigen应该是要加的,我算它的另一个途径的过渡态时不加就会出错,应该不是这个问题。
来 是 空 言 去 绝 踪
3楼2009-03-24 23:20:09
已阅   回复此楼   关注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的回帖

sshadow

铁虫 (初入文坛)

★ ★
yuhuobuku(金币+2,VIP+0):欢迎参加讨论 4-8 09:12
noeigen这个按照说明书上的说法是说不要程序检查负本征值的数量,如果数目不正确就终止,如果要用就最好把maxcycle设置成最小5,我对这个关键字一直也很疑惑,有高手指点下就好了。
4楼2009-04-04 22:28:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程调剂289 +41 yang婷 2026-04-07 47/2350 2026-04-09 23:23 by wolf97
[考研] 301求调剂 +5 静静想想 2026-04-05 5/250 2026-04-09 20:58 by 运气yunqi
[考研] 化工求调剂! +34 RichLi_ 2026-04-06 34/1700 2026-04-09 20:39 by zhouxiaoyu
[考研] 085400 328分 求调剂 +8 喂你一个大橙子 2026-04-09 12/600 2026-04-09 19:20 by zl8213662
[考研] 292求调剂 +8 笑笑袁 2026-04-09 8/400 2026-04-09 18:14 by 1753564080
[考研] 一志愿中科大070300化学,314分求调剂 +11 wakeluofu 2026-04-09 11/550 2026-04-09 18:03 by lijunpoly
[考研] 262求调剂 +10 天下第一文 2026-04-04 13/650 2026-04-09 15:16 by 探123
[硕博家园] 新一代电子信息294求调剂 不挑学校 +5 Ytyt11 2026-04-09 6/300 2026-04-09 14:40 by Ytyt11
[考研] 一志愿0807 数一英一 313 有没有二轮调剂 +11 emokidd 2026-04-08 12/600 2026-04-09 09:24 by wyf236
[考研] 一志愿北京2,材料与化工308求调剂 +17 熊二想上岸 2026-04-04 18/900 2026-04-09 09:13 by leyan1127
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 一志愿郑州大学085600求调剂 +21 吃的不少 2026-04-05 24/1200 2026-04-08 16:47 by sunhuadong
[考研] 297分083200求助 +9 aekx 2026-04-05 9/450 2026-04-06 20:57 by flysky1234
[考研] 求调剂 +10 chenxrlkx 2026-04-05 10/500 2026-04-06 11:31 by 猪会飞
[考研] 308求调剂 +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by 蓝云思雨
[考研] 348求调剂 +6 wukira 2026-04-04 6/300 2026-04-05 18:11 by 猪会飞
[考研] 288求调剂,一志愿华南理工大学071005 +6 ioodiiij 2026-04-04 6/300 2026-04-05 10:09 by guoweigw
[考研] 考研调剂 +11 小sun要好运 2026-04-04 11/550 2026-04-05 08:02 by qlm5820
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-04 20:37 by 蓝云思雨
信息提示
请填处理意见