计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 高斯算ESP电荷，同种原子的电荷不唯一，求助

高斯算ESP电荷，同种原子的电荷不唯一，求助

查看数: 20257 | 评论数: 7 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

小瓶子

发布时间: 2016-11-6 10:42

正文摘要:

根据文献内容，利用高斯软件求算分子的ESP电荷。先用Gaussianview限制分子的对称性为D3h（见下图），然后生成gjf文件进行计算。按理说，分子中的6个Cu应该处于同一种化学环境中，可是结果中的ESP电荷却是三种不同的 ...

回复 Reply

小瓶子 发表于 Post on 2016-11-9 21:40:35

非常感谢sob老师，我再好好读一下文章，谢谢啦！

sobereva 发表于 Post on 2016-11-9 20:50:27

小瓶子发表于 2016-11-9 20:00
还想请问sob老师，在上面那个分子中，对于苯环上连接H的C原子，文献得到的ESP电荷是负值（-0.014），而我计 ...

拟合静电势电荷用很多不同计算方法，看这里
http://bbs.keinsci.com/forum.php?mod=viewthread&tid=221
即便是同种拟合静电势电荷计算方法，格点设置不同、静电势计算级别和方法的不同、几何结构的细微不同等等，都会造成数值的不同。你的数值和文献的数值虽然符号不同，但差异很小，都很接近0，所以姑且可当成一样。

小瓶子 发表于 Post on 2016-11-9 20:00:35

还想请问sob老师，在上面那个分子中，对于苯环上连接H的C原子，文献得到的ESP电荷是负值（-0.014），而我计算所得的ESP电荷为正值（0.011），请问老师可能是什么原因呢？谢谢老师啦！

小瓶子 发表于 Post on 2016-11-7 10:23:22

谢谢sob老师指点

sobereva 发表于 Post on 2016-11-6 20:57:47

小瓶子发表于 2016-11-6 16:12
Cu半径本来只写了一次，但是优化结果提示，缺少Cu的范德华半径，故此又加了一次半径，然后就可以顺利优化得 ...

没有那种参数，至多只能取个平均

小瓶子 发表于 Post on 2016-11-6 16:12:42

Cu半径本来只写了一次，但是优化结果提示，缺少Cu的范德华半径，故此又加了一次半径，然后就可以顺利优化得到结果了。
听从sob老师的建议，查看了mulliken电荷，发现结果分子中的O原子有四种电荷（但我认为O应该只有一种，即一种电荷），而Cu只有一种电荷（Cu有一种，这个没问题）。请问在这种情况下，可以调整参数使得同种原子的ESP电荷保持一致吗？谢谢！

sobereva 发表于 Post on 2016-11-6 12:39:42

Cu 1.32没必要写两次

要想得到相同的拟合静电势电荷，必须原子的化学环境严格相同（不同的Cu与其它各个原子距离一致），而且电荷分布也一致（若波函数低于几何结构的对称性，此时显然没戏）。而且一般的拟合静电势方法由于格点分布原因，哪怕上两条都满足，等价的原子结果也可能因为数值因素而不同（G09里支持将格点连续化的拟合方法倒是应该没这个问题）。

建议你先看看这些Cu的mulliken电荷是否相同。这都不同的话，CHELPG电荷注定不可能相同。

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