计算化学公社

标题: Benzene-Cu2+优化时分离的问题 [打印本页]

作者
Author: Augustan 时间: 2015-10-30 21:50
标题: Benzene-Cu2+优化时分离的问题

我把Cu2+放在了苯平面上（在C-H键的延长线上），而不是苯环中心的正上方。当我在用b3lyp/avtz优化该体系（+2, 2）时，cu2+竟然跑开了。
有没有谁知道原因啊
先谢谢了！

作者
Author: liyuanhe211 时间: 2015-10-30 23:10
放在侧面你想得到什么结构？跑开了是还在苯平面上还是跑到了苯环上面

我个人理解可能会这么做：
考虑到这个结构的朝向“分解”的能垒可能很小，可尝试CalcFC甚至（考虑到体系较小）CalcAll，并尝试减小步长。
计算ion-pi弱相互作用可以考虑M06-2X-D3、wB97XD。

作者
Author: sobereva 时间: 2015-10-31 05:36
本来就应该跑开。H带正电荷，Cu2+也带正电荷，静电互斥很强，肯定呆不住。

作者
Author: Augustan 时间: 2015-10-31 10:10
还有就是，优化到距离非常远（100 angstrom）也没有收敛

Maximum Force          0.000352    0.000450    YES
RMS    Force          0.000083    0.000300    YES
Maximum Displacement    0.117222    0.001800    NO
RMS    Displacement    0.019233    0.001200    NO

另外就是，在最后的结构（未收敛，Cu2+ -- Benzene距离为100A左右）的能量与单独Cu2+ ,单独benzene的能量之和差别非常大有
200多千卡（dimer能量低）

作者
Author: helpme 时间: 2015-10-31 11:14
距离很远的时候，计算是很不准的。
你以为是 A + （Cu2+），结果可能算的是 A+ + Cu+, 或者分数占据。

作者
Author: Augustan 时间: 2015-11-2 09:02
谢谢各位的回答

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