计算化学公社

标题: 显式溶剂模型计算电荷转移 [打印本页]

作者
Author: wbn 时间: 2017-8-15 10:03
标题: 显式溶剂模型计算电荷转移
我想定性地了解一下离子液体阴阳离子间电荷转移的情况。请大家帮忙看一下下面的方法可不可行？
我已经有纯离子液体动力学模拟的轨迹，我随机在里面取一个阴离子，然后取出所有离这个阴离子一定距离以内的离子，适当加减离子使电荷为0，用PM6-D3优化整个团簇的结构，然后用B3LYP-D3/6-31+g**计算MK电荷。然后重复这个过程多次，取最中间阴离子的电荷平均值，来比较电子转移的程度。请问这个计算可行吗？估计最后可能有200-300个原子，都是C,H,N,O,我描述的DFT方法计算量现实不现实？

谢谢大家

作者
Author: sobereva 时间: 2017-8-15 12:14
两三百个原子单点算得动
不建议用MK电荷，因为拟合静电势电荷的通病就是被包埋的原子的拟合质量很低，因此没法合理体现电荷转移问题。用NPA、ADCH之类都比这合适
其实没必要加弥散，时间能省很多，对于合适的原子电荷计算方法而言，弥散对原子电荷计算结果也不会有太大影响

作者
Author: wbn 时间: 2017-8-15 22:14

sobereva 发表于 2017-8-15 12:14
两三百个原子单点算得动
不建议用MK电荷，因为拟合静电势电荷的通病就是被包埋的原子的拟合质量很低，因此 ...

好，谢谢SOB老师，那我就用ADCH电荷，然后把弥散去了

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