计算化学公社

标题: 在构建聚合物时重复单元的电荷如何限定才算合理 [打印本页]

作者
Author: Weijie 时间: 2023-3-13 09:51
标题: 在构建聚合物时重复单元的电荷如何限定才算合理
请教一下各位前辈，我之前在构建聚合物时向人请教时得出需要在B3LYP/6-31G**级别下进行RESP得到的电荷才算能接受，随后又问询得知在得到电荷后需要进行电荷限定，不可自己手动调节，所以根据Multiwfn中氨基酸残基的教程进行了电荷限定的限制，但是在限制后，所限制的原子序号范围内的原子电荷加和并没有实现完全匹配到原本限定的电荷数，总会多一点或者少一点，不知这种情况是否正常，我随后是否还需要进一步进行电荷调整，以实现完全限定呢？

我注意到一点的是，官网给出的教程进行限定后所得的原子电荷总和也是所构建基元的所有原子电荷总和符合限定，而不是限定的原子范围内的电荷符合限定条件

(, 下载次数 Times of downloads: 18)
如图，我所限定的是1-32号原子，限定电荷用chgcons.txt总和为0（限定内容为 1-32 0）

作者
Author: sobereva 时间: 2023-3-13 21:03

官网给出的教程进行限定后所得的原子电荷总和也是所构建基元的所有原子电荷总和符合限定，而不是限定的原子范围内的电荷符合限定条件

没有这回事

RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）
里面天冬氨酸的例子里，1-12号原子电荷加和就是用户自己设的-1

只可能是你自己弄错了，或者选择选的是两步式RESP电荷拟合，此时没法用电荷约束。博文里该说的都明确说了

在标准两步式RESP电荷计算时也照样可以设置电荷约束、等价性约束，但只是对第一步起效（默认情况下，第一步的时候既没有用任何等价性约束也没有使用电荷约束）。对于此例，载入上述eqvcons.txt后选择选项0做两步式RESP拟合，会发现最终结果里O5和O6的电荷相同了，但是C7、C11的电荷不同，不同甲基上的氢的电荷也不同。这是因为自定义约束对于第二步不生效，而且按照定义，第二步RESP拟合时甲基氢、甲基碳会被重新拟合。

作者
Author: Weijie 时间: 2023-3-14 10:18

sobereva 发表于 2023-3-13 21:03
没有这回事

RESP拟合静电势电荷的原理以及在Multiwfn中的计算

老师，是我的这边关注的问题，我只关注了博文中总电荷，没有细去核算1-12的电荷之和，这是我的疏忽，非常抱歉。
我确实是做的标准两步RESP拟合，那是做一步的ESP就可以用了嘛

作者
Author: sobereva 时间: 2023-3-14 21:35

Weijie 发表于 2023-3-14 10:18
老师，是我的这边关注的问题，我只关注了博文中总电荷，没有细去核算1-12的电荷之和，这是我的疏忽，非常 ...

是

作者
Author: Weijie 时间: 2023-3-15 09:25

sobereva 发表于 2023-3-14 21:35
是

好的老师，谢谢老师啦

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