24小时热门版块排行榜    

查看: 1229  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yangmou

木虫 (小有名气)

[求助] 一个四原子自由基COONa怎么也优化不出来了!!!!! 已有2人参与

输入文件
%mem=4GB
%nprocshared=8
%chk=COONa-2.chk
#p opt=calcall freq=noraman ub3lyp/6-311+g(d,p) nosymm
scf=(qc,maxcycle=365,tight)

COONa

0 2
C              
O                  1            B1
O                  1            B2    2            A1
Na                 3            B3    1            A2    2            D1    0

   B1             1.24377491
   B2             1.24386992
   B3             2.42551046
   A1           133.73088482
   A2            81.25155496
   D1            21.67057927


最后输出:
Maximum Force            0.111497     0.000450     NO
RMS     Force            0.045425     0.000300     NO
Maximum Displacement          NaN     0.001800     YES
RMS     Displacement          NaN     0.001200     YES
Predicted change in Energy=          NaN
GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad

                          Input orientation:                          
---------------------------------------------------------------------
Center     Atomic      Atomic             Coordinates (Angstroms)
Number     Number       Type             X           Y           Z
---------------------------------------------------------------------
      1          6           0             NaN         NaN         NaN
      2          8           0             NaN         NaN         NaN
      3          8           0             NaN         NaN         NaN
      4         11           0             NaN         NaN         NaN
---------------------------------------------------------------------
                    Distance matrix (angstroms):
                    1          2          3          4
     1  C    0.000000
     2  O    0.000000   0.000000
     3  O    0.000000   0.000000   0.000000
     4  Na   0.000000   0.000000   0.000000   0.000000
Small interatomic distances encountered:      2     1     3     1     3     2     4     1     4     2
Small interatomic distances encountered:      4     3
Problem with the distance matrix.
Error termination via Lnk1e in /public1/software/gaussian09/g09/l202.exe at Sat Sep  6 11:09:40 2014.
Job cpu time:  0 days  0 hours  0 minutes 55.5 seconds.
File lengths (MBytes):  RWF=     83 Int=      0 D2E=      0 Chk=      2 Scr=      2



我用了各种基组和加不同关键词就是总遇到这个错误,构型也调整过了,也换过不同机器。
这是我投稿的一篇文章里计算出来的裂解产物COONa自由基。我原来算出来的是这种构型:
一个四原子自由基COONa怎么也优化不出来了!!!!!
所以我认为COONa会解离为CO2和Na原子,不能稳定存在。审稿人说可以出现这种稳定构型:
一个四原子自由基COONa怎么也优化不出来了!!!!!-1
并且给我列出了详细的图和几何参数,计算出来的频率(上图)。我就是怎么也优化不出审稿人说的构型,为什么啊,求救啊 啊;啊; !!!!!!

那位牛人帮我用 cbs-qb3 geom=connectivity nosym scf=(tight,maxcycle=5000)方法算一下啊,多谢了啊
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
坐标是对的
6楼2014-09-10 16:11:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

对不起,我看错了,坐标是对的。既然有好的对称性,不必加nosymm。
9楼2014-09-11 08:40:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

我觉得这个体系在电子结构上比较特殊,需要考虑电子相关。不能考虑电子相关的ub3lyp,优化时原子会四散奔逃。分子如此之小,建议用ccsd法。
11楼2014-09-11 15:31:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangmou 的主题更新
信息提示
请填处理意见