24小时热门版块排行榜    

查看: 3315  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 ningpingm 的主题更新
信息提示
请填处理意见