24小时热门版块排行榜    

查看: 2501  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境工程调剂 +6 大可digkids 2026-03-16 6/300 2026-03-16 17:16 by barlinike
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] 材料工程,326分,求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 285化工学硕求调剂(081700) +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
信息提示
请填处理意见