d电子重费米子的发现范式Cr掺杂CsFe2As2中的轨道选择性关联增强

PHYS. REV. LETT. 134, 076504 (2025)

Paradigm for Finding d-Electron Heavy Fermions

导读导读：重费米子行为几乎全部出现在f电子体系中，d电子重费米子极为罕见。本文提出并实验验证了寻找d电子重费米子的通用策略：在Hund金属CsFe2As2中通过Cr掺杂向half-filling推进，利用轨道选择性关联增强，迫使部分Fe-3d轨道局域化，形成d轨道间的Kondo杂化。实验获得铁基超导体中最高的Sommerfeld系数（gamma=270 mJ/mol K^2），并构建了完整的Doniach型相图。这一工作为d电子重费米子材料的设计开辟了新范式。

一、前言背景

重费米子：从f电子到d电子

重费米子（heavy fermion）材料是强关联电子体系的经典代表，其特征是低温下有效质量m*可达自由电子质量的100-1000倍，Sommerfeld系数gamma极大（>400 mJ/mol K^2）。传统重费米子体系几乎全部基于f电子（Ce、Yb、U等稀土/锕系元素），其中局域f电子通过Kondo效应与巡游电子杂化，形成重准粒子。

然而，d电子体系的重费米子行为极为罕见。d电子比f电子更扩展，通常不易局域化，但Hund耦合金属（Hund metal）中的轨道选择性Mott相变（OSMT）为d电子重费米子提供了新机制。

本文提出并验证了寻找d电子重费米子的通用策略：在Hund金属中通过掺杂向half-filling推进，利用轨道选择性关联增强，迫使部分d轨道局域化，形成"d电子提供的重费米子态"。

Cr掺杂CsFe2As2：Hund金属中的重费米子

CsFe2As2是典型的Hund金属，Fe-3d电子处于半满的half-filling状态（~5.5 electrons/Fe）。通过Cr掺杂引入空穴，将体系推向half-filling，增强轨道选择性关联效应。

核心发现：(1) Cr掺杂将Sommerfeld系数gamma从~30 mJ/mol K^2提升至~270 mJ/mol K^2--铁基超导体中的最高值；(2) 输运性质和磁化率符合重费米子Fermi液体行为；(3) 在更高掺杂量下出现反铁磁序，形成Doniach型相图。

二、研究方法

理论框架：Hund金属+轨道选择性Mott转变

Hund金属是一类多轨道关联体系，其中Hund耦合J_H增强了轨道自由度的关联效应。在Hund金属中，不同轨道可以具有截然不同的关联强度--部分轨道局域化（轨道选择性Mott相），部分轨道保持巡游。

当体系向half-filling推进时，Hund耦合驱动的轨道选择性关联增强，局域轨道与巡游轨道之间的杂化产生Kondo型重费米子行为。与f电子重费米子的关键区别：这里的"局域"和"巡游"电子都来自d轨道，只是不同的轨道自由度。

DFT+DMFT（动力学平均场）是处理这一问题的标准方法。DFT提供能带结构，DMFT处理局域关联效应，两者自洽耦合。

Sommerfeld系数：正比于Fermi面态密度，有效质量m*的增强直接反映在gamma中

Wilson比率：磁化率与比热的比值，重费米子Fermi液体中R_W~2

Doniach相图：Kondo温度T_K与RKKY温度T_RKKY的竞争，决定重费米子态vs磁有序

三、核心结果

图 1：Sommerfeld系数gamma随空穴掺杂的变化。Cr掺杂CsFe2As2将gamma从~30 mJ/mol K^2提升至~270 mJ/mol K^2，为铁基超导体中最高值。虚线为理论预测的轨道选择性Mott转变趋势。

Sommerfeld系数的巨大增强

gamma从母体CsFe2As2的~30 mJ/mol K^2增加到Cr掺杂后的~270 mJ/mol K^2，增强近10倍。这一数值已进入重费米子区间（通常gamma>200 mJ/mol K^2即被视为重费米子）。

理论分析表明，gamma的增强主要来源于dxz/dyz轨道（Fe-3d的t2g轨道）的有效质量增加。这些轨道在half-filling附近趋于局域化，而dxy轨道保持巡游，形成了d轨道间的Kondo杂化。

