| 查看: 530 | 回复: 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 |
回复此楼» 猜你喜欢
国自然面上复盘~欢迎讨论
已经有18人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有8人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
20