24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ly67255776

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 850.6
帖子: 48
在线: 60.4小时
虫号: 732402

[交流] 【求助】边界能量处理

求教，我在一个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 ------------------------"

回复此楼