24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3792  |  回复: 3

幽夜流星雨

银虫 (初入文坛)

[求助] LAMMPS Tersoff提示出错 已有1人参与

我在用LAMMPS编写铜在硅上沉积的过程中,使用了Tersoff style,然后提示Tersoff使用的不正确,找了好久问题没找到,因此求教各位大神。错误如图,代码如下。其中加下划线部分为提示出错的地方,但是我认为不应该是这里的问题。
附in文件:
#input file for copper deposition

clear
units metal
dimension 3
atom_style atomic
boundary p p f

#Definition of the simulation box
lattice fcc 5.43
region box block 0 20 0 20 0 22
create_box 4 box

#Definition of mass
mass 1 28
mass 2 28
mass 3 28
mass 4 63.5

#Definition of pair style and pair coefficient
newton on
pair_style tersoff
pair_coeff * * Si.tersoff Si
pair_coeff 3 4 Si.tersoff.modc Si Cu
pair_style eam/alloy
pair_coeff 4 4 Cu_u3.eam Cu

#Definition of substrate atoms
region substrate_down block 0 20 0 20 0 1
create_atoms 1 region substrate_down
group fix_atoms region substrate_down

region substrate_middle block 0 20 0 20 1 4
create_atoms 2 region substrate_middle
group temperature_atoms region substrate_middle

region substrate_up block 0 20 0 20 4 6
create_atoms 3 region substrate_up
group mobile_atoms region substrate_middle

neigh_modify delay 0

#Definition of the deposition atoms
group Cu_atoms type 4
region atoms_region block 0 20 0 20 6 22
group atoms_region region atoms_region

#Definition of the computation information
compute atoms_temperature atoms_region temp
compute_modify        atoms_temperature dynamic yes extra 0

#Setting of other conditions
fix                1 fix_atoms nve
fix                2 temperature_atoms langevin 220 220 1 587283
fix                3 temperature_atoms nve
fix                4 mobile_atoms langevin 220 220 1 587283
fix                5 mobile_atoms nve
fix                6 Cu_atoms nve

#Position of deposition atoms
region          slab block 0 20 0 20 16 22

#
fix                7 Cu_atoms deposit 2200 4 100 12345 region slab near 4.0 &
                         vz -0.0055 -0.0055
fix                8 Cu_atoms wall/reflect zhi EDGE
thermo_style        custom step atoms temp epair etotal press
thermo          100
thermo_modify        temp atoms_temperature lost/bond ignore lost warn

dump                1 all atom 50 dump.deposit.atom

#dump                2 all image 50 image.*.jpg type type &
#                axes yes 0.8 0.02 view 80 -30
#dump_modify        2 pad 5

dump                3 all movie 50 movie.mpg type type &
                axes yes 0.8 0.02 view 80 -30
dump_modify        3 pad 5

run             300000

LAMMPS Tersoff提示出错
错误.png
回复此楼

» 猜你喜欢

自己不坚强,懦弱给谁看?
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zx108s221

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
幽夜流星雨: 金币+36, ★★★很有帮助 2017-12-15 22:50:58
你好几个pair style  定义的时候要用hybrid
回去好好看手册吧 手册上写的非常清楚明白
I'maslowlywalker,butIneverwalkback.
2楼2017-12-14 18:18:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

幽夜流星雨

银虫 (初入文坛)

找到问题的原因了,一部分如楼上所述,多个pair_style时,应该用hybrid,而不是单独列出每一个
此外,多少种原子类型,tersoff后面则需要出现对应数目的参数。
自己不坚强,懦弱给谁看?
3楼2017-12-15 22:55:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
4楼2018-01-14 18:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 幽夜流星雨 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工349分求调剂 +14 李木子啊哈哈 2026-04-12 15/750 2026-04-14 00:14 by Equinoxhua
[考研] 085600材料与化工,求调剂 +14 won_qii 2026-04-07 14/700 2026-04-13 22:21 by pies112
[考研] 食品与营养(0955)271求调剂 +9 升格阿达 2026-04-12 9/450 2026-04-13 14:12 by maddjdld
[考研] 368化学求调剂 +14 wwwwabcde 2026-04-07 15/750 2026-04-13 08:36 by lhj2009
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 求调剂 +6 电气300求调剂不 2026-04-08 6/300 2026-04-11 20:14 by 逆水乘风
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 0854调剂 +5 音像店听花鼓戏 2026-04-10 5/250 2026-04-11 10:49 by qingpingzhu
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 材料与化工调剂 +12 否极泰来2026 2026-04-10 13/650 2026-04-11 00:28 by wangjihu
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 一志愿中科大070300化学,314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 332,085601求调剂 +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[考研] 337求调剂 +4 Gky09300550, 2026-04-09 4/200 2026-04-09 17:18 by 帕尔马拉特
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 287求调剂 +6 Fnhc 2026-04-07 6/300 2026-04-08 10:05 by xingguangj
信息提示
请填处理意见