使用BoltzTraP2计算Seebeck系数

详细安装步骤及教程可参考网站：

tutorial · Wiki · Jesús Carrete Montaña / BoltzTraP2 · GitLab

下面以硅烯为例讲解步骤

（1）使用vasp进行高精度的自洽计算，输入文件如下

POSCAR：

POSCAR\(1\1\1)
1.0
        3.8657999039         0.0000000000         0.0000000000
       -1.9328999519         3.3478809227         0.0000000000
        0.0000000000         0.0000000000        20.7891006470
   Si
    2
Direct
     0.000000000         0.000000000         0.037960000
     0.333330005         0.666670024         0.000000000

KPOINTS：

Automatic mesh
0
Gamma
101 101 1
0. 0. 0.

INCAR：

PREC = High
IBRION = 2
NSW = 0
ISMEAR = 0
SIGMA = 0.01
LREAL = .FALSE.
ADDGRID = .TRUE.
LWAVE = .FALSE.
LCHARG = .TRUE.
LOPTICS = .TRUE.
LVEL=.TRUE.

（2）输入命令，产生interpolation.bt2文件，后面的计算都会基于这个文件，其中/public/home/xiaqian/AlSb_2__NodleLine/seebeck/2dSi.vasp是自洽计算的路径

btp2 -vv interpolate -m 3 /public/home/xiaqian/AlSb_2__NodleLine/seebeck/2dSi.vasp

（3）自洽算完之后，可以使用BoltzTraP2拟合能带，看是否能和VASP算的能带匹配，输入以下命令：

btp2 plotbands interpolation.bt2 ["[-0.50000,   0.5000,    0.00000], [0.00000 ,  0.00000 ,   0.00000], [-0.33333 ,  0.66667  ,  0.00000], [-0.50000 ,  0.50000 ,   0.00000]"]

输入完之后会产生能带图，如下所示：

注意要使用具有可视化功能的集群连接软件，不然无法显示图片，当然也可以在自己电脑上用python运行。

（4）分别输入以下命令：

btp2 -vv integrate interpolation.bt2 300:500:1
btp2 plot -u -c  '["xx","yy","xy", "zz"]' -s 50 interpolation.btj S

最终产生eebeck关于能量的图片，如下所示：