| 查看: 2042 | 回复: 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都没怎么变化，这是哪里出问题了呢？不都应该随着温度升高而变大的么？ |
回复此楼» 猜你喜欢
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有9人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
国自然面上复盘~欢迎讨论
已经有19人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 为什么烧结温度越高，密度越低，而微波介电性能却先升后降
已经有21人回复
- 酶类是不是保存温度越低越好
已经有8人回复
- 为什么温度越高，三氧化二锰越多？
已经有18人回复
- 【求助】柱压越冲越高怎么回事啊
已经有22人回复
- 【求助】gulp中改变温度，晶格常数没变？
已经有27人回复
- 【求助】H2-TPR,H2浓度一般多少，为什么说H2浓度越低，还原温度越高？
已经有6人回复
- 【讨论】环境温度越低 开路电压应该越低还是越高...？
已经有10人回复
- 【求助】怎么计算镁合金的晶格常数？
已经有7人回复
- 【求助】知道空间群和晶格常数，如何知道原子坐标！！
已经有6人回复
- 【讨论】重结晶温度不一定是越低越好吧？
已经有6人回复
huai
木虫 (职业作家)
- 应助: 4 (幼儿园)
- 金币: 1159.7
- 散金: 118
- 红花: 6
- 帖子: 4542
- 在线: 554小时
- 虫号: 420388
- 注册: 2007-07-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2011-11-24 18:45:30
5