24小时热门版块排行榜

>论坛更新日志 (3681)
>文献求助 (444)
>虫友互识 (322)
>导师招生 (315)
>考博 (174)
>论文投稿 (129)
>基金申请 (111)
>硕博家园 (90)
>博后之家 (80)
>休闲灌水 (80)
>招聘信息布告栏 (79)
>公派出国 (61)
>教师之家 (60)
>找工作 (58)
>考研 (48)
>论文道贺祈福 (36)

返回列表

dingfengbo

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 79
帖子: 70
在线: 38.4小时
虫号: 939022
注册: 2010-01-08
专业: 凝聚态物性 II ：电子结构

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

我用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都没怎么变化，这是哪里出问题了呢？不都应该随着温度升高而变大的么？

回复此楼