计算化学公社

标题: 请问如何将计算得到的RESP电荷分配到构建的聚合物上 [打印本页]

作者
Author:
高处裹棉被    时间: 2022-3-14 20:16
标题: 请问如何将计算得到的RESP电荷分配到构建的聚合物上
本帖最后由 高处裹棉被 于 2022-3-14 20:16 编辑

各位老师好,我使用Avogadro建立cml文件后,采用PolyParGen脚本后在OPLS-AA力场下创建了3个重复单元的PET聚合物,并选择了默认的mulliken电荷,得到了PET的itp文件

然后我使用了sob老师给的高斯懒人脚本在gas环境下计算3个重复单元PET的RESP电荷,并得到chg文件,现在我想用UltraEdit将其charge复制黏贴过去,但请问我该怎么一一对应上?
谢谢各位老师了




作者
Author:
sobereva    时间: 2022-3-15 14:25
确保3个重复单元的寡聚物模型里头部、中部、尾部的原子顺序和你实际聚合物里三类单元的原子顺序一样。末端单元的电荷直接用ultraedit列模式复制过去,中间的自己写个脚本,重复N输出N次,N是中间的重复单元数,然后挪过去就完了




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3