24小时热门版块排行榜    

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

noiradoris

新虫 (初入文坛)

[求助] Co的计算问题!新手求助

本人计算方面的新手,最近想做一下Co的羰基化合物的计算,但是对于Co(CO)3一直无法得到收敛结果,请各位大神帮忙看下问题处在哪里?

gif:# opt  ub3lyp  genecp

    0 2

Co
C                  1              B1
C                  1              B1    2             120.0
C                  1              B1    2             120.0    3              180.0
O                  4              B2   1               180.0   2               0.0
O                  3              B2    1              180.0    2              0.0
O                  2              B2    1              180.0    4             0.0

B1=1.84
B2=1.14

2 3 4 5 6 7 0
6-311G(d)
****
1 0
LANL2DZ
****

1 0
LANL2DZ

因为一直无法收敛 也有尝试过去增大循环步数,减小步长,添加gdiis等语句 但是始终无法收敛 问题出在哪里?
回复此楼

» 猜你喜欢

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

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

noiradoris

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 一起来化学 at 2012-04-28 14:48:09:
可以看看你的输出文件不?要不没法解决。。。。。

Error termination request processed by link 9999.
Error termination via Lnk1e in

Maximum Force            0.058108     0.000450     NO
RMS     Force            0.033504     0.000300     NO
Maximum Displacement     0.000923     0.001800     YES
RMS     Displacement     0.000545     0.001200     YES
Predicted change in Energy=-1.203569D-02
Optimization stopped.
    -- Number of steps exceeded,  NStep=  25
    -- Flag reset to prevent archiving.

之前有试过添加maxclcye=300 opt=gdiis语句 但是也不能正常收敛
3楼2012-04-28 15:02:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

noiradoris

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by jiewei at 2012-04-28 15:59:53:
没有看你的结构,但是有一个方法是可以尝试的, 先用较小的基组优化一遍,在优化结果上 在用你现在的基组优化

我最开始有用3-21G去优化 无法收敛 最后的键长大致停在1.9

然后用6-311G优化 也是无法收敛 最后的键长停留在1.7?

我是按照文献上给出的键长去优化的 最开始没有固定键角和二面角 无法收敛 后来固定了也无法收敛

真心不知道问题出在哪里
5楼2012-04-28 19:56:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

noiradoris

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by yjcmwgk at 2012-04-28 20:02:35:
换基组咯
最近我用SDD基组算Co蛮好的

啊我去试试

其实我是在试着做做Fe Co Ni的羰基络合物的比较 用的混合基组一直是6-311G(d)/LANL2DZ Ni和Fe都没什么问题 Co的四配位也被我用小基组到大基组的方法算下来了 现在就是一直卡在Co的三配位上 为什么会出这种问题
7楼2012-04-28 20:09:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

noiradoris

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by yjcmwgk at 2012-04-28 20:02:35:
换基组咯
最近我用SDD基组算Co蛮好的

啊我去试试

其实我是在试着做做Fe Co Ni的羰基络合物的比较 用的混合基组一直是6-311G(d)/LANL2DZ Ni和Fe都没什么问题 Co的四配位也被我用小基组到大基组的方法算下来了 现在就是一直卡在Co的三配位上 为什么会出这种问题
8楼2012-04-28 20:09:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

noiradoris

新虫 (初入文坛)

引用回帖:
9楼: Originally posted by jiewei at 2012-04-28 20:18:00:
你的O原子是否有问题? 我怎么复制你给我的文件 老是提示第一个氧原子有问题..

什么意思?

用gview打开没问题啊 文献上最后结果也是这种平面三角形结构
10楼2012-04-28 20:21:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

noiradoris

新虫 (初入文坛)

引用回帖:
9楼: Originally posted by jiewei at 2012-04-28 20:18:00:
你的O原子是否有问题? 我怎么复制你给我的文件 老是提示第一个氧原子有问题..

什么意思?

用gview打开没问题啊 文献上最后结果也是这种平面三角形结构
11楼2012-04-28 20:25:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 noiradoris 的主题更新
信息提示
请填处理意见