24小时热门版块排行榜    

查看: 2199  |  回复: 1

山涧清泉HC

新虫 (初入文坛)

[求助] BN如何建模

如何用lammps建立六方结构的氮化硼(h-BN)的模型?顺便帮忙看一下我这样建模的问题。附件是h-BN的晶体结构。
# ------------------------ INITIALIZATION ----------------------------
units                 metal#单位?
dimension        3
boundary        p        p        p
atom_style        atomic

#neighbor      2.0 bin
#neigh_modify  one 5
lattice       hcp 2.5



# ----------------------- ATOM DEFINITION ----------------------------
lattice custom 2.5 a1 0.23 0.0 0.0 a2 0.0 0.39 0.0 a3 0.0 0.0 0.38  &
             basis 0.5 0.0 0.0 basis 0.0 0.5 0.0 basis 0.5 0.0 0.5  &
             basis 0.0 0.5 0.5
region                whole block 0 5 0 5 0 5

create_box        2 whole
create_atoms  2 box basis 1 1 basis 2 2 basis 3 1 basis 4 2

# ------------------------ FORCE FIELDS ------------------------------
# Tersoff Potentials

pair_style    tersoff
pair_coeff    * * BNC.tersoff   B   N
mass         1 10
mass         2 14
######################################
write_data data.BN
# EQUILIBRATION
reset_timestep        0
timestep 0.001
velocity all create 300 12345 mom yes rot no
fix 1 all npt temp 300 300 1 iso 0 0 1 drag 1

# Use cfg for Ovito
#dump 1 all custom 10 dump.comp.* id type xs ys zs c_csym c_eng c_strs[1]&
      c_strs[2] c_strs[3]
#dump_modify 1 element B  N

# Set thermo output
thermo 1000
thermo_style custom step lx ly lz press pxx pyy pzz pe temp

# Run for at least 10 picosecond (assuming 1 fs timestep)
run 20000
unfix 1

# Store final cell length for strain calculations
variable tmp equal "lx"
variable L0 equal ${tmp}
print "Initial Length, L0: ${L0}"

######################################
# SIMULATION DONE
print "All done"

BN如何建模
75X]CK@5]Y8QCUUYZ)W@KBO.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

山涧清泉HC

新虫 (初入文坛)

嗯,对了,这个应该是hex
2楼2015-03-31 10:06:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 山涧清泉HC 的主题更新
信息提示
请填处理意见