计算化学公社

标题: MOF的top文件如何获取电荷 [打印本页]

作者
Author: 邱美佳 时间: 2020-7-17 19:19
标题: MOF的top文件如何获取电荷
我得到了MOF的UFF力场参数后，想获取它的电荷。体系有四百多个原子，用高斯算感觉很吃力，不太算的动。请问能不能用vasp计算bader电荷的方法获取这个体系的电荷分布然后写进去top里呢？或者有没有其他更好的方法？

作者
Author: ggdh 时间: 2020-7-17 21:06
算单体啊，难道单体就400个原子？

作者
Author: sobereva 时间: 2020-7-18 10:33
用bader（AIM）是算这种问题最最最不行的电荷，拿这个电荷做模拟碰见稍微懂行的人立刻打回去
AIM对静电势重现性是所有原子电荷里几乎最垃圾的，显然绝对不能用于模拟计算，仔细看《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）里关于静电势重现性讨论。

抠出一个单体，边界恰当饱和，然后用Gaussian或ORCA等量化程序结合Multiwfn算RESP电荷，再把原子电荷填到整个体系所有原子。
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）

作者
Author: 邱美佳 时间: 2020-9-1 12:04

sobereva 发表于 2020-7-18 10:33
用bader（AIM）是算这种问题最最最不行的电荷，拿这个电荷做模拟碰见稍微懂行的人立刻打回去
AIM对静电势 ...

sob老师，我可不可以直接把MOF这四百个原子一次性丢进Multiwfn里面算EEM电荷？感觉算这个电荷挺快的

作者
Author: sobereva 时间: 2020-9-1 19:19

邱美佳发表于 2020-9-1 12:04
sob老师，我可不可以直接把MOF这四百个原子一次性丢进Multiwfn里面算EEM电荷？感觉算这个电荷挺快的

完全可以。但前提是内置的EEM参数有你体系里的所有元素（看手册3.9.15节的说明）。另外由于EEM经验性强，算完之后最好人工大致检查一下合理性

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