24小时热门版块排行榜    

查看: 1939  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 081700 调剂 267分 +3 迷人的哈哈 2026-03-23 3/150 2026-03-23 11:16 by barlinike
[考研] 336化工调剂 +3 王大坦1 2026-03-23 4/200 2026-03-23 11:07 by laoshidan
[考研] 北科281学硕材料求调剂 +8 tcxiaoxx 2026-03-20 8/400 2026-03-23 10:16 by 冠c哥
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 求调剂 +6 十三加油 2026-03-21 6/300 2026-03-22 17:00 by i_cooler
[考研] 材料与化工(0856)304求 B区 调剂 +3 邱gl 2026-03-21 3/150 2026-03-21 13:47 by lature00
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 332求调剂 +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +4 solanXXX 2026-03-20 4/200 2026-03-20 23:49 by alg094825
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 广西大学家禽遗传育种课题组2026年硕士招生(接收计算机专业调剂) +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
信息提示
请填处理意见