24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1633  |  回复: 23
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yongning88

铜虫 (小有名气)

[求助] 势函数的构件 已有1人参与

我按照文献中的公式和参数计算了一下pair-potentiial, 算出的值和文献中给出的图一样(如图3所示),可为什么当 r 小于1.5 A 时,对势就很大(如图5所示),文献中没给出这部分的值,而LAMMPS中给出的势文件却是有个峰值的(如图4所示)。跪求大牛们的帮助。

势函数的构件
function.PNG


势函数的构件-1
parmeters.PNG


势函数的构件-2
pair_potential_in-paper.PNG


势函数的构件-3
pair_potential_for_LAMMPS.PNG


势函数的构件-4
pair_potenial_in_paper.PNG
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

yongning88

铜虫 (小有名气)

引用回帖:
7楼: Originally posted by 老虎大王 at 2014-05-05 09:13:09
原子离得很近的情况,比如r<1.0ang,那基本上是没有什么实际意义的,原子不会靠得这么近,如果你模拟中原子真要是靠得这么近,那一定是弄错了。这一部分的势函数不需要拟合,因为实际上用不到。只要在r 比较小的 ...

谢谢您的解释,既然r很小是可以忽略不用拟合,那为什么在有些已有的势函数文件中,对势值是从零开始增长,直到出现峰值,他这有事怎么处理的呢??
8楼2014-05-05 09:42:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答

lsloneil

专家顾问 (正式写手)

你是如何用lammps得到图4和图5的?
2楼2014-04-29 10:31:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongning88

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by lsloneil at 2014-04-29 10:31:08
你是如何用lammps得到图4和图5的?

这个不是用LAMMPS得到的,是用MATLAB 算的,我是想自己构架那个可以在LAMMPS中用的势函数
3楼2014-04-29 10:37:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)


月只蓝: 金币+1, 鼓励交流! 2014-04-30 13:22:50
mark一下。
楼主的那张图跟文献一样?我觉得从这个尺度看哪有一样的,或者楼主标清楚一下图几。。。
最好把文献贴出来,帖子里面只放自己的图。
4楼2014-04-30 09:37:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 2026硕士调剂_能动_河南农业大学 +3 河南农业大学-能 2026-04-12 3/150 2026-04-13 12:44 by 木鱼.
[考研] 求调剂 +7 何气正 2026-04-13 8/400 2026-04-13 09:33 by 零零二
[考研] 295分求调剂 +13 ?要上岸? 2026-04-10 13/650 2026-04-12 15:37 by laoshidan
[考研] 调剂 +10 月@163.com 2026-04-11 10/500 2026-04-12 09:14 by zhouyuwinner
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 0860004 求调剂 309分 +9 Yin DY 2026-04-08 9/450 2026-04-11 22:55 by dongdian1
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[考研] 283求调剂,工科! +12 苏打水7777 2026-04-08 12/600 2026-04-11 10:28 by 逆水乘风
[考研] 302分求调剂 +9 凡语祈愿 2026-04-08 10/500 2026-04-10 23:26 by 314126402
[考研] 0858求调剂 5+5 Gky09300550, 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 机械专368 有去处吗 +4 种大树 2026-04-10 4/200 2026-04-10 15:31 by jiajinhpu
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 08工学 309分求调剂 +6 Yin DY 2026-04-08 6/300 2026-04-10 09:18 by Delta2012
[考研] 070300化学 求调剂 +13 73372112 2026-04-08 13/650 2026-04-09 20:22 by maddjdld
[考研] 367求调剂 +10 hffQAQ 2026-04-09 10/500 2026-04-09 18:06 by lijunpoly
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 338求调剂 +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by 无际的草原
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
信息提示
请填处理意见