24小时热门版块排行榜

返回列表

小白小孟

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 4
帖子: 5
在线: 6.4小时
虫号: 4451455
注册: 2016-03-01
专业: 凝聚态物性I:结构、力学和

[求助] 求救。。。烈日炎炎~有劳各位大神啦。小弟我愿奉上我所有金币已有1人参与

本人lammps新手，最近在计算钨的表面能，不知为什么一直提示 displace-box有错，检查了好多遍也没找到原因？？？还有就是（本人看很多例子里都要先删掉一部分原子，制造真空层。为什么啊？求赐教，不删可以吗？）是希望坛里大神动个手帮助小弟我一下。。拜谢拜谢
in文件
units          metal
atom_style    atomic
boundary       p p p
lattice       bcc 3.14
region       box block 0 20 0 20 0 40
create_box    3 box
create_atoms 1 box
group       W type 1
group       H type 2
group       He type 3
group       all region box
include pot.eam
compute 1 all pe/atom
compute 2 all ke/atom
variable E equal pe
timestep  0.005
min_style sd
minimize 1e-10 1.0e-10 10000 10000
print @@@@@=$E

region       up1 block INF INF INF INF 19.9 INF
region       down1 block INF INF INF INF INF 19.8
group       up1 region up1
group       down1 region down1
dump 1 all custom 1 all_minimize_before.lammpstrj id type x y z c_1 c_2
displace_box all z delta 0.0 40.0 remap none
displace_atoms up1 move 0 0 40
undump 1
min_style sd
minimize 1e-10 1.0e-10 10000 10000
print @@@@@-----final=$E

回复此楼