📚 课程目录第3章宽禁带材料晶格动力学

3. 宽禁带材料晶格动力学

宽禁带半导体（SiC、GaN、AlN、β-Ga₂O₃等）的热输运行为由晶格振动主导，而晶格振动的完整描述需要四个相互关联的物理量：声子色散关系 ω(q)、声子态密度 g(ω)、声子群速度 v_g(q)=∂ω/∂q，以及表征非谐性的 Grüneisen 参数 γ。四者并非并列的独立指标，而是构成一条从"力常数→色散→态密度/群速度→非谐散射强度"的因果链：色散是晶格动力学矩阵对角化的直接输出，态密度是色散在布里渊区的积分，群速度是色散的一阶导数，Grüneisen 参数则来自声子频率对体积的响应。理解这条链条，才能判断第一性原理计算结果与实验测量为何在有些材料上吻合良好、在另一些材料上出现系统性偏差。

一、第一性原理计算流程与关键参数

当前主流做法是密度泛函理论（DFT）结合密度泛函微扰理论（DFPT）或有限位移法。以 4H-SiC 为例，典型流程如下：

# 以 Quantum ESPRESSO 为例的声子计算关键输入（示例参数）
&control
  calculation = 'scf'
/
&system
  ibrav = 4, celldm(1) = 6.09, celldm(3) = 1.63   # 4H-SiC 六方胞近似
  nat = 8, ntyp = 2
/
&electrons
  ecutwfc = 60, ecutrho = 480
/
K_POINTS automatic
  8 8 4 0 0 0

# 随后进行 ph.x 声子计算，q 点网格取 4×4×2
# 再经 q2r.x 与 matdyn.x 得到完整色散

关键收敛参数包括：截断能（SiC 通常 60–80 Ry）、k 点网格（需保证力常数收敛到 1 cm⁻¹ 量级）、q 点网格（决定傅里叶插值精度）。对宽禁带材料，交换关联泛函的选择会显著影响声子频率：LDA 通常高估频率约 3%–5%，PBE 低估约 2%–4%，而杂化泛函 HSE06 可将误差压缩到 1% 以内，但计算成本增加一个数量级。因此，若目标是定量对比实验，建议至少用 PBE 加经验修正，或对关键声子支采用 HSE06 验证。

二、声子色散与态密度的解读

4H-SiC 的声子色散沿 Γ–M–K–Γ–A 路径呈现明显的各向异性：光学支在 Γ 点分裂为横声学（TA）、纵声学（LA）、横光学（TO）和纵光学（LO）四类，其中 LO–TO 分裂由长程库仑相互作用引起，在极性宽禁带材料中可达 100–200 cm⁻¹。态密度 g(ω) 则把色散在布里渊区积分，其峰位对应声子支的平坦区域——平坦支意味着群速度小、态密度高，对热导率贡献反而低。这一点常被误解：态密度高不等于热导率高，必须结合群速度一起看。

三、群速度与 Grüneisen 参数的物理意义

群速度直接决定声子输运效率。在 4H-SiC 中，声学支在 Γ 点附近 v_g 可达 8000–10000 m/s，而光学支通常只有 1000–3000 m/s。这意味着热导率主要由声学支贡献，光学支更多充当散射通道。计算 v_g 时需注意：色散曲线在 q 点网格上的数值微分会引入误差，建议用傅里叶插值后的密集网格（如 20×20×10）求导，或直接由动力学矩阵的解析导数获得。

Grüneisen 参数 γ 定义为 γ_i = −∂lnω_i/∂lnV，表征声子频率对体积变化的敏感度。γ 越大，非谐性越强，声子-声子散射（Umklapp 过程）越剧烈，热导率越低。对宽禁带材料，γ 的典型值在 0.5–1.5 之间：SiC 约 0.8–1.0，GaN 约 0.7–1.2，AlN 约 0.6–0.9。计算 γ 有两种途径：一是准谐近似（QHA），在不同体积下分别计算声子频率再拟合；二是 DFPT 直接计算 ∂ω/∂V。QHA 实现简单但忽略高阶非谐，DFPT 更准确但需修改代码。实际工作中，常用 QHA 得到 γ 的温度依赖，再与实验热膨胀系数对比验证。

四、与实验对比的判据与常见偏差

实验上，声子色散主要通过非弹性中子散射（INS）或非弹性 X 射线散射（IXS）测量，态密度可由比热、拉曼/红外光谱间接推断，群速度可由热导率拟合或皮秒声学测量获得，Grüneisen 参数则来自热膨胀系数与比热的关系：γ = αV_B/(C_V·κ_T)。对比时需注意以下判据：

色散对比：光学支频率误差应小于 3%，声学支在 Γ 点附近应趋近于零（声学支求和规则）。若声学支在 Γ 点出现明显非零截距，说明力常数收敛不足或存在虚频。
态密度对比：比热 C_V 在低温区（T < 100 K）对 g(ω) 的低频部分敏感，若计算 C_V 低于实验值，通常意味着低频声学支被低估。
群速度对比：由热导率反推的 v_g 是平均值，与第一性原理的支分辨 v_g 对比时需按态密度加权。
Grüneisen 对比：实验 γ 来自热膨胀与比热，是各声子支的加权平均；计算 γ 应做同样加权后再比较，否则会因光学支 γ 偏大而高估整体非谐性。

常见偏差来源：（1）DFT 泛函选择导致频率整体偏移；（2）q 点网格不足使色散插值失真；（3）忽略 LO-TO 分裂的极性修正；（4）QHA 在高温区失效，γ 被高估；（5）实验样品中的同位素无序和缺陷展宽使谱峰变宽，与理想计算谱直接对比会误判。建议在对比前先做频率缩放因子校准，并明确标注计算与实验的温度条件。

五、实操建议与结果判读

若以 4H-SiC 为对象开展一次完整的计算-对比练习，可按以下顺序推进：先做 SCF 收敛测试，确定 ecutwfc 和 k 网格；再用 DFPT 在 4×4×2 q 网格上计算力常数；经傅里叶插值得到沿 Γ–M–K–Γ–A 的色散；由色散积分得 g(ω)，数值微分得 v_g(q)；最后在 0.98V₀、V₀、1.02V₀ 三个体积下重复声子计算，拟合 γ。结果判读时，优先检查声学支求和规则和 Γ 点光学支频率，这两项通过后再看整体色散和态密度。若与实验偏差超过 5%，应优先排查泛函和 q 网格，而非直接调整经验参数。

需要强调的是，第一性原理计算与实验对比的目的不是追求完全重合，而是建立可解释的偏差来源。宽禁带材料的强极性、各向异性和高温非谐性使得任何单一近似都有局限，只有把色散、态密度、群速度和 Grüneisen 参数作为一个整体来审视，才能对热输运行为做出可靠判断。