| 查看: 2845 | 回复: 2 | ||||
[求助]
lammps模拟压缩过程 已有1人参与
|
|
各位朋友大家好，我是一名自学lammps的新手，有个问题想请教一下各位，先谢谢各位。 是这样的，我是用lammps做拉伸压缩的，之前在一篇文章上有看到有人先在Z方向进行了一个预压缩，然后再在Y方向进行压缩模拟。 他的条件是这样的：先将模拟物体分为上中下三部分，固定下端，移动上端刚性原子，用的是Berendsen控温50K，Y和Z方向是周期性边界条件，X方向是自由表面，10的-8次方应变速率进行小变形的预压缩。然后在Y方向进行一个压缩模拟。 我的问题是：1、他用的是Berendsen控温，是不是说明他用的不是NPT系综，而是NVE系综+Berendsen控温组成NVT系综？2、他将模拟的物体分为了三段，据我所知这应该是用velocity进行压缩的一个前设置，但是如果他用的是velocity的话，他YZ方向都是周期性边界条件，不是会出现压不动的情况么，从一边出去就会从另一边进来，这块我学的不是很明白，望老师能够指教一下。如果不是用velocity的话，那就是用的deform，可是用deform的话，那么何必将模拟物体分为三段呢？ 新手请教一下，谢谢各位帮忙解答一下或者互相交流学习一下啦。 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有9人回复
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
【答案】应助回帖
|
按你所说，该文献陈述的确是令人不解。建议按照自己的理解来进行模拟。 1、他用的是Berendsen控温，是不是说明他用的不是NPT系综，而是NVE系综+Berendsen控温组成NVT系综？ 差不多可以这样理解，只是说法不专业。其实这种非平衡过程，谈不上系综，就是控温而已。 2、他将模拟的物体分为了三段，据我所知这应该是用velocity进行压缩的一个前设置，但是如果他用的是velocity的话，他YZ方向都是周期性边界条件，不是会出现压不动的情况么，从一边出去就会从另一边进来，这块我学的不是很明白，望老师能够指教一下。如果不是用velocity的话，那就是用的deform，可是用deform的话，那么何必将模拟物体分为三段呢？ 我也觉得不好理解。可能文章描述不清。 |
2楼2017-12-30 02:20:43
3楼2018-01-01 15:18:42
