电声耦合计算详解：DFPT声子谱与超导转变温度预测

电声耦合计算最让项目组兴奋的一刻是：MgB2的McMillan公式给出Tc≈39K，实验值也是39K。这种精度在DFT预测中极为罕见。但这个结果的背后是三周的参数收敛测试和四轮迭代，中间走过两次弯路。电声耦合不是”跑通计算”就行的任务，每一步参数都可能让最终Tc偏差10K以上。

电声耦合的物理图像

电声耦合描述的是电子与晶格振动（声子）之间的能量交换。在常规超导体中，这个耦合是配对机制的来源：电子通过交换虚声子形成Cooper对。耦合强度用电声耦合常数λ表征，λ越大配对越强，Tc越高。

电声耦合的计算需要三个要素：

1. 电子能带结构（从DFT获得）
2. 声子谱（从DFPT获得）
3. 电子-声子矩阵元（描述特定声子模式对特定电子态的散射能力）

三者通过Eliashberg谱函数α²F(ω)联系在一起：

α²F(ω) = Σ_{qν} δ(ω - ω_{qν}) × Σ_{k} |g_{k,qν}|² δ(ε_k) δ(ε_{k+q})

这个公式的核心是电子-声子矩阵元g，它描述了一个电子从k态散射到k+q态、同时产生或湮灭一个qν声子的概率幅度。

DFPT计算声子谱：参数选择

项目组使用VASP的DFPT（密度泛函微扰理论）方法计算声子谱，配合PHONOPY做后处理。

声子计算的关键参数：

• ENCUT=600 eV：比常规结构优化高20%。声子频率对力常数极敏感，而力常数对平面波截断的收敛比能量慢——500 eV时声子频率可能还差2-3 cm⁻¹。
• EDIFF=1E-8：电子步收敛标准比常规计算严格两个数量级。力常数的精度取决于力对电子密度的微分，电子密度不收敛则力常数噪声大。
• PREC=Accurate：避免网格平移引起的系统误差
• IBRION=6或7/8：6是有限位移法，7/8是DFPT。项目组用IBRION=8（DFPT，考虑非解析项修正）
• 超胞：2×2×2或3×3×3。MgB2用的是4×4×3超胞（96原子），因为B原子的声子频率高（~800-1000 cm⁻¹），需要密集的q点采样

踩过一个坑：最初用2×2×2超胞，声子谱在Γ点附近的声学支出现虚频——不是结构不稳定，是超胞太小导致力常数矩阵截断误差大。换到4×4×3后虚频消失。

电子-声子矩阵元与Eliashberg函数

有了声子谱后，计算电子-声子矩阵元需要密集的k点和q点网格。项目组使用的方案：

1. 先在密集k点网格（24×24×16）上计算电子能带
2. 在q点网格（6×6×4）上计算声子频率和本征矢
3. 在每个(k, q)对上计算电子-声子矩阵元g

这个计算量极大——MgB2的24×24×16 k网格对应约9000个k点，6×6×4 q网格对应144个q点，每个q点计算所有k→k+q的矩阵元。实际使用的是Wannier插值方法：先在粗网格上做DFPT计算，再通过Wannier函数插值到密集网格。项目组使用EPW软件完成这一步。

EPW的计算流程：

1. VASP计算粗网格上的电子能带和DFPT力常数
2. Wannier90构建MLWF（Mg的s、p轨道 + B的sp²/sp_z轨道）
3. EPW在Wannier基下插值电子-声子矩阵元到密集网格
4. 计算α²F(ω)和λ

λ的收敛测试

电声耦合常数λ的计算公式：

λ = 2 ∫ α²F(ω)/ω dω

λ的收敛对q点密度极其敏感。项目组做了系统收敛测试：

从粗网格到极密网格，λ从约0.5涨到约0.8，Tc从约28K涨到约39K——差距超过10K。如果不做收敛测试，直接用粗网格的结果，会严重低估MgB2的超导转变温度。这就是为什么Wannier插值必不可少——直接在24×24×16 q网格上做DFPT计算在算力上不可行。

McMillan公式与Tc预测

有了λ，用Allen-Dynes修正的McMillan公式预测Tc：

Tc = ω_log / 1.2 × exp[-1.04(1+λ) / (λ - μ*(1+0.62λ))]

这里有两个关键参数：ω_log（对数平均声子频率）和μ*（有效库仑赝势）。

ω_log从α²F(ω)直接积分获得，不引入额外误差。MgB2的ω_log约60 meV（约700K）。

μ是最大的不确定性来源。它描述电子之间的库仑排斥对配对的抑制，通常取0.10-0.15。项目组的做法是：先取μ=0.10计算Tc，再取μ*=0.15计算，给出Tc的范围。MgB2的结果：

• μ*=0.10 → Tc≈42K
• μ*=0.12 → Tc≈39K（与实验最吻合）
• μ*=0.15 → Tc≈35K

μ*=0.12是经验上对金属间化合物最常用的值。但这个选择本身就有0.5左右的偏差空间——这是整个电声耦合计算中最大的系统性误差来源。

反思：理论预测的边界

MgB2的电声耦合计算能预测到40K附近，部分原因是MgB2是一个”干净”的常规超导体——电声耦合机制明确，没有强关联效应，没有磁性竞争。对于铁基超导体或铜基超导体，电声耦合不是主要配对机制，上述方法给出的Tc远低于实验值——这时候需要考虑自旋涨落或其他非常规配对机制。

电声耦合计算的适用边界：常规BCS超导体（MgB2、Nb3Sn、Pb等）效果好；非常规超导体效果差甚至给出定性错误结论。更多VASP高级计算的实战经验，可以参考VASP/第一性原理栏目，或返回科研学术网首页。