应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » LAMMPS里模拟纳米团簇熔化时关于Lindemann指数的计算问题
52/1返回列表
查看: 3561  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

维特的日记

木虫 (小有名气)

[求助] LAMMPS里模拟纳米团簇熔化时关于Lindemann指数的计算问题 已有1人参与

我在做一个金属纳米团簇熔化的模拟，由于需要Lidemann指数来判断熔点，lammps里没有现成的命令来计算把？感觉是不是要自己编程来算。。。
我对<Rij^2>和<Rij>^2不是很理解，在计算的时候怎么区别他们呢？怎样是先取了系综平均<>，我不知道用程序怎么实现他，希望研究过的各位前辈指点一下，谢谢哈！
提前祝大家新年快乐！

LAMMPS里模拟纳米团簇熔化时关于Lindemann指数的计算问题
0.png
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2015-02-16 12:38:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)
引用回帖:
15楼: Originally posted by xjniudong at 2015-04-07 22:59:48
对于符号"<>"系纵平均来说，不是要经历一个对时间的平均么？我对这个系纵平均的概念不是太懂，看版主的意思是通过一个特定时刻下各个粒子的坐标信息求得林得曼参数，不用考虑时间平均么？？多谢 ...

我的理解，对于一个md过程，除非是在特别的研究驰誉过程，否则我们都会以体系稳定后才开始研究，如何体系稳定，说明体系对时间上的变化已趋于固定值 ，因此，不再考虑时间上的平均(当然了，也可以取几个时间进行平均一下，但我认为差别不大)。 这是我的理解，不知正确与否。不敢赐教，共同讨论。
赞一下 回复此楼
16楼2015-04-08 07:52:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答

abinitio

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
维特的日记: 金币+10, ★★★很有帮助 2015-02-16 18:41:56
ben_ladeng: 金币+2, 谢谢指教 2015-02-16 19:33:20
维特的日记: 金币+20, ★★★★★最佳答案 2015-02-17 18:09:36
<Rij^2>是两点距离平方的平均值，<Rij>^2是两点距离平均的平方值，系综平均<>就是将所有任何两个原子之间遍历后求总和再除原子数，
编程思路是：
如果有N个原子，先做i=1：N的循环，，再j=1：N-1
计算i原子为中心与原子j 的距离Rij=，此时要注意一下周期性边界条件
把所有的Rij累加，把所有Rij平方累加，然后除以总原子数
赞一下 回复此楼
2楼2015-02-16 15:09:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

维特的日记

木虫 (小有名气)
引用回帖:
2楼: Originally posted by abinitio at 2015-02-16 15:09:36
<Rij^2>是两点距离平方的平均值，<Rij>^2是两点距离平均的平方值，系综平均<>就是将所有任何两个原子之间遍历后求总和再除原子数，
编程思路是：
如果有N个原子，先做i=1：N的循环，，再j=1：N ...

谢谢你的回答！大体已经明白了，但是那个周期性边界要怎么处理呢，这个地方没怎么想明白
赞一下 回复此楼
3楼2015-02-16 18:38:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 维特的日记 at 2015-02-16 18:38:50
谢谢你的回答！大体已经明白了，但是那个周期性边界要怎么处理呢，这个地方没怎么想明白...

如果在某一方向的位移值绝对值超过了晶格参数的一半，就将该方向的位移减或加上一个晶格参数，你想一下周期性边界条件的定义，就会理解了。
赞一下 回复此楼
4楼2015-02-17 09:00:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 25 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 272求调剂 +9 材紫有化 2026-02-28 9/450 2026-03-02 20:22 by hypershenger
[考研] 290分材料工程085601求调剂 数二英一 +5 llx0610 2026-03-02 5/250 2026-03-02 20:20 by hypershenger
[考研] 调剂 +3 13853210211 2026-03-02 3/150 2026-03-02 19:59 by hypershenger
[考研] 化学，材料，环境类求调剂 +7 考研版棒棒 2026-03-02 7/350 2026-03-02 19:56 by hypershenger
[考研] 课题组接收材料类调剂研究生 +4 gaoxiaoniuma 2026-02-28 5/250 2026-03-02 19:33 by Xinyuu
[考研] 275求调剂 +7 明远求学 2026-03-01 7/350 2026-03-02 19:22 by zhukairuo
[考研] 289求调剂 +8 yang婷 2026-03-02 9/450 2026-03-02 19:08 by zhukairuo
[考研] 268求调剂 +3 好运连绵不绝 2026-03-02 3/150 2026-03-02 19:05 by zhukairuo
[考研] 一志愿中石油（华东）本科齐鲁工业大学 +3 石能伟 2026-03-02 3/150 2026-03-02 18:54 by caszguilin
[考研] 材料调剂 +3 恒顺自然 2026-03-02 3/150 2026-03-02 18:49 by L135790
[考研] 303求调剂 +5 今夏不夏 2026-03-01 5/250 2026-03-02 15:01 by 向上的胖东
[考研] 化工专硕348，一志愿985求调剂 +6 弗格个 2026-02-28 9/450 2026-03-02 14:09 by liyongv
[考研] 338求调剂 +3 18162027187 2026-03-02 3/150 2026-03-02 13:12 by houyaoxu
[基金申请] 此成果不能导入原因：元数据必填信息不完整，可 进行补充。 +4 Kittylucky 2026-03-02 5/250 2026-03-02 11:07 by jurkat.1640
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 279求调剂 +3 dua1 2026-03-01 4/200 2026-03-02 00:23 by 大脸蛋子
[考研] 299求调剂 +3 Y墨明棋妙Y 2026-02-28 5/250 2026-03-01 21:01 by tangxiaotian
[考研] 求调剂 +6 repeatt?t 2026-02-28 6/300 2026-03-01 14:37 by Sakura绘
[考研] 317一志愿华南理工电气工程求调剂 +6 Soliloquy_Q 2026-02-28 11/550 2026-03-01 11:14 by 歌liekkas
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
信息提示
关闭
请填处理意见
关闭