应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 【求助】过渡态优化的问题
41/1返回列表
查看: 691  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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 ]
回复此楼

» 猜你喜欢
来 是 空 言 去 绝 踪
1楼 2009-03-24 15:38:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jasminer

铁杆木虫 (著名写手)
★ ★
xuefei06(金币+2,VIP+0):3Q! 3-24 23:58
你把中间优化方向没有错误之前的结构取出来，按同样的关键字再接着算吧。
如果还不行，观察一下离收敛还有多远时开始偏离你想要的方向，如果离收敛不远，同样把结构再取出来，加calcall参数优化，应该就行了。
另外我对noeigen参数没多少把握，你确定你这个优化加上它合适吗？
赞一下(10人) 回复此楼
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，我对这个关键字一直也很疑惑，有高手指点下就好了。
赞一下(6人) 回复此楼
4楼2009-04-04 22:28:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dralexandero 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