| 查看: 1821 | 回复: 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 |
回复此楼» 猜你喜欢
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
ChemiAndy
木虫 (正式写手)
- 模拟EPI: 40
- 应助: 52 (初中生)
- 金币: 2573.3
- 散金: 18
- 红花: 104
- 帖子: 590
- 在线: 553.8小时
- 虫号: 836907
- 注册: 2009-08-30
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★
月只蓝: 金币+3, 鼓励交流 2014-05-16 20:11:27
月只蓝: 金币+3, 鼓励交流 2014-05-16 20:11:27
|
貌似楼主的主要困惑是他得到的势能曲线在 r< 0.67 下与文献结果不一致。特别是，paper给出的曲线是无穷大，但LAMMPS得到的却是向下的。各位回答强调6<0.67A下没有意义，因为实际模拟当中原子间距不可能这么小，但是，仍然没有解决楼主的困惑。 我的意见是你的计算中肯定有哪个地方出错了，或者没有完全理解文献的意思。比如，你给出的图3和图5，势能最低点明显不一致。图5很可能只是用来作者用来示意势能函数中排斥部分的。如果加上吸引部分，总势能很有可能得到类似你图4的曲线，但我也不确信你的matlab计算是否完全follow了文献的公式。图2参数表给出的参数，有几个明显没有出现在图1的公式当中，你是否遗漏了什么？ |
11楼2014-05-12 08:11:56
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 (小学生)
- 金币: 2024.2
- 红花: 23
- 帖子: 516
- 在线: 659.8小时
- 虫号: 1069215
- 注册: 2010-08-04
- 专业: 金属结构材料
4楼2014-04-30 09:37:56
5