24小时热门版块排行榜    

查看: 2405  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 WqtG0Aem8Ufl 2026-10-24 5/250 2026-10-26 02:55 by cLOojpOxkX0G
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 5/250 2026-10-26 02:51 by cLOojpOxkX0G
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3xN9SualEUPt 2026-10-24 6/300 2026-10-26 02:42 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:25 by cLOojpOxkX0G
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-26 02:16 by cLOojpOxkX0G
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-26 02:03 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 5/250 2026-10-26 01:26 by cLOojpOxkX0G
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-26 01:07 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:54 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:45 by cLOojpOxkX0G
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 7xHjyB3bAiwg 2026-10-23 9/450 2026-10-26 00:34 by cLOojpOxkX0G
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:32 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:19 by cLOojpOxkX0G
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 3/150 2026-10-25 14:23 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:18 by XiYP04It3Evb
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-25 00:15 by XcOYG53d3x4p
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 18:23 by WqtG0Aem8Ufl
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 15:41 by Equinoxhua
[文学芳草园] 亲密关系 +7 宋昙 2026-10-21 13/650 2026-10-23 17:07 by 小牧11
[考博] 郑州大学申博 +5 stars~ 2026-10-19 8/400 2026-10-22 15:06 by stars~
信息提示
请填处理意见