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年的指南还没有出来？已经有4人回复
纳米粒子粒径的测量已经有8人回复
疑惑？已经有5人回复
国自然申请面上模板最新2026版出了吗？已经有14人回复
计算机、0854电子信息（085401-058412）调剂已经有5人回复
Materials Today Chemistry审稿周期已经有5人回复
溴的反应液脱色已经有7人回复
推荐一本书已经有12人回复
基金申报已经有4人回复
常年博士招收(双一流，工科) 已经有4人回复

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

Si原子，O原子，C原子，H原子的L-J势函数的参数已经有6人回复
非能量最小点的频率计算(任意结构的频率计算) 已经有14人回复
想用GULP拟合纯钼的morse势函数参数已经有13人回复
Mg金屬結構分子模擬已经有5人回复
关于势函数选取的问题。。。已经有9人回复
分子模拟的关于势函数以及初始构型的几个小问题已经有15人回复
求助：哪里能找到各种物质的勒纳德琼斯势函数已经有5人回复
对于MD，其中的势函数参数选取，如何自己计算已经有5人回复
如何在lammps里加入新的势函数形式？已经有6人回复
Gulp 拟合势函数初值设置的问题已经有12人回复
国内有开发分子动力学势函数的课题组吗？已经有12人回复
关于势函数的求助已经有23人回复
日本开发出在磁场中变硬的橡胶新材料已经有4人回复
Si的S-W势函数已经有3人回复
【求助】优化的结构非常奇怪，Au10 簇被优化成了4组小簇，为什么，输入及结构已贴出已经有12人回复
【求助】ms中如何假如势函数，比如EAM，比如自己修改过的势函数。已经有12人回复
【讨论】铜和水之间的LJ势函数的参数是多少？已经有19人回复

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

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

yongning88

铜虫 (小有名气)

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

引用回帖:

9楼: Originally posted by 老虎大王 at 2014-05-06 09:36:17
他既然提出了一个函数，这个函数的参数拟合出来之后，自然会有一个函数值随r变化的过程，r从零开始增长，不久出现峰值，这是这个函数的性质所决定的，是一个数学问题。而这个函数在r很小的时候的性质，是不影响它的 ...

对于我问题中给出的对势公式和拟合出来的参数值，我用MATLAB从r=0开始算了下，如图5所示，这个关于r的函数就没有峰值，也就是从无穷大趋向于零。如果我要用这个公式来建立一个势函数文件，需要从r=0，dr=0.0134A，截断半径取为6.68A，也就是500个对势值。选取r=0.67A时，对势值为910eV，这个为最小的原子间距，也就是最大的对势值，那对于r<0.67A的对势值怎么选取？