使用RESP电荷计算静电势的一些问题

hzz_xizhe · 发表于 Post on 2025-5-9 10:34:18

本帖最后由 hzz_xizhe 于 2025-5-9 10:33 编辑

请教一下老师，我研究的是HI-(H2O)n团簇，属于弱相互作用体系。

近期阅读了一篇文献，文献作者计算了团簇中所有的水分子在H和I两个原子位置之间产生的电场的值，之后对计算得到的电场进行拟合，发现解离的构型的电场都大于某个阈值，文献作者将此阈值称为临界电场。具体做法是在H和I两个位置设置鬼原子（不考虑H和I两个原子产生的影响），之后通过Gaussian计算H和I两个位置的静电势，之后对这两个位置的电势做差再除以它们之间的距离，求出了它们之间的电场。

我想复现这个结果，在sob老师的谈谈怎么计算化学体系中“原子的静电势” - 思想家公社的门口：量子化学·分子模拟·二次元帖子中，通过第三个方法“计算原子核位置的静电势”成功复现了临界电场的数值，具体做法也是在H和I两个位置设置鬼原子，然后用Multiwfn分析Gaussian计算产生的fchk文件，得到H和I两个位置的原子核位置的静电势，之后对这两个位置的电势做差再除以它们之间的距离，求出了它们之间的电场。

之后，我想用sob老师提到的第一种方法“原子的拟合静电势电荷”，也是在H和I两个位置设置了鬼原子，然后计算RESP拟合静电势电荷，得到了关于电荷的.chg文件，之后我用Multiwfn的主功能3去计算基于原子电荷产生的静电势，通过设置H和I的坐标为两个端点的坐标来绘制曲线图，但是我发现在I位置计算出来的静电势非常大有1000a.u.，不符合实际情况，想请教一下老师这是什么原因导致的？还是说这个做法不合理？

这是我用Multiwfn画出来的图像，是计算HI-(H2O)3团簇的某一个构型，具体的文件都上传了。

sobereva · 发表于 Post on 2025-5-9 16:48:57

在Multiwfn里计算任意一个点的静电势完全不需要设鬼原子，设鬼原子完全多此一举

基于拟合静电势电荷算原子核位置的静电势毫无意义，数值是无穷大的，因为1/r的分母为0。为了避免数值无穷大，Multiwfn刻意把数值设成了1000。

Stardust0831 · 发表于 Post on 2025-5-9 20:34:10

感觉你会对这个感兴趣《使用Multiwfn计算分子片段的偶极矩和复合物中单体的偶极矩》（http://sobereva.com/558）

hzz_xizhe · 发表于 Post on 2025-5-10 14:41:11

sobereva 发表于 2025-5-9 16:48
在Multiwfn里计算任意一个点的静电势完全不需要设鬼原子，设鬼原子完全多此一举

基于拟合静电势电荷算原 ...

谢谢sob老师的回答

hzz_xizhe · 发表于 Post on 2025-5-10 14:41:32

Stardust0831 发表于 2025-5-9 20:34
感觉你会对这个感兴趣《使用Multiwfn计算分子片段的偶极矩和复合物中单体的偶极矩》（http://sobereva.com/ ...

谢谢您

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

[波函数分析求助] 使用RESP电荷计算静电势的一些问题

浏览过的版块