注:对计算结果的准确性不做任何保证。对于参数的设置没怎么验证,可能会有些错误。
主要使用ShengBTE、thirdorder(两者在ShengBTE网站下载) 和phonopy软件。
ShengBTE计算热导率前所需要准备的文件
FORCE_CONSTANTS_2ND #见二阶力常数计算
FORCE_CONSTANTS_3RD #见三阶力常数计算
注:基于玻尔兹曼输运方程计算晶格热导率的软件一般都需要计算二阶力和三阶力常数。
二阶力常数计算
计算流程与声子谱计算流程相同,用DFPT和有限位移法都可以。将得到的FORCE_CONSTANTS重命名为FORCE_CONSTANTS_2ND
注:对于二阶力常数计算采用的超胞,个人觉得应该起码保证结果是晶格动力学稳定的(声子谱基本没有虚频)。
三阶力常数计算
必要的文件:
create.sh #生成计算三阶力常数的POSCAR和文件夹
##thirdorder里面的文件复制到自己的文件夹中。
thirdorder_common.py
thirdorder_vasp.py #基于vasp
run_3rd.sh #运行计算并提取三阶力常数
POTCAR #赝势
INCAR KPOINTS POSCAR可单独给出,也可在creat.sh中设置
create.sh脚本设置如下:
#!/bin/bash
#优化后