| 查看: 523 | 回复: 4 | ||
| 当前主题已经存档。 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
cxz314木虫 (正式写手)
|
[交流]
【求助】NVT系综下液态金属
|
|
|
小弟用MC模拟NVT系综下液态金属（以金属铝为例）的性质，不知道液态金属的初始密度应该怎么设置？是不是可以随意的设置，比如就设置成固态金属的密度，还是要设成其他的密度。 请各位指教！ [ Last edited by zdhlover on 2009-11-13 at 02:56 ] |
回复此楼
cxz314
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 贵宾: 0.2
- 金币: 3815.1
- 散金: 993
- 红花: 1
- 帖子: 997
- 在线: 189小时
- 虫号: 475745
- 注册: 2007-12-07
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2009-04-02 16:57:03
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢 4-2 11:54
lei0736(金币+4,VIP+0):谢谢 4-2 11:54
|
如果不涉及相变，可以用晶体密度，忽略固液体积差别。 也可以先用NPT充分弛豫，平衡之后得到一个平均体积（密度），将之用于NVT模拟。 也可以使用实验的液体密度。Al这样的常见金属，是有液态密度数据的。 以上三种方法，以第二种较为常用，因为它提供的是与你的势函数有关的密度参数。虽然这样得到的数据通常比实验值为小（体积偏大），但它保证了计算的自洽性，并使得计算中的原子应力（或体系的应力张量）降低。不过，除非你专门研究密度或体积问题，一般不会有人认真关心模拟中的密度参数，只要它与实验值的偏离不是很离谱。 |
2楼2009-04-02 11:18:15
cxz314
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 贵宾: 0.2
- 金币: 3815.1
- 散金: 993
- 红花: 1
- 帖子: 997
- 在线: 189小时
- 虫号: 475745
- 注册: 2007-12-07
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2009-04-02 15:00:51
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★
yuhuobuku(金币+5,VIP+0):欢迎参加讨论 4-8 09:25
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
