共聚物RESP电荷计算求助

lts！ · 发表于 Post on 4 day ago

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

sobereva · 发表于 Post on 3 day ago

lts！ · 发表于 Post on 前天 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运行输出文件放在这个帖子下面了，您可以看看我是哪里没有弄对吗

chgcons.txt (61 Bytes, 下载次数 Times of downloads: 1)

multiwfn_resp_out.txt (37.18 KB, 下载次数 Times of downloads: 4)

multiwfnresp.pbs (1.28 KB, 下载次数 Times of downloads: 0)

sobereva · 发表于 Post on yesterday 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拟合时甲基氢、甲基碳会被重新拟合。

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

[Multiwfn使用咨询] 共聚物RESP电荷计算求助