24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yongning88

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 683.8
散金: 18
红花: 1
帖子: 81
在线: 40.6小时
虫号: 2140323
注册: 2012-11-21
专业: 金属基复合材料

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

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

function.PNG

parmeters.PNG

pair_potential_in-paper.PNG

pair_potential_for_LAMMPS.PNG

pair_potenial_in_paper.PNG

回复此楼

» 猜你喜欢

1楼 2014-04-29 09:51:37

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

老虎大王

木虫 (著名写手)

模拟EPI: 17
应助: 26 (小学生)
贵宾: 0.17
金币: 4774.1
散金: 8
红花: 42
帖子: 1361
在线: 215.2小时
虫号: 659094
注册: 2008-11-21
专业: 金属结构材料

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
月只蓝: 金币+10, 感谢指导！ 2014-05-16 20:11:55

引用回帖:

17楼: Originally posted by yongning88 at 2014-05-14 21:28:31
谢谢您的回复，其他的都明白了，可现在还是不明白的是，如果势函数不能取很大，也就是 r 的值不能太小，但构件势函数文件时，r得以零为起点，那对于这个部分的势函数值怎么得到？...

这个，我想，我已经说得比较清楚了。可是你好像还是不大清楚。现在我给你一个例子，比如说，我要对某种立方晶体拟合L-J势函数的参数： E=(A/r)^6-(B/r)^12，要拟合的就是A和B两个参数。我选择晶体的结合能Ec，晶胞体积V0，体弹性模量B，以及弹性系数C11和C12等来做为拟合目标。拟合之后：
1）我获得了A和B的值。现在我的势函数就完整了。
2）现在我做E-r图。做为一个数学函数，r 可以从零取到无穷大，都有对应的E值。
3）但是对于我的模拟工作来说，r很小的时候，虽然有对应的E值，但这些值在模拟上是用不到的。r很大的时候，同样如此。
4）我发表这个势函数的曲线图的时候，我不需要从r=0开始画图，并且我要避免从r=0开始画图。因为r很小的时候，势能非常大，比如在30000以上，这样我就需要把纵坐标范围定得很大，这就把后面实际上有用的部分掩盖掉了，后面的E值的波动在图上就反映不出来了。
5）所以我画图的时候，要从某个r值开始画，这个时候，E值还是比较大，但是能够保证我后面的E值的波动可以清楚地显示。
6）这就是你看到的文献的图上，r不从零画的原因。

总结一下，
我没有说势函数不能取很大，我是说两点：
A) 作图的时候，如果把r 很小的时候的势函数曲线作出来，纵坐标范围就要定得很大，后面的E值的波动在图上就反映不出来了，而只有后面的E值才对于模拟是真正有用的。
B） r很小的时候，模拟上用不到，所以你根本不要去纠结r很小时候的势能的情况，你计算势函数值的时候，要把精力放在后面的部分。后面的部分和文献一致了，你就算对了，后面的部分和文献不一致，那么说明你和文献作者之间至少有一方出错了。但是前面的r很小的部分的一致性是不需要考虑的。