| 查看: 2242 | 回复: 5 | ||
| 【悬赏金币】回答本帖问题，作者lasergaga将赠送您 150 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
lasergaga新虫 (初入文坛)
|
[求助]
都进来看看啊，求教Interatomic Potentials Database参数如何引用，我快走投无路了 已有1人参与
|
|
 如题，想在https://apot.mgedata.cn/Interato ... E4CA9E390754E3DC2F7这里引用Al/MgO的势参数，却看不懂这是用的什么势函数，读论文说是一种修正的对势，也就是在morse的基础上有改进，他给的参数也和morse不一样，所以即使有参数，我也不清楚在lammps中该用哪种pair_style，pair_coeff又该如何写？请各位大神帮帮忙，卡这里半个多月没进展了，被老板劈死了   (AlMgO(110)).png  (AlMgO(110)).png  微信图片_20220712184227.png |
回复此楼
田淮安
新虫 (初入文坛)
- 应助: 4 (幼儿园)
- 金币: 179.1
- 帖子: 23
- 在线: 15.6小时
- 虫号: 2895705
- 注册: 2013-12-26
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2022-08-07 22:11:05
lasergaga
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 56.1
- 帖子: 14
- 在线: 5.4小时
- 虫号: 5111338
- 注册: 2016-10-12
- 性别: MM
- 专业: 金属基复合材料
2楼2022-07-21 00:40:19
田淮安
新虫 (初入文坛)
- 应助: 4 (幼儿园)
- 金币: 179.1
- 帖子: 23
- 在线: 15.6小时
- 虫号: 2895705
- 注册: 2013-12-26
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
【答案】应助回帖
| 你好，你用的势能公式叫做Rahman–Stillinger– Lemberg (RSL2) potential，这一名称就写在你给的文献的Eq(19)附近。我查了一下，这个势能在LAMMPS中好像没有专门的pair_style。不过因为它是pairwise potential，所以很容易用pair_style table实现；也就是说，你可以直接算出这个公式在一系列离散R值下的势能U和力F，然后写成一个table让LAMMPS读入即可。这个势能的形式比较简单，F的公式直接对U的公式求导就能得到。另外，我发现比较近的一篇文章用LAMMPS实现了这个势能，见 Atomistic simulation study on the shear behavior of Ag/MgO interface, Computational Materials Science 155 (2018) 116–128。你可以直接发邮件给这篇文章的作者，问他们是怎么在LAMMPS中实现这个势能的。作者是中科院和北大的老师，应该很容易获得帮助。祝你顺利。 |
3楼2022-08-04 05:24:24
lasergaga
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 56.1
- 帖子: 14
- 在线: 5.4小时
- 虫号: 5111338
- 注册: 2016-10-12
- 性别: MM
- 专业: 金属基复合材料
4楼2022-08-06 13:32:10
