关于Bi2Te3的势函数，文献中用的是GULP，我想把其改写为LAMMPS的。
势函数类型是morseelectrostatic
short-range interactions: morse potentiallong-range electrostatic interactions: Ewald summation methodwith uniform cutoff radius.
morse势参数如下：
电荷量是：0.38, −0.26, and −0.24 for Bi, Te1, Te2 atoms
我感觉文献中似乎少了
long-range electrostatic interactions的截断，我随便给了一个。
我个人编写的脚本如下，原子类型对应方式：
1-Bi
2-Te1
3-Te2pair_style hybrid/overlay morse 5.5 coul/long 12pair_coeff 1 1 morse 0.085 2.212 4.203 5.5
pair_coeff 1 2 morse 0.975 1.285 3.089 4.0
pair_coeff 1 3 morse 0.582 1.257 3.251 4.0
pair_coeff 2 2 morse 0.076 1.675 3.642 5.0
pair_coeff 2 3 morse 0.807 0.731 4.497 5.5
pair_coeff 3 3 morse 0.066 2.876 4.312 5.0
pair_coeff * * coul/long
mass 1 208.98
mass 2 127.60
mass 3 127.60
set type 1 charge 0.38
set type 2 charge -0.26
set type 3 charge -0.24
kspace_style ewald 1.0e-4复制代码但在使用lammps的例子ELASTIC计算弹性常数时，跟文献中计算的总不一样。
不知道哪里设置有问题，求指教呀！先谢谢啦！

楼主你好，请问下，你这样设置势模型参数跟直接应用lammps potentials 里面的势文件有什么区别呀？是不是你这样设置后，lammps利用这些数据产生产生势文件？

私聊