计算化学公社

标题: 开发water-hBN相互作用力场中电荷计算的求助 [打印本页]

作者
Author: hfdsakjhf 时间: 2022-11-23 14:51
标题: 开发water-hBN相互作用力场中电荷计算的求助
前辈们好，我目前在做水和二维材料氮化硼的力场拟合，结合能和滑移势垒面使用VASP的SCAN泛函计算，数据已计算完毕，在计算电荷是遇到了问题。问题是：
1.拟合力场需要添加water和hBN的库仑相互作用，但是B和N的电荷量该选用哪种电荷来计算？
2.目前我学习了使用VASP计算Bader电荷（AIM）、使用FHI-aims计算Hirshfeld电荷，发现两者差别很大。仔细阅读过卢天老师的文章“原子电荷计算方法的对比”，文章里老师提到“基于原子电荷的分子模拟任务需要原子电荷能够较好地描述静电相互作用. 通常来说以MK为代表的拟合静电势电荷是最适合用于分子模拟的.”，但不知MK电荷是否适合描述water-hBN的库仑相互作用，如果合适的话该如何计算？

作者
Author: hfdsakjhf 时间: 2022-11-23 15:01
计算Mulliken电荷算出来N原子大约带+0.01电荷，B原子大约带-0.01电荷

作者
Author: sobereva 时间: 2022-11-24 02:18
二维材料体系用REPEAT电荷是很好的选择，用CP2K可以直接计算，Multiwfn可以直接产生输入文件，看
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）
AIM和Hirshfeld电荷都绝对不能用

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