24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1993  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +11 1234567wang 2026-05-17 13/650 2026-05-20 18:55 by tanfancy
[基金申请] 评审有感 +14 popular289 2026-05-18 23/1150 2026-05-20 17:39 by lincyb
[基金申请] 提交了我也来说说感想 +6 fummck 2026-05-20 6/300 2026-05-20 17:29 by 姜太公烤鱼
[基金申请] 国自然评分 +3 无名者登山 2026-05-20 4/200 2026-05-20 13:29 by vito刘
[考博] 如果工作了想读博,可以边工作边读全日制嘛? 30+3 铁达火车 2026-05-18 5/250 2026-05-20 09:33 by tfang
[基金申请] 河北省自然科学基金 +3 Peterchao 2026-05-18 3/150 2026-05-20 08:57 by 霸_霸
[教师之家] 上海大学实验技术岗位非升即走 +12 嘻嘻哈哈乐呵呵 2026-05-15 13/650 2026-05-20 08:34 by xli1984
[考博] 云南大学材料与能源学院解琳课题组钙钛矿博士招生 +3 光伏爱好者 2026-05-17 5/250 2026-05-19 19:13 by 光伏爱好者
[基金申请] 2026山东省优青 +3 luckinging 2026-05-14 4/200 2026-05-19 16:15 by 52wls
[考博] 找博士生导师 +6 小代想上岸 2026-05-15 7/350 2026-05-19 10:22 by free_fisher
[考博] 26/27申博自荐-锂/钠电池方向 5+3 狗头军师. 2026-05-15 4/200 2026-05-19 09:10 by moonboat
[基金申请] 同样的基金本子,换个专家直接从C变A! (金币-10) +3 国自然国社科中 2026-05-19 3/150 2026-05-19 08:50 by Equinoxhua
[基金申请] 别被青基扩招骗了!26年科研内卷才刚刚开始 +3 国自然国社科中 2026-05-14 4/200 2026-05-19 08:48 by archvillain
[考博] 2026博士还有哪些学校有名额 +7 小王求读研 2026-05-15 8/400 2026-05-19 08:27 by zhyzzh
[考博] 博士申请 +5 星…… 2026-05-18 6/300 2026-05-18 23:49 by 糊糊涂涂好
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[硕博家园] 我在等一个没有答案的答案 +3 Love_MH 2026-05-17 3/150 2026-05-18 02:22 by 竹林孤影
[文学芳草园] 半夜喝咖啡 +3 myrtle 2026-05-15 5/250 2026-05-18 01:03 by 小沈2018
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
信息提示
请填处理意见