24小时热门版块排行榜    

查看: 3894  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ningpingm

银虫 (正式写手)

[求助] 混合势函数的问题,求教!感激不尽!问题如果解决了再继续追加50金币! 已有2人参与

今天自己倒腾了一个in文件,想把硅的SW势函数和ZBL势函数都用上,试了好几种写法都不行,有大神能指教一二吗?
一直在报错:ERROR: Pair coeff for hybrid has invalid style<../pair_hybrid_overlay.cpp:68>

请教大神们,下面的语句到底哪里错了?
pair_style hybrid/overlay sw zbl 1.1 1.8
pair_coeff * * Si.sw Si
pair_coeff * * zbl 14.0

如果帮我解决了,再追加金币,谢谢啦!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

影子政

铜虫 (初入文坛)

谢谢楼主和几个大神的解答,也解了我的燃眉之急,我还有个疑问:

@王雨辰phy 所提到的:
   pair_coeff * * SiC.sw Si Si Si C
那么pair_coeff * * zbl 14.0 怎么改变?
在 LAMMPS manual 中提到ZBL时:
The following coefficients must be defined for each pair of atom types via the pair_coeff command as in the examples above, or in the LAMMPS data file.
Z_i (atomic number for first atom type, e.g. 13.0 for aluminum)
Z_j (ditto for second atom type)
我想问的是在pair_coeff * * zbl  怎么定义Si 、C 的 原子序数。

看LAMMPS中说
When used with hybrid/overlay and pairs are assigned to more than one sub-style, the mixing rule is not used and each pair of types interacting with the ZBL sub-style must be included in a pair_coeff command.
我试了一下,就只用了两种原子类型:
pair_style   hybrid/overlay eam/cd zbl 1 2
pair_coeff   * * eam/cd FeCr.cdeam Fe Cr
pair_coeff   1 1 zbl 26
pair_coeff   2  2 zbl 24
这样可以,但是pair_coeff 1 2 又怎么定义呢? 求问大神们!!
if I rest,I rust!
9楼2016-02-29 11:42:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

xbb20121991

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
ningpingm: 金币+100, ★★★★★最佳答案, 你的回答完美地解决了我的问题!谢谢! 2015-08-03 12:25:32
pair_style hybrid/overlay sw zbl 1.1 1.8
pair_coeff * * sw Si.sw Si (我感觉这行应该填上sw,楼主原句子中第二行没有声明势的style)
pair_coeff * * zbl 14.0

我不做硅,用这类势比较少。你先试试。
我觉得sw的例子可以和tersoff的类比一下
2楼2015-08-03 07:10:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王雨辰phy

木虫 (小有名气)

老菜鸟

【答案】应助回帖

★
感谢参与,应助指数 +1
ben_ladeng: 金币+1, 辛苦了 2015-08-04 06:06:51
As an example, imagine a file SiC.sw has Stillinger-Weber values for Si and C. If your LAMMPS simulation has 4 atoms types and you want the 1st 3 to be Si, and the 4th to be C, you would use the following pair_coeff command:

pair_coeff * * SiC.sw Si Si Si C

你在in文件里定义了几种原子类型,sw后面就要分别罗列出。
楼上说的没有声明style也是正确的

» 本帖已获得的红花(最新10朵)

但行好事,莫问前程
3楼2015-08-03 09:32:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangdu98

铁虫 (小有名气)

我今天上午试了一下,确实是这个问题~ 我用lammps不熟,多谢你的帮忙,帮我解决大问题了,谢谢!金币给你啦!
4楼2015-08-03 12:04:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 BMTdnp4OJO9f 2026-11-06 5/250 2026-11-07 17:56 by Xc2FxHsHkBNc
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 2rKfVQg2jQww 2026-11-06 5/250 2026-11-07 17:34 by Xc2FxHsHkBNc
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 2rKfVQg2jQww 2026-11-06 5/250 2026-11-07 17:23 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 17:01 by Xc2FxHsHkBNc
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 4/200 2026-11-07 16:28 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 5sJFDKp2kwuI 2026-11-06 6/300 2026-11-07 15:44 by Xc2FxHsHkBNc
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
请填处理意见