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

回复此楼

» 猜你喜欢

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

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:10:35

引用回帖:

6楼: Originally posted by yongning88 at 2014-05-03 19:02:45
不好意思这几天放假，没及时回复，那对于原子很近的情况是怎么处理的呢？...

原子离得很近的情况，比如r<1.0ang，那基本上是没有什么实际意义的，原子不会靠得这么近，如果你模拟中原子真要是靠得这么近，那一定是弄错了。这一部分的势函数不需要拟合，因为实际上用不到。只要在r 比较小的的时候势能足够大，能够阻止原子继续靠近就可以了。你忽略掉这一部分就行了，不必纠结。主要是看后面的是否一致。