| 查看: 1067 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
急求，1.lammps中计算分子线的问题；2.是否以在计算过程中，改变原子的属性，如电荷。
|
||
|
请问，1.用lammps模拟一个分子晶体，boundary p p p 用npt; 然后，试图在x,y,方向上增加真空层，形成纳米线，用nvt 再算，可是总是提示“ERROR on proc 0: Bond atoms 1457 1479 missing on proc 0 at step 84 (../neigh_bond.cpp:65)”，这是不是因为change_box 之后，原子坐标不是unwrap的,导致键过长的原因。请问高人该如何实现这个过模拟。 2.在计算的过程中可否改变原子的属性，如电荷等。 请高手指点。 |
回复此楼» 猜你喜欢
国自然面上和省基金B类撒花
已经有17人回复
-
有没有学校收留
已经有3人回复
-
312求调剂
已经有3人回复
-
华师大读博
已经有5人回复
-
又一批高校组建人工智能学院 师资行吗 不是骗人吗
已经有5人回复
-
急需审稿人！！！
已经有3人回复
-
申博/考博
已经有8人回复
-
295分求调剂
已经有6人回复
-
085600材料与化工调剂
已经有6人回复
-
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
4楼2016-09-12 23:11:57
|
问题1中，换一种表述是这样的：如果在npt->nvt->npt->nvt……的过程中，在npt之后，用change_box 命令改变了box的大小，那接下来的nvt就会出错： ERROR on proc 0: Bond atoms 1457 1479 missing on proc 0 at step 84 (../neigh_bond.cpp:65) …… 这怎么解决？ 似乎是因为npt时部分原子穿过边界，由于周期性边界条件，原子跑去另一边了。如果这里加大了box然后nvt的话，穿越边界的键就相当于被拉长了，所以导致系统崩溃。然后，有没有一个类似于unwrap的命令，把穿过边界的原子放回来，然后再change_box？请高手指点 |
2楼2015-11-04 12:59:31
809011916
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 316.4
- 帖子: 55
- 在线: 42.3小时
- 虫号: 1501800
- 注册: 2011-11-20
- 性别: GG
- 专业: 化工热力学和基础数据
3楼2016-02-05 03:13:49
5