计算化学公社

标题: 原子电荷和静电势的关系 [打印本页]

作者
Author: 小菜鸟求关照 时间: 2015-11-3 20:47
标题: 原子电荷和静电势的关系
请问老师：在分子ABC中，A处的电荷和B处的电荷都是正电荷且A处的电荷<B处的电荷，但是A处的静电势要比B处的静电势更正，请问这是为什么呢？谢谢老师！

作者
Author: sobereva 时间: 2015-11-3 21:39
用对静电势重现性最好的CHELPG或MK电荷，看看是否与静电势观测结果相一致。

作者
Author: 小菜鸟求关照 时间: 2015-11-3 21:45

sobereva 发表于 2015-11-3 21:39
用对静电势重现性最好的CHELPG或MK电荷，看看是否与静电势观测结果相一致。

请问老师：我用的MP2方法，在用multiwfn结合VMD做静电势分析时，需要用到.wfx文件，对于用gauss生成.wfx文件时需要加什么特殊的关键词吗？ gauss版本是 Gaussian 09: ES64L-G09RevD.01

作者
Author: sobereva 时间: 2015-11-3 21:50

小菜鸟求关照发表于 2015-11-3 21:45
请问老师：我用的MP2方法，在用multiwfn结合VMD做静电势分析时，需要用到.wfx文件，对于用gauss生成.wfx ...

不建议用MP2，用DFT最划算

DFT计算直接out=wfn就完了
MP2需要写density out=wfn

wfx的效果和wfn一样，用它的话out=wfx

作者
Author: 小菜鸟求关照 时间: 2015-11-3 21:51

sobereva 发表于 2015-11-3 21:50
不建议用MP2，用DFT最划算

DFT计算直接out=wfn就完了

为什么不建议用mp2呢？我这个版本默认生成.wfx时就是含有density关键词的吧

作者
Author: sobereva 时间: 2015-11-3 21:52

小菜鸟求关照发表于 2015-11-3 21:51
为什么不建议用mp2呢？

耗时比DFT多啊
而且也并不比DFT结果有明显优势

对于考察静电势，B3LYP/6-311G**就够了，如果有的原子带显著负电荷，再给相应原子加上弥散就足够了。

作者
Author: 小菜鸟求关照 时间: 2015-11-3 21:59

sobereva 发表于 2015-11-3 21:52
耗时比DFT多啊
而且也并不比DFT结果有明显优势

我那个gauss版本默认生成.wfn时用的是mp2density吧，不需要再加density关键词了吧

作者
Author: sobereva 时间: 2015-11-3 22:14

小菜鸟求关照发表于 2015-11-3 21:59
我那个gauss版本默认生成.wfn时用的是mp2density吧，不需要再加density关键词了吧

是的。

作者
Author: 虎王 时间: 2015-11-4 13:26

sobereva 发表于 2015-11-3 21:39
用对静电势重现性最好的CHELPG或MK电荷，看看是否与静电势观测结果相一致。

也就是说尽管不是同一个概念，电荷反映的规律与静电势反映的规律原则上应该一致？

作者
Author: sobereva 时间: 2015-11-4 16:06

虎王发表于 2015-11-4 13:26
也就是说尽管不是同一个概念，电荷反映的规律与静电势反映的规律原则上应该一致？

对

但注意毕竟原子电荷的模型过于简单，没法把分子所有性质特征都照顾到，完全为了使静电势重现性好而提出的原子电荷，在其它方面就可能没那么好。

作者
Author: 小菜鸟求关照 时间: 2015-11-10 19:36

sobereva 发表于 2015-11-4 16:06
对

但注意毕竟原子电荷的模型过于简单，没法把分子所有性质特征都照顾到，完全为了使静电势重现性好而 ...

老师：我计算了CHELPG电荷之后，和静电势重现性无法一致，电荷仍然是上述情况，这个原因是什么呢？谢谢老师

作者
Author: 小菜鸟求关照 时间: 2015-11-10 20:52

sobereva 发表于 2015-11-4 16:06
对

但注意毕竟原子电荷的模型过于简单，没法把分子所有性质特征都照顾到，完全为了使静电势重现性好而 ...

老师：您在这个帖子的回复说的好像就是http://bbs.keinsci.com/forum.php ... =%BE%B2%B5%E7%CA%C6
原子电荷和静电势是两个概念吧，它们之间的关系好像不是一致的

作者
Author: sobereva 时间: 2015-11-10 21:05

小菜鸟求关照发表于 2015-11-10 20:52
老师：您在这个帖子的回复说的好像就是http://bbs.keinsci.com/forum.php?mod=viewthread&tid=451&highl ...

对。

不一致就说明你对比原子电荷和静电势的方式不对。

作者
Author: 小菜鸟求关照 时间: 2015-11-10 21:26

sobereva 发表于 2015-11-10 21:05
对。

不一致就说明你对比原子电荷和静电势的方式不对。

对比的方式不对是什么意思呢？

作者
Author: sobereva 时间: 2015-11-10 23:22

小菜鸟求关照发表于 2015-11-10 21:26
对比的方式不对是什么意思呢？

一定要弄清楚拟合静电势电荷的原理，静电势精确计算（基于密度计算）时的表达式，以及基于原子电荷计算时的表达式，自然而然就能理解体现问题的异同

作者
Author: 小菜鸟求关照 时间: 2015-11-14 16:10

sobereva 发表于 2015-11-10 23:22
一定要弄清楚拟合静电势电荷的原理，静电势精确计算（基于密度计算）时的表达式，以及基于原子电荷计算 ...

老师您好，原子电荷计算方法的对比我看到您这篇文献中，在静电势重现性的计算中有一个这样的公式，这个公式中每一个量的数值是怎么确定的呢？谢谢老师

作者
Author: sobereva 时间: 2015-11-14 17:25
这篇你也看看
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=221

Vreal是通过量化算出来的电子密度得到的静电势，这是当前级别下精确的静电势。Vmodel是通过原子电荷计算的静电势，是近似的静电势。通过最小化这个，就能让原子电荷尽可能好地重现当前级别下精确的静电势。

作者
Author: 小菜鸟求关照 时间: 2015-11-14 18:53

sobereva 发表于 2015-11-14 17:25
这篇你也看看
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=221

好的谢谢老师

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