手机版

VASP 载流子迁移率计算

发布时间：2026-08-31 来源：科研学术网

小中大

字号：

一、项目背景与科学问题

载流子迁移率（μ）是半导体材料的核心输运参数，它决定晶体管开关速度、太阳能电池载流子收集效率与发光器件的注入特性。对二维材料（如单层 MoS₂），载流子迁移率不仅影响器件性能，还与其柔性电子、光电探测等应用前景直接相关。实验上测量二维材料迁移率受接触电阻、衬底散射等因素干扰，而第一性原理计算可从原子层面获得本征迁移率的可靠估计。

本项目以单层 MoS₂ 为对象，利用 VASP 计算其电子结构（有效质量）与形变势，结合形变势理论（deformation potential theory）估算室温下载流子迁移率，分析电子与空穴的输运差异。

二、计算方法与理论背景

2.1 形变势理论（Bardeen-Shockley）

在长波声学声子散射主导的输运机制下，载流子迁移率可表示为：

μ = (2√2π e ħ³ C₂D) / [3(k_B T)^{3/2} (m*)⁵ᐟ² E₁²]

其中 C₂D 为二维弹性模量（面内刚度），m* 为载流子有效质量，E₁ 为形变势常数（能带边相对应变的变化率），T 为温度。该公式显示：迁移率对有效质量（m* 的 2.5 次方）与形变势（平方）极为敏感，轻有效质量与弱电子-声子耦合的材料具有高迁移率。

2.2 有效质量

有效质量由能带边附近的色散曲率决定：m* = ħ²(d²E/dk²)⁻¹。对单层 MoS₂，导带底（CBM）与价带顶（VBM）均位于 K 点（直接带隙半导体），需沿 K→Γ 与 K→M 两个方向拟合抛物线提取有效质量，取平均获得各向同性估计。

2.3 形变势与弹性模量

形变势 E₁ = ΔE_edge/Δε 为单轴应变下能带边能量位移的斜率；二维弹性模量 C₂D 由总能-应变关系二次拟合得到，反映材料抵抗面内拉伸的刚度。两者均由 VASP 对应变结构的能量计算获得。

三、计算设置与模型构建

软件：VASP；泛函：GGA-PBE（含 DFT-D3 色散修正）。
模型：单层 MoS₂ 原胞，真空层 18 Å 抑制层间作用。
k 点：高密度路径采样提取能带，K 点附近加密（精细网格）。
应变：在 −2% 至 +2% 范围内以 0.5% 步长施加双轴应变，计算能带边位移与总能。
温度：300 K（载流子迁移率计算）。

四、结果分析与案例图解读

图 3 的左面板给出单层 MoS₂ 在 K 点附近的能带结构示意。蓝色曲线为导带（CBM），红色曲线为价带（VBM），两带在 K 点（横轴中心）形成直接带隙 Eg = 1.7 eV。CBM 的抛物线开口较宽、曲率小，对应较轻的电子有效质量 m_e* = 0.48 m₀；VBM 抛物线曲率稍大，对应空穴有效质量 m_h* = 0.55 m₀。由带隙与带边有效质量可直接读取：单层 MoS₂ 为直接带隙半导体，电子与空穴有效质量相近但电子略轻，预示电子迁移率略高于空穴。该能带结构为迁移率计算提供了核心输入参数。

图 3 的右面板给出形变势拟合与迁移率结果。散点为 CBM/VBM 能量随双轴应变的位移，线性拟合斜率即为形变势常数：CBM 形变势 E₁ᵉ = −8.2 eV（蓝色，斜率陡峭），VBM 形变势 E₁ʰ = −2.1 eV（红色，斜率平缓）。结合面内弹性模量 C₂D = 128 N/m 与上述有效质量，代入形变势公式计算得室温（300 K）电子迁移率 μ_e = 210 cm²/(V·s)、空穴迁移率 μ_h = 165 cm²/(V·s)（图中文本框标注）。电子迁移率高于空穴约 27%，主要源于电子更轻的有效质量（尽管其形变势更大）。计算值与实验报告的单层 MoS₂ 场效应迁移率（约 100–200 cm²/(V·s)）处于同一量级，说明本征声子散射机制正确描述了室温输运。

两面板结合，左面板提供能带结构与有效质量，右面板提供形变势与迁移率结果，共同完成”从电子结构到宏观输运”的完整计算链条，为二维半导体器件设计提供定量迁移率数据。