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

回复此楼

» 猜你喜欢

英文综述是否需要润色及查重已经有5人回复
救命帖已经有5人回复
限项规定已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
招博士已经有3人回复
存款400万可以在学校里躺平吗已经有35人回复
最失望的一年已经有18人回复
求推荐英文EI期刊已经有5人回复

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

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

yongning88

铜虫 (小有名气)

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

引用回帖:

11楼: Originally posted by ChemiAndy at 2014-05-12 08:11:56
貌似楼主的主要困惑是他得到的势能曲线在 r< 0.67 下与文献结果不一致。特别是，paper给出的曲线是无穷大，但LAMMPS得到的却是向下的。各位回答强调6<0.67A下没有意义，因为实际模拟当中原子间距不可能这么小 ...

谢谢您的回答，我只列出了对势求解的公式，对于电子密度和嵌入式能的公式没有给出，这就是为什么有的参数在公式上没有出现，
matlab计算没有问题，我特意取了个值手算了一下，结果是一样的。
我把文献帖出来，希望您能给点意见

Structural stability and lattice defects in copper: Ab initio, tight-binding,
and embedded-atom calculations

PDF文件怎么上传？？？