计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 隐式溶剂模型剧烈影响RESP结果

隐式溶剂模型剧烈影响RESP结果

查看数: 1203 | 评论数: 3 | 收藏 Add to favorites 0

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

帖子模式

ZetaFunction

发布时间: 2023-10-14 22:46

正文摘要:

对图中的分子拟合RESP电荷：先用B3LYP-D3/6-311g+(2d,p)以乙醚为溶剂优化了几何结构，然后计算单点能用以拟合RESP。按照RESP0.5的方法，应该分别计算气相和液相下的电荷然后取平均，我选用SMD的水模型，但是实 ...

回复 Reply

sobereva 发表于 Post on 2023-10-15 15:49:15

ZetaFunction 发表于 2023-10-15 03:31
感谢老师解答。

关于这两个甲基在气相下的巨大电荷差异是否可以如此理解：即拟合的目标是通过点电荷尽 ...

对于刚性区域，不用管拟合出来的原子电荷是多少、等价性如何，只要原子电荷能把周围的静电势描述好就够了、就达到实际目的了、就可以用于动力学模拟了

基组问题，到了ma-def2-TZVP程度就充分收敛了，能算得动就用这个就完了，不用再顾虑什么。和溶剂环境下差异大是溶剂环境的极化造成的，模拟一个分子在某种溶剂下的情况用RESP2(0.5)比RESP更强，避免高估等效的极化程度。对于模拟配体在蛋白质内部的情况，用隐式溶剂模型表现较低介电常数（如eps=4）的方式比较常见

ZetaFunction 发表于 Post on 2023-10-15 03:31:07

sobereva 发表于 2023-10-14 23:55
气相和溶剂环境下原子电荷差异大是正常的事情
化学等价的原子的电荷等价性差这是拟合静电势电荷的通病，和 ...

感谢老师解答。

关于这两个甲基在气相下的巨大电荷差异是否可以如此理解：即拟合的目标是通过点电荷尽可能的重现分子的静电势，所以这两个甲基即使中心碳原子电荷差异很大但是对外的整体静电势差异却很小？这两个甲基并不是完全等价，直接添加等价性约束似乎不太合适。

鉴于气相下电荷随着基组波动过大，和液相下的结果相差也太大，此时使用乙醚溶剂下的RESP是否是比RESP0.5更好的选择？有研究表明蛋白内部的介电常数和乙醚比较接近。

sobereva 发表于 Post on 2023-10-14 23:55:04

气相和溶剂环境下原子电荷差异大是正常的事情
化学等价的原子的电荷等价性差这是拟合静电势电荷的通病，和溶剂没直接关系，解决办法是在Multiwfn计算RESP电荷的时候设置等价性约束，或者用它们空间等价性更好的构象用Multiwfn算RESP电荷

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