24小时热门版块排行榜    

查看: 3955  |  回复: 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-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 04:45 by 4wMiSEwB6436
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 7/350 2026-08-15 00:40 by 4wMiSEwB6436
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +5 工厂打螺丝 2026-08-08 5/250 2026-08-14 11:03 by 玖戈弋
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见