Benzene-Cu2+优化时分离的问题

Augustan · 发表于 Post on 2015-10-30 21:50:19

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

Augustan · 发表于 Post on 2015-11-2 09:02:35

谢谢各位的回答

helpme · 发表于 Post on 2015-10-31 11:14:12

距离很远的时候，计算是很不准的。
你以为是 A + （Cu2+），结果可能算的是 A+ + Cu+, 或者分数占据。

Augustan · 发表于 Post on 2015-10-31 10:10:46

还有就是，优化到距离非常远（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能量低）

sobereva · 发表于 Post on 2015-10-31 05:36:02

本来就应该跑开。H带正电荷，Cu2+也带正电荷，静电互斥很强，肯定呆不住。

liyuanhe211 · 发表于 Post on 2015-10-30 23:10:58

放在侧面你想得到什么结构？跑开了是还在苯平面上还是跑到了苯环上面

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] Benzene-Cu2+优化时分离的问题

评分 Rate