| 查看: 488 | 回复: 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 |
回复此楼» 猜你喜欢
献血感触
已经有6人回复
-
2026山东省优青
已经有4人回复
-
评审有感
已经有16人回复
-
反应很差，大量原料没有反应
已经有3人回复
-
Sci. Bull. 悲剧经验
已经有4人回复
-
找博士生导师
已经有7人回复
-
上海大学实验技术岗位非升即走
已经有11人回复
-
26/27申博自荐-锂/钠电池方向
已经有4人回复
-
同样的基金本子，换个专家直接从C变A！
已经有3人回复
-
别被青基扩招骗了！26年科研内卷才刚刚开始
已经有4人回复
5