24小时热门版块排行榜    

查看: 461  |  回复: 0

上官七七

金虫 (小有名气)

[求助] 表面弛豫

刚学习lammps不久
做表面吸附原子扩散前 先对基底进行弛豫
结果中间的原子都不知道跑到哪去了  被掏空了
而且温度也不能稳定在设定值
请大家帮忙分析错误
in文件如下
# 3d EAM Ag surface diffusion

units                metal
boundary        p p s
atom_style        atomic

lattice                fcc 4.08 orient x 1 -1 0 orient y 1 1 -2 orient z 2 2 2
region                box block 0 5 0 5 0 10     #To define 5 lattice
create_box        1 box
create_atoms        1 box

region          TopAtom block INF INF INF INF 4.9 INF    #To define a delete region
group           Top region TopAtom
delete_atoms    group Top

region          Bottom block INF INF INF INF INF 2.9    #To define a fix region
group           Bottom region Bottom
###create_atoms    1 single 2.5 2.5 5.5
###create_atoms    1 single 2.5 2.5 -0.5

###variable        T index 450 500 550 600 650

###print           "--------Temperature=$T K-----------"

timestep        0.001

pair_style        eam
pair_coeff        * * Ag_u3.eam
###mass            1 107

fix             freeze Bottom setforce 0.0 0.0 0.0      #To fix the 2 layers of bottom atoms
velocity        Bottom set 0.0 0.0 0.0

neighbor        2 bin
neigh_modify        every 20 delay 0 check yes

fix             myRelax all nvt temp 300.0 300.0 100.0
compute         3 all pe/atom
compute         4 all ke/atom
compute         5 all coord/atom 1.0
dump            id all xyz 50 dump_sur2.atom

thermo_style    custom step temp etotal press
thermo          100
run                5000
回复此楼

» 猜你喜欢

德胜才即君子
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 上官七七 的主题更新
信息提示
请填处理意见