| 查看: 695 | 回复: 1 | |||
[交流]
升温熔化模拟求助，没报错，但也不在计算 已有1人参与
|
|
我这两天做了一个金属化合物升温的模拟，但是遇到一些问题，下面是in文件，还望高人看看 # in.melt # --------------------- INITIALIZAITION --------------------- clear units metal dimension 3 boundary p p p atom_style atomic variable a equal 3.572 variable N equal 88353 # --------------------- ATOM DEFINITION --------------------- lattice custom ${a} a1 1 0 0 a2 0 1 0 a3 0 0 1 basis 0.5 0.5 0 basis 0.5 0 0.5 basis 0 0.5 0.5 basis 0 0 0 & orient x 1 -1 0 orient y 1 1 -2 orient z 1 1 1 region mybox sphere 0 0 0 50 units box create_box 2 mybox create_atoms 2 box basis 1 1 basis 2 1 basis 3 1 basis 4 2 mass 1 58.69 #Ni mass 2 26.98 #Al # --------------------- FORCE FIELDS ------------------------ pair_style eam/alloy pair_coeff * * NiAl.eam.alloy Ni Al # --------------------- SETTINGS ---------------------------- neighbor 0.5 bin neigh_modify every 5 delay 0 check yes timestep 0.01 min_style cg minimize 1.0e-12 1.0e-12 1000 1000 velocity all create 275 1234567 dist gaussian fix 1 all nvt temp 275 275 0.1 thermo_style custom step temp etotal pe ke press thermo 1000 compute 1 all pe/atom compute 2 all reduce sum c_1 run 15000 unfix 1 fix 1 all npt temp 275 2000 4 iso 0 0 6 drag 0.2 variable p equal c_2/${N} variable v equal vol/${N} fix def1 all print 1000 "${temp} ${p} ${v}" file potential_peratom.txt screen no dump 1 all atom 1000 dump.melt run 2000000 |
回复此楼» 猜你喜欢
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
jdl272
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 766.8
- 帖子: 67
- 在线: 129.7小时
- 虫号: 1464573
- 注册: 2011-10-28
- 专业: 凝聚态物性I:结构、力学和
2楼2018-02-16 22:39:16
20