24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3641  |  回复: 24
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

维特的日记

木虫 (小有名气)

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

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

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

» 收录本帖的淘帖专辑推荐

熔化理论

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

维特的日记

木虫 (小有名气)

引用回帖:
8楼: Originally posted by abinitio at 2015-02-17 11:33:42
你说得对,我讲错了,是盒子大小。
比如,一个2D的,盒子大小是8,如果一个原子原始坐标是(8,8),由于运动和周期性边界条件,它的坐标可能会变为(1,1),如果按它xy各运动了7个单位,肯定是错误的,因为它在x和y ...

嗯嗯 明白了 就是说如果两个原子之间的距离大于盒子的一半 肯定会存在一个更近的原子 所以应该考虑更近的原子与它之间的距离
9楼2015-02-17 11:41:51
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +11 Xxjc1107. 2026-04-13 11/550 2026-04-13 15:28 by chw1980_0
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 53/2650 2026-04-13 14:20 by 张zhihao
[教师之家] 转长聘了 +3 简单化xn 2026-04-13 3/150 2026-04-13 14:18 by wwwkkk83
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 280求调剂 +7 兮兮夜夜 2026-04-09 10/500 2026-04-12 00:33 by 蓝云思雨
[考研] 材料工程日语考生求调剂 +7 0856?调剂 2026-04-10 7/350 2026-04-11 21:33 by 蓝云思雨
[考研] 352 求调剂 +6 yzion 2026-04-11 8/400 2026-04-11 16:24 by 明月此时有
[考研] 化学308分求调剂 +22 你好明天你好 2026-04-07 24/1200 2026-04-11 11:14 by ChemPharm
[考研] 农学0904 312求调剂 +6 Say Never 2026-04-10 6/300 2026-04-11 10:33 by wwj2530616
[考研] 283求调剂,工科! +12 苏打水7777 2026-04-08 12/600 2026-04-11 10:28 by 逆水乘风
[考研] 346,工科0854求调剂,专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 一志愿华南理工大学331分材料求调剂 +9 天下ww 2026-04-09 9/450 2026-04-10 22:58 by Ftglcn90
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 一志愿中科院105500专业总分315求调剂 +6 lallalh 2026-04-09 7/350 2026-04-09 17:51 by lallalh
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 085501机械英二77总分294求调剂,接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 349学科化学045106求调剂,化学类都可以 +8 保好懂懂 2026-04-08 8/400 2026-04-09 14:03 by xulei3024
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
信息提示
请填处理意见