24小时热门版块排行榜    

查看: 3775  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +12 wangmingqi 2026-08-28 17/850 2026-08-29 08:43 by appleren
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 8/400 2026-08-29 08:32 by foolishmani
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 05:45 by gy1nBQXYQJqL
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见