24小时热门版块排行榜    

查看: 3478  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085700资源与环境308求调剂 +7 墨墨漠 2026-03-18 7/350 2026-03-20 05:55 by laoshidan
[考研] 294求调剂材料与化工专硕 +14 陌の森林 2026-03-18 14/700 2026-03-19 22:38 by 学员8dgXkO
[考研] 一志愿武汉理工材料工程专硕调剂 +5 Doleres 2026-03-19 5/250 2026-03-19 20:14 by 制度的
[考研] 321求调剂 +8 何润采123 2026-03-18 10/500 2026-03-19 16:46 by 何润采123
[考研] 本人考085602 化学工程 专硕 +17 不知道叫什么! 2026-03-15 19/950 2026-03-19 15:06 by 尽舜尧1
[考研] 324分 085600材料化工求调剂 +3 llllkkkhh 2026-03-18 3/150 2026-03-19 14:22 by houyaoxu
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 274求调剂 +6 S.H1 2026-03-18 6/300 2026-03-19 09:34 by 花店相见
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
信息提示
请填处理意见