24小时热门版块排行榜    

查看: 3957  |  回复: 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的回帖
4楼2018-01-14 18:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 04:21 by 4wMiSEwB6436
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 奇怪,两个人的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 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 读博的好处 +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能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 为什么网上很多人说本周 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有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
信息提示
请填处理意见