24小时热门版块排行榜    

查看: 2478  |  回复: 9

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
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
pair_coeff        * * eam Cu_u3.eam NULL NULL Cu

» 本帖已获得的红花(最新10朵)

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的回帖
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 的主题更新
信息提示
请填处理意见