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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 怎么使用pair_style hybrid以及怎么正确使用EAM势能？
51/1返回列表
查看: 2597  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cl1994861

金虫 (小有名气)

[求助] 怎么使用pair_style hybrid以及怎么正确使用EAM势能？ 已有1人参与

我想使用pair_style hybrid来使用EAM势能，代码如下，但是模拟得到的结果是水盒子中的铜颗粒完全发散成一个个原子了，也就是说Cu与Cu之间的力小于O与Cu之间的作用力，导致这样的结果，但是EAM势能不是更准确么，不应该出现这个问题啊，所以我怀疑是这个代码的问题，有没有大神指导下？

pair_style        hybrid lj/cut/coul/long 9 lj/cut 9 eam     
kspace_style    pppm 1.0e-5
pair_modify     tail yes

#interaction destruction
pair_coeff        1 1 lj/cut/coul/long 0.1553 3.166
pair_coeff        1 2 lj/cut/coul/long 0.0000 0.0000
pair_coeff        1 3 lj/cut 1.211 2.752
pair_coeff        2 2 lj/cut/coul/long 0.0000 0.0000
pair_coeff        2 3 lj/cut 0 0
pair_coeff        3 3 eam Cu_u3.eam
回复此楼

» 猜你喜欢
1楼 2017-11-20 18:40:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cl1994861

金虫 (小有名气)
引用回帖:
6楼: Originally posted by lx_PICO at 2017-11-24 10:46:01
你这个势函数组合是文章里给的吗？因为这个EAM应该是LAMMPS自带的那个，是通过纯金属体系物理性质拟合的，可以直接算和离子（H、O）的相互作用吗？ 我原来试过铁和 Y2O3的，直接这样加上去好像是不行的。

这个EAM势能是用的lammps自带的，但是我只是用它算铜铜的相互作用，其他铜与氧，以及水分子之间的作用我用的是LJ势能，我看到过有文献这样处理。
赞一下 回复此楼
7楼2017-11-24 11:32:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
pair_coeff        * * eam Cu_u3.eam NULL NULL Cu
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

cl1994861
2楼2017-11-22 10:09:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cl1994861

金虫 (小有名气)
送红花一朵
引用回帖:
2楼: Originally posted by lx_PICO at 2017-11-22 10:09:06
pair_coeff        * * eam Cu_u3.eam NULL NULL Cu

您好，我刚刚试了一下，出现这样的error：ERROR: Incorrect args for pair coefficients (../pair_eam.cpp:359)，能不能详细说一下这个怎么用
赞一下 回复此楼
3楼2017-11-22 10:55:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by cl1994861 at 2017-11-22 10:55:40
您好，我刚刚试了一下，出现这样的error：ERROR: Incorrect args for pair coefficients (../pair_eam.cpp:359)，能不能详细说一下这个怎么用...

我错了，这个纯的EAM势还是你原来那样写才对。你可以用pair_write把势能写出来看看对不对。
赞一下 回复此楼
4楼2017-11-22 15:12:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿A区211，22408 321求调剂 +7 随心所欲☆ 2026-04-15 8/400 2026-04-21 08:22 by Equinoxhua
[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +11 否极泰来2026 2026-04-15 13/650 2026-04-20 22:31 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考研] 297，工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 297，工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 一志愿华中农业071010，320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9，生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
信息提示
关闭
请填处理意见
关闭