24小时热门版块排行榜    

查看: 3780  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 幽夜流星雨 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 一志愿东华大学化学070300,求调剂 +6 2117205181 2026-03-21 7/350 2026-03-22 14:59 by allen-yin
[考博] 招收博士1-2人 +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
信息提示
请填处理意见