应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » pair_style hybrid/overlay指令用法
311/1返回列表
查看: 4098  |  回复: 30
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

2009云飞扬

铁杆木虫 (著名写手)

[交流] 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

19861226ya

木虫 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
2009云飞扬(金币+1): 谢谢参与
2009云飞扬: 金币+10, 很有帮助，谢谢 2018-12-30 14:01:14
NULL表示的意思是不对这个类型原子/分子施加那种类型的作用势，比如说你现在有两种类型原子，Fe的原子序号是1，C的原子序号是2。在分配作用势的时候比如只是对于C的作用势时就必须忽略Fe。

发自小木虫Android客户端
赞一下(1人) 回复此楼
12楼2018-12-27 10:28:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2009云飞扬

铁杆木虫 (著名写手)
请大家指点指点吧

发自小木虫Android客户端
赞一下 回复此楼
6楼2018-12-26 22:41:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cumt2010meng

木虫 (著名写手)
★
2009云飞扬(金币+1): 谢谢参与
我是来领金币的～当然首先要送上祝福啦～～

发自小木虫Android客户端
赞一下 回复此楼
8楼2018-12-27 04:04:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuqiang68

至尊木虫 (文坛精英)
★
2009云飞扬(金币+1): 谢谢参与
加油！
回复此楼
10楼2018-12-27 06:14:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Nater1ver

木虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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朵）

2009云飞扬
11楼2018-12-27 09:18:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2009云飞扬

铁杆木虫 (著名写手)
送红花一朵
引用回帖:
11楼: Originally posted by Nater1ver at 2018-12-27 09:18:05
随便简单说一说。势函数文件自己打开看一下就知道，是一个索引，你引用的元素名是文件里面引的，而in文件中，但凡要定义和使用这种势函数，必须遍历in文件中使用的所有原子种类，因此会出现 pair_xxx ** xxx.eam ty ...

谢谢朋友 在咨询咨询，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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pl911

禁虫 (小有名气)
★
2009云飞扬(金币+1): 谢谢参与
本帖内容被屏蔽
29楼2020-01-03 15:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pl958

禁虫 (小有名气)
★
2009云飞扬(金币+1): 谢谢参与
本帖内容被屏蔽
30楼2020-01-06 11:40:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pl937

禁虫 (小有名气)
★
2009云飞扬(金币+1): 谢谢参与
本帖内容被屏蔽
31楼2020-01-09 16:08:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
youngen2楼
2018-12-26 19:16   回复  
2009云飞扬(金币+1): 谢谢参与
tzynew3楼
2018-12-26 19:42   回复  
2009云飞扬(金币+1): 谢谢参与
dongmings4楼
2018-12-26 19:52   回复  
2009云飞扬(金币+1): 谢谢参与
发自小木虫IOS客户端
slytjiaofei5楼
2018-12-26 22:24   回复  
2009云飞扬(金币+1): 谢谢参与
psylhh7楼
2018-12-26 22:51   回复  
2009云飞扬(金币+1): 谢谢参与
xiaoqiuqiu209楼
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   回复  
2009云飞扬(金币+1): 谢谢参与
发自小木虫Android客户端
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   回复  
2009云飞扬(金币+1): 谢谢参与
发自小木虫IOS客户端
相关版块跳转 我要订阅楼主 2009云飞扬 的主题更新
311/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！ +8 zju2000 2026-09-26 9/450 2026-09-29 13:47 by ahsoarli
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
关闭
请填处理意见
关闭