24小时热门版块排行榜    

查看: 4031  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cxz314

木虫 (正式写手)

[交流] 【讨论】lammps 用sw势计算硅的熔点怎么这么高(2400k),附in文件 已有7人参与

lammps 用sw势计算硅的熔点怎么这么高(2400k)?请各位指教
附in文件

# density=2.3292 pho=0.049941614

# bulk Sillicon lattice



units         metal

boundary         p p p

atom_style         atomic



lattice          diamond   5.43095

region                 box block 0 10.0 0 10.0 0 10.0

create_box         1 box

create_atoms         1 box



pair_style         sw

pair_coeff         * * Si.sw Si

mass        1 28.0855

timestep         0.002
thermo                 1000



neighbor         0.5 bin

neigh_modify         every 1 delay 5 check yes



velocity         all create 100.0 4928459 rot yes dist gaussian

fix                   1 all nvt 100.0 100.0 0.1 drag 0.7



thermo_style        custom step temp press vol etotal
log                      log100.sw



run                    50000


unfix       1


fix                2 all npt 100.0 2800.0 0.1 xyz 0.0 0.0 0.1  drag 0.7

thermo_style        custom step temp press vol etotal


log                log100_2800.sw



run      500000

unfix     2

fix                3 all npt 2800.0 2800.0 0.1 xyz 0.0 0.0 0.1  drag 0.7

thermo_style        custom step temp press vol etotal

log                log2800.sw

run      50000

dump                1 all custom 600 dump_2800.sw id type x y z xu yu zu
run      12000
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛排小黑娇

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
12楼: Originally posted by xy515258 at 2014-04-17 19:27:04
你是想直接从100K升温到2800K然后看能量突变的点?   那当然高得飞起啦。加热速度那么高   当然会出现过热很多才熔化啊       2400K根本就不是你测出来的熔点   只是过热度而已吧       建议有两个方法   一个是你跑 ...

亲 有木有这个两相法的程序呢 我也遇到这个问题
没走过的,是路;走过的,才是人生。
14楼2014-12-23 16:42:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cxz314 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见