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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 怎么使用pair_style hybrid以及怎么正确使用EAM势能？
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 2562  |  回复: 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的回帖

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的回帖

cl1994861

金虫 (小有名气)
引用回帖:
4楼: Originally posted by lx_PICO at 2017-11-22 15:12:51
我错了，这个纯的EAM势还是你原来那样写才对。你可以用pair_write把势能写出来看看对不对。...

我问了其他人，他们也说没问题，可是就是计算过程中，一个几十个原子的铜颗粒完全发散成一个个原子到水盒子去了，这明显与实际不符，而且我用lj势没有出现发散的问题，所以才问怎么正确用这个势能
赞一下 回复此楼
5楼2017-11-22 18:15:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
cl1994861: 金币+5, ★★★很有帮助 2017-11-24 11:32:42
你这个势函数组合是文章里给的吗？因为这个EAM应该是LAMMPS自带的那个，是通过纯金属体系物理性质拟合的，可以直接算和离子（H、O）的相互作用吗？ 我原来试过铁和 Y2O3的，直接这样加上去好像是不行的。
赞一下 回复此楼
6楼2017-11-24 10:46:01
已阅   回复此楼   关注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的回帖

18370993420
感谢分享

发自小木虫Android客户端
8楼2017-11-24 20:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zx108s221

铜虫 (小有名气)
你把所有的水分子都删掉 只留铜 再用EAM 看会怎么样？
赞一下 回复此楼
I'maslowlywalker,butIneverwalkback.
9楼2017-12-05 09:30:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蔷薇妮妮

新虫 (初入文坛)
楼主你好，请问怎么解决的呢？遇到了同样的问题，想请教下，谢谢
赞一下 回复此楼
10楼2020-12-11 21:29:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cl1994861 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 总分293求调剂 +8 加一一九 2026-03-25 11/550 2026-03-29 19:53 by 无际的草原
[考研] 求调剂，一志愿 南京航空航天大学 ，080500材料科学与工程学硕，总分289分 +7 @taotao 2026-03-29 7/350 2026-03-29 12:03 by longlotian
[考研] 086000生物与医药调剂 +5 Feisty。 2026-03-28 9/450 2026-03-29 12:02 by longlotian
[考研] 085701环境工程，267求调剂 +16 minht 2026-03-26 16/800 2026-03-28 12:16 by zllcz
[考研] 286求调剂 +12 PolarBear11 2026-03-26 12/600 2026-03-28 12:14 by zllcz
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 291求调剂 +6 HanBeiNingZC 2026-03-24 6/300 2026-03-28 07:55 by baoball
[考研] 308求调剂 +7 墨墨漠 2026-03-27 7/350 2026-03-28 07:43 by 热情沙漠
[考研] 07化学280分求调剂 +10 722865 2026-03-23 10/500 2026-03-27 15:51 by Plutoqq
[考研] 286求调剂 +4 lim0922 2026-03-26 4/200 2026-03-27 10:28 by guoweigw
[考研] 359求调剂 +4 王了个楠 2026-03-25 4/200 2026-03-27 08:43 by 不吃魚的貓
[考研] 求调剂，一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +4 @taotao 2026-03-26 5/250 2026-03-27 08:10 by hypershenger
[考研] 304材料求调剂 +4 钟llll 2026-03-26 4/200 2026-03-27 03:42 by wxiongid
[考研] 342求调剂 +3 加油a李zs 2026-03-26 3/150 2026-03-27 00:29 by wxiongid
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +5 星沉uu 2026-03-26 6/300 2026-03-26 19:02 by macy2011
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 材料/农业专业，07/08开头均可，过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪！？！
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
信息提示
关闭
请填处理意见
关闭