拟合静电势电荷的计算方法

sobereva · 发表于 Post on 2014-11-9 08:15:00

这篇文章专门介绍拟合静电势电荷的原理、特点，对于MK、CHELPG、RESP、AM1-BCC等各种拟合静电势电荷都专门进行了介绍和讨论。

拟合静电势电荷的计算方法.pdf (50.53 KB, 下载次数 Times of downloads: 2178)

补充：后来笔者在Multiwfn程序（http://sobereva.com/multiwfn）中实现了RESP电荷的计算（MK、CHELPG也都支持），并撰写了下文，强烈建议一看
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441
计算RESP原子电荷的超级懒人脚本（一行命令就算出结果）
http://sobereva.com/476（http://bbs.keinsci.com/thread-12858-1-1.html）
RESP2原子电荷的思想以及在Multiwfn中的计算
http://sobereva.com/531（http://bbs.keinsci.com/thread-16190-1-1.html）
Multiwfn在计算RESP电荷上是所有程序中最灵活、最强大、使用最方便的。

数据挖掘 · 发表于 Post on 2015-3-3 09:19:28

非常感谢。
npa电荷 mullekin电荷和半径无关，设定半径计算chelpg电荷不，会改变npa和mullekin电荷

肖富贵 · 发表于 Post on 2016-1-7 00:00:51

Sob前辈，我会用到您的文章里提到的CHELP-BOW方法。此方法的课题组主页上有，但是链接失效，不能下载；和开发者联系，也无果。因此想自己写程序，请问这会不会是一个复杂的程序？

sobereva · 发表于 Post on 2016-1-7 00:07:49

肖富贵发表于 2016-1-7 00:00
Sob前辈，我会用到您的文章里提到的CHELP-BOW方法。此方法的课题组主页上有，但是链接失效，不能下载；和开 ...

勤快点的话估计也不难，主要就是引入了拟合点的权重函数，没法像一般拟合静电势电荷那样通过矩阵运算一下子就得到，而需要通过非线性优化算法来拟合，这有一些现成的库可以用。

肖富贵 · 发表于 Post on 2016-1-7 00:11:52

sobereva 发表于 2016-1-7 00:07
勤快点的话估计也不难，主要就是引入了拟合点的权重函数，没法像一般拟合静电势电荷那样通过矩阵运算一下 ...

多谢多谢多谢Sob前辈.

lty子于 · 发表于 Post on 2016-1-18 13:46:01

感谢老师的共享

zhiweigu123 · 发表于 Post on 2017-5-23 13:38:43

感谢您的分享！

meatball1982 · 发表于 Post on 2017-9-12 15:30:11

sob：
看了文章，学了很多。
现在，想看Ref73.Ref78.Ref86进一步阅读一下。

查了一下，请帮忙看一下是不是对应的73，86两个文献。
另外，能不能把Ref78的题目给一下。

Ref73
@article{doi:10.1021/j100494a019,
author = {Momany, Frank A.},
title = {Determination of partial atomic charges from ab initio molecular electrostatic potentials. Application to formamide, methanol, and formic acid},
journal = {The Journal of Physical Chemistry},
volume = {82},
number = {5},
pages = {592-601},
year = {1978},
doi = {10.1021/j100494a019},
}

Ref86

@article{Horn2006,
  doi = {10.1007/s00894-006-0137-8},
  url = {https://doi.org/10.1007/s00894-006-0137-8},
  year  = {2006},
  month = {aug},
  publisher = {Springer Nature},
  volume = {13},
  number = {2},
  pages = {381--392},
  author = {Anselm H. C. Horn and Timothy Clark},
  title = {Multipole electrostatic potential derived atomic charges in {NDDO}-methods with spd-basis sets},
  journal = {Journal of Molecular Modeling}
}

sobereva · 发表于 Post on 2017-9-12 16:29:44

meatball1982 发表于 2017-9-12 15:30
sob：
看了文章，学了很多。
现在，想看Ref73.Ref78.Ref86进一步阅读一下。

是。

[78] Franc, M. M., Chirlian, L. E. Reviews in Computational Chemistry. Lipkowitz, K. B.; Boyd, D. B., Eds. VCH Publishers: New York, 2000; Vol. 14, pp 1-31.

meatball1982 · 发表于 Post on 2017-9-12 20:35:32

sobereva 发表于 2017-9-12 16:29
是。

[78] Franc, M. M., Chirlian, L. E. Reviews in Computational Chemistry. Lipkowitz, K. B.; B ...

谢sob.

熊猫黑又白 · 发表于 Post on 2020-10-27 19:28:24

请问您的这篇文章发表了吗？想学习一下DMA方法

sobereva · 发表于 Post on 2020-10-29 13:13:03

熊猫黑又白发表于 2020-10-27 19:28
请问您的这篇文章发表了吗？想学习一下DMA方法

没
相关信息看MOLECULAR PHYSICS, 1985, VOL. 56, No. 5, 1047-1064

熊猫黑又白 · 发表于 Post on 2020-10-29 16:08:27

sobereva 发表于 2020-10-29 13:13
没
相关信息看MOLECULAR PHYSICS, 1985, VOL. 56, No. 5, 1047-1064

好的，非常感谢！

dgwfly · 发表于 Post on 2022-2-20 22:29:02

谢谢前辈的帖子。想请教一下，我想通过B48N48H24分子得到氮化硼的表面电荷。首先我采用B3LYP/631G(d)对B48N48H24分子进行结构优化，然后用Gaussian的pop=chelpg进行原子电荷的计算，分别测试了HF/jun-cc-pvdz 和wB97X-D/6-33++G(d,p)这两种方法。得到的位于分子最内环的B和N的平均电荷大概是1.6e，而大多文献中计算的氮化硼平板的电荷一般是小于1e的。我的结果有点偏大，不知道我的计算是不是存在问题。附件是我的输入文件和最后的电荷分布图。

dircumboronene.out (107.62 KB, 下载次数 Times of downloads: 1)

dircumboronene.gjf (7.17 KB, 下载次数 Times of downloads: 1)

wzkchem5 · 发表于 Post on 2022-2-21 01:04:51

dgwfly 发表于 2022-2-20 15:29
谢谢前辈的帖子。想请教一下，我想通过B48N48H24分子得到氮化硼的表面电荷。首先我采用B3LYP/631G(d)对B48N ...

文献也是CHELPG电荷吗？
把文献也贴一下吧

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

[波函数分析交流] 拟合静电势电荷的计算方法

评分 Rate