二维材料输运计算为什么参数多？

引言

二维材料输运计算最容易让人困惑的地方，是公式里同时出现有效质量、形变势常数、弹性常数和温度项。只看能带很轻，不能保证迁移率一定高；只看材料很硬，也不能说明载流子散射弱。迁移率本质上是在问载流子加速、声子散射和晶格响应三件事怎样同时发生。

对 VASP 结果来说，这三个量来自不同后处理：有效质量来自带边曲率，形变势来自小应变下 VBM 或 CBM 的移动，二维弹性常数来自应变能量曲线。它们任何一个取错方向、载流子类型或应变范围，最后的 μ 都会偏离真实比较对象。

图1 二维半导体高迁移率筛选框架及关键因子贡献。迁移率不只由带隙决定，还受到有效质量、结构特征和局域成键环境共同影响。DOI: 10.1038/s41467-024-49686-z。

迁移率公式里到底有哪几个入口？

常用二维形变势理论把声学声子散射写成一个可计算近似。对某一输运方向，迁移率通常写成下面这种形式，其中每个量都对应一次独立计算，而不是从同一张能带图里直接读完。

μ_2D = eℏ³C_2D / (k_BT m* m_d E₁²)　式(1)

式(1)中，C₂D 是二维弹性常数，m* 是输运方向的有效质量，m_d 是态密度有效质量，E₁ 是带边相对真空能级或平均静电势随应变变化的斜率。温度 T 通常取目标工作温度或 300 K 参考值。

有效质量决定载流子在外场下是否容易加速，形变势决定带边对应变有多敏感，弹性常数决定晶格涨落的能量代价。因此迁移率不是一个单参数排序；三者相乘后，某个看起来优势明显的材料可能被另一个参数抵消。

m_d 也不能忽略。二维各向异性体系常把它写成 m_d = (m_xm_y)^1/2，这意味着 x 方向迁移率仍会受到 y 方向带边曲率影响。若只把输运方向的 m* 放进公式，而没有处理垂直方向的态密度有效质量，结果会系统性偏高或偏低。

有效质量为什么先限制速度响应？

有效质量来自带边附近 E-k 曲线的二阶导数。带边越平，二阶导数越小，对应的 m* 越大，载流子在同样外场下获得的速度响应越弱；带边越弯，m* 越小，输运方向上更容易形成高迁移率。

二维材料还要区分 x 与 y 方向，因为许多单层结构不是各向同性。黑磷、Janus 结构、矩形晶格和低对称单层都可能出现一个方向轻、另一个方向重的情况。此时不能把 Γ-X、Γ-Y 或 K-M 的曲率混成一个平均值。

图2 二维 tetrahex-GeC₂ 的能带、带边电荷密度和方向相关有效质量。能带曲率越明显，载流子对外场的响应越强。DOI: 10.1021/acsami.0c23017。

图2中的方向相关有效质量说明，同一种材料里电子和空穴也可以给出不同趋势。计算迁移率时必须先确定载流子类型、带边位置和输运方向；否则后面的形变势和弹性常数即使算得再细，也是在为错误对象服务。

拟合有效质量时还要控制 k 点窗口。窗口太宽会把非抛物线区间纳入拟合，窗口太窄又容易受数值噪声影响。对带边附近存在多个近简并谷的二维半导体，还应分别确认真正的 VBM、CBM 和次级谷位置，避免把局域曲率当成整体输运性质。

形变势为什么比带隙变化更敏感？

形变势常数 E₁ 描述的是带边能级对晶格微小拉伸或压缩的响应，而不是简单看带隙变大还是变小。实际计算时，应在小应变区间内分别施加 +ε 和 -ε，读取 VBM 或 CBM 相对稳定参照能级的位置，再做线性拟合。

这里的参照很关键。二维 slab 有真空层，绝对 Kohn-Sham 能级会随计算设置漂移；更可靠的做法是使用真空平台或平面平均静电势作为能量零点，再比较带边移动。图3展示的就是带边能量随单轴应变的线性变化。

