开发water-hBN相互作用力场中电荷计算的求助

hfdsakjhf · 发表于 Post on 2022-11-23 14:51:58

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

hfdsakjhf · 发表于 Post on 2022-11-23 15:01:44

计算Mulliken电荷算出来N原子大约带+0.01电荷，B原子大约带-0.01电荷

sobereva · 发表于 Post on 2022-11-24 02:18:20

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

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

[综合交流] 开发water-hBN相互作用力场中电荷计算的求助

浏览过的版块