24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yangmou

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2220.5
红花: 2
帖子: 254
在线: 369小时
虫号: 555886
注册: 2008-05-09
专业: 工业心理学

[求助] 一个四原子自由基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会解离为CO2和Na原子，不能稳定存在。审稿人说可以出现这种稳定构型：

并且给我列出了详细的图和几何参数，计算出来的频率（上图）。我就是怎么也优化不出审稿人说的构型，为什么啊，求救啊啊；啊；！！！！！！

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

回复此楼

1楼2014-09-09 10:24:51

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

yangmou

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 2220.5
红花: 2
帖子: 254
在线: 369小时
虫号: 555886
注册: 2008-05-09
专业: 工业心理学

引用回帖:

2楼: Originally posted by gauss98 at 2014-09-09 12:48:00
用B3LYP-D3  6-311+G** 优化你的第二个结果
E(UB3LYP) =  -350.934411163
! R1 R(1,2)                1.2359       -DE/DX = -0.0001             !
! R2 R(1,3)                1.2359 ...