24小时热门版块排行榜

返回列表

shiling1813

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 38
红花: 1
帖子: 128
在线: 88.2小时
虫号: 725356
注册: 2009-03-18
专业: 金属材料的微观结构

[求助] 势函数数值列表时的嵌入能和�?势能的数值是根据什么表达式生成的？

问大神们个几个问题啊，困惑很久了，期待解答。
Quantum Sutton-Chen势函数公式生成lammps势函数数值列表时的嵌入能和�?势能的数值是根据什么表达式生成的？
另外，生成lammps势函数数值列表时，截断半径一般取多大？电子密度rho的取值间隔根据什么原则取，单位是什么？
我想根据Sutton-Chen势函数的解析表达式生成关于元素Ir和Rh的数值列表，用在lammps软件中。
根据Sutton-Chen势函数公式生成元素Ir和Rh的lammps势函数数值列表的过程是不是下面这样的？问题比较长，详细的内容写成word文件放在附件中了。

1.首先，Quantum Sutton-Chen势函数的表达式为（参考文献在附件中：Kimura Y, Qi Y, Agin T C, et al. The Quantum Sutton-Chen Many-Body Potential for Properties of fcc Metals[J]. Csdrm.caltech.edu, 1998, 8(2-3):151-159.）：

2根据lammps使用手册，数据在势函数数值列表文件中出现的顺序应该是下面这样的吧
一万个关于Ir的嵌入能F(rho)数据值（假如取一万个数据）
一万个关于IrIr的电子密度rho(r)数据值
一万个关于IrRh的电子密度rho(r)数据值
一万个关于Rh的嵌入能F(rho)数据值
一万个关于RhRh的电子密度rho(r)数据值
一万个关于RhIr的电子密度rho(r)数据值
一万个关于Ir-Ir的�?势能和距离的乘积phi(r)*r数据值
一万个关于Rh-Rh的�?势能和距离的乘积phi(r)*r数据值
一万个关于Ir-Rh的�?势能和距离的乘积phi(r)*r数据值

3.产生数值列表时的第j个原子对第i个原子的电子密度项表达式是不是下面这个表达式？

计算Ir原子对Ir原子的电子密度项时，（文献中已给）a=3.8344,m=6，从而得到一个列表（即第2个一万行电子密度项数值）；
计算Ir原子对Rh原子的电子密度项时，（文献中已给）a=3.8164,m=5.5，从而得到一个列表（即第3个一万行电子密度项数值）;

4.产生数值列表时的嵌入能表达式是下面这个表达式吗？是不是将一万个不同的距离r代入下面这个表达式，就可以得到一万个嵌入能F数据值？

那么怎么根据这个嵌入能列表中的数值得到实际的第i个原子的嵌入能呢？

5.产生数值列表时的第j个原子对第i个原子的对势能表达式是下面这个表达式吗？

那么怎么根据这个对势能列表中的数值得到实际的第i个原子的对势能呢？

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : The_Quantum_Sutton-Chen_Many-Body_Potential_for_Properties_of_fcc_Metals.pdf

2017-12-18 19:17:29, 242.93 K