24小时热门版块排行榜    

查看: 2172  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 3/150 2026-08-15 12:41 by KxMI1BYBxWX1
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 3/150 2026-08-15 12:33 by KxMI1BYBxWX1
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 HFw0lei2R37i 2026-08-14 9/450 2026-08-15 11:33 by KxMI1BYBxWX1
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 09:28 by sunzitan
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 7/350 2026-08-15 00:40 by 4wMiSEwB6436
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见