24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1954  |  回复: 1

dingfengbo

新虫 (小有名气)

[求助] 温度越高晶格常数却越低,是怎么回事?

我用lammps做SiO2的分子动力学模拟,用的in文件如下:
units        metal
atom_style        atomic
boundary        p p p
variable        x index 0.01 500 1000 1500 2000 2500 3000 3500

print        "----------------temperature=$x K---------"

read_data      data.SiO2

pair_style        tersoff
pair_coeff        * * SiO.tersoff O Si

neighbor           0.5 bin
neigh_modify  every 4 delay 0 check yes

group              SiO2 type 1 2         

timestep            0.0001
thermo            1000

velocity        SiO2 create $x 825577 dist gaussian

fix             3 SiO2 npt temp $x $x 0.1 iso 0.0 0.0 500 drag 0.2
run              1200000
unfix                3
fix             3 SiO2 nvt temp $x $x 0.1 drag 0.2
dump            2 all atom 1000 OSi_$x.lammpstrj
run                10000
clear
next                x
jump                in.temp
用的data文件开头如下:
192 atoms

   2 atom types

     0.000000000    14.32 xlo xhi
     0.000000000    14.32 ylo yhi
     0.000000000    14.32 zlo zhi

Masses

    1 15.9989995956
    2 28.0849990845

Atoms

1       1       0.67875 0.67875 0.67875
2       1       0.67875 3.39375 3.39375
3       1       3.39375 0.67875 3.39375
4       1       3.39375 3.39375 0.67875
结果在0.01K下得到的xlo与xhi间的距离是17.8,在500、1000、1500、2000、2500、3000K下得到的差不多都是14.32,随着温度升高晶格常数反而变小了,或者说从500到3000K都没怎么变化,这是哪里出问题了呢?不都应该随着温度升高而变大的么?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huai

木虫 (职业作家)

★ ★ ★
ghcacj(金币+3): 谢谢 2011-11-25 10:37:58
我利用一种势计算钨的晶格常数。发现它的晶格常数是随着温度的升高而减小。也就是说热膨胀系数是负的。这仅仅是这种势描述的时候导致的。换另外一种势就对了。势函数不是万能的,它有可能算某些性质是对的,而另外一些性质是错的或不好。希望对你排查错误有帮助。
2楼2011-11-24 18:45:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dingfengbo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 人工智能320调剂08工类还有机会吗 +12 振—TZ 2026-04-10 12/600 2026-04-12 07:54 by wj165256
[考研] 299求调剂 +8 ZVVZ13 2026-04-08 8/400 2026-04-12 00:40 by 蓝云思雨
[考研] 296求调剂 +8 汪!?! 2026-04-09 8/400 2026-04-11 21:02 by 逆水乘风
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 270求调剂 +14 杨乐369 2026-04-11 14/700 2026-04-11 20:16 by 蓝云思雨
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[考研] 吉大计算机技术331分,英语六级,求调剂 +3 峰峰021116 2026-04-09 3/150 2026-04-10 20:01 by chemisry
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 273求调剂 +51 麦小叮当 2026-04-06 58/2900 2026-04-10 15:54 by jiajinhpu
[基金申请] 有爆料,一个青年教师卖房得400万,然后换了一个四青帽子 +9 babu2015 2026-04-08 9/450 2026-04-10 11:43 by 苏东坡二世
[考研] 材料复试求调剂 +20 xhhdjdjsjks 2026-04-09 20/1000 2026-04-10 10:25 by 孙小小12457
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 297求调剂 +27 GENJIOW 2026-04-07 30/1500 2026-04-09 23:20 by wolf97
[考研] 初试分332,一志愿报考西北工业大学, +11 故人?? 2026-04-09 11/550 2026-04-09 21:54 by JineShine
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH(等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 本科郑州大学,一志愿华东师范大学282求调剂 +23 熊哥xtk 2026-04-07 26/1300 2026-04-09 17:17 by 18446523
[考研] 085501机械英二77总分294求调剂,接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
信息提示
请填处理意见