这一机制与f电子重费米子类似：局域磁矩（d轨道）通过杂化与巡游电子（另一d轨道）耦合，产生重准粒子。

图 2：热力学和输运性质。(a) 电子比热Ce/T vs T^2；(b) 磁化率chi(T)；(c) 电阻率rho(T)。所有数据符合重费米子Fermi液体行为。

Fermi液体行为的实验验证

比热：Ce/T vs T^2在低温下为线性，外推得到gamma值。在最低温下（<2 K），Ce/T趋于常数，是Fermi液体的标志。

磁化率：chi(T)在低温下趋于常数（Pauli顺磁），符合Fermi液体预期。Wilson比率R_W通过chi_0和gamma的比值计算，落在重费米子Fermi液体区间。

电阻率：rho(T)~T^2行为（Fermi液体特征），系数A与gamma^2成正比（Kadowaki-Woods关系），进一步确认重费米子态。

图 3：(a) 理论Doniach相图，展示Kondo温度T_K和RKKY温度T_RKKY的竞争；(b) 实验(H,T)磁相图，显示反铁磁序在更高掺杂下的出现。

Doniach相图与磁有序

Doniach相图描述了重费米子体系中Kondo屏蔽（T_K~exp(-1/rho J)）与RKKY磁交换（T_RKKY~J^2 rho）的竞争。当T_K > T_RKKY时，体系为重费米子Fermi液体；当T_RKKY > T_K时，体系为磁有序态。

在Cr掺杂CsFe2As2中，低掺杂区（x<0.3）为重费米子态，高掺杂区出现反铁磁序。这符合Doniach相图的预期：向half-filling推进时，RKKY相互作用增强，最终超过Kondo屏蔽。

实验(H,T)相图显示反铁磁转变温度T_N~40 K，且随磁场增加而降低，与巡游反铁磁体的预期一致。

DFT Tips

【DFT Tip 1】DFT+DMFT：处理强关联体系的正确姿势

Hund金属和重费米子体系不能仅用DFT描述。DFT会严重低估d/f电子的关联效应，导致错误的能带结构和有效质量。DFT+DMFT是处理这类体系的标准方法。

基本流程：(1) DFT计算能带结构（VASP/WIEN2k）；(2) 用Wannier90构建紧束缚模型；(3) 将TB模型输入DMFT求解器（如TRIQS/EDMFT/w2dynamics）；(4) 自洽更新自能和化学势。

关键参数：Hubbard U=3-5 eV for Fe-3d, Hund J_H=0.7-0.9 eV。U和J_H的值可以通过cRPA方法从第一性原理计算得到，比经验值更可靠。

【DFT Tip 2】轨道选择性Mott转变的DFT特征：轨道分辨能带

轨道选择性Mott转变（OSMT）是Hund金属的核心特征。在OSMT中，部分轨道（如dxz/dyz）局域化，形成Mott绝缘体行为，而其他轨道（如dxy）保持金属性。

DFT+DMFT诊断OSMT：(1) 计算轨道分辨的谱函数A(omega)；(2) 检查各轨道在Fermi面处的谱权重--局域化轨道的谱权重在Fermi面处显著降低；(3) 计算轨道分辨的自能虚部Im Sigma(omega=0)--大虚部表示强关联。

常见错误：仅看DFT能带就声称OSMT--DFT无法正确描述轨道选择性的关联效应，必须用DMFT。

【DFT Tip 3】CsFe2As2的DFT计算：Fe-3d的Hubbard U

CsFe2As2是铁基超导体家族成员，Fe-3d电子需要适当的Hubbard U修正。U值过小-->低估关联效应，U值过大-->虚假Mott绝缘体。

推荐U值：PBE+U中Ueff=3-4 eV for Fe-3d。建议通过cRPA计算（VASP的LDAUTYPE=3）自洽确定U值，或与实验XPS/ARPES能带对比校准。

注意：铁基超导体中Fe-3d的U值通常在3-5 eV范围，但具体值依赖于泛函（PBE vs PBEsol）、赝势和结构。不同铁基材料（如BaFe2As2、FeSe）的U值可能不同。

【DFT Tip 4】Cr掺杂DFT计算：超胞与虚拟晶体近似

