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角度部分参数不变,但是对于角的两个边的参数会有问题。所以需要检查。