lammps问题求助

ywang

各位大神好，我有一个关于lammps eam势的问题，求各位大神指点我想做一个FS type的potential ，其计算能量的公式如下

公式里面并没有embeding function和电子密度的项，而是用binding和repulsion来描述能量。其中binding项文献里说是电子密度开根号，但是我自己算的时候发现并不是电子密度。lammps的说明书里面列的fs type列举的能量方程是用电子云密度和embeding function算的，如下图



这两种能量的表达方式应该如何转化？ EAM势文件需要对embeding function， 电子密度和相互作用列表，如何从FS势的表达式中获得embeding function和电子密度？
谢谢！

ywang

bluewhale 发表于 2016-1-27 14:00
你的理解是不是走入误区了，上面的公式是一目了然的事。

公式是一目了然，但是我按照文献里的公式做了tabulation以后用lammps跑了个NPT equilibration，结果报错跳出了，压强过大。然后我用公式直接列了势能和力的table用table style再算了一遍equilibration结果模型一直在膨胀，无法平衡。求指导啊

bluewhale

你的理解是不是走入误区了，上面的公式是一目了然的事。