24小时热门版块排行榜    

查看: 3954  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 幽夜流星雨 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 7lpolszZVXgi 2026-08-14 4/200 2026-08-15 01:08 by 4wMiSEwB6436
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 00:52 by 4wMiSEwB6436
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] filecode +6 cratir 2026-08-14 10/500 2026-08-14 18:29 by 笑叹辞穷
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +5 工厂打螺丝 2026-08-08 5/250 2026-08-14 11:03 by 玖戈弋
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见