GULP不同版本导出ShengBTE用力常数文件，存在一些问题

GULP是基于经验势能的分子静力学计算软件，可以导出ShengBTE计算热导率用2阶和3阶力常数文件，相较于用LAMMPS做分动计算，可以先快速判断势能3阶非谐下的热导率理论值。

但是，大家要注意GULP不同版本里输出的力常数文件，有些存在问题：

1. GULP-5版本，导出的2阶力常数文件，存在原子扩胞后超胞内原子排序和phononpy等的不一致，从而导致热导率计算问题，这个在GULP-6版本中已经纠正；

2. 但是，GULP-6版本中导出的2阶力常数文件正确，3阶力常数文件存在错误，表现在三阶力常数文件中对超胞内相邻胞的位置x和y存在颠倒，反倒GULP-5版本中的是正确的。

3. 可以自己通过thirdorder.py来计算，这时候，gulp输出中力的单位为eV, 对正交晶格，需要有一个转换系数，-1*输出力（eV）/相应方向的晶格模量，需要自己编脚本做转换。

4.使用SW势能时，对于SW3部分，gulp输入和gulp输出中对三个原子的次序可能会有变化，需要检查，比如P2 P2 P1，在输出中认为的势能次序是P2 P1 P2,虽然SW3角度部分参数不变，但是对于角的两个边的参数会有问题。所以需要检查。