24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2589  |  回复: 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的回帖

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

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

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的回帖
相关版块跳转 我要订阅楼主 cl1994861 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 322求调剂 +7 123安康 2026-04-12 14/700 2026-04-15 21:48 by noqvsozv
[考研] 290调剂生物0860 +37 哇哈哈,。 2026-04-11 43/2150 2026-04-15 18:55 by zhuwenxu
[考博] 26年电池方向博士申请 10+3 15615939259 2026-04-09 4/200 2026-04-15 14:22 by kwansei
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +26 慕绝cc 2026-04-09 30/1500 2026-04-14 18:50 by 蔡苏阳
[考研] 一志愿中国科学院上海有机所,有机化学356分找调剂 +12 Nadiums 2026-04-09 13/650 2026-04-14 17:54 by lhj2009
[考研] 农学0904 312求调剂 +4 Say Never 2026-04-11 4/200 2026-04-14 09:10 by zs92450
[考研] 一志愿华南理工大学331分材料求调剂 +10 天下ww 2026-04-09 11/550 2026-04-13 23:25 by pies112
[考研] 一志愿双非085400电子信息344 求调剂,对材料和化学方向也感兴趣 +12 无情的小羊 2026-04-09 13/650 2026-04-13 14:17 by 张zhihao
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 材料考研调剂 +29 云木达达 2026-04-11 31/1550 2026-04-13 13:32 by lyh鲁老师
[考研] 生物学调剂 +11 小冉要努力 2026-04-10 13/650 2026-04-13 11:46 by 电化学及催化
[考研] 电气工程专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 10:47 by zhouyuwinner
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 求调剂 +11 翩翩一书生 2026-04-09 11/550 2026-04-11 19:57 by 逆水乘风
[考研] 263能源动力专硕求调剂 +3 加大号饭盒袋 2026-04-10 3/150 2026-04-10 22:23 by 286640313
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
[考研] 吉大计算机技术331分,英语六级,求调剂 +3 峰峰021116 2026-04-09 3/150 2026-04-10 20:01 by chemisry
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 314求调剂 +18 xhhdjdjsjks 2026-04-09 19/950 2026-04-10 18:53 by HPUCZ
[考研] 一志愿中科大070300化学,314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
信息提示
请填处理意见