MgFe6Ge6原型双层Kagome家族的高通量筛选
CHINESE PHYSICS LETTERS 2026
MgFe₆Ge₆原型双层Kagome家族的高通量筛选
High-Throughput Screening MgFe₆Ge₆-Prototype Bilayer Kagome Family
导读：对 168 种 AM₆X₆ 双层 Kagome 化合物进行高通量第一性原理筛选，系统评估热力学/动力学稳定性和磁性，识别出 31 种弱耦合 BCS 超导体（Tc = 0.65–3.97 K，λ = 0.37–0.62）。EPC 主要由 Kagome 层振动驱动，含 Sn 材料的低频面外软模进一步增强 λ。Nb 基 Kagome 层展现出最优超导性能，HNb₆Sb₆ 的 Tc 最高（约 3.97 K）。本研究建立了双层 Kagome 超导体的全局图谱，为实验探索提供了理论数据库和设计原则。
一、前言背景
1.1 Kagome 晶格：几何阻挫的典范
Kagome 晶格由共角三角形网络构成，是凝聚态物理中几何阻挫体系的原型。这种独特的二维结构天然支持 Dirac 点、van Hove 奇点和几何诱导平带等非平庸拓扑电子态。
早期理论工作（Syôzi, 1951）奠定了 Kagome 晶格的统计力学基础，随后 Ohgushi 等人（2000）预言了 Kagome 晶格中自旋-轨道耦合驱动的反常霍尔效应，开启了 Kagome 材料电子性质研究的新纪元。
Kagome 材料的核心魅力在于：晶格几何与电子结构的深度纠缠——共角三角形的排列方式不仅决定了能带拓扑，还催生了阻挫磁性、电荷密度波（CDW）和非常规超导电性等丰富量子现象。
1.2 AV₃Sb₅ 家族：Kagome 超导研究的里程碑
AV₃Sb₅（A=K, Rb, Cs）家族是 Kagome 材料研究中最具代表性的体系。Ortiz 等人（2019–2021）的系列工作系统揭示了该家族中 CDW 与超导电性的复杂关联——CDW 转变温度约 78–103 K，超导 Tc 约 0.9–2.5 K。
大量实验手段——包括角分辨光电子能谱（ARPES）、扫描隧道显微镜（STM）、μ子自旋旋转（μSR）和核磁共振（NMR）——证实了 AV₃Sb₅ 中存在拓扑非平庸的电子结构、手性 CDW 序和非常规超导电性配对对称性。
理论方面，第一性原理计算揭示了 van Hove 奇点附近的电子不稳定性是 CDW 的驱动力，而 CDW 与超导电性的竞争/共存关系仍是该领域最活跃的研究方向之一。
1.3 "166" 家族的发现与演化
早在 1990 年代，RMn₆X₆（R=稀土，X=Sn, Ge）材料已被发现结晶于 HfFe₆Ge₆ 型结构（P6/mmm），其 Mn 原子构成 Kagome 层。Venturini 等人（1991–2004）系统研究了该家族的复杂磁有序——包括螺旋磁结构、变磁转变和各向异性磁相图。
与 AV₃Sb₅ 的 "135" 化学计量比不同，"166" 结构（AM₆X₆）具有双层 Kagome 特征：两个 Kagome 层被三角形 AX₂ 和蜂窝状 X₁X₃ 子晶格隔开，层间耦合更强，化学自由度更大。
近年来，非磁性变体如 ScV₆Sn₆ 的发现（Pokharel et al., 2021）标志着 "166" 研究从磁性主导转向电子序主导的新阶段——ScV₆Sn₆ 在 92 K 处表现出 CDW 转变，且能带结构中存在 Dirac 点、鞍点和平带等 Kagome 特征电子态。
1.4 RMn₆X₆ 磁性 Kagome 材料的启示
RMn₆X₆ 家族为理解 Kagome 层中磁性与电子结构的相互作用提供了丰富案例。实验上，该家族展现出多样的磁输运性质：包括反常霍尔效应（AHE）、拓扑霍尔效应（THE）和巨大的磁阻效应。
在 GdMn₆Ge₆ 中，非共面自旋织构产生的标量自旋手性引发了拓扑霍尔效应，揭示了 Kagome 晶格中实空间 Berry 相位与动量空间 Berry 曲率之间的深层联系。
这些磁性 "166" 材料的研究为后续非磁性变体的探索奠定了重要的结构-性质关联基础：Mn 在 Kagome 层中的磁性行为暗示了过渡金属 d 轨道填充对 Kagome 层电子态的关键调控作用。
1.5 非磁性 "166" Kagome 金属的崛起
Ortiz 等人（2025）系统考察了 AM₆X₆ Kagome 金属，引入 LnNb₆Sn₆ 系列，展示了 (1/3, 1/3, 1/3) CDW 有序。Yang 等人（2025）报道了 LuNb₆Sn₆ 中费米面驱动的电荷有序阻挫——观察到平光学声子支的显著软化和竞争电荷序关联，揭示了费米面拓扑、声子和电荷序之间在 166 型 Kagome 金属中的强耦合。
ARPES 测量结合 DFT 计算表明，RV₆Sn₆（R=Gd, Tb, Lu）的能带在不引入重整化的情况下即可重现主要色散特征，暗示 V-Kagome 家族中的电子关联较弱（Ding et al., 2023）。
这些发现为非磁性 "166" Kagome 金属的热力学和动力学可行性提供了有力证据，也激发了本工作对该家族进行系统性高通量筛选的决心。
1.6 本文的研究动机、目标与创新点
研究动机：尽管 "166" Kagome 家族具有丰富的化学组成和多样的电子序，但对该家族中弱耦合超导电性的系统性高通量筛选仍然缺失。
核心目标：(1) 对 168 种 AM₆X₆ 化合物进行高通量第一性原理筛选，系统评估其热力学稳定性、动力学稳定性和磁性；(2) 对满足双重稳定性判据的化合物进行 EPC 计算和 Tc 估算；(3) 建立双层 Kagome AM₆X₆ 家族中弱耦合 BCS 超导电性的全局图谱。
创新点：首次将热力学、动力学和磁性稳定性评估与完整的声子和 EPC 分析相结合，系统绘制 "166" 双层 Kagome 家族的超导电性全景图，为实验探索提供理论数据库和设计原则。
二、研究方法
2.1 高通量 DFT 计算框架总览
所有计算在密度泛函理论（DFT）水平上进行，使用 QUANTUM ESPRESSO 软件包（Giannozzi et al., 2009）。采用超软赝势（Vanderbilt, 1990）描述电子-离子相互作用，PBE 泛函（Perdew, Burke & Ernzerhof, 1996）处理交换关联效应，Grimme DFT-D3 修正（Grimme et al., 2010）考虑范德华相互作用。
该计算方案已在对同构 Kagome 家族材料的描述中得到验证（Tan et al., 2021），确保了对层状 Kagome 结构中范德华力和层间耦合的合理描述。
高通量工作流涵盖：结构优化 → 形成能计算 → 声子色散（DFPT）→ EPC 计算 → Tc 估算，全程使用统一参数以保证可比性。
2.2 赝势选择与计算参数
波函数截断能：80 Ry；电荷密度截断能：800 Ry。高斯展宽方法（σ=0.01 Ry）用于电荷密度计算。结构优化采用 16×16×12 的 Γ 中心 Monkhorst-Pack k 点网格，力收敛标准 10⁻⁵ Ry/Å，总能收敛阈值 10⁻⁶ Ry。
声子色散通过密度泛函微扰理论（DFPT, Baroni et al., 2001）计算。EPC 计算采用更密集的 20×20×12 k 点网格和 5×5×3 q 点网格。
平面波截断能和 BZ 采样通过代表性化合物的收敛性测试验证——总能量、声子频率和 EPC 相关量在高通量筛选所需的精度范围内均收敛。SOC 在筛选阶段被排除以减少计算成本。
2.3 热力学稳定性评估方法
形成能 Ef 定义为化合物总能量与组成元素单质总能量之差除以原子数：Ef = (E(AM₆X₆) − E(A) − 6E(M) − 6E(X)) / 13。负的 Ef 表示放热过程，化合物相对于元素单质热力学稳定。
元素参考能量取自 Materials Project 数据库（Jain et al., 2013），确保参考态的一致性和可重复性。所有计算使用统一参数。
需注意：单独的 Ef 不足以完全判断可合成性——部分亚稳材料（如稍正的 Ef）可在非平衡或高压条件下实验合成（Sun et al., 2016）。因此，Ef 筛选应结合动力学稳定性综合判断。
2.4 动力学稳定性评估方法
动力学稳定性通过 DFPT 声子色散计算评估——声子谱中无虚频（所有声子频率 ω > 0）即表示动力学稳定。虚频的出现意味着存在晶格畸变以降低能量的方向，对应 CDW 或其他结构相变。
声子计算采用 DFPT 线性响应方法，在 5×5×3 的 q 点网格上计算动力学矩阵，通过傅里叶插值获得沿高对称路径的声子色散。
对于部分填充 f 壳层的稀土体系，f-in-core 处理使声子计算的定性趋势有意义，但定量讨论需谨慎——因此本文的定量声子、EPC 和超导性质讨论仅限于不含部分填充 f 壳层的化合物。
2.5 EPC 计算与超导 Tc 估算
超导临界温度 Tc 基于 Allen-Dynes 修正的 McMillan 公式计算。该公式以对数平均声子频率 ω_log 和 EPC 常数 λ 为输入，引入 Coulomb 赝势 μ* 描述屏蔽库仑排斥。
EPC 常数 λ(ω) 由 Eliashberg 谱函数 α²F(ω) 的频率积分定义。α²F(ω) 反映声子模式对电子散射的贡献权重，是连接晶格动力学和超导电性的核心物理量。
为评估 μ* 的影响，Tc 值同时给出 μ*=0.10 和 μ*=0.13 两种取值的结果。μ* 的典型取值范围为 0.10–0.15，本文选择这两个值以覆盖常见金属体系的 Coulomb 赝势范围。
2.6 f-in-core 处理策略与适用范围
稀土元素采用 f-in-core 赝势——4f 电子被纳入冻结核心，不作为价电子处理。这一策略在高通量筛选阶段保证了计算效率和一致性。
f-in-core 处理的物理合理性：对于 "166" 家族中的非磁性金属化合物，过渡金属 d 态在费米面附近形成宽能带，电子关联较弱（Ding et al., 2023 对 RV₆Sn₆ 的 ARPES 研究证实了这一判断）。
在此框架下，引入在位 Hubbard U 不会定性改变稳定性判据或 EPC 排序——因此在工作流中排除了 DFT+U 以保持一致性。
适用范围：定量的 EPC 和 Tc 讨论仅限于不含部分填充 f 壳层的化合物；含稀土体系的结果主要用于展示定性的化学和结构趋势。
三、晶体结构与稳定性筛选
图 1：AM₆X₆ 材料的晶体结构。(a) 元素选择——蓝色、橙色和深绿色方块分别代表 A、M 和 X 原子；(b) 侧视图，展示沿 c 轴的四层堆叠——X₁、X₂、M 和 AX₃；(c) 俯视图，展示 Kagome 网络（M 原子共角三角形）、三角层（A 和 X₂）和蜂窝状子晶格（X₁/X₃）。
3.1 晶体结构精析
AM₆X₆ 属于六方晶系 P6/mmm 空间群（No. 191），为 HfFe₆Ge₆ 同构型，由单晶 X 射线衍射验证（Pokharel et al., 2021）。沿 c 轴堆叠四层：X₁ 层、X₂ 层、过渡金属 M 层（构成 Kagome 网络）和 AX₃ 层。
A 元素选自氢、碱金属（Li–Cs）、碱土金属（Be–Ba）和稀土金属（Sc, Y, La–Lu），共 28 种。M 元素选自 V、Nb、Ta（3 种）。X 元素为 Sb 和 Sn（2 种）。理论组合空间为 28×3×2=168 种化合物。
与 AV₃Sb₅ "135" 家族的单层 Kagome 结构相比，"166" 结构的双层 Kagome 特征导致更复杂的层间耦合和更丰富的配位环境——这为调控电子不稳定性（CDW 形成、轨道杂化、超导电性）提供了更多自由度。
3.2 与 "135" 家族的结构对比
135 结构（AV₃Sb₅）：单层 Kagome 层 + 简单层间填充，层间相互作用最小化（Ortiz et al., 2020）。166 结构：双层 Kagome 层 + 三角形（AX₂）和蜂窝状（X₁X₃）子晶格，层间耦合显著增强。
尽管两者共享相同的六方晶格对称性（P6/mmm），但 166 结构的单位胞包含更多不等价原子位置（4 个不同 Wyckoff 位置），导致更复杂的声子模式和更多的 EPC 通道。
结构复杂度的增加同时带来了化学可调性的提升：A 位、M 位和 X 位各自独立调控，为筛选最优超导组成提供了广阔的元素空间。
图 2：AM₆X₆ 化合物的高通量稳定性与磁性筛选。(a) 形成能热图——深绿色（ΔE<−30 meV/atom）、绿色（−30≤ΔE<−10 meV/atom）和黄色（ΔE≥−10 meV/atom）；(b) 动力学稳定性图谱——蓝色（无虚频）vs. 红色（有虚频）；(c) 磁性筛选——深色（非磁性，零磁矩）vs. 浅色（磁性基态）。
3.3 热力学稳定性的化学趋势
在 168 种候选材料中，157 种（93.5%）具有负形成能，表明热力学稳定性是该家族材料的普遍特征。形成能热图揭示了清晰的化学趋势：
A 位效应：碱土金属（Mg, Ca, Sr, Ba）A 位阳离子通常比碱金属（Li, Na, K, Rb, Cs）具有更低的形成能，表明更高的离子电荷有利于增强与 Kagome 框架的静电相互作用和结构匹配。
X 位效应：Sb 基化合物的形成能系统性地低于 Sn 基类似物——这与 Sb 更强的电负性和更小的原子半径导致的更紧密的键合有关。这些趋势与 Cai 等人（2023）对相关 Kagome 候选材料的高通量筛选结果一致。
3.4 动力学稳定性的筛选结果
94 种（56.0%）化合物动力学稳定——显著少于热力学稳定的 157 种，突显了两种稳定性判据之间的非平凡关系。动力学稳定性表现出强烈的 M 位和 X 位依赖性：
V₆Sb₆ 基元材料几乎全部动力学不稳定（SmV₆Sb₆ 是唯一例外），而 Nb₆Sb₆ 和 Ta₆Sb₆ 基元材料大多稳定（除 LiNb₆Sb₆、BeNb₆Sb₆、MgNb₆Sb₆ 和对应的 Ta 类似物外）。
含 Sn 材料中，碱金属和碱土金属基结构大多动力学不稳定，仅 HNb₆Sn₆、KNb₆Sn₆、CsNb₆Sn₆、HTa₆Sn₆ 和 NaTa₆Sn₆ 例外。总体而言，Nb 和 Ta 基 Kagome 层比 V 基 Kagome 层具有更好的动力学稳定性。
3.5 磁性筛选与综合统计
自旋极化 DFT 计算（Zeller, 2006）揭示：97 种（57.7%）材料具有零总磁矩，为非磁性材料。含 V₆Sn₆ 以及 Gd-、Dy-、Ho-、Er-等重稀土基元的材料表现出非零磁矩，具有本征磁性。
综合筛选漏斗：168 → 157（热力学稳定）→ 94（动力学稳定）→ 97（非磁性）→ 31（同时满足热力学+动力学稳定 + 非磁性 + 金属性）。
这 31 种材料构成了后续 EPC 分析和超导电性评估的核心候选集。值得注意的是，实验已知的 YV₆Sn₆、(Gd–Tm)V₆Sn₆、YbV₆Sn₆ 和 LuV₆Sn₆ 在筛选流程中均被预测为动力学稳定，验证了 DFPT 工作流的定性可靠性。
四、声子振动与超导电性
图 3：非磁性 AM₆X₆ 化合物的超导临界温度 Tc 分布。Tc 通过 Allen-Dynes 修正 McMillan 公式计算（μ*=0.10），色标从紫色（低 Tc）到红色（高 Tc）。空白区域表示未纳入定量 Tc 分析的组成。仅包含不含部分填充 f 壳层的化合物。
4.1 Tc 分布的整体特征
31 种候选材料的 EPC 常数 λ 范围为 0.37–0.62，Tc 范围为 0.65–3.97 K，属于典型的弱耦合 BCS 超导体。λ 和 Tc 呈正相关关系——高 λ 值对应高 Tc。
化学趋势：具有相同 M₆X₆ 单元的材料的 Tc 值相近，表明超导电性主要由 Kagome 层的电子和声子特性决定。(1) 含 Sb 材料通常比含 Sn 材料具有更高的 Tc；(2) 含 Nb 的 Kagome 层材料表现出最大的 Tc 值；(3) A 位原子对 Tc 的影响有限。
四种代表性高 Tc 材料——HNb₆Sb₆、ScNb₆Sb₆、LuNb₆Sb₆ 和 CsNb₆Sn₆——被选为进一步深入分析的对象，以阐明双层 Kagome 家族的晶格动力学和 EPC 特征。
4.2 Tc 的化学趋势解析
Sb vs. Sn：Sb 基材料 Tc 较高的原因可能与 Sb 的更强电负性导致的费米面附近电子态密度（DOS）增强有关，同时也与 Sb 原子质量较大导致的更低特征声子频率有关。
Nb vs. V/Ta：Nb 基材料 Tc 最高的现象可归因于 Nb 的 d 轨道能级与费米面的最优匹配——Nb 的 4d 轨道在能量上比 V 的 3d 轨道更接近费米面，而比 Ta 的 5d 轨道具有更少的自旋-轨道耦合效应。
A 位弱效应：除 H 外，A 位原子对费米面 DOS 的贡献极小，其振动模式与费米面附近电子态的耦合也较弱——这解释了为何 A 位替换对 Tc 的调控作用有限。
图 4：四种代表性化合物的声子色散（加权原子振动模式）和原子振动模式分辨的声子态密度。(a) HNb₆Sb₆；(b) ScNb₆Sb₆；(c) LuNb₆Sb₆；(d) CsNb₆Sn₆。不同颜色线表示不同原子振动模式的贡献。
4.3 四种代表性材料的声子振动特征
HNb₆Sb₆：H 原子振动集中在高频区域（100–700 cm⁻¹），面内-面外振动明显分离。Nb 原子振动主要分布在 100–200 cm⁻¹ 的中频区域，面内和面外振动耦合。Sb 原子振动在 0–200 cm⁻¹ 的低频区域，软模通常来自面外（z 方向）振动。
ScNb₆Sb₆ 和 LuNb₆Sb₆：声子模式相似——Nb 原子在中高频区振动，Sb 原子在中低频区振动。Lu 原子的面外振动几乎可忽略，反映了重原子质量对振动模式的抑制效应。
CsNb₆Sn₆：Cs 原子的面内和面外振动部分耦合。Nb 原子面外振动主导高频区（175–225 cm⁻¹），面内振动主要分布在中频区（120–230 cm⁻¹）。Sn 原子的 z 方向面外振动产生声学软模——这是增强 EPC 的关键因素。
4.4 软模与 CDW 不稳定性的关联
若干材料的声子谱在高对称点 A→L 路径上出现声子频率趋于零甚至变为虚频的现象，包括 KV₆Sb₆、KTa₆Sn₆、CsTa₆Sn₆、GdV₆Sb₆、DyV₆Sb₆、HoV₆Sb₆、ScV₆Sn₆、LaV₆Sn₆ 和 CeV₆Sn₆。
这种声子不稳定性通常与 CDW 相变相关——声子特征向量方向的晶格畸变可以打开能隙，降低电子能量。ScV₆Sn₆ 的 CDW 转变已被实验广泛证实和讨论（Arachchige et al., 2022; 等），进一步支持了本文声子计算的可靠性。
软模的存在对超导电性具有双重影响：一方面，软模增强 EPC（低频率 → 大的 λ 贡献），有利于提高 Tc；另一方面，软模可能驱动 CDW 转变，与超导电性竞争费米面态密度。
4.5 声子数据的系统性总结
所有 31 种超导候选材料的完整声子色散曲线和原子振动模式分辨的声子态密度见补充材料 Figs. S2–S19。这些数据提供了双层 Kagome 家族中声子-电子耦合的微观视角。
总体规律：A 位原子振动分布范围广泛（取决于原子质量），但除 H 外，对超导性质的贡献可忽略。EPC 主要贡献来自过渡金属原子形成的 Kagome 层振动。
含 Sn 材料中，Sn 原子的低频面外（z 方向）软模提供了额外的 EPC 通道——这一独特振动行为类似于 CDW 体系中软模与电子序的耦合机制，对 Tc 的增强具有重要贡献。
五、电子-声子耦合机制分析
声子线宽 γ_{qν}：描述波矢 q、模式 ν 的声子对电子散射的贡献权重，由 EPC 矩阵元 |g_{qν}|² 和费米面态密度共同决定。
EPC 常数 λ(ω)：Eliashberg 谱函数 α²F(ω) 的频率积分，反映所有声子模式对电子配对的净耦合强度。
Allen-Dynes Tc 公式：以 ω_log、λ 和 μ* 为输入，是弱耦合 BCS 超导体 Tc 的标准估算方法。
对数平均声子频率 ω_log：α²F(ω) 加权的对数平均频率，是 Tc 估算中的关键参数。
图 5：电子-声子耦合分析。(a) HNb₆Sb₆ 的声子线宽加权色散（橙色圆点），右侧面板为 Eliashberg 谱函数 α²F(ω)（灰色阴影）和 EPC 积分强度 λ(ω)（红线）；(b) ScNb₆Sb₆；(c) LuNb₆Sb₆；(d) CsNb₆Sn₆。
5.1 HNb₆Sb₆ 的 EPC 特征
HNb₆Sb₆ 的声子线宽分析表明，EPC 的主要贡献来自 H 原子的面外（z 方向）振动，Nb 原子有少量贡献。这是所有 31 种材料中唯一一个 A 位原子对 EPC 有显著贡献的案例。
H 原子的独特作用源于其极小的原子质量——高频振动模式（100–700 cm⁻¹）与费米面附近的电子态发生有效耦合。Eliashberg 谱函数 α²F(ω) 在高频区域显示出明显的峰值，与 H 原子的振动频率一致。
λ(ω) 的累积曲线显示，H 原子的高频振动贡献了 λ 总值的约 40–50%，Nb 原子振动贡献了剩余的 50–60%。
5.2 ScNb₆Sb₆ 和 LuNb₆Sb₆ 的 EPC 特征
对于 ScNb₆Sb₆ 和 LuNb₆Sb₆，声子线宽的主要贡献来自 Nb 原子振动。Sc 和 Lu 的振动模式对 EPC 的贡献几乎为零——这进一步证实了 A 位原子（除 H 外）对超导电性的影响可忽略。
两种材料的 EPC 机制高度相似：Nb 原子在 Kagome 层中的面内和面外振动模式与费米面附近的 Nb-d 轨道电子态发生有效耦合。α²F(ω) 在中低频区域（50–200 cm⁻¹）显示出主要峰值，与 Nb 原子振动频率一致。
这一相似性表明，当 A 位为重原子时，超导电性几乎完全由 Kagome 层的本征电子-声子耦合决定，A 位原子的主要作用是提供结构稳定性和电荷平衡。
5.3 CsNb₆Sn₆ 的软模增强 EPC
CsNb₆Sn₆ 的 EPC 机制最为复杂：声子线宽的贡献同时来自 Kagome 层（Nb 原子）振动和 Sn 原子的 z 方向振动，尤其是低频区域的软模。这是含 Sn 材料中软模增强 EPC 的典型案例。
α²F(ω) 在低频区域（< 50 cm⁻¹）显示出额外的峰值，对应 Sn 原子的面外软模。这些软模的频率极低，根据 λ(ω) 的定义（被 1/ω 加权），低频声子对 λ 的贡献被显著放大。
这一机制与 LuNb₆Sn₆ 中观察到的平光学声子支软化现象（Yang et al., 2025）类似——Sn 原子在 z 方向的振动自由度与费米面拓扑之间存在强耦合，是 166 型 Kagome 金属中 CDW 和超导电性竞争的关键微观机制。
5.4 Kagome 层振动主导 EPC 的普遍规律
对所有 31 种材料的系统分析揭示了普适规律：EPC 的主要贡献源自过渡金属原子形成的 Kagome 层振动。这一规律在 Nb 基材料中最为显著——Nb 原子的 d 轨道与费米面附近的电子态具有最优的能级匹配和空间重叠。
M 位选择对 Tc 的深远影响：当 M=Nb 时，Tc 显著高于 V 或 Ta。这可以从三个方面理解：(1) Nb 的 BCC 结构元素单质本身就是中等 Tc 超导体（Tc≈9.2 K），其电子-声子耦合本征较强；(2) Nb 的 4d 轨道在 Kagome 配位场下的能级分裂与费米面匹配最佳；(3) V 的 3d 轨道更局域化，Ta 的 5d 轨道受 SOC 影响更大。
含 Sn 材料中，Sn 原子的低频面外软模提供了额外的 EPC 通道。软模的频率越低，对 λ 的贡献越大——这解释了为何尽管含 Sn 材料整体 Tc 低于含 Sb 材料，但 CsNb₆Sn₆ 仍具有较高的 Tc（约 2.4 K）。
5.5 A 位原子与 Sn 软模的角色分工
A 位原子的角色：主要提供结构稳定性（离子半径匹配 Kagome 框架）和电荷平衡（价电子数调控费米面位置）。除 H 外，A 位振动对 EPC 的贡献可忽略——这一结论对所有 31 种材料普遍成立。
Sn 软模的角色：Sn 原子在 z 方向的低频面外振动是增强 λ 的 "第二通道"。软模的产生与 Sn 原子较重的质量、较弱的层间键合以及费米面嵌套效应有关。
设计原则：若目标是最大化 Tc，应优先选择 Nb 为 Kagome 层金属，Sb 为 X 位元素，并考虑引入轻质量 A 位原子（如 H 或 Li）以提供额外的高频 EPC 通道。若目标是研究 CDW-超导竞争，则应选择含 Sn 材料。
六、讨论与展望
6.1 弱耦合超导电性的微观起源
所有 31 种材料的 λ 值均小于 0.62，远低于强耦合阈值（λ ≈ 1.0–2.0），表明它们属于弱耦合 BCS 超导体。电子-声子耦合不足以驱动强耦合配对或非常规超导电性——超导配对机制应为传统的声子媒介 BCS 机制。
Kagome 平带效应：尽管理论上 Kagome 晶格的平带可以提供大的 DOS 以增强 Tc，但在本文研究的 "166" 材料中，平带主要位于费米面以下较深能量处，对费米面 DOS 的贡献有限——这可能是 Tc 偏低的一个重要原因。
低 Tc 并不意味着这些材料不重要——它们作为弱耦合 BCS 超导体的模型体系，为理解 Kagome 晶格中电子-声子耦合的基本规律提供了理想平台。
6.2 与已有研究的系统对比
Yi 等人（2022, Phys. Rev. B）对有限数量的 AM₆X₆ 材料进行了 DFT 计算，预测了若干非磁性化合物可能具有超导电性，但未进行详细的声子计算或 EPC 评估。
Cai 等人（2023, Chin. Phys. Lett.）应用高通量方法对更广泛的 Kagome 候选材料进行了筛选，但主要关注拓扑相和 CDW 相关不稳定性，未系统性地处理超导电性和 EPC。
本工作填补了关键空白：首次将热力学、动力学和磁性稳定性评估与完整的声子和 EPC 分析相结合，系统性地绘制了双层 Kagome AM₆X₆ 家族中弱耦合 BCS 超导电性的全景图。识别出 59 种稳定非磁性金属材料，并确定了 31 种弱耦合超导体的 Tc、λ 和声子线宽分辨贡献。
6.3 与 MgB₂ 及其他 Kagome 超导体的对比
MgB₂ 是经典的声子媒介 BCS 超导体（Tc=39 K），其高 Tc 来源于 B-B 键的高频拉伸模式（E₂g 声子）与 σ 能带空穴型费米面的强耦合。与之相比，本文 "166" 材料的 Tc 偏低（0.65–3.97 K），主要因为：(1) 过渡金属 d 轨道的 EPC 固有强度弱于轻元素的共价键振动；(2) Kagome 层的费米面 DOS 相对较低。
与 AV₃Sb₅ 家族（Tc≈0.9–2.5 K）相比，本文 "166" 材料的 Tc 范围（0.65–3.97 K）相当甚至略高——特别是 HNb₆Sb₆ 的 Tc 约 3.97 K，接近 AV₃Sb₅ 家族中 Tc 最高的 CsV₃Sb₅（约 2.5 K）的两倍。
这一对比突显了 "166" 双层 Kagome 结构在超导电性方面的独特优势：双 Kagome 层提供了更多的 EPC 通道，且化学组成的灵活性允许对 Tc 进行更精细的调控。
6.4 Nb 元素的特殊作用
Nb 在 "166" 家族中展现出最优的超导电性表现——当 M=Nb 时，Tc 系统性地高于 M=V 或 M=Ta。这一现象可以从三个层面理解：
(1) 本征 EPC 强度：Nb 元素单质是中等 Tc 超导体（Tc≈9.2 K, λ≈0.8），其费米面附近的电子对声子扰动具有较大的响应函数。
(2) 能级匹配：在 Kagome 配位场下，Nb 的 4d 轨道晶体场分裂后的能级与费米面的能量匹配最优——V 的 3d 轨道更局域（更小的径向扩展），Ta 的 5d 轨道受 SOC 影响更大（Z 更大）。
(3) 动力学稳定性：Nb 基 Kagome 层比 V 基 Kagome 层具有更好的动力学稳定性，这为实验合成提供了更有利的条件。
6.5 CDW 与超导电性的竞争与共存
若干 "166" 材料（如 ScV₆Sn₆）同时表现出 CDW 转变和潜在的超导电性，引出了该家族中 CDW 与超导电性竞争/共存关系的问题。这与 AV₃Sb₅ 家族中观察到的现象类似，但 "166" 结构的化学自由度更大，为理解这一竞争提供了更多调控维度。
从 EPC 角度看，驱动 CDW 的软模同时也增强 λ——但 CDW 会打开能隙，减少费米面 DOS，不利于超导电性。因此，最优 Tc 可能出现在 CDW 被抑制但仍保留部分软模增强效应的组成区域。
Yang 等人（2025）在 LuNb₆Sn₆ 中观察到的竞争电荷序关联和声子软化，强烈暗示费米面拓扑在 CDW-超导竞争中扮演核心角色——这为未来研究提供了重要方向。
6.6 未来展望
实验验证：31 种候选材料（特别是 HNb₆Sb₆、ScNb₆Sb₆、LuNb₆Sb₆ 和 CsNb₆Sn₆ 四种高 Tc 材料）等待实验合成和超导电性测量。其中 Nb-Sb 基材料因不含稀土元素，合成难度相对较低。
CDW 与超导电性的竞争：含 Sn 材料中软模驱动的 CDW 不稳定性与超导电性的相互作用值得深入研究——特别是通过压力、掺杂或化学替代调控 CDW 转变温度，以探索 Tc 的最优值。
强关联效应：本文排除了 DFT+U 以保持高通量计算的一致性，但对于含部分填充 d 或 f 轨道的材料，关联效应可能需要更精确的处理——如 DMFT 或 GW 方法。
扩展元素空间：除 V、Nb、Ta 和 Sb、Sn 外，探索其他过渡金属（如 Cr、Mo、W）和主族元素（如 Bi、Pb）可能发现新的 Kagome 超导体。机器学习辅助的高通量筛选有望加速这一探索。
七、总结
7.1 核心结论
(1) 对 168 种 AM₆X₆ 双层 Kagome 化合物进行了高通量第一性原理筛选——157 种热力学稳定，94 种动力学稳定，97 种非磁性。
(2) 识别出 31 种同时满足热力学和动力学稳定性判据的非磁性金属弱耦合 BCS 超导体，Tc=0.65–3.97 K，λ=0.37–0.62。
(3) EPC 主要贡献来自过渡金属原子形成的 Kagome 层振动。含 Sn 材料中，Sn 原子的低频面外（z 方向）软模提供额外的 EPC 通道，增强 λ。
(4) 建立了双层 Kagome AM₆X₆ 家族中弱耦合超导电性的全局图谱，为实验探索提供了理论数据库和设计原则。
7.2 方法学贡献
(1) 展示了高通量 DFT+DFPT 工作流在 Kagome 材料超导电性筛选中的有效性——从 168 种候选材料中高效筛选出 31 种有前景的弱耦合超导体。
(2) f-in-core 处理策略使得含稀土元素体系的定性趋势分析成为可能，同时将定量 EPC 讨论限制在物理上可靠的范围。
(3) 声子线宽分辨的 EPC 分析揭示了不同原子振动模式对超导电性的贡献，为理解 Kagome 超导体中的电子-声子耦合机制提供了微观视角。
(4) 本文建立的高通量筛选工作流可迁移至其他结构类型的 Kagome 材料体系，为未来更大规模的 Kagome 超导体搜索提供了方法论基础。
7.3 研究意义与影响
本研究是首个系统性绘制双层 Kagome AM₆X₆ 家族中弱耦合 BCS 超导电性全景图的工作，填补了该家族超导电性理论研究的空白。
通过整合晶格动力学和电子结构分析，本文揭示了 Kagome 层振动主导 EPC 的普适规律以及 Sn 软模增强 λ 的独特机制，为 Kagome 超导体的理性设计提供了微观物理基础。
随着实验合成技术的进步，本文预测的 31 种超导候选材料有望在未来被实验验证，推动 Kagome 超导体从理论研究走向实际应用。
九、主要参考文献
核心引用（本工作及直接相关文献）
[1] Shi LT, Si JG, Zhou HB, et al. Chin. Phys. Lett. 43, 030711 (2026) — 本工作。
[2] Ortiz BR, et al. Phys. Rev. Mater. 3, 094407 (2019) — AV₃Sb₅ 的发现。
[3] Ortiz BR, et al. Phys. Rev. Lett. 125, 247002 (2020) — CsV₃Sb₅ 中 CDW 与超导。
[4] Pokharel G, et al. Phys. Rev. B 104, 235139 (2021) — ScV₆Sn₆ 结构与电子性质。
[5] Cai G, et al. Chin. Phys. Lett. 40, 117101 (2023) — Kagome 候选材料高通量筛选。
[6] Allen PB and Dynes RC, Phys. Rev. B 12, 905 (1975) — McMillan 公式的 Allen-Dynes 修正。
[7] Yang FZ, et al. Phys. Rev. B 112, 245113 (2025) — LuNb₆Sn₆ 费米面驱动电荷有序。
[8] Ortiz BR, et al. J. Am. Chem. Soc. 147, 5279 (2025) — LnNb₆Sn₆ 系列。
补充引用（方法与背景）
[9] Giannozzi P, et al. J. Phys.: Condens. Matter 21, 395502 (2009) — QUANTUM ESPRESSO。
[10] Baroni S, et al. Rev. Mod. Phys. 73, 515 (2001) — DFPT 方法。
[11] Perdew JP, Burke K, Ernzerhof M. Phys. Rev. Lett. 77, 3865 (1996) — PBE 泛函。
[12] Jain A, et al. APL Mater. 1, 011002 (2013) — Materials Project。
[13] Yi XW, et al. Phys. Rev. B 106, L220505 (2022) — AM₆X₆ 超导预测。
[14] Ding J, et al. J. Phys.: Condens. Matter 35, 405502 (2023) — RV₆Sn₆ ARPES。
[15] Arachchige HWS, et al. Phys. Rev. Lett. 129, 216402 (2022) — ScV₆Sn₆ CDW。
[16] Venturini G, et al. J. Magn. Magn. Mater. 94, 35 (1991) — RMn₆X₆ 磁结构。
L.-T. Shi, J.-G. Si, H.-B. Zhou, et al. | Chinese Physics Letters 43, 030711 (2026) | 双层Kagome · 高通量筛选 · 电子-声子耦合 · 弱耦合超导
