24小时热门版块排行榜    

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

维特的日记

木虫 (小有名气)

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

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

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注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 ...

后来我编程的时候发现那个公式里的sigma求和就是指的系综平均吧 不然我们做系综平均后那个sigma无法求和呀 因为在系综平均的时候已经求和过一次了

发自小木虫IOS客户端
20楼2015-09-12 22:58:25
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856材料求调剂 +12 hyf hyf hyf 2026-02-28 13/650 2026-03-02 20:19 by hypershenger
[考研] 材料270求调剂 6+4 Eiiiio 2026-03-01 7/350 2026-03-02 20:17 by hypershenger
[考研] 0856求调剂285 +11 吕仔龙 2026-02-28 11/550 2026-03-02 20:15 by hypershenger
[考研] 材料284求调剂,一志愿郑州大学英一数二专硕 +15 想上岸的土拨鼠 2026-02-28 15/750 2026-03-02 20:13 by hypershenger
[考研] 化学,材料,环境类求调剂 +7 考研版棒棒 2026-03-02 7/350 2026-03-02 19:56 by hypershenger
[考研] 0805总分292,求调剂 +10 幻想之殇 2026-03-01 10/500 2026-03-02 19:51 by zhukairuo
[考研] 求调剂 +10 yunziaaaaa 2026-03-01 11/550 2026-03-02 19:17 by caszguilin
[考研] 一志愿中石油(华东)本科齐鲁工业大学 +3 石能伟 2026-03-02 3/150 2026-03-02 18:54 by caszguilin
[考研] 化工京区271求调剂 +6 11ing 2026-03-02 6/300 2026-03-02 18:52 by caszguilin
[考研] 材料调剂 +3 恒顺自然 2026-03-02 3/150 2026-03-02 18:49 by L135790
[考研] 一志愿山东大学材料与化工325求调剂 +5 半截的诗0927 2026-03-02 5/250 2026-03-02 18:37 by 明亮9527
[考研] 282求调剂 +4 2103240126 2026-03-02 6/300 2026-03-02 18:07 by 2103240126
[考研] 化工专硕348,一志愿985求调剂 +6 弗格个 2026-02-28 9/450 2026-03-02 14:09 by liyongv
[考研] 材料与化工328求调剂 +3 。,。,。,。i 2026-03-02 3/150 2026-03-02 13:09 by houyaoxu
[考研] 26考研报考西工大材料308分求调剂 +4 weizhong123 2026-03-01 4/200 2026-03-02 12:46 by 无际的草原
[考研] 材料学调剂 +10 提神豆沙包 2026-02-28 12/600 2026-03-02 09:26 by 李老师!
[考研] 材料工程269求调剂 +3 白刺玫 2026-03-02 3/150 2026-03-02 09:25 by 一休哥FU
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 290求调剂 +9 材料专硕调剂; 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
信息提示
请填处理意见