24小时热门版块排行榜    

查看: 1250  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

superbuff

新虫 (初入文坛)

[求助] 急求,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.在计算的过程中可否改变原子的属性,如电荷等。
请高手指点。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

时间就是用来浪费的,不是浪在这儿,就得费在那儿!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

809011916

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by superbuff at 2015-11-04 12:59:31
问题1中,换一种表述是这样的:如果在npt->nvt->npt->nvt……的过程中,在npt之后,用change_box 命令改变了box的大小,那接下来的nvt就会出错:
ERROR on proc 0: Bond atoms 1457 1479 missing on proc ...

楼主,请问下你解决了这个问题吗?我也遇到了这种问题,也很着急,不知道你是否分享下你的解决办法。非常感谢!!!
选择最喜欢的,专注到无情。一个,就够了。
3楼2016-02-05 03:13:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

superbuff

新虫 (初入文坛)

问题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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fyz_ok

金虫 (正式写手)

引用回帖:
2楼: Originally posted by superbuff at 2015-11-04 12:59:31
问题1中,换一种表述是这样的:如果在npt->nvt->npt->nvt……的过程中,在npt之后,用change_box 命令改变了box的大小,那接下来的nvt就会出错:
ERROR on proc 0: Bond atoms 1457 1479 missing on proc ...

你的问题解决了吗?
4楼2016-09-12 23:11:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 9/450 2026-11-15 02:04 by xSrQNSSjcgoX
[文学芳草园] 厉娜 +4 张丹枫1991 2026-11-09 6/300 2026-11-15 00:48 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 FChKX8meu1ds 2026-11-09 8/400 2026-11-15 00:35 by Equinoxhua
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 G2rhSvqvhNzl 2026-11-10 7/350 2026-11-14 21:58 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-13 11:08 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 Z7qfCsdIiD1T 2026-11-11 5/250 2026-11-12 11:44 by HajghKOMynmz
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:33 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 11:25 by HajghKOMynmz
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 11:14 by HajghKOMynmz
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:11 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:09 by HajghKOMynmz
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:00 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 10:49 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:38 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 09:41 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 12/600 2026-11-12 09:21 by HajghKOMynmz
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 cqv1zKnzd7hr 2026-11-09 8/400 2026-11-12 05:08 by 9fOTqj074tIg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:03 by TgvkbZZl9tqB
信息提示
请填处理意见