VASP载流子迁移率计算：变形势理论与BoltzTraP方法解析

半导体材料的载流子迁移率是连接电子结构和器件性能的关键桥梁。计算迁移率的需求在二维材料、有机半导体和钙钛矿光伏研究中尤其集中——实验测量受样品质量影响大，理论计算反而能给出更”干净”的本征值。VASP提供了计算所需的电子结构数据，具体迁移率的提取则依赖不同的后处理框架。

变形势近似：适用于声学声子散射主导的体系

变形势（Deformation Potential，DP）理论是计算声学声子散射限制的载流子迁移率最常用的近似方法，适用于室温下非极性半导体中以声学声子为主要散射机制的体系。

迁移率的表达式为：

μ = e·ħ³·C / (k_B·T·m*²·E₁²)

其中C是弹性常数（描述材料对应变的刚度），m*是有效质量，E₁是变形势常数（能带边缘随晶格应变的变化率）。

VASP计算流程：

第一步：结构优化。充分弛豫到能量和力收敛，获得平衡态构型。

第二步：弹性常数计算。INCAR中设置IBRION=6，ISIF=3，VASP会自动计算弹性张量，从OUTCAR中提取所需的弹性常数C（对二维材料需要做厚度归一化）。

第三步：能带极值有效质量提取。通过高精度k点路径的能带计算，在导带底（CBM）或价带顶（VBM）附近拟合抛物线，得到有效质量m*。对于各向异性体系，需要分别提取沿不同晶向的有效质量。

第四步：变形势常数E₁。对晶格施加一系列小应变（通常±0.5%、±1.0%、±1.5%），分别计算每个应变态下的CBM或VBM能量位置，以应变量为横轴、能带边缘移动量为纵轴线性拟合，斜率即E₁。

这一步最容易犯的错误是应变范围选取不当。应变太小（<0.2%），数值噪声淹没信号；应变太大（>2%），进入非线性区，拟合结果不可靠。此外，每个应变态下的VASP计算需要使用同一套k点和截断能，确保比较的一致性。

BoltzTraP方法：适用于多带输运与温度依赖性分析

BoltzTraP（Boltzmann Transport Properties）是基于玻尔兹曼输运方程的后处理工具，能够从VASP的电子结构（能带、DOS）出发，计算塞贝克系数、电导率（含迁移率贡献）和热导率的电子部分，特别擅长处理多带输运和温度依赖行为。

VASP输入准备：BoltzTraP需要稠密k点网格的本征值数据，通常比普通能带计算密一倍以上（如普通用12×12×12，BoltzTraP输入用24×24×24或更密）。在INCAR中设置ISMEAR=0（高斯展宽），SIGMA=0.01 eV（小展宽保证精度）。计算完成后，BoltzTraP读取EIGENVAL、DOSCAR和结构信息文件。

弛豫时间近似：BoltzTraP本身不包含散射时间τ的计算，输出的是τ×σ（电导率×弛豫时间）。要得到绝对迁移率，需要引入τ的估计值——常见做法是通过变形势理论或实验测量值反推τ，或者直接计算τ×σ/τ×n得到μ=σ/(n·e)。对于比较不同材料或不同掺杂浓度的相对迁移率，τ相同的假设通常是合理的。

BoltzTraP2是更新的版本，支持Python接口，与pymatgen的集成使批量计算和数据处理更加方便。

两种方法的适用边界

对于二维半导体（如MoS₂、黑磷、MXene）的初步迁移率估算，变形势近似因其流程清晰、计算量小被广泛使用；对于热电材料需要全温区塞贝克系数和功率因子分析，BoltzTraP是更合适的工具。

科研学术网提供半导体电输运性质全套计算，包括载流子迁移率、塞贝克系数及热电优值ZT，欢迎访问 https://www.keyanxueshu.com 了解详情。