| 查看: 3954 | 回复: 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”是代表缺省吗? |
回复此楼» 猜你喜欢
全人类对于AI的整体认知，尚且停留在幼儿园启蒙阶段
已经有7人回复
-
自然基金已经彻底沦为学术资源的交换工具
已经有29人回复
-
周末思考-婚姻长久底色
已经有12人回复
-
27考博求助
已经有4人回复
-
深圳大学材料学院“智能分子材料团队”招聘博士后
已经有4人回复
-
面上的水很深
已经有32人回复
-
国社科通讯评审
已经有8人回复
» 抢金币啦！回帖就可以得到:
-
浙江理工大学能源催化团队诚聘专职教师/博士后
+2/282
-
中国石油大学（北京）未来能源学院“电催化能源工程”课题组博士后招聘
+1/171
-
不用懂统计，用初中数学锤爆这篇造假学位论文：10+处铁证全解析
+1/76
-
钙钛矿团队诚聘青年博士、博士后 2–4 人
+1/76
-
浙江大学侯阳团队(杰青)高薪招聘博士后
+1/68
-
哈尔滨工业大学物理学院光电材料和器件团队：招收2027级秋季和2028级春季博士生
+1/28
-
帮杭州程序员同学发帖
+1/28
-
海南大学教职 入职要求 考核要求
+1/25
-
国家自然科学基金申报者发榜前心理状态调查，麻烦帮忙填一下哇！
+3/12
-
东北大学招收2027级硕士、博士研究生（推免、报考）凝聚态物理/微电子/电子等相关方向
+1/7
-
东北大学招收2027级硕士、博士研究生（推免、报考）凝聚态物理/微电子/电子等相关方向
+1/7
-
中国地质大学（北京）李梅教授团队招生（含推免生，人工智能方向）
+1/6
-
东北大学招收2027级硕士、博士研究生（推免、报考）凝聚态物理/微电子/电子等相关方向
+1/6
-
东北大学招收2027级硕士、博士研究生（推免、报考）凝聚态物理/微电子/电子等相关方向
+1/6
-
国自然重大研究计划培育项目啥时候公布结果
+1/5
-
2027年直招入伍攻读硕士研究生招生-军事医学研究院生物技术药物代谢研究方向
+1/5
-
中国地质大学（北京）李梅教授团队博士招生（人工智能方向）
+1/5
-
中科院东莞材料所招聘科研助理（长期有效）
+1/4
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/2
-
高级智能体工作流架构师
+1/1
送红花一朵 |
谢谢朋友 在咨询咨询，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
6楼2018-12-26 22:41:13
8楼2018-12-27 04:04:16
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
简单回复
2018-12-26 19:52
回复
