24小时热门版块排行榜    

查看: 1786  |  回复: 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)方法算一下啊,多谢了啊
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangmou

木虫 (小有名气)

引用回帖:
11楼: Originally posted by zhou2009 at 2014-09-11 15:31:15
我觉得这个体系在电子结构上比较特殊,需要考虑电子相关。不能考虑电子相关的ub3lyp,优化时原子会四散奔逃。分子如此之小,建议用ccsd法。

我其他的分子都用的b3lyp,这个用ccsd做,相互比较的话有意义么,怎么解释呢,还请指教啊。
整个论文有几十个反应,审稿人让我用cbs方法做里面几个小分子的和b3lyp对比一下,我做了。就这一个两个方法都无法算。我本来也想把另外几十个原子数较多的反应一起用cbs做,但是实在算不动。
一个审稿人提了33个问题,呃呃呃呃。

[ 发自手机版 http://muchong.com/3g ]
12楼2014-09-12 00:10:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

gauss98

禁虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhou2009: 金币+8 2014-09-09 14:56:05
本帖内容被屏蔽

2楼2014-09-09 12:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangmou

木虫 (小有名气)

引用回帖:
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    ...

,b3lyp/6-311+g(d,p) 和 cbs-qb3这两个方法能优化出来么,我其他分子都是用的这种方法,所以想一致起来。
3楼2014-09-09 17:24:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gauss98

禁虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yangmou: 金币+10, ★★★很有帮助 2014-09-10 18:38:06
本帖内容被屏蔽

4楼2014-09-10 11:13:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 3/150 2026-11-11 17:47 by TgvkbZZl9tqB
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:14 by TgvkbZZl9tqB
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 16:47 by TgvkbZZl9tqB
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-11 15:24 by TgvkbZZl9tqB
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 C6CRn66RYYXS 2026-11-10 4/200 2026-11-11 14:24 by cdqGD6ftHmbl
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 YHiw1gou1Eda 2026-11-10 4/200 2026-11-11 14:13 by cdqGD6ftHmbl
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 FChKX8meu1ds 2026-11-09 4/200 2026-11-11 14:04 by cdqGD6ftHmbl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 14:02 by cdqGD6ftHmbl
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 13:56 by cdqGD6ftHmbl
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 13:53 by cdqGD6ftHmbl
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-11 13:39 by cdqGD6ftHmbl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 3/150 2026-11-11 12:05 by gRND4NVAZ8VG
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-11 09:20 by Equinoxhua
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
信息提示
请填处理意见