24小时热门版块排行榜    

查看: 483  |  回复: 4
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cxz314

木虫 (正式写手)

[交流] 【求助】NVT系综下液态金属

小弟用MC模拟NVT系综下液态金属(以金属铝为例)的性质,不知道液态金属的初始密度应该怎么设置?是不是可以随意的设置,比如就设置成固态金属的密度,还是要设成其他的密度。
请各位指教!

[ Last edited by zdhlover on 2009-11-13 at 02:56 ]
回复此楼

» 猜你喜欢

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

cxz314

木虫 (正式写手)

引用回帖:
Originally posted by 老虎大王 at 2009-4-2 11:18:
如果不涉及相变,可以用晶体密度,忽略固液体积差别。

也可以先用NPT充分弛豫,平衡之后得到一个平均体积(密度),将之用于NVT模拟。

也可以使用实验的液体密度。Al这样的常见金属,是有液态密度数据的。
...

非常感谢!
那我就用第二种方法试一下。另外如果方便的话能不能把常见金属的液态密度告诉我,我找了一下没能找到(哈哈哈,要求有点过分了)
再次表示感谢!
3楼2009-04-02 15:00:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

老虎大王

木虫 (著名写手)

★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢 4-2 11:54
如果不涉及相变,可以用晶体密度,忽略固液体积差别。

也可以先用NPT充分弛豫,平衡之后得到一个平均体积(密度),将之用于NVT模拟。

也可以使用实验的液体密度。Al这样的常见金属,是有液态密度数据的。

以上三种方法,以第二种较为常用,因为它提供的是与你的势函数有关的密度参数。虽然这样得到的数据通常比实验值为小(体积偏大),但它保证了计算的自洽性,并使得计算中的原子应力(或体系的应力张量)降低。不过,除非你专门研究密度或体积问题,一般不会有人认真关心模拟中的密度参数,只要它与实验值的偏离不是很离谱。
2楼2009-04-02 11:18:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★ ★ ★
yuhuobuku(金币+5,VIP+0):欢迎参加讨论 4-8 09:25
可以去图书馆找CRC handbook of chemiistry and physics。

对于液态金属密度,有关联公式: rho=rhom-k(t-tm)
tm是熔点(℃),对Al为660.32
rhom是熔点时密度(g/cm^-3 ), 对Al 为2.375
k是参数(g/cm^-3 ℃^-1 ),对Al为0.000233
此式可用的最大温度是1340℃以内。

参考文献:D.R.Lide, CRC handbook of chemistry and physics, 84th ed., CRC press, Boca Raton. 2004

这个手册每一两年就出一本,如果找不到这一本,找到其它年份的也可以,内容不会有太大变化。
4楼2009-04-02 16:37:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cxz314

木虫 (正式写手)

引用回帖:
Originally posted by 老虎大王 at 2009-4-2 16:37:
可以去图书馆找CRC handbook of chemiistry and physics。

对于液态金属密度,有关联公式: rho=rhom-k(t-tm)
tm是熔点(℃),对Al为660.32
rhom是熔点时密度(g/cm^-3 ), 对Al 为2.375
k是参数(g/cm^-3 ...

太感谢你了!
5楼2009-04-02 16:57:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见