求助：阴离子聚合物的RESP2(0.5)电荷计算问题

陈欣羽

我的阴离子聚合物结构为AM-AM-AA-AM-AM-AM-AM-AM-AA-AA-AM-AM-AA-AM-AM-AA-AM-AM-AM-AM，带有5个负电荷的20聚物，想计算它的RESP2(0.5)电荷。AM为丙烯酰胺，AA为丙烯酸。
具体分为三部分：
1、AM-AM-AM:计算取中间-AM-的电荷和链首尾AM的电荷（结构优化：# opt b3lyp/6-311g(d,p) geom=connectivity em=gd3bj；单点计算：# b3lyp/6-311+g(d,p) em=gd3bj pop=mk iop(6/33=2,6/42=6)）
2、AM-AM-AA-AM-AM:计算取中间-AM-AA-AM-的电荷（结构优化：# opt=modredundant b3lyp/6-311g(d,p) scrf=(iefpcm,solvent=water) geom=connectivity em=gd3bj （gjf文件末尾：2 41 F）；单点计算：# b3lyp/6-311+g(d,p) em=gd3bj pop=mk  iop(6/33=2,6/42=6)）
3、AM-AM-AA-AA-AM-AM：计算取中间-AM-AA-AA-AM-的电荷（结构优化：# opt b3lyp/6-311g(d,p) scrf=(iefpcm,solvent=water) geom=connectivity em=gd3bj nosymm（gjf文件：首尾的AM坐标前加上-1冻死）；单点计算：# b3lyp/6-311+g(d,p) pop=mk geom=connectivity em=gd3bj iop(6/33=2,6/42=6)）
以上单点计算为气相单点的使用基组，若是计算液相单点会加入scrf=(iefpcm,solvent=water)
因为是带负电荷的聚合物，所以单点计算统一加上弥散函数。
之后导出的fch文件，在multiwfn中7-18-5-1-6-1-2（加载chgcons.txt和eqvcons.txt，并进行单步计算电荷）
问题：1、不带电荷的寡聚物AM-AM-AM的结构优化是否也要统一在溶剂中进行结构优化？
2、因为在计算带电荷的AM-AM-AA-AA-AM-AM用普通的结构优化总是发生部分蜷曲，试了很多种限制方法，最终只有直接冻结首尾的AM单体才能防止对中间的羧基基团产生影响。那是否所有都要统一，AM-AM-AM和AM-AM-AA-AM-AM都要将两端的AM给冻结？
3、如果冻结了两侧，计算两端的AM是不是不用冻结？
4、以上的步骤是否还有其他不合理的地方？
感谢各位老师的指导

sobereva

1 是
2 更优雅的做法是把链最两端的一对原子之间加上距离限制，从而避免优化时链收缩起来

陈欣羽

老师我之前AM-AM-AA-AM-AM利用给两端原子进行距离限制，我在末尾加了2 49 F，确实是成功了，但是我的AM-AM-AA-AA-AM-AM不管怎么限制都在链的左侧会有一个小的U型，虽然主链是直的，但支链互相吸引，下面的图是我今天做的结构优化前后的图，明显结构优化后链首支链聚集，外层的AM对内层的AA电荷产生影响。
所以我采用了直接冻结两端单体的方法，不知道老师是否还有更好的办法，如果只能采用冻结是否之前所有寡聚物的结构优化都得冻结收尾单体