计算化学公社

标题: 优化配合物的时候配体总是没配上 [打印本页]

作者
Author: 喵帕斯 时间: 2015-12-29 15:37
标题: 优化配合物的时候配体总是没配上
本帖最后由喵帕斯于 2015-12-29 20:33 编辑

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

作者
Author: liyuanhe211 时间: 2015-12-29 19:40
可能解离所需的势垒很低，建议参考已知体系，确定一个近似的M-C(O2)键长，然后锁定该键长进行优化，优化后放开限制，使用较小的优化步长进行优化。

作者
Author: 喵帕斯 时间: 2015-12-29 20:19

liyuanhe211 发表于 2015-12-29 19:40
可能解离所需的势垒很低，建议参考已知体系，确定一个近似的M-C(O2)键长，然后锁定该键长进行优化，优化后 ...

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

作者
Author: root 时间: 2015-12-29 22:13
你把CO2竖着放上去试试，让O原子靠近钴，OCO键角稍微弯一下，150°吧，用这样的初始猜想优化看看会怎么样

作者
Author: liyuanhe211 时间: 2015-12-29 22:23

喵帕斯发表于 2015-12-29 20:19
冻结键长的方法我试过了，结果还是那样。不过倒是没有缩小步长，这种我再去试试，多谢啦

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

作者
Author: 喵帕斯 时间: 2015-12-30 01:13

root 发表于 2015-12-29 22:13
你把CO2竖着放上去试试，让O原子靠近钴，OCO键角稍微弯一下，150°吧，用这样的初始猜想优化看看会怎么样

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

作者
Author: 喵帕斯 时间: 2015-12-30 01:13

liyuanhe211 发表于 2015-12-29 22:23
冻结键长先优化只是为了减少小步长下的优化循环数，减少计算时间。如果有其他提供好的初猜的办法也一样。

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)