计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12205|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 优化配合物的时候配体总是没配上

[复制链接 Copy URL]

29

帖子

0

威望

1870

eV
积分
1899

Level 5 (御坂)

本帖最后由 喵帕斯 于 2015-12-29 20:33 编辑

最近在算一个卟啉钴的体系,其中一个中间产物是M-CO2-的形式,结果优化完成后,CO2总是配不上。另外有意思的是,另一个中间产物M-COOH倒是配得很好,和文献给出的结构一致。现在想请问一下这种情况该怎么处理?忘了说了,用的是b3lyp/6-31+g(d,p),钴原子用的是lanl2dz基组。

1.png (26.58 KB, 下载次数 Times of downloads: 48)

优化前

优化前

2.png (32.6 KB, 下载次数 Times of downloads: 45)

优化后

优化后

3098

帖子

28

威望

1万

eV
积分
16892

Level 6 (一方通行)

2#
发表于 Post on 2015-12-29 19:40:32 | 只看该作者 Only view this author
可能解离所需的势垒很低,建议参考已知体系,确定一个近似的M-C(O2)键长,然后锁定该键长进行优化,优化后放开限制,使用较小的优化步长进行优化。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

29

帖子

0

威望

1870

eV
积分
1899

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2015-12-29 20:19:33 | 只看该作者 Only view this author
liyuanhe211 发表于 2015-12-29 19:40
可能解离所需的势垒很低,建议参考已知体系,确定一个近似的M-C(O2)键长,然后锁定该键长进行优化,优化后 ...

冻结键长的方法我试过了,结果还是那样。不过倒是没有缩小步长,这种我再去试试,多谢啦

24

帖子

0

威望

358

eV
积分
382

Level 3 能力者

4#
发表于 Post on 2015-12-29 22:13:47 | 只看该作者 Only view this author
你把CO2竖着放上去试试,让O原子靠近钴,OCO键角稍微弯一下,150°吧,用这样的初始猜想优化看看会怎么样

3098

帖子

28

威望

1万

eV
积分
16892

Level 6 (一方通行)

5#
发表于 Post on 2015-12-29 22:23:07 | 只看该作者 Only view this author
喵帕斯 发表于 2015-12-29 20:19
冻结键长的方法我试过了,结果还是那样。不过倒是没有缩小步长,这种我再去试试,多谢啦

冻结键长先优化只是为了减少小步长下的优化循环数,减少计算时间。如果有其他提供好的初猜的办法也一样。

29

帖子

0

威望

1870

eV
积分
1899

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2015-12-30 01:13:32 | 只看该作者 Only view this author
root 发表于 2015-12-29 22:13
你把CO2竖着放上去试试,让O原子靠近钴,OCO键角稍微弯一下,150°吧,用这样的初始猜想优化看看会怎么样

不过这个结构和文献给出的结构不一致,所以还是很矛盾

29

帖子

0

威望

1870

eV
积分
1899

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2015-12-30 01:13:53 | 只看该作者 Only view this author
liyuanhe211 发表于 2015-12-29 22:23
冻结键长先优化只是为了减少小步长下的优化循环数,减少计算时间。如果有其他提供好的初猜的办法也一样。

嗯,已经提交上去了,感谢帮助

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-27 18:27 , Processed in 0.169121 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list