计算化学公社

标题: 共聚物RESP电荷计算求助 [打印本页]

作者
Author: lts！ 时间: 2026-5-28 20:30
标题: 共聚物RESP电荷计算求助
我需要构建一个20个由DA构成的聚合物链，现在已经对DADADA三个单元的寡聚物DFT计算完成，想问一下如何对这个寡聚物计算RESP，电荷算三次是直接写chgcons.txt对头部中部尾部的DA原子约束电荷总和为0然后运行一次multiwfn吗？求助

作者
Author: sobereva 时间: 2026-5-29 00:58
是

作者
Author: lts！ 时间: 2026-5-30 17:24
本帖最后由 lts！于 2026-5-30 17:26 编辑

sobereva 发表于 2026-5-29 00:58
是

老师，我就是按照电荷约束写了chgcons.txt，然后在里面设置了头部尾部中部电荷分别约束为0，在multiwfn输出时也输出了
Loading charge constraint setting from /home/ningbocheng/pm6topgro/PM6multiwfn/chgcons.txt
Charge constraint 1: 156 atoms, charge: 0.000000
Charge constraint 2: 157 atoms, charge: 0.000000
Charge constraint 3: 157 atoms, charge: 0.000000
No atom equivalence constraint is imposed in this fitting stage
然后最后得到的结果总电荷为0，但是三个部分的电荷总和都不为0，导致我在生成多聚体的时候电荷总量也不为0，我把multiwfn运行输出文件放在这个帖子下面了，您可以看看我是哪里没有弄对吗
(, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 1)

作者
Author: sobereva 时间: 2026-5-31 02:39

lts！发表于 2026-5-30 17:24
老师，我就是按照电荷约束写了chgcons.txt，然后在里面设置了头部尾部中部电荷分别约束为0，在multiwfn输 ...

进选项6时提示得非常清楚： Note: For standard two-stage RESP fitting, below options only take effect for the first stage

下文也明确说了
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441（http://bbs.keinsci.com/thread-10880-1-1.html）

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

作者
Author: lts！ 时间: 2026-6-1 10:05

sobereva 发表于 2026-5-31 02:39
进选项6时提示得非常清楚： Note: For standard two-stage RESP fitting, below options only take effec ...

是我没看仔细！非常感谢老师！

作者
Author: WLJFYQ 时间: 2026-7-1 21:27

sobereva 发表于 2026-5-29 00:58
是

请问sob老师，该帖子里面电荷算三次（下文粗体字所述内容），是否可以分成单独地三次Multiwfn计算RESP2操作呢？（分别约束头部，中间，尾部电荷为0）。我自己的体系，如果一次性同时约束三部分电荷均分别为0的话，multiwfn会报错不收敛。

我需要构建一个20个由DA构成的聚合物链，现在已经对DADADA三个单元的寡聚物DFT计算完成，想问一下如何对这个寡聚物计算RESP，电荷算三次是直接写chgcons.txt对头部中部尾部的DA原子约束电荷总和为0然后运行一次multiwfn吗？

作者
Author: sobereva 时间: 2026-7-2 02:07

WLJFYQ 发表于 2026-7-1 21:27
请问sob老师，该帖子里面电荷算三次（下文粗体字所述内容），是否可以分成单独地三次Multiwfn计算RESP2操 ...

可以。是

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