计算化学公社

标题: lammps问题求助 [打印本页]

作者
Author: ywang 时间: 2016-1-25 16:31
标题: lammps问题求助
本帖最后由 ywang 于 2016-1-25 16:34 编辑

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

(, 下载次数 Times of downloads: 43)

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

(, 下载次数 Times of downloads: 32)

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

作者
Author: bluewhale 时间: 2016-1-27 14:00
你的理解是不是走入误区了，上面的公式是一目了然的事。

作者
Author: ywang 时间: 2016-2-9 07:14

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)