24小时热门版块排行榜    

查看: 1506  |  回复: 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

铜虫 (小有名气)

引用回帖:
21楼: Originally posted by 老虎大王 at 2014-05-15 14:29:48
我真的是要败给你了。我没有说不能很大!我是说作图的时候不要把它作出来!因为前面的数值太大了影响后面数值的显示。你看看你的图,后面都成一条水平线了。我希望你把你的数据和文献上取同样的r 范围,再作图比较 ...

好吧,您的意思我现在了解了,为了更好的反映后面数值的效果,在文献中就只能取这一部分的值来做图。当r值很小时,其势函数值对模拟的结果作用不是很大,所以这一部分的值可以完全是按照公式算出的值,也可以是自己给出的势函数值,这就是为什么有些现有的势函数文件中,当r=0时,phi=0 eV了,就像附件中给出的Al-Mg对势的结果一样;而不是严格按照公式趋向于无穷大,这样理解对吗?

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : AnisotropicsurfacesegregationinA1-Mgalloys.pdf
  • 2014-05-15 18:37:32, 1.09 M
  • 附件 2 : Al_Mg.eam.alloy
  • 2014-05-15 18:37:50, 83.82 K
22楼2014-05-15 18:40:04
已阅   回复此楼   关注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的回帖
信息提示
请填处理意见