图3 ZnSbX 单层在 armchair 与 zigzag 方向应变下的带边移动。带边能量对应变的线性斜率可用于提取形变势常数。DOI: 10.1039/d2ra04782a。

在式(1)里，E₁ 以平方形式出现在分母。形变势稍微增大，就会显著压低迁移率。这也是为什么只看有效质量常常会误判：一个带边很轻的材料，如果带边对应变极其敏感，声学声子散射仍然会很强。

应变幅度建议控制在小形变区间，例如 ±0.5%、±1.0% 或 ±1.5%，并检查带边移动是否保持近似线性。若大应变已经改变带边位置、诱导能谷切换或改变直接/间接带隙性质，那么得到的斜率就不再代表声学声子微扰下的形变势。

弹性常数为什么要和形变势一起看？

二维弹性常数 C₂D 来自应变能量曲线，常用能量相对平衡面积的二阶拟合得到。它反映材料在该方向上抵抗拉伸或压缩的能力。晶格越容易被热涨落扰动，声学声子散射背景通常越强。

但 C₂D 不能单独解释迁移率。较大的 C₂D 会提高式(1)中的分子项，看起来有利于迁移率；如果同一方向上的 E₁ 也很大，带边会对微小形变高度敏感，最终结果仍可能降低。图4里的弹性能曲线正适合用来理解这种力学入口。

图4 二维 tetrahex-GeC₂ 在不同应变方式下的弹性能和面内力学响应。弹性常数描述晶格承受形变时的能量代价。DOI: 10.1021/acsami.0c23017。

弹性常数和形变势必须来自同一个晶格方向、同一应变范围和同一结构弛豫口径。例如 x 方向输运不能随手搭配 y 方向弹性常数；电子迁移率也不能直接套用空穴带边的形变势。

二维 C₂D 的单位常写成 N m⁻¹，也可由 eV Å⁻² 转换得到。计算时需要用面内面积归一化，而不是把真空层厚度当作真实体积参与三维弹性常数换算；否则真空层越厚，表观弹性量越小。

二维输运为什么必须拆成多组结果？

一套完整迁移率数据至少要拆成电子 x、电子 y、空穴 x、空穴 y 四组。低对称结构还可能需要沿任意角度拟合有效质量和弹性响应。若只给一个平均 μ，容易掩盖真实的方向选择性。

图5给出了二维 III-V 材料中迁移率及形变势贡献的差异。它提醒我们，迁移率高不是因为某个变量单独漂亮，而是 m*、m_d、E₁、C₂D 和温度条件共同组合后仍保持优势。

图5 二维 III-V 材料的载流子迁移率及形变势贡献。不同方向、不同载流子类型会给出明显不同的迁移率结果。DOI: 10.1126/sciadv.adu5294。

落到 VASP 任务安排上，建议先做高精度结构优化，再做带边曲率拟合、单轴小应变序列和应变能拟合，最后把所有量按方向合并。迁移率表格的核心不是列越多越好，而是每一列都能回到同一个输运方向和同一个带边。

如果要比较不同材料，计算精度也要保持一致。ENCUT、k 点密度、泛函、真空层、偶极修正和结构弛豫标准都会影响带边曲率、应变能和能级参照。迁移率数值通常跨度很大，不能把不同精度下的小数点差异当作可靠排序。

本章要点总结

要点：二维迁移率由有效质量、态密度有效质量、形变势常数、二维弹性常数和温度共同决定。

要点：有效质量来自带边曲率，必须区分电子、空穴和输运方向。

要点：形变势来自带边随小应变的线性移动，能量参照应保持一致。

要点：弹性常数来自应变能量曲线，需要和形变势使用同一方向。

要点：二维材料输运结果应按载流子类型和方向分组，不能只给一个平均值。