24小时热门版块排行榜    

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

维特的日记

木虫 (小有名气)

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

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

LAMMPS里模拟纳米团簇熔化时关于Lindemann指数的计算问题
0.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fuga

铜虫 (初入文坛)

这跟lammps里面自带的msd(mean squared distance)计算有什么区别吗? 我怎么觉得不需要自己编程
18楼2015-08-20 23:39:26
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-26 21:15 by xxniao123
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 基金不中,共勉 +3 eulota 2026-08-26 3/150 2026-08-26 16:17 by 路yyyy
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见