| 查看: 699 | 回复: 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 |
回复此楼» 猜你喜欢
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有11人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
面上函评意见出来了，像什么等级？
已经有14人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
-
面上意见出来了
已经有23人回复
jdl272
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 766.8
- 帖子: 67
- 在线: 129.7小时
- 虫号: 1464573
- 注册: 2011-10-28
- 专业: 凝聚态物性I:结构、力学和
2楼2018-02-16 22:39:16
5
投票:
