24小时热门版块排行榜    

查看: 2792  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 言之凿凿,8月21号(明天)放榜 +13 医学老男孩 2026-08-20 15/750 2026-08-20 15:19 by chengyan1220
[基金申请] 时间戳变了,能看出什么问题? +17 基诺咪客 2026-08-17 21/1050 2026-08-20 14:54 by 昝占全2025
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 21/1050 2026-08-20 10:37 by Vivilian
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见