计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 关于计算RESP2电荷的一个疑问

关于计算RESP2电荷的一个疑问

查看数: 3804 | 评论数: 4 | 收藏 Add to favorites 1

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

帖子模式

a9471163

发布时间: 2021-12-24 18:09

正文摘要:

本帖最后由 a9471163 于 2021-12-24 18:18 编辑各位老师，我想做氧化石墨烯（GO）在OPLS-AA力场下的分子动力学模拟，由于羧基会去质子化形成带电体系，我选择RESP2（0.5）电荷作为原子电荷。首先，我用xtb在ex ...

回复 Reply

a9471163 发表于 Post on 2021-12-24 19:18:40

eagletyr 发表于 2021-12-24 19:03
你在MD的时候是去质子化的，拟合RESP电荷的时候肯定就要用去质子化的模型来计算

感谢大佬解惑，计算真空里RESP电荷时模型选错了，应该改为与水溶液里使用的模型完全一致，并且施加-2的总电荷

eagletyr 发表于 Post on 2021-12-24 19:03:32

你在MD的时候是去质子化的，拟合RESP电荷的时候肯定就要用去质子化的模型来计算

a9471163 发表于 Post on 2021-12-24 18:26:45

本帖最后由 a9471163 于 2021-12-24 18:30 编辑

喵星大佬发表于 2021-12-24 18:24
一定哪里操作出问题了，不管用不用溶剂模型，RESP电荷总数肯定是你的净电荷

您好。我没加溶剂模型时，是带着羧基上的H原子去计算了RESP电荷，这时候羧基的H原子化学价是+1，但实际的RESP电荷肯定是比+1小的，算出来就是+0.36左右，这里不是整数，然后整个体系的RESP电荷是0，是整数，减去羧基H原子的RESP电荷以后，整个体系的RESP电荷就不是整数了，问题就是在这一步出现的。

我知道只需加个电荷约束，让羧基的H原子带的电量是+1就能解决问题，但合理吗，羧基的H原子化学价是+1，但实际的RESP电荷肯定是比+1小的。我规定成+1，感觉确实有点勉强

喵星大佬 发表于 Post on 2021-12-24 18:24:07

一定哪里操作出问题了，不管用不用溶剂模型，RESP电荷总数肯定是你的净电荷

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