来源论文: https://arxiv.org/abs/2603.06966v1 生成时间: Mar 10, 2026 09:39

0. 执行摘要

拓扑物态与强关联电子系统的交叉研究是当前凝聚态物理的最前沿领域之一。重费米子（Heavy-Fermion, HF）体系由于其强烈的电子-电子相互作用，展现出诸如近藤效应（Kondo effect）、量子临界性以及非常规超导等奇特现象。近年来，Weyl-Kondo 半金属的发现证明了拓扑序可以与强关联效应共存。

本研究聚焦于非中心对称的重费米子化合物 $CeCoGe_3$。通过结合密度泛函理论（DFT）与动力学平均场理论（DMFT），研究团队揭示了该体系从高温不相干状态到低温相干重准粒子的演变过程。关键发现包括：在 $T=25$ K 时，体系表现出高达 52.6 倍的质量增强（$m^*/m_{DFT}$）；由于强关联、自旋-轨道耦合（SOC）与 $I4mm$ 晶格对称性的协同作用，费米能级附近（10 meV 以内）形成了高度平坦的拓扑节线（Nodal Lines）。这些平坦能带带来的极高态密度（DOS）可能在压力诱导的非常规超导中起到中介作用。本工作正式确立了 $CeCoGe_3$ 作为拓扑节线近藤半金属（Topological Nodal Line Kondo Semimetal）的原型地位。

1. 核心科学问题，理论基础，技术难点，方法细节

1.1 核心科学问题

在重费米子体系中，拓扑特征通常由 $4f$ 电子与传导电子（$c$-$f$ 杂化）形成的重准粒子能带决定。然而，如何定量描述强关联环境下拓扑能带的演变，特别是当体系缺乏空间反演对称性时，对称性保护的能带交叉点（如节线）如何保持稳定，是该领域的关键科学难题。对于 $CeCoGe_3$，虽然实验已观测到其压力诱导超导，但其底层拓扑物态是否参与了超导配对仍不明确。

1.2 理论基础：Kondo 物理与拓扑节线

  • 近藤晶格模型：$CeCoGe_3$ 中的 $Ce^{3+}$ 离子贡献了局部化的 $4f$ 矩，在低于近藤温度 $T_K$ 时，这些局部矩被传导电子屏蔽，形成重准粒子能带。
  • 非中心对称超导：$CeCoGe_3$ 结晶于 $BaNiSn_3$ 型四方结构（空间群 $I4mm$），缺乏反演中心。这意味着自旋不再是好量子数，能带会发生反对称自旋-轨道耦合分裂（ASOC），从而允许 Weyl 点或节线的存在。
  • 拓扑节线半金属（NLS）：能带交叉在布里渊区（BZ）形成闭合曲线。在本项目中，研究者区分了“本质（Essential）”节线和“偶然（Accidental）”节线，前者由空间群对称性强制锁定，后者由能带反转诱导。

1.3 技术难点

  • 强关联处理：传统的 DFT 在处理 $Ce$- $4f$ 轨道时会严重低估关联效应，无法捕捉质量增强和近藤共振峰。
  • 多尺度计算：需要同时处理原子尺度的库仑排斥（$U$）和宏观尺度的布里渊区积分，且需在不同温度下保持电荷自洽。
  • 能带平坦化：当质量增强因子 $1/Z$ 极大时，能带极其细窄，对数值精度提出了严苛要求。

1.4 方法细节：DFT+DMFT 流程

研究采用了电荷自洽的 DFT+DMFT 方法：

  1. DFT 部分：基于 WIEN2k 软件包，使用全电势线性增强平面波（FP-LAPW）方法。交换关联势采用 PBE-GGA。
  2. DMFT 部分：使用 EDMFTF 插件。将 $Ce$-$4f$ 轨道视为关联轨道，自能 $\Sigma(i\omega)$ 通过连续时间量子蒙特卡洛（CTQMC）杂质求解器获得。参数设定为 $U=6.0$ eV, $J_H=0.7$ eV。
  3. 准粒子哈密顿量构建：为了分析拓扑性质，构建了有效哈密顿量： $$H_{QP} = \sqrt{Z} [H_0 + \Re\Sigma(0) - \mu] \sqrt{Z}$$ 其中 $Z$ 是准粒子权重，$H_0$ 是从 Wannier 函数提取的单体项。

2. 关键 benchmark 体系，计算所得数据，性能数据

2.1 质量增强与相干性演变

计算结果显示了显著的温度依赖性：

  • 高温区 (T = 500 K)：虚轴自能 $-\text{Im}\Sigma(i0^+)$ 约为 2.29 eV，体系表现出典型的 Mott 式不相干行为，光谱函数呈现弥散状。
  • 低温区 (T = 25 K)：自能骤降至 0.125 eV，标志着近藤相干性的建立。计算得到的近藤温度 $T_K$ 约为 46 K，与实验值（约 50-80 K）吻合良好。
  • 质量增强因子：在 $T=25$ K，准粒子权重 $Z \approx 0.019$，对应有效质量增强 $m^*/m_{DFT} \approx 52.6$。这解释了实验观察到的巨大比热系数 $\gamma = 111 \text{ mJ/(K}^2\text{ mol)}$。

2.2 费米面重构 (Fermi Surface Reconstruction)

  • 在 $T > T_K$ 时，费米面主要由 $Ge$-$4p$ 和 $Co$-$3d$ 轨道的空穴口袋组成。
  • 在 $T < T_K$ 时，由于 $4f$ 电子参与杂化，费米面发生显著变化：$\Gamma$ 点附近的空穴口袋由于杂化间隙的产生而消失，取而代之的是高度重整化的 $f$ 性质能带交叉点。这一演变与 ARPES 实验结果一致。