Cr掺杂CsFe2As2的DFT计算面临掺杂建模的挑战。两种主要方法：(1) 超胞法--构建大超胞，用Cr原子替换特定位置的Fe原子；(2) 虚拟晶体近似（VCA）--在赝势中混合Fe和Cr的原子序数。

超胞法更准确但，VCA更简便但可能丢失局域无序效应。对于重费米子行为，电子结构对局域环境敏感，推荐使用超胞法（至少2x2x2超胞）。

建议：比较不同Cr掺杂浓度（x=0.1, 0.2, 0.3, 0.5）的超胞计算结果，验证gamma增强趋势的一致性。

【DFT Tip 5】Sommerfeld系数的DFT计算：态密度+有效质量

Sommerfeld系数gamma=(pi^2/3)k_B^2 N(E_F)，其中N(E_F)为Fermi面态密度。DFT可以直接计算N(E_F)并估算gamma。

但DFT的gamma通常被严重低估（因为关联效应未包含），DFT+DMFT的gamma更接近实验。在DMFT中，gamma可以从自能修正后的有效质量m*/m=1-dRe Sigma/domega|_EF计算。

验证：将DFT+DMFT预测的gamma与实验比热数据对比，是验证理论框架正确性的关键步骤。

【DFT Tip 6】Kadowaki-Woods关系的DFT验证

Kadowaki-Woods关系A/gamma^2=常数是重费米子Fermi液体的标志。A为电阻率rho~AT^2的系数，gamma为Sommerfeld系数。

DFT+DMFT可以同时计算A和gamma：(1) gamma从有效质量增强得到；(2) A从DMFT的输运计算得到（Kubo公式）。

注意：A/gamma^2的比值与轨道简并度N有关--对于N=1的Kondo晶格，A/gamma^2~1e-5 muOhm cm (mol K/mJ)^2；对于N=2，比值约为其1/4。

【DFT Tip 7】Hund耦合J_H对轨道选择性的影响

Hund耦合J_H是Hund金属中轨道选择性的驱动力。J_H越大，轨道选择性越强，部分轨道越容易局域化。

在DFT+DMFT中，J_H通常取0.7-0.9 eV for Fe-3d（与原子物理值一致）。但需注意：J_H对结果的影响与U相当，需要做J_H的收敛性测试。

建议：在U=3-5 eV范围内，J_H=0.7-0.9 eV范围内，做U-J_H参数扫描，验证轨道选择性和重费米子行为的鲁棒性。

【DFT Tip 8】Wannier90构建Fe-3d轨道模型：投影选择

从DFT到DMFT的关键步骤是构建Wannier紧束缚模型。对于CsFe2As2，需要投影Fe-3d轨道（五个d轨道：dxy, dyz, dxz, dx^2-y^2, dz^2）和As-4p轨道。

Wannier90设置：(1) 初始投影选择Fe的d轨道和As的p轨道；(2) 纠缠能窗覆盖Fermi面上下的2-3 eV能带；(3) 验证TB模型与DFT能带的一致性（偏差<10 meV）。

常见错误：只投影Fe-3d而忽略As-4p，导致Fe-As杂化信息丢失，DMFT结果不准确。

【DFT Tip 9】铁基超导体的磁性：条纹反铁磁序

铁基超导体中常见的磁序是条纹反铁磁（stripe AFM），即沿一个方向铁磁排列、垂直方向反铁磁排列。在DFT中，需要构建相应的磁超胞。

对于CsFe2As2，磁序对电子结构有显著影响--条纹AFM序打开Fermi面nesting能隙，显著改变态密度和gamma值。

建议：(1) 比较非磁和条纹AFM态的总能量，确认磁基态；(2) 在磁基态上计算DMFT，得到更准确的gamma和有效质量。

【DFT Tip 10】重费米子态的ARPES特征：DFT辅助谱学分析

重费米子态在ARPES中有独特的谱学特征：Fermi面附近出现重准粒子峰（窄峰）+Hubbard带（宽峰）。DFT+DMFT可以预测这些谱学特征。

计算流程：(1) DMFT自洽计算得到自能Sigma(omega)；(2) 用自能修正构造谱函数A(k,omega)；(3) 模拟ARPES谱（包含矩阵元效应和分辨率展宽）。

