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

回复此楼

» 猜你喜欢

西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复
救命帖已经有6人回复
存款400万可以在学校里躺平吗已经有37人回复
英文综述是否需要润色及查重已经有5人回复
限项规定已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
招博士已经有3人回复
最失望的一年已经有18人回复

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

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

ChemiAndy

木虫 (正式写手)

模拟EPI: 40
应助: 52 (初中生)
金币: 2573.3
散金: 18
红花: 104
帖子: 590
在线: 553.8小时
虫号: 836907
注册: 2009-08-30
性别: GG
专业: 理论和计算化学

★ ★ ★
月只蓝: 金币+3, 鼓励交流 2014-05-16 20:11:27

引用回帖:

10楼: Originally posted by yongning88 at 2014-05-05 15:14:56
对于我问题中给出的对势公式和拟合出来的参数值，我用MATLAB从r=0开始算了下，如图5所示，这个关于r的函数就没有峰值，也就是从无穷大趋向于零。如果我要用这个公式来建立一个势函数文件，需要从r=0，dr=0.0134A， ...

貌似楼主的主要困惑是他得到的势能曲线在 r< 0.67 下与文献结果不一致。特别是，paper给出的曲线是无穷大，但LAMMPS得到的却是向下的。各位回答强调6<0.67A下没有意义，因为实际模拟当中原子间距不可能这么小，但是，仍然没有解决楼主的困惑。

我的意见是你的计算中肯定有哪个地方出错了，或者没有完全理解文献的意思。比如，你给出的图3和图5，势能最低点明显不一致。图5很可能只是用来作者用来示意势能函数中排斥部分的。如果加上吸引部分，总势能很有可能得到类似你图4的曲线，但我也不确信你的matlab计算是否完全follow了文献的公式。图2参数表给出的参数，有几个明显没有出现在图1的公式当中，你是否遗漏了什么？