计算化学公社

标题: sobtop+GROMACS模拟MOF、COF体系时可以用的快速计算原子电荷的方法：PACMAN [打印本页]

作者
Author: sobereva 时间: 2024-6-6 23:03
标题: sobtop+GROMACS模拟MOF、COF体系时可以用的快速计算原子电荷的方法：PACMAN
sobtop（http://sobereva.com/soft/Sobtop）创建MOF和COF体系的GROMACS拓扑文件用的原子电荷我最推荐CP2K算的REPEAT电荷，但需要研究者会CP2K基本使用、安装CP2K。今天刚看到JCTC上有一篇新文章https://doi.org/10.1021/acs.jctc.4c00434，介绍了称为PACMAN的基于机器学习快速预测MOF和COF体系的原子电荷的方法，其中提供了在线预测程序，只要提供cif文件，就可以快速预测REPEAT、AIM、CM5等原子电荷。文中测试表现得不错，绝大多数预测的原子电荷和严格方式计算的相应类型的原子电荷的误差都在零点零几的程度。因此，不会用CP2K者也可以靠这个计算sobtop+GROMACS模拟MOF和COF的原子电荷了（但有条件的话还是争取以严格方式计算，免得这个预测方法在个别体系上得到误差较大的结果）。

网站是https://pacman-charge-mtap.streamlit.app/。上传cif即可在半分钟内得到含有原子电荷的cif文件。比如我上传了 (, 下载次数 Times of downloads: 58) ，选择重现REPEAT电荷，花了十几秒就返回了结果 (, 下载次数 Times of downloads: 42) 。这个文件里面可以看到这样的部分

loop_
_atom_site_type_symbol
_atom_site_label
_atom_site_symmetry_multiplicity
_atom_site_fract_x
_atom_site_fract_y
_atom_site_fract_z
_atom_site_occupancy
_atom_site_charge
Zn Zn1 1.0 0.29337 0.29337 0.20663 1.0000 1.210338
Zn Zn2 1.0 0.70663 0.70663 0.79337 1.0000 1.210338
Zn Zn3 1.0 0.70663 0.70663 0.20663 1.0000 1.201583
Zn Zn4 1.0 0.29337 0.29337 0.79337 1.0000 1.217656
Zn Zn5 1.0 0.70663 0.29337 0.79337 1.0000 1.210338
Zn Zn6 1.0 0.29337 0.70663 0.20663 1.0000 1.210338

复制代码

上面最后一列就是PACMAN预测的原子电荷。可以手动删除多余的列，整理成chg格式，即下面这样的文件。
(, 下载次数 Times of downloads: 172)

然后sobtop创建MOF或COF拓扑文件的时候，要求sobtop从这个chg文件里读取原子电荷即可。注意cif文件里的是分数坐标因此这个chg文件里的也是，所以sobtop目录下的sobtop.ini里的ichggeom应该设为0，从而要求sobtop读取chg文件的时候不把chg文件里的坐标替换当前坐标。

这个方法很新，我也尚没时间仔细研究PACMAN的原文，大家觉得结果靠谱就用，若某些情况得到的结果明显违背化学直觉就算了。

作者
Author: qmlearner 时间: 2026-3-6 18:28
老师请问下，我要研究MOF与聚合物的界面性质，需要建立MOF slab表面模型，就是沿着某个晶面切一个表面模型，另外我还想考察下在MOF 有机配体上接上一些官能团的情况，请问下这种slab模型的拓扑也可以用sobtop来创建吗？谢谢老师！

作者
Author: sobereva 时间: 2026-3-7 02:14

qmlearner 发表于 2026-3-6 18:28
老师请问下，我要研究MOF与聚合物的界面性质，需要建立MOF slab表面模型，就是沿着某个晶面切一个表面模型 ...

可以

作者
Author: qmlearner 时间: 2026-3-7 15:31

sobereva 发表于 2026-3-7 02:14
可以

谢谢老师：）

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