| 查看: 6463 | 回复: 12 | ||||
feng0603新虫 (小有名气)
|
[求助]
lammps非周期性边界怎么实现压力控制已有3人参与
|
|||
| lammps中周期性边界可以通过fix npt、fix press/berendsen、fix box/relax等命令实现（上述命令只针对周期性边界条件）。故对于非周期性边界控压问题怎么实现？请大神指点！谢谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 分子动力学模拟 |
» 猜你喜欢
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有9人回复
-
心脉受损
已经有6人回复
-
博士申请都是内定的吗？
已经有8人回复
-
博士读完未来一定会好吗
已经有32人回复
-
读博
已经有5人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
Springer期刊投稿求助
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
请教：在lammps中非周期性边界下怎么进行结构弛豫
已经有4人回复- 周期性边界条件的影响
已经有9人回复
- lammps报错
已经有2人回复
- VMD由于周期性边界造成的显示问题
已经有6人回复
- 求一个LAMMPS里面聚合物适用的拉伸in文件
已经有2人回复
- Lammps运行结果问题！无数据，出现nan！
已经有4人回复
- Lammps计算油水界面，油相密度、平衡等问题求指导
已经有6人回复
- 关于用lammps建立无限延伸聚合物晶体模型的问题
已经有2人回复
- 关于石墨加压的lammps模拟
已经有1人回复
- 在LAMMPS中怎么设置冷却速度以及在非周期性边界条件方向上添加stress
已经有3人回复
- 一个关于非周期性边界条件的问题
已经有5人回复
- 【分享】LAMMPS的语法中文解释
已经有137人回复
- 【求助】Lammps 中加外力的问题
已经有1人回复
- 【求助】LAMMPS中非周期性边界条件和Born-Mayer-Huggins势函数不会设置
已经有4人回复
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
feng0603: 金币+2 2016-01-08 19:25:48
感谢参与，应助指数 +1
feng0603: 金币+2 2016-01-08 19:25:48
|
请仔细阅读lammps手册 http://lammps.sandia.gov/doc/fix_nh.html#fix-npt-command X, y, z cannot be barostatted if the associated dimension is not periodic. 当你的盒子在x,y,z三个方向上都是非周期性边界条件时，不能用fix npt，因为体积不可变，只能使用fix nvt |
3楼2016-01-07 12:02:46
feng0603
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2.9
- 散金: 12
- 帖子: 108
- 在线: 61.4小时
- 虫号: 2973885
- 注册: 2014-02-17
- 专业: 微/纳机械系统
2楼2016-01-06 19:17:26
bluewhale
铁杆木虫 (正式写手)
- 应助: 58 (初中生)
- 金币: 8438.3
- 散金: 50
- 红花: 10
- 帖子: 890
- 在线: 362.2小时
- 虫号: 278289
- 注册: 2006-09-10
- 专业: 凝聚态物性I:结构、力学和
4楼2016-01-07 19:13:45
wolfshow
铜虫 (初入文坛)
- 应助: 2 (幼儿园)
- 金币: 126.4
- 散金: 200
- 帖子: 28
- 在线: 17.9小时
- 虫号: 575281
- 注册: 2008-06-17
- 性别: GG
- 专业: 固体力学
5楼2016-01-08 15:36:37
feng0603
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2.9
- 散金: 12
- 帖子: 108
- 在线: 61.4小时
- 虫号: 2973885
- 注册: 2014-02-17
- 专业: 微/纳机械系统
6楼2016-01-08 19:12:30
feng0603
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2.9
- 散金: 12
- 帖子: 108
- 在线: 61.4小时
- 虫号: 2973885
- 注册: 2014-02-17
- 专业: 微/纳机械系统
7楼2016-01-08 19:21:31
feng0603
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2.9
- 散金: 12
- 帖子: 108
- 在线: 61.4小时
- 虫号: 2973885
- 注册: 2014-02-17
- 专业: 微/纳机械系统
8楼2016-01-08 19:24:43
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
|
假定你的表面是固定的，可以参照这个帖子里的做法 http://lammps.sandia.gov/threads/msg29869.html 请仔细阅读这个帖子，我只是粗略看了一下。 你需要定义两个group,一个叫immobile，代表固定的表面；一个叫mobile，代表可以自由运动的原子。 根据这个帖子，主要是根据这两行语句实现的 fix Ffreezestuff immobile rigid single force * off off off torque * off off off 这句话相当于把表面固定住。 fix Fmovestuff mobile npt temp 300.0 300.0 100.0 z 10000.0 10000.0 1000.0 dilate mobile 这一部分就是控压 如果你的自由表面内所有原子是活动的，这个会相对简单一些 boundary p p f fix 1 all npt 300.0 300.0 100.0 aniso 0.5 0.5 1000.0 这里使用aniso选项对x,y,z方向分别进行控压。 |
» 本帖已获得的红花（最新10朵）
feng06039楼2016-01-09 11:41:08
feng0603
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2.9
- 散金: 12
- 帖子: 108
- 在线: 61.4小时
- 虫号: 2973885
- 注册: 2014-02-17
- 专业: 微/纳机械系统
10楼2016-01-09 12:58:47