注意：重费米子态的特征温度（相干温度T_coh）通常在10-100 K量级，ARPES测量需要在足够低的温度下进行。

知识扩展

【知识扩展 1】重费米子：Kondo晶格与Doniach相图

【理论解释】重费米子体系的核心物理是Kondo晶格问题：周期性排列的局域磁矩通过Kondo效应与巡游电子杂化，形成重准粒子。Kondo温度T_K标志了局域磁矩被巡游电子屏蔽的能标，Fermi液体相干温度T_coh标志了重准粒子形成的温度。

【Doniach相图】Doniach（1977）提出了描述Kondo晶格中Kondo屏蔽与RKKY磁交换竞争的经典相图。当交换耦合J很小时，T_RKKY > T_K，体系为磁有序态；当J增大到临界值时，T_K > T_RKKY，体系为顺磁重费米子态。两者之间的量子临界点（QCP）是重费米子研究的核心。

【经典体系】CeCu6（gamma~1600 mJ/mol K^2）、CeAl3（gamma~1620）、YbRh2Si2（QCP附近的非Fermi液体）、URu2Si2（隐藏序）。

【经典参考】Hewson, "The Kondo Problem to Heavy Fermions" (1993)；Coleman, "Heavy Fermions: Electrons at the Edge of Magnetism" (2007)；Stewart, RMP 56, 755 (1984)。

【知识扩展 2】Hund金属与轨道选择性Mott转变

【理论解释】Hund金属是多轨道关联电子体系，其中Hund耦合J_H（而非单纯的Hubbard U）是关联效应的主要驱动力。J_H将不同轨道的电子自旋对齐，降低动能，增强关联效应。在Hund金属中，不同轨道可以具有截然不同的关联强度--轨道选择性Mott转变（OSMT）。

【与Mott绝缘体的区别】传统Mott绝缘体（如VO2、NiO）由Hubbard U驱动，所有轨道同时局域化。Hund金属的OSMT是部分轨道局域化，部分轨道保持巡游--这是"坏金属"行为（rho(T)在高温下饱和）的微观起源。

【经典参考】de Medici et al., PRL 107, 256401 (2011); Georges et al., RMP 85, 421 (2013)--Hund金属综述; Yin et al., PRL 107, 216401 (2011)--铁基超导体中OSMT的DMFT证据。

【迁移能力】Hund金属+OSMT框架适用于所有多轨道关联体系：铁基超导体、钌氧化物（Sr2RuO4）、铬氧化物（CrO2）、铁氧体等。

科研经验

【科研经验 1】DFT+DMFT收敛性：自洽循环的陷阱

问题：DFT+DMFT自洽计算不收敛，自能和化学势在迭代中振荡。

原因：(1) 轨道间电荷转移的过度响应--DMFT更新自能后，DFT计算给出大幅变化的电荷密度，导致下一步DMFT的自能又大幅变化；(2) 化学势调整步长太大。

解决方案：(1) 使用混合（mixing）技术--新自能=alpha*新算+(1-alpha)*旧自能，alpha=0.3-0.5；(2) 逐步减小化学势调整步长；(3) 在DMFT中增加Monte Carlo采样数，减小统计噪声。

建议：DFT+DMFT的收敛是计算中最耗时的步骤。建议从较小的U值开始，逐步增大U，用前一个U的收敛解作为初始猜测。

【科研经验 2】实验与理论的对比：如何判断重费米子态

问题：实验中观察到gamma增强，如何判断是否真的进入了重费米子态？

判断标准（至少满足3条）：(1) gamma > 200 mJ/mol K^2；(2) Wilson比率R_W在1-3之间；(3) 电阻率rho~AT^2，且A/gamma^2满足Kadowaki-Woods关系；(4) 磁化率chi在低温下趋于常数（Pauli顺磁）；(5) 低温比热中出现Schottky异常（来自局域磁矩的激发）。

排除效应：gamma增强也可能来自磁涨落（非Fermi液体）或剩余Schottky贡献。需要结合多种测量手段，排除这些可能性。

建议：如果只有gamma增强但其他判据不满足，可能是"近重费米子"而非真正的重费米子态。