| 查看: 755 | 回复: 2 | |||
[交流]
【求助】边界能量处理
|
|
求教，我在一个box里面设置上下两个region不同的晶格取向，左边的晶格是右边绕Z方向旋转45度所得，我本意是研究界面的能量，但是结果边界处的能量异常高，而界面处的能量却无法突出，求教如何解决这个问题！！！！in文件如下 units metal boundary p p p atom_style atomic lattice fcc 3.61 region box block 0 6 0 6 0 36 create_box 1 box create_atoms 1 box region up block INF INF INF INF 18 INF delete_atoms region up lattice none lattice fcc 3.61 orient x 1 1 0 orient y -1 1 0 create_atoms 1 region up pair_style eam/alloy pair_coeff * * jin_copper_lammps.setfl Cu neighbor 0.5 bin neigh_modify every 5 delay 0 check yes timestep 0.01 thermo 10 variable N equal atoms variable E equal pe compute pe all pe/atom compute ke all ke/atom compute coord all coord/atom 3.0 dump 1 all cfg 1 a.*.cfg id type xs ys zs c_pe c_ke c_coord dump_modify 1 element Cu run 0 undump 1 min_style cg minimize 1.0e-30 1.0e-15 1000 1000 print "^---------fcc_boundary------------- N=$N E=$E ------------------------"  |
回复此楼» 猜你喜欢
面上意见出来了
已经有16人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
国自然面上复盘~欢迎讨论
已经有16人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
有没有仍没收到信息的
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 振动与波动理论 讨论
已经有6人回复
- 求 摘要 引言 结论翻译 （我已经自己翻译了一遍，希望高人帮忙修改一下）
已经有1人回复
- 极性表面的表面能怎么计算呢
已经有9人回复
- 【求助完毕】vasp-slab-边界附近原子顶位吸附小分子的结构优化很难收敛吗？
已经有5人回复
- 【求助】vasp为什么采用周期性边界条件（或超原胞模型）来处理原子分子表面体系等，有
已经有8人回复
- 【求助】纳米线拉伸性时，截面的面积大小
已经有20人回复
2楼2011-01-09 16:24:03
3楼2011-01-13 10:06:45
30