| 查看: 1572 | 回复: 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 |
回复此楼» 猜你喜欢
323求调剂
已经有6人回复
-
一志愿北京化工大学 070300 学硕 336分 求调剂
已经有4人回复
-
352求调剂
已经有3人回复
-
一志愿东华大学化学070300，求调剂
已经有8人回复
-
277材料科学与工程080500求调剂
已经有7人回复
-
317求调剂
已经有18人回复
-
293求调剂
已经有5人回复
-
280分求调剂 一志愿085802
已经有7人回复
-
0854电子信息求调剂
已经有3人回复
-
263求调剂
已经有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人回复
老虎大王
木虫 (著名写手)
- 模拟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
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
5楼2014-04-30 10:20:36
yongning88
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 683.8
- 散金: 18
- 红花: 1
- 帖子: 81
- 在线: 40.6小时
- 虫号: 2140323
- 注册: 2012-11-21
- 专业: 金属基复合材料
6楼2014-05-03 19:02:45
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
7楼2014-05-05 09:13:09
yongning88
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 683.8
- 散金: 18
- 红花: 1
- 帖子: 81
- 在线: 40.6小时
- 虫号: 2140323
- 注册: 2012-11-21
- 专业: 金属基复合材料
8楼2014-05-05 09:42:33
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
9楼2014-05-06 09:36:17
yongning88
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 683.8
- 散金: 18
- 红花: 1
- 帖子: 81
- 在线: 40.6小时
- 虫号: 2140323
- 注册: 2012-11-21
- 专业: 金属基复合材料
10楼2014-05-06 11:14:56
40