手机版

电声耦合计算

发布时间：2026-09-02 来源：科研学术网

小中大

字号：

一、项目背景与科学问题

电声耦合（electron-phonon coupling, EPC）描述了电子与晶格振动的相互作用，是解释金属超导电性、电荷密度波、电输运与热输运等现象的核心物理量。二硼化镁（MgB₂）以 39 K 的临界温度打破了传统 BCS 理论的麦克米兰极限，其超导电性源于金属硼层中强耦合的 E₂g 声子模与 σ 带电子之间的相互作用，是研究强电声耦合的典型体系。

本项目的科学问题是：MgB₂ 的 Eliashberg 谱函数 α²F(ω) 具有怎样的能量分布？电声耦合常数 λ 有多大？由 McMillan-Allen-Dynes 公式预测的临界温度 Tc 能否重现实验的 39 K？通过第一性原理电声耦合计算可以对这些量给出定量回答。

二、计算方法与理论背景

2.1 Eliashberg 谱函数

电声耦合的核心量为 Eliashberg 谱函数 α²F(ω)，它描述了能量为 ω 的声子对电子散射的加权态密度。由 α²F(ω) 可积分得到电声耦合常数 λ = 2∫α²F(ω)/ω dω。λ 越大表示电子-声子耦合越强。对 MgB₂，α²F(ω) 在约 35 meV 处出现源于 B 原子 E₂g 面内伸缩模的强峰，这一特征峰是 λ 的主要来源。

2.2 麦克米兰临界温度

基于强耦合理论，McMillan-Allen-Dynes 公式给出 Tc = (ω_log/1.2)exp[−1.04(1+λ)/(λ−μ*(1+0.62λ))]，其中 ω_log 为对数平均声子频率，μ* 为库仑赝势（通常取 0.1–0.15）。对 MgB₂，λ ≈ 0.87、ω_log ≈ 67 meV、μ* = 0.12 时，Tc ≈ 39 K，与实验一致。该公式是连接第一性原理电声耦合计算与宏观超导转变温度的桥梁。

2.3 计算方法

电声耦合矩阵元可用密度泛函微扰理论（DFPT）或 EPW 方法计算：先在均匀 q 网格上计算声子色散与电声矩阵元，再通过 Wannier 插值获得密集网格上的 α²F(ω)。

三、计算设置与模型构建

软件与方法：量子力学第一性原理（DFPT/EPW 框架），平面波赝势。
结构：MgB₂ 六方结构，晶格常数 a = 3.086 Å、c = 3.524 Å，包含 Mg 与 B 的 AB 层状排列。
电子计算：均匀 k 点网格 24×24×18，平面波截断 60 Ry。
声子计算：q 点网格 6×6×4，DFPT 求解动力矩阵与电声矩阵元。
后处理：Wannier 插值得到密集 α²F(ω)，积分求 λ 与 ω_log，代入 McMillan 公式求 Tc。

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

图 5 的左面板（a）给出 MgB₂ 的 Eliashberg 谱函数 α²F(ω) 与声子态密度（DOS）对比。α²F(ω) 呈现三个主要峰：约 35 meV 处的最强峰对应 B 原子面内 E₂g 声子模，是电声耦合的主导来源；约 62 meV 处的次强峰对应 B 面内光学模；约 80 meV 处还有较弱的贡献。声子 DOS 与 α²F(ω) 的差异清晰表明：并非所有声子都对耦合有同等贡献，只有与 σ 带电子强烈耦合的 E₂g 模在 α²F(ω) 中占据主导，这正是 MgB₂”声子选择性强耦合”的特征。

图 5 的右面板（b）给出累计电声耦合常数 λ(ω) 随能量的增长曲线。λ 从低能端快速上升，在约 40 meV 处已贡献主要部分，最终收敛于 λ = 0.87。图中同时标注了 McMillan 公式的输入参数：ω_log = 67 meV、μ* = 0.12，由此预测 Tc ≈ 39 K，与实验测得的 39 K 临界温度高度吻合。这一结果定量说明了 MgB₂ 的高温超导性源于强电声耦合与高频率声子的协同作用：E₂g 模既提供大的 λ，又因频率适中保证了 ω_log 不至于过低，从而突破传统金属氢化物中超导转变温度的上限。

两面板结合，左面板揭示了 α²F(ω) 的能量分布与声子选择机制，右面板展示了 λ 的累积过程与 Tc 的预测路径，共同完成了 MgB₂ 电声耦合的第一性原理分析，验证了 BCS 强耦合框架对 MgB₂ 的适用性。