24小时热门版块排行榜

返回列表

» 猜你喜欢

1楼 2015-08-10 11:29:42

已阅关注TA 给TA发消息送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

专家经验: +49
模拟EPI: 2
应助: 148 (高中生)
金币: 8640.9
散金: 507
红花: 30
帖子: 505
在线: 598.4小时
虫号: 556471
注册: 2008-05-11
性别: GG
专业: 理论和计算化学
管辖: 分子模拟

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
303103263(abinitio代发): 金币+5, 鼓励交流，正解 2015-08-10 15:04:05

http://lammps.sandia.gov/doc/pair_sw.html
lammps手册上说不是说你要包含一个sw文件来指定你的sw势的参数吗？

“For a single-element simulation, only a single entry is required (e.g. SiSiSi). For a two-element simulation, the file must contain 8 entries (for SiSiSi, SiSiC, SiCSi, SiCC, CSiSi, CSiC, CCSi, CCC), that specify SW parameters for all permutations of the two elements interacting in three-body configurations. Thus for 3 elements, 27 entries would be required, etc.”

如果你有三种元素，在sw文件里要指定27行参数。怎么能说lammps中sw参数的输入只针对两种元素的情况呢？

赞一下

回复此楼

2楼2015-08-10 14:01:52

已阅关注TA 给TA发消息送TA红花 TA的回帖

303103263

禁虫 (小有名气)

本帖内容被屏蔽

3楼2015-08-10 16:27:32

已阅关注TA 给TA发消息送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

专家经验: +49
模拟EPI: 2
应助: 148 (高中生)
金币: 8640.9
散金: 507
红花: 30
帖子: 505
在线: 598.4小时
虫号: 556471
注册: 2008-05-11
性别: GG
专业: 理论和计算化学
管辖: 分子模拟

引用回帖:

3楼: Originally posted by 303103263 at 2015-08-09 20:27:32
我不是这个意思。因为SW是个三体势，在它的函数表达式里有σij和σik，αij和αik这样的参数，对于c-si-c来说，j和k代表的都是c原子，因此参数σij和σik一样，αij和αik一样，对c-si-o来说，σij和σik是不一样的 ...

这的确是个问题，不过这个问题之前已经有人在lammps mailing list里面提出了
http://lammps.sandia.gov/threads/msg16418.html

为了计算C-Si-O的三体作用，需要C-Si和C-O的两体作用参数。这些两体作用参数的sigma和epsilon可以从C-O-O和C-Si-Si所在行的sigma和epsilon获得。也就是说，在sw文件中，C-Si-O所在行的sigma和epsilon其实并不参与到C-Si-O三体作用的计算中。

如同提问者所说的那样，这可能是lammps在处理Stillinger-Weber势的一个缺陷，但也只好这样了。

赞一下

回复此楼

4楼2015-08-11 09:09:26

已阅关注TA 给TA发消息送TA红花 TA的回帖

303103263

禁虫 (小有名气)

本帖内容被屏蔽

5楼2015-08-12 15:44:15

已阅关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 303103263 的主题更新

返回列表