| 查看: 3602 | 回复: 2 | |||
1207150929金虫 (小有名气)
|
[求助]
关于用lammps建立无限延伸聚合物晶体模型的问题
|
|
我现在用lammps（17 Oct 2012 version）用非平衡分子动力学的方法做聚合物晶体聚合物热导率的计算。模拟建模时，通过拓展polymer的unit cell， 生成supercell，这个supercell即为模拟的box，力场用GAFF，周期性边界条件加在box的三个方向。 我的具体建模做法如下： 1、用antechamber处理gaussian的输出文件，生成polymer单胞内单体的mol2文件和frcmod文件； 2、用tleap工具将单体的力场信息load到超胞中，生成超胞的prmtop和inpcrd文件； 3、用amber2lammps将超胞的prmtop和inpcrd文件转化成lammps可读的data文件。 但由于我需要模拟无限延伸的polymer体系，故在box的边界总会有bond穿越box的边界。所以我在使用tleap工具将单体的力场信息load到超胞中时，就有意将 polymer链最两端的原子bond到一起。这样生成的data文件并无错误信息。但之后我测试用lammps模拟能量最小化时，在log输出文件中总会出现错误信息： Bond/angle/dihedral extent > half of periodic box length，模拟过程终止。我不知道是不是我建模的过程有问题；我也不确定lammps在模拟的过程中 是否遵循最小镜像原理，如果遵循的话，lammps应该就不会误认为将链最两端的原子bond到一起，而是一个box中的最边缘原子和它近邻box中的镜像bond到一起。 |
回复此楼» 猜你喜欢
27博士申请
已经有3人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
厉娜
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有12人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有15人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有16人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有15人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 分子动力学模拟的时间、体系大小和温度的争议问题！！
已经有23人回复
- 哪位大神能告诉我LAMMPS可以模拟聚合物的聚合过程吗？
已经有4人回复
- 能否给我解答下什么是周期性系统和非周期性系统啊？
已经有3人回复
- MS建模求助！！！
已经有3人回复
- 如图，以下模型怎么建立，详细步骤
已经有8人回复
- 晶体材料中的1D 2D 3D是什么意思？
已经有8人回复
- LAMMPS氯化钠晶体结构构建
已经有11人回复
- 关于分子模拟中建模的初级问题
已经有9人回复
- 活性炭怎么用MS建模？
已经有9人回复
- 可视化软件怎么选择
已经有2人回复
- 【讨论】周期性边界条件什么时候用
已经有11人回复
809011916
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 316.4
- 帖子: 55
- 在线: 42.3小时
- 虫号: 1501800
- 注册: 2011-11-20
- 性别: GG
- 专业: 化工热力学和基础数据
2楼2016-02-05 07:48:10
gaolimei55
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 415.2
- 散金: 20
- 帖子: 174
- 在线: 34小时
- 虫号: 3569470
- 注册: 2014-12-01
- 专业: 理论和计算化学
3楼2017-08-24 11:04:30
