| 查看: 3837 | 回复: 9 | ||
[求助]
混合势函数的问题，求教！感激不尽！问题如果解决了再继续追加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 如果帮我解决了，再追加金币，谢谢啦！！！ |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有14人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有5人回复
-
申请科研项目纯靠实力能成功吗？
已经有12人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
极限拉扯，如何抉择
已经有15人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Multiwfn支持的弱相互作用的分析方法概览
已经有10人回复
- 杂谈量子化学中的简化招数，和SCC-DFTB模型
已经有24人回复
xbb20121991
银虫 (小有名气)
- 应助: 18 (小学生)
- 金币: 443.6
- 红花: 7
- 帖子: 140
- 在线: 114.1小时
- 虫号: 3435358
- 注册: 2014-09-23
- 专业: 金属材料的微观结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
ningpingm: 金币+100, ★★★★★最佳答案, 你的回答完美地解决了我的问题！谢谢！ 2015-08-03 12:25:32
感谢参与，应助指数 +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
王雨辰phy
木虫 (小有名气)
老菜鸟
- 应助: 2 (幼儿园)
- 金币: 1830.5
- 散金: 410
- 红花: 69
- 帖子: 103
- 在线: 141.5小时
- 虫号: 2326670
- 注册: 2013-03-07
- 性别: GG
- 专业: 半导体材料
【答案】应助回帖
★
感谢参与，应助指数 +1
ben_ladeng: 金币+1, 辛苦了 2015-08-04 06:06:51
感谢参与，应助指数 +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朵）
ningpingm
3楼2015-08-03 09:32:18
4楼2015-08-03 12:04:39
5楼2015-08-03 12:05:50
xbb20121991
银虫 (小有名气)
- 应助: 18 (小学生)
- 金币: 443.6
- 红花: 7
- 帖子: 140
- 在线: 114.1小时
- 虫号: 3435358
- 注册: 2014-09-23
- 专业: 金属材料的微观结构
6楼2015-08-03 12:12:04
7楼2015-08-03 12:26:47
8楼2015-08-03 12:27:45
|
谢谢楼主和几个大神的解答，也解了我的燃眉之急，我还有个疑问： @王雨辰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 的 原子序数。 看ＬＡＭＭＰＳ中说 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 又怎么定义呢？ 求问大神们！！ |
9楼2016-02-29 11:42:37
努力学习i
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 25.5
- 散金: 10
- 帖子: 10
- 在线: 42.3小时
- 虫号: 11128011
- 注册: 2018-11-02
- 专业: 力学中的基本问题和方法
10楼2019-07-09 20:11:34
10