| 查看: 2595 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
有用lammps模拟过聚合物的相关性能吗？
|
|||
|
各位虫友，有人用lammps做过聚合物的相关性能吗？有几个问题请指教： 1.我在in里边定义了一个group， region 1 block INF 2 INF INF INF INF units box group left region 1 在运行之前，查看data文件，x坐标小于2的atom有好多，但是，运行的时候，log文件里显示0 atoms in group left。（从log文件里看，在这条语句之前，没有任何模拟的成分，只是读取数据。） 我可以用其他方法定义group，但是我不懂为什么还没有运行，原子的坐标就变化了。 2.在运行的过程中，经常会报 warning：dihedral problem： 查过手册，说Conformation of the 4 listed dihedral atoms is extreme; you may want to check your simulation geometry。 四个原子，总有坐标很大或者很小的，但是我不知道怎么才能避免这种情况，大家有没有什么经验，可以告诉我一下。 也许问题很不上道，但是这个对我真的很困扰，希望得到你的帮助，先谢谢了！ |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
国自然面上复盘~欢迎讨论
已经有15人回复
-
中青基了要发朋友圈吗？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【讨论】LAMMPS 在高性能计算机上测试，这样的结果正常吗
已经有5人回复
luxujie3
金虫 (著名写手)
- 应助: 54 (初中生)
- 金币: 4974.1
- 散金: 27
- 红花: 4
- 帖子: 1221
- 在线: 266.5小时
- 虫号: 756922
- 注册: 2009-04-24
- 性别: GG
- 专业: 高分子合成化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-10-25 09:00:56
xdr: 金币+5, ★★★很有帮助, 非常感谢！真心的。谢谢给我很有用的提供建议…… 2012-10-25 11:01:08
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-10-25 09:00:56
xdr: 金币+5, ★★★很有帮助, 非常感谢！真心的。谢谢给我很有用的提供建议…… 2012-10-25 11:01:08
| 我是使用lammps的，模拟高聚物一般都需要经过粗粒化的过程处理，且要找到合适的力场参数才能模拟得到比较理想的结果，看楼主的问题是说没有运行就发生坐标改变肯定是没有的事情我个人的经验就说明了，即使是使用的那个minimize命令初始构型的坐标也不会发生什么变化的。着说明一个问题是你定义group时候命令使用的不够准确 |
10楼2012-10-25 07:58:55
2楼2012-10-24 19:47:50
chaizhm
木虫 (著名写手)
- 模拟EPI: 2
- 应助: 123 (高中生)
- 贵宾: 0.25
- 金币: 4805.4
- 散金: 139
- 红花: 13
- 帖子: 1472
- 在线: 267.9小时
- 虫号: 1253209
- 注册: 2011-04-02
- 性别: GG
- 专业: 半导体材料
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-10-24 20:13:50
jiaoyixiong: 金币+2, 鼓励交流 2012-10-24 20:14:55
xdr: 金币+10, ★★★很有帮助, 谢谢，希望通过改结构解决，力场对于目前的我来说有点困难…… 2012-10-24 21:07:20
感谢参与，应助指数 +1
jiaoyixiong: 金币+2, 鼓励交流 2012-10-24 20:13:50
jiaoyixiong: 金币+2, 鼓励交流 2012-10-24 20:14:55
xdr: 金币+10, ★★★很有帮助, 谢谢，希望通过改结构解决，力场对于目前的我来说有点困难…… 2012-10-24 21:07:20
|
第二个问题加周期性边界条件了吗？ 警告是说二面角的问题，二面角是有四个不共面的原子构成的。我感觉是力场原则或者初始结构选择不恰当造成的 第一个问题还真不知道因为我不用lammps做计算 |
3楼2012-10-24 20:03:53
chaizhm
木虫 (著名写手)
- 模拟EPI: 2
- 应助: 123 (高中生)
- 贵宾: 0.25
- 金币: 4805.4
- 散金: 139
- 红花: 13
- 帖子: 1472
- 在线: 267.9小时
- 虫号: 1253209
- 注册: 2011-04-02
- 性别: GG
- 专业: 半导体材料
xdr4楼2012-10-24 20:05:54
投票:
30