24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1952  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 366求调剂 +8 不知名的小卅 2026-04-11 8/400 2026-04-12 19:49 by zhouxiaoyu
[考研] 286求调剂 +26 Faune 2026-04-06 26/1300 2026-04-12 10:24 by lhj2009
[考研] 297求调剂 +17 ORCHID1 2026-04-10 18/900 2026-04-12 09:59 by zhouwenxian
[考研] 调剂 +10 月@163.com 2026-04-11 10/500 2026-04-12 09:14 by zhouyuwinner
[考研] 材料与化工300求调剂 +39 肖开文 2026-04-09 43/2150 2026-04-12 01:30 by 秋豆菜芽
[考研] 288求调剂,一志愿华南理工大学071005 +18 ioodiiij 2026-04-08 18/900 2026-04-11 20:25 by liyun12321
[考研] 求调剂 +10 璃茉一定上岸 2026-04-10 10/500 2026-04-11 13:31 by 1005715100
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[材料工程] 材料调剂推荐 +8 蛋糕x2 2026-04-07 8/400 2026-04-10 23:13 by Ftglcn90
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +10 努力奋斗112 2026-04-07 10/500 2026-04-08 15:01 by screening
[考研] 求调剂 +15 熊二想上岸 2026-04-06 15/750 2026-04-08 04:53 by 无际的草原
[考研] 机械工程264学硕求调剂 +3 qiushangxian 2026-04-06 3/150 2026-04-08 01:53 by Linzejun
[考研] 求考研材料调剂 +3 材化李可 2026-04-07 3/150 2026-04-08 00:21 by JourneyLucky
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
[考研] 327考研调剂推荐 +6 呜呜呜呜呢 2026-04-06 6/300 2026-04-06 21:39 by 啵啵啵0119
信息提示
请填处理意见