| 查看: 3612 | 回复: 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”是代表缺省吗? |
回复此楼» 猜你喜欢
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
磺酰氟产物，毕不了业了！
已经有6人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有9人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 抢金币啦！回帖就可以得到:
-
北京-89175-事业单位-诚征女友
+1/89
-
福建师范大学柔性电子学院招收2026年博士（储能材料与柔性电子器件）
+1/84
-
北京—征老婆
+1/82
-
中国石油大学（华东）电气工程专业博士研究生招生
+1/81
-
湖南师范大学医工交叉科研团队招收计算机博士生
+1/80
-
上海大学昝鹏教授、军事医学研究院伯晓晨研究员/倪铭副研究员 课题组招聘博士生
+1/77
-
美国密歇根州立大学林学系杜海顺课题组招收全奖博士生及联合培养博士生
+1/75
-
山东征女友，坐标济南
+1/60
-
SCI，计算机相关可以写
+1/26
-
SCI，计算机相关可以写
+1/23
-
太原理工大学电工部招聘老师-偏电类专业的博士们快来看啊
+1/23
-
SCI，计算机相关可以写
+1/20
-
福建师范大学柔性电子学院招收2026年博士（储能材料与柔性电子器件）
+2/18
-
中国科学院大学-杨晗课题组-诚聘-博士后、副研究员
+1/15
-
中国科学院大学-杨晗课题组-诚聘-博士后、副研究员
+1/14
-
上海理工顾敏院士/李蔚团队招收2026级博士研究生 （集成光学、量子信息方向）
+1/9
-
博士/硕士招生
+1/7
-
中科院动物所招收2026年博士生（优先少干专项计划、化学或生命科学背景）
+1/3
-
土木、交通工程专业博士后站有吗？（无博士毕业3年要求+可接受兼职博后）
+1/1
-
电子科技大学激光与光子制造团队招硕士博士
+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
回复