24小时热门版块排行榜    

查看: 1916  |  回复: 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 的主题更新
信息提示
请填处理意见