24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

2009云飞扬

铁杆木虫 (著名写手)

应助: 1 (幼儿园)
金币: 6399.2
帖子: 1796
在线: 364.5小时
虫号: 1745038

[交流] pair_style hybrid/overlay指令用法

各位大神关于pair_style hybrid/overlay指令用法小弟有点不是太明白，能否详细说说呢？如这个例子
using, EAM, Tersoff and LJ for Fe-Fe,C-C and Fe-C interactions respectively. Fe is atom type 1, C is atom type 2.
pair_style hybrid eam/fs lj/cut 10.0 tersoff
pair_coeff * * eam/fs ../potentials/Fe_mm.eam.fs Fe NULL
pair_coeff 1 2 lj/cut 0.00759200 3.5 9.7625
pair_coeff * * tersoff ../potentials/Si.tersoff NULL Si
它包含了三个势函数类型，两种原子类型Fe、C，在系数设置上面，第一条 ** 不是代表所有原子的相互作用吗？而原本是想铁铁的作用势选用EAM势，铁碳作用势力用的是LJ对势，这个我能看懂，可是最后一条不是太懂。“NULL”是代表缺省吗?

回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2018-12-26 18:44:30

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

Nater1ver

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3481.2
帖子: 768
在线: 47.9小时
虫号: 1526557

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
2009云飞扬(金币+1): 谢谢参与
2009云飞扬: 金币+10, 很有帮助，谢谢 2018-12-30 14:01:01

随便简单说一说。势函数文件自己打开看一下就知道，是一个索引，你引用的元素名是文件里面引的，而in文件中，但凡要定义和使用这种势函数，必须遍历in文件中使用的所有原子种类，因此会出现 pair_xxx ** xxx.eam type1 type2 type3 type4这种定义方法(type1到type4表示每种原子在势函数索引中对应的名称，并不代表原子自己，即使type1是氟，你完全可以在这里填上氯，用氯原子的势函数来计算，即使这种操作是不正常的)，你可以看到，它是按照type顺序声明的。即使有的原子没有这种势函数的相互作用，这些原子也要包含进去(这显然与常见的势函数比如lj不同，lj中只用定义参与相互作用的原子类型)，当然假如没有参与相互作用就用NULL代替了，表示不参与作用

发自小木虫Android客户端