24小时热门版块排行榜    

查看: 2881  |  回复: 10

gulubaozi

银虫 (小有名气)

[求助] LAMMPS跟踪原子 已有2人参与

LAMMPS手册上介绍说,原子的ID不变,但是就输出的文件而言,同一个序号的原子,真的是同一个原子吗?
如果不是如何实现跟踪原子。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

gulubaozi

银虫 (小有名气)

★
abinitio: 金币+1, 鼓励交流 2014-11-12 09:56:38
引用回帖:
2楼: Originally posted by abinitio at 2014-11-11 08:50:18
这一点没有疑问吧,肯定就是同一原子。原子的ID不确定,如果会通过势的梯度来算加速度呢?楼主可以自己测试一下,建一个固体结构,在熔点之下驰誉,你会发现,边界处的原子由于周期性条件会转移到相应位置外,其它的 ...

我刚试了一下,对于同一种格式的输出数据文件,ID是不变的。
但是不同格式文件中同一个原子的ID不一样。
因为电脑很渣,所以只算了2步。so 这个实验结果可信度不高。
期待更准确更详细的说明。
4楼2014-11-11 09:38:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

abinitio

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gulubaozi: 金币+10, ★★★很有帮助 2014-11-11 09:35:27
这一点没有疑问吧,肯定就是同一原子。原子的ID不确定,如何 会通过势的梯度来算加速度呢?楼主可以自己测试一下,建一个固体结构,在熔点之下驰誉,你会发现,边界处的原子由于周期性条件会转移到相应位置外,其它的基本在原地振动
2楼2014-11-11 08:50:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gulubaozi

银虫 (小有名气)

引用回帖:
2楼: Originally posted by abinitio at 2014-11-11 08:50:18
这一点没有疑问吧,肯定就是同一原子。原子的ID不确定,如果会通过势的梯度来算加速度呢?楼主可以自己测试一下,建一个固体结构,在熔点之下驰誉,你会发现,边界处的原子由于周期性条件会转移到相应位置外,其它的 ...

我在用不同格式输出时比如cfg,xyz,甚至把compute得到的结果输出,对于每一个输出文件,同一个序号的参数组,都是同一个原子的性质参数,你是这个意思吗?
3楼2014-11-11 09:04:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huai

木虫 (职业作家)

★ ★ ★ ★ ★
gulubaozi(abinitio代发): 金币+5, 感谢互助。应助+1 2014-11-12 17:13:16
原子的ID是不变的。但是两个不时刻构型的文件,原子的排列顺序不一样。所以最好用lammps里面的排序命令,输出的时候,每一次都是1,2,3,。。。,N。这样算起来方便一些,也不易出错。
5楼2014-11-12 12:56:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gulubaozi

银虫 (小有名气)

引用回帖:
5楼: Originally posted by huai at 2014-11-12 12:56:31
原子的ID是不变的。但是两个不时刻构型的文件,原子的排列顺序不一样。所以最好用lammps里面的排序命令,输出的时候,每一次都是1,2,3,。。。,N。这样算起来方便一些,也不易出错。

每次输出前使用一下LAMMPS的排序命令,是这样吗?
6楼2014-11-12 16:19:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huai

木虫 (职业作家)

引用回帖:
6楼: Originally posted by gulubaozi at 2014-11-12 16:19:02
每次输出前使用一下LAMMPS的排序命令,是这样吗?...

是的,可以把原子的ID输出来一看就明白了。
7楼2014-11-12 16:33:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gulubaozi

银虫 (小有名气)

abinitio: 我代劳了。:) 2014-11-12 17:13:43
引用回帖:
7楼: Originally posted by huai at 2014-11-12 16:33:15
是的,可以把原子的ID输出来一看就明白了。...

非常感谢!
为什么不选择应助呢?楼主很愿意给你金币。
8楼2014-11-12 16:39:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huai

木虫 (职业作家)

引用回帖:
8楼: Originally posted by gulubaozi at 2014-11-12 16:39:32
非常感谢!
为什么不选择应助呢?楼主很愿意给你金币。...

呵呵,有点帮助就行。
9楼2014-11-12 17:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
这似乎与你的输出格式有关,我刚看了一下我输出的CFG格式,原子ID并不是固定的,但我以前输出的CUSTOM我好像有印象是固定的
10楼2014-11-13 21:51:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gulubaozi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:45 by Z9YWQ5EAO3qp
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:26 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
请填处理意见