| 查看: 1608 | 回复: 23 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
yongning88铜虫 (小有名气)
|
[求助]
势函数的构件 已有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 |
回复此楼» 猜你喜欢
调剂310
已经有11人回复
-
300求调剂，材料科学英一数二
已经有9人回复
-
296求调剂
已经有3人回复
-
一志愿南昌大学324求调剂
已经有3人回复
-
299求调剂
已经有9人回复
-
22408 359分调剂
已经有6人回复
-
085600 材料与化工 329分求调剂
已经有12人回复
-
一志愿北京化工大学材料与化工（085600）296求调剂
已经有14人回复
-
085601材料工程找调剂
已经有11人回复
-
总分293求调剂
已经有9人回复
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
月只蓝: 金币+10, 感谢指导！ 2014-05-16 20:11:55
月只蓝: 金币+10, 感谢指导！ 2014-05-16 20:11:55
|
这个，我想，我已经说得比较清楚了。可是你好像还是不大清楚。现在我给你一个例子，比如说，我要对某种立方晶体拟合L-J势函数的参数： E=(A/r)^6-(B/r)^12，要拟合的就是A和B两个参数。我选择晶体的结合能Ec，晶胞体积V0， 体弹性模量B，以及弹性系数C11和C12等来做为拟合目标。拟合之后： 1）我获得了A和B的值。现在我的势函数就完整了。 2）现在我做E-r图。做为一个数学函数，r 可以从零取到无穷大，都有对应的E值。 3）但是对于我的模拟工作来说，r很小的时候，虽然有对应的E值，但这些值在模拟上是用不到的。r很大的时候，同样如此。 4）我发表这个势函数的曲线图的时候，我不需要从r=0开始画图，并且我要避免从r=0开始画图。因为r很小的时候，势能非常大，比如在30000以上，这样我就需要把纵坐标范围定得很大，这就把后面实际上有用的部分掩盖掉了，后面的E值的波动在图上就反映不出来了。 5）所以我画图的时候，要从某个r值开始画，这个时候，E值还是比较大，但是能够保证我后面的E值的波动可以清楚地显示。 6）这就是你看到的文献的图上，r不从零画的原因。 总结一下， 我没有说势函数不能取很大，我是说两点： A) 作图的时候，如果把r 很小的时候的势函数曲线作出来，纵坐标范围就要定得很大，后面的E值的波动在图上就反映不出来了，而只有后面的E值才对于模拟是真正有用的。 B） r很小的时候，模拟上用不到，所以你根本不要去纠结r很小时候的势能的情况，你计算势函数值的时候，要把精力放在后面的部分。后面的部分和文献一致了，你就算对了，后面的部分和文献不一致，那么说明你和文献作者之间至少有一方出错了。但是前面的r很小的部分的一致性是不需要考虑的。 |
18楼2014-05-15 08:25:12
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2014-04-29 10:31:08
yongning88
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 683.8
- 散金: 18
- 红花: 1
- 帖子: 81
- 在线: 40.6小时
- 虫号: 2140323
- 注册: 2012-11-21
- 专业: 金属基复合材料
3楼2014-04-29 10:37:08
fanarsenal
木虫 (正式写手)
- 应助: 47 (小学生)
- 金币: 2023.2
- 红花: 23
- 帖子: 515
- 在线: 659.5小时
- 虫号: 1069215
- 注册: 2010-08-04
- 专业: 金属结构材料
4楼2014-04-30 09:37:56
10