24小时热门版块排行榜

>论坛更新日志 (3512)
>文献求助 (342)
>导师招生 (263)
>虫友互识 (262)
>考博 (140)
>休闲灌水 (96)
>招聘信息布告栏 (83)
>硕博家园 (80)
>博后之家 (73)
>基金申请 (59)
>论文投稿 (56)
>论文道贺祈福 (48)
>考研 (48)
>公派出国 (42)
>找工作 (38)
>教师之家 (33)

返回列表

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

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

回复此楼

» 猜你喜欢

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

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的回帖

老虎大王

木虫 (著名写手)

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

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

引用回帖:

8楼: Originally posted by yongning88 at 2014-05-05 09:42:33
谢谢您的解释，既然r很小是可以忽略不用拟合，那为什么在有些已有的势函数文件中，对势值是从零开始增长，直到出现峰值，他这有事怎么处理的呢？？...

他既然提出了一个函数，这个函数的参数拟合出来之后，自然会有一个函数值随r变化的过程，r从零开始增长，不久出现峰值，这是这个函数的性质所决定的，是一个数学问题。而这个函数在r很小的时候的性质，是不影响它的使用的。一个势函数的可用性和势函数的好坏，不是看这个函数在r很小的时候的情况，而是看这个函数在原子最近邻距离的附近开始之后的情况，比如r＝1.0ang 或者1.5ang以后的情况。