| 查看: 2229 | 回复: 2 | ||
[求助]
计算W-Mo合金的熔点，程序lost atoms，请教能人解答 已有1人参与
|
|
参照上海交大那个求Cu或Al熔点的例子，把lattice部分换成了我Mo-0.25W成分合金的原子分布，in文件如下： ++++++++++++++++++++++++++++++++++++++++++++++++++++++ in.melt_temp ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ # LAMMPS Melt_Mo-W_temp units metal boundary p p p atom_style atomic variable x index 10 500 1000 1800 2000 variable y equal 3.1466/15.733 variable z universe 15.733 lattice custom $z a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 $y basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 $y basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.2 basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 region box block -4.0 4.0 -3.0 3.0 0.1 45.0 create_box 2 box region Mo-0.25W block -4.0 4.0 -3.0 3.0 0.1 45.0 create_atoms 2 region Mo-0.25W basis 1 1 basis 2 1 basis 3 1 basis 4 1 basis 5 1 basis 6 1 basis 7 1 basis 8 1 basis 9 1 basis 10 1 basis 11 1 basis 12 1 basis 13 1 basis 14 1 basis 15 1 basis 16 2 basis 17 2 basis 18 2 basis 19 2 basis 20 2 group Mo-0.25W region Mo-0.25W region Mo-0.25W11 block -4.0 4.0 -3.0 3.0 43.0 45.0 region Mo-0.25W22 block -4.0 4.0 -3.0 3.0 0.1 3.0 group Mo-0.25W11 region Mo-0.25W11 group Mo-0.25W22 region Mo-0.25W22 fix region Mo-0.25W11 setforce 0 0 0 fix region Mo-0.25W22 setforce 0 0 0 mass 1 95.941 mass 2 183.841 #fix wall all wall/reflect zlo zhi timestep 0.01 thermo 100 pair_style hybrid morse 10.0 pair_coeff 1 1 morse 0.77019 1.434 3.012 pair_coeff 1 2 morse 13.68932 1.82699 2.30154 pair_coeff 2 2 morse 0.969483 1.385 3.053 neighbor 0.6 bin neigh_modify every 5 delay 0 check yes velocity all create $x 825577 dist gaussian fix 1 all npt $x $x 2.0 xyz 0.0 0.0 45.0 drag 0.2 #在x温度下保持外压为零弛豫 compute 3 all pe/atom compute 4 all ke/atom compute 5 all coord/atom 3.0 run 10000 unfix 1 fix 1 all nvt $x $x 2.0 drag 0.2 #NVT控温 dump 1 all custom 50 dump_$x.atom id xs ys zs c_3 c_4 c_5 fix 2 all msd 1 msd_Mo-W_$x.dat #输出MSD文件 run 10000 clear next x jump in.melt_Mo-W_temp +++++++++++++++++++++++++++++++++++++++++++++++++++= 以下是log文件 +++++++++++++++++++++++++++++++++++++++++++++++++++++ LAMMPS (7 Jul 2009) jump in.melt_temp # LAMMPS Melt_Mo-W_temp units metal boundary p p p atom_style atomic variable x index 10 500 1000 1800 2000 variable y equal 3.1466/15.733 variable z universe 15.733 Initial ${z} setting: value 1 on partition 0 lattice custom $z a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 $y basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 $y basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.2 basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 Lattice spacing in x,y,z = 6.2932 15.733 3.1466 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 $y basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.2 basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 Lattice spacing in x,y,z = 6.2932 15.733 3.1466 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.2 basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 Lattice spacing in x,y,z = 6.2932 15.733 3.1466 region box block -4.0 4.0 -3.0 3.0 0.1 45.0 create_box 2 box Created orthogonal box = (-25.1728 -47.199 0.31466) to (25.1728 47.199 141.597) 1 by 1 by 1 processor grid region Mo-0.25W block -4.0 4.0 -3.0 3.0 0.1 45.0 create_atoms 2 region Mo-0.25W basis 1 1 basis 2 1 basis 3 1 basis 4 1 basis 5 1 basis 6 1 basis 7 1 basis 8 1 basis 9 1 basis 10 1 basis 11 1 basis 12 1 basis 13 1 basis 14 1 basis 15 1 basis 16 2 basis 17 2 basis 18 2 basis 19 2 basis 20 2 Created 42720 atoms group Mo-0.25W region Mo-0.25W 42720 atoms in group Mo-0.25W region Mo-0.25W11 block -4.0 4.0 -3.0 3.0 43.0 45.0 region Mo-0.25W22 block -4.0 4.0 -3.0 3.0 0.1 3.0 group Mo-0.25W11 region Mo-0.25W11 1920 atoms in group Mo-0.25W11 group Mo-0.25W22 region Mo-0.25W22 2880 atoms in group Mo-0.25W22 fix region Mo-0.25W11 setforce 0 0 0 fix region Mo-0.25W22 setforce 0 0 0 mass 1 95.941 mass 2 183.841 #fix wall all wall/reflect zlo zhi timestep 0.01 thermo 100 pair_style hybrid morse 10.0 pair_coeff 1 1 morse 0.77019 1.434 3.012 pair_coeff 1 2 morse 13.68932 1.82699 2.30154 pair_coeff 2 2 morse 0.969483 1.385 3.053 neighbor 0.6 bin neigh_modify every 5 delay 0 check yes velocity all create $x 825577 dist gaussian velocity all create 10 825577 dist gaussian fix 1 all npt $x $x 2.0 xyz 0.0 0.0 45.0 drag 0.2 #在x温度下保持外压为零弛豫 fix 1 all npt 10 $x 2.0 xyz 0.0 0.0 45.0 drag 0.2 fix 1 all npt 10 10 2.0 xyz 0.0 0.0 45.0 drag 0.2 compute 3 all pe/atom compute 4 all ke/atom compute 5 all coord/atom 3.0 run 10000 Memory usage per processor = 40.9131 Mbytes Step Temp E_pair E_mol TotEng Press Volume 0 10 -1187595.7 0 -1187540.5 -3305841.1 671447.7 ERROR: Lost atoms: original 42720 current 41050 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 请教一下大神们，这个问题要怎么解决啊？我试图加一段能量最小化的代码，可是运行后是各种数据都输出为0……同样不知道怎么解决~ 请教大神，望各位能积极援助，不胜感激~~~ @老虎大王 @月只蓝 |
回复此楼» 猜你喜欢
346求调剂[0856]
已经有7人回复
-
一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！
已经有3人回复
-
310求调剂
已经有3人回复
-
机械专硕299求调剂至材料
已经有4人回复
-
070300化学319求调剂
已经有7人回复
-
08工科 320总分 求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有12人回复
-
307求调剂
已经有3人回复
-
265求调剂
已经有3人回复
-
301求调剂
已经有10人回复
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2014-09-27 19:40:33
3楼2018-10-25 10:33:05