2.3 拓扑节线的表征数据

研究识别出两类节线：

  1. 本质 Weyl 节线：位于 $\Gamma-Z$ 和 $X-P$ 路径上。由于 $C_{2z}$ 旋转对称性与垂直镜像对称性 $M_v$ 在 SOC 存在下的反对易关系，这些路径上的能带强制简并。能带在费米能级 $\pm 10$ meV 范围内极其平坦。
  2. 偶然拓扑节线：位于四个垂直镜像面上，由能带反转引起，受镜像特征值保护。计算得到的拓扑不变量 $|\xi|=1$，证明其具有鲁棒的拓扑保护。

3.1 核心软件栈

  • WIEN2k：用于高精度全电势 DFT 计算。 官网链接
  • EDMFTF (Haule’s DMFT)：用于执行电荷自洽的 DMFT 循环。 开源链接/教程
  • Wannier90：用于从 Bloch 态构造最大局域化 Wannier 函数（MLWFs）。 GitHub Repo
  • WannierTools：用于扫描布里渊区中的拓扑节点、节线及计算表面态。 GitHub Repo
  • VASP：用于结构优化和压力下的晶格弛豫。 官网链接

3.2 复现步骤指南

  1. 结构准备：使用 VASP 在 PBE 泛函下优化 $CeCoGe_3$ 的原子坐标。注意保持 $I4mm$ 对称性。
  2. DFT 初始化：在 WIEN2k 中运行常规自洽过程。建议 $RK_{max} = 7.5$，k 网格至少 $12\times 12\times 12$。
  3. DMFT 关联设置：
    • 在 case.indmft 文件中指定 $Ce$-$4f$ 轨道。
    • 设置 $U=6.0, J=0.7$，选择 CTQMC 求解器。
    • 运行 run_dmft0 进行电荷自洽。由于 $Z$ 非常小，收敛可能较慢，建议调整混合参数。
  4. 拓扑分析：
    • 使用 dmft_to_wannier 脚本导出重整化后的哈密顿量。
    • 将 $\sqrt{Z}$ 缩放后的参数输入 WannierTools。
    • 设置 SearchNodalLine = T 扫描镜像面上的节线位置。

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献

  1. [32] H.-H. Lai, et al., PNAS 115, 93 (2018)：提出了 Weyl-Kondo 半金属的理论框架。
  2. [39] A. Georges, et al., Rev. Mod. Phys. 68, 13 (1996)：DMFT 领域的奠基性综述。
  3. [40] G. Kotliar, et al., Rev. Mod. Phys. 78, 865 (2006)：DFT+DMFT 方法论的权威参考。
  4. [37] P. Li, et al., Phys. Rev. B 107, L201104 (2023)：关于 $CeCoGe_3$ 的关键 ARPES 实验证据。

4.2 局限性评论

  • 局部关联近似：DMFT 本质上是一种局域理论，忽略了动量相关的自能纠正（$\Sigma(\mathbf{k}, \omega)$）。在量子临界点附近，非局域关联（如反铁磁波动）可能对拓扑节线的精细结构产生影响。
  • $U$ 值的敏感性：虽然 $U=6$ eV 是 Cerium 体系的常用参数，但不同的双计数（Double-counting）方案可能导致费米能级附近 $f$ 带位置的微小漂移，进而改变“偶然”节线的精确位置。
  • 压力诱导的结构变化：尽管计算表明 $I4mm$ 结构在压力下稳定，但极高压下可能存在的微小畸变可能进一步破坏对称性，导致节线分裂为 Weyl 点。

5. 其他必要补充：对称性分析与超导意义

5.1 对称性保护机制的数学表达

在 $\Gamma-Z$ 路径上，波向量 $\mathbf{k} = (0, 0, k_z)$ 具有 $C_{4v}$ 小群。考虑二倍旋转 $C_{2z}$ 和镜像 $M_x$： 在存在 SOC 时，$C_{2z} = -i\sigma_z$，$M_x = -i\sigma_x$。由于 $\{\sigma_z, \sigma_x\} = 0$，对于 $C_{2z}$ 的本征态 $|c\rangle$，其镜像变换态 $M_x|c\rangle$ 具有相反的本征值且能量简并。这解释了为何红色的“本质”节线会被锁定在轴线上而不会被杂化打开能隙。

5.2 对非常规超导的启示

$CeCoGe_3$ 的平坦节线产生了一个巨大的态密度峰值。根据 BCS 或 Eliashberg 理论，超导转变温度 $T_c$ 指数依赖于费米能级的 DOS。在压力下，虽然有效质量下降（从 52.6 降至 22.6），但相干性的增强使得准粒子寿命延长。这种“强关联+高 DOS+非中心对称”的组合是实现拓扑超导（Topological Superconductivity）的理想温床。动量锁定的自旋织构（Spin Texture）可能诱导三重态配对分量，使 $CeCoGe_3$ 成为研究 Majorana 费米子的潜在平台。

5.3 结论

本研究通过数值严谨的 DFT+DMFT 计算，不仅定性解释了 $CeCoGe_3$ 的重费米子行为，更定量的预测了其平坦拓扑节线的存在。这为后续的高压 ARPES 实验和输运测量提供了明确的理论指引。