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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 一个四原子自由基COONa怎么也优化不出来了！！！！！
181/2返回列表12下一页
查看: 1341  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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朵）

yangmou

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-09-09 10:24:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

yangmou

木虫 (小有名气)
引用回帖:
9楼: Originally posted by zhou2009 at 2014-09-11 08:40:29
对不起，我看错了，坐标是对的。既然有好的对称性，不必加nosymm。

不加nosymm也计算不出来呢，优化后原子都四散奔逃了，求解救啊，给出ub3lyp/6-311+g(d,p)的能量和零点能矫正值
赞一下(1人) 回复此楼
10楼2014-09-11 12:05:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

gauss98

禁虫 (正式写手)
本帖内容被屏蔽
5楼2014-09-10 11:14:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

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

yangmou

木虫 (小有名气)
引用回帖:
6楼: Originally posted by zhou2009 at 2014-09-10 16:11:48
你提问时列出的内坐标全是错的，将坐标在GV中打开就可看见不是你想要的。

为什么是错误的呢，我用GV打开没发现呢？检查键长和键角没发现不太合理的地方呢，求指教，确实卡在这里了。
赞一下 回复此楼
7楼2014-09-10 22:49:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangmou

木虫 (小有名气)
引用回帖:
5楼: Originally posted by gauss98 at 2014-09-10 11:14:59
Temperature=               298.150000 Pressure=                       1.000000
E(ZPE)=                      0.010225 E(Thermal)=                     0.014536
E(SCF)=                   -349.537821 ...

%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
赞一下 回复此楼
8楼2014-09-10 22:50:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)
对不起，我看错了，坐标是对的。既然有好的对称性，不必加nosymm。
赞一下 回复此楼
9楼2014-09-11 08:40:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangmou 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +4 @taotao 2026-03-20 4/200 2026-03-20 10:37 by EBSD
[考研] 085600材料与化工 求调剂 +17 enenenhui 2026-03-13 18/900 2026-03-20 08:55 by xingguangj
[考研] 295材料求调剂，一志愿武汉理工085601专硕 +3 Charlieyq 2026-03-19 3/150 2026-03-20 08:53 by xingguangj
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考博] 东华理工大学化材专业26届硕士博士申请 +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by 轻松不少随
[考研] 一志愿中海洋材料工程专硕330分求调剂 +7 小材化本科 2026-03-18 7/350 2026-03-19 10:46 by Linda Hu
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 一志愿华中科技大学，080502，354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 11408 一志愿西电，277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
信息提示
关闭
请填处理意见
关闭