计算化学公社

标题: 求助：阴离子聚合物的RESP2(0.5)电荷计算问题 [打印本页]

作者
Author: 陈欣羽 时间: yesterday 13:49
标题: 求助：阴离子聚合物的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、以上的步骤是否还有其他不合理的地方？
感谢各位老师的指导

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

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