“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register
Views: 57|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 求助:阴离子聚合物的RESP2(0.5)电荷计算问题

[复制链接 Copy URL]

13

帖子

0

威望

99

eV
积分
112

Level 2 能力者

我的阴离子聚合物结构为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、以上的步骤是否还有其他不合理的地方?
感谢各位老师的指导


本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-4-15 08:29 , Processed in 0.178433 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list