24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3211  |  回复: 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的回帖

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的回帖

tangdu98

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by xbb20121991 at 2015-08-03 07:10:08
pair_style hybrid/overlay sw zbl 1.1 1.8
pair_coeff * * sw Si.sw Si (我感觉这行应该填上sw,楼主原句子中第二行没有声明势的style)
pair_coeff * * zbl 14.0

我不做硅,用这类势比较少。你先试试。
我觉 ...

我今天上午试了一下,确实是这个问题~ 我用lammps不熟,多谢你的帮忙,帮我解决大问题了,谢谢!金币给你啦
5楼2015-08-03 12:05:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbb20121991

银虫 (小有名气)

引用回帖:
5楼: Originally posted by tangdu98 at 2015-08-03 12:05:50
我今天上午试了一下,确实是这个问题~ 我用lammps不熟,多谢你的帮忙,帮我解决大问题了,谢谢!金币给你啦...

多谢楼主,最近没关注新版的lammps。都没注意到这个ZBL势,有时间好好研究一下怎么可以搭配。
这样也算是互相学习交流了
6楼2015-08-03 12:12:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ningpingm

银虫 (正式写手)

送红花一朵
引用回帖:
3楼: Originally posted by 王雨辰phy at 2015-08-03 09:32:18
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 foll ...

你好,你的答案也是正确的,但是我只能把金币给最先回答正确的人,不好意思。。。
不过我还是很感谢你的回复哦,真的很感激,谢谢!送你一朵小红花~
7楼2015-08-03 12:26:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ningpingm

银虫 (正式写手)

引用回帖:
6楼: Originally posted by xbb20121991 at 2015-08-03 12:12:04
多谢楼主,最近没关注新版的lammps。都没注意到这个ZBL势,有时间好好研究一下怎么可以搭配。
这样也算是互相学习交流了...

我这个账号是借来的,说好的追加50个金币,我应该怎么操作啊?sorry,对小木虫的这个追加机制不熟悉。。。
8楼2015-08-03 12:27:45
已阅   回复此楼   关注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的回帖

努力学习i

新虫 (初入文坛)

同类问题已解决非常感谢
10楼2019-07-09 20:11:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ningpingm 的主题更新
信息提示
请填处理意见