| 查看: 3906 | 回复: 30 | |||
[交流]
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”是代表缺省吗? |
回复此楼» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
有谁可曾问过你过的还好吗？
已经有21人回复
-
E0414, 我的本子有没有希望？
已经有7人回复
-
一篇论文同时出现在两个期刊，一模一样，这算不算学术不端，请各位老师斧正。
已经有12人回复
-
希望面上有个好结果
已经有7人回复
-
三区计算机方向期刊推荐
已经有5人回复
-
sci论文二审求助
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
郑州大学急招1名2026级博士生
+1/471
-
硫化物 / 氧化物 / 卤化物全拆解：固态电池手套箱选型核心指南
+1/84
-
多聚组氨酸，大分子量
+1/80
-
Cu2O纳米线
+5/80
-
坐标北京，诚征男友，找到会删帖。试试运气！说不定就有了！
+1/78
-
青年拔尖人才还有指标 25w*8年
+1/74
-
哈工大深圳-材料学院-招收申请审核制2027年春季/秋季入学博士生（2026年9月报名）
+1/36
-
东南大学有机多孔功能材料团队（国家杰青团队）招收2027级推夏令营免硕士/直博生
+1/33
-
上海交通大学刘祖明课题组招聘博士后
+1/31
-
博士毕业之后申请发明专利
+1/28
-
重庆大学金属塑性成形方向招收2027年博士研究生（学博）
+2/22
-
墨尔本大学(26年QS19)招全奖博士/CSC博士(补齐全奖)等-材料/生物电子/器官芯片等
+1/11
-
墨尔本大学(26年QS19)招全奖博士/CSC博士(补齐全奖)/访问学者等-材料/生物医学等
+1/8
-
同济大学环境学院 肖倩研究员课题组 招聘2027级硕士/博士
+1/8
-
天津大学高分子方向博士招生
+1/6
-
密苏里大学生物材料合成生物学博士后招聘
+1/6
-
美国普渡大学(Purdue University)生物化学系招聘博士后
+1/6
-
上海大学微电子学院杨军教授团队招聘带编专任教师
+1/3
-
上海大学微电子学院杨军教授团队招聘带编专任教师
+1/2
-
中科院深圳先进技术研究院成会明院士/唐永炳杰青团队博士后招聘
+1/1
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
2009云飞扬(金币+1): 谢谢参与
2009云飞扬: 金币+10, 很有帮助，谢谢 2018-12-30 14:01:14
2009云飞扬(金币+1): 谢谢参与
2009云飞扬: 金币+10, 很有帮助，谢谢 2018-12-30 14:01:14
|
NULL表示的意思是不对这个类型原子/分子施加那种类型的作用势，比如说你现在有两种类型原子，Fe的原子序号是1，C的原子序号是2。在分配作用势的时候比如只是对于C的作用势时就必须忽略Fe。 发自小木虫Android客户端 |
12楼2018-12-27 10:28:12
6楼2018-12-26 22:41:13
8楼2018-12-27 04:04:16
10楼2018-12-27 06:14:21
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
2009云飞扬(金币+1): 谢谢参与
2009云飞扬: 金币+10, 很有帮助，谢谢 2018-12-30 14:01:01
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客户端 |
» 本帖已获得的红花（最新10朵）
11楼2018-12-27 09:18:05
送红花一朵 |
谢谢朋友 在咨询咨询，hybrid/overlay我理解的是把几个作用势集合起来，在咨询咨询“pair_xxx ** xxx.eam type1 type2 type3 type4”我原来理解的是在eam势函数文件里面，包含的原子种类，就是后面的type1-4。比如Fe_mm.eam.fs ，这个势函数里面就只有铁的势函数，因而写成pair_coeff 1 1 eam/fs ../potentials/Fe_mm.eam.fs Fe （铁在in文件中代表原子类型为1），。你的意思是说in文件里面假如有两种原子类型，则在定义势函数系数的时候，都必须遍历取值，因而在这条指令pair_coeff * * eam/fs ../potentials/Fe_mm.eam.fs Fe NULL 是** 表示in文件中原子类型（铁 、碳）所有原子都参与了（铁铁、铁碳、碳碳三种作用势都采用Fe_mm.eam.fs 势函数，但是事实上它只是铁铁的势函数，而我们又只需要表示铁铁的作用势力，因而后面type类型中标注为Fe NULL ，“NULL ”的作用就是表示碳原子没有参与到这个势函数作用）。不知道我理解的是否正确，请你指点指点，实在是有点迷糊 了。 |
13楼2018-12-27 16:27:25
★
2009云飞扬(金币+1): 谢谢参与
2009云飞扬(金币+1): 谢谢参与
|
本帖内容被屏蔽 |
29楼2020-01-03 15:51:30
★
2009云飞扬(金币+1): 谢谢参与
2009云飞扬(金币+1): 谢谢参与
|
本帖内容被屏蔽 |
30楼2020-01-06 11:40:15
★
2009云飞扬(金币+1): 谢谢参与
2009云飞扬(金币+1): 谢谢参与
|
本帖内容被屏蔽 |
31楼2020-01-09 16:08:51
简单回复
youngen2楼
2018-12-26 19:16
回复
2009云飞扬(金币+1): 谢谢参与
tzynew3楼
2018-12-26 19:42
回复
2009云飞扬(金币+1): 谢谢参与
2018-12-26 19:52
回复
2018-12-26 22:24
回复
2009云飞扬(金币+1): 谢谢参与
psylhh7楼
2018-12-26 22:51
回复
2009云飞扬(金币+1): 谢谢参与
2018-12-27 05:11
回复
2009云飞扬(金币+1): 谢谢参与
syhorchid14楼
2018-12-27 16:56
回复
2009云飞扬(金币+1): 谢谢参与
nono200915楼
2018-12-27 21:03
回复
2009云飞扬(金币+1): 谢谢参与
·
Biotage16楼
2018-12-27 23:13
回复
2009云飞扬(金币+1): 谢谢参与
baowahao17楼
2018-12-27 23:19
回复
2009云飞扬(金币+1): 谢谢参与
sz_kebaoyuan18楼
2018-12-27 23:20
回复
2009云飞扬(金币+1): 谢谢参与
syhorchid19楼
2018-12-27 23:57
回复
miaojiabing20楼
2018-12-28 08:30
回复
2009云飞扬(金币+1): 谢谢参与
miaojiabing21楼
2018-12-29 08:29
回复
youngen22楼
2018-12-31 13:04
回复
xhmaohan23楼
2018-12-31 13:24
回复
slytjiaofei24楼
2018-12-31 15:30
回复
guanlianwu25楼
2019-01-01 17:08
回复
2009云飞扬(金币+1): 谢谢参与
xautldhldh26楼
2019-01-02 14:30
回复
2009云飞扬(金币+1): 谢谢参与
瑞驰雨林27楼
2019-07-27 20:17
回复
doublerobust28楼
2019-12-31 09:45
回复