| 查看: 520 | 回复: 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 |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
100