超导Kagome金属CsV3Sb5中电荷有序的软模起源

《Nature Communications》：Soft mode origin of charge ordering in superconducting kagome CsV3Sb5

　　

摘要：Kagome金属AV3Sb5（A＝K, Cs, Rb）中电荷密度波（CDW）序与超导共存，引发了关于其交织相驱动机制的基本问题。研究人员结合高分辨非弹性X射线散射（IXS）与第一性原理计算，揭示了CsV3Sb5中CDW形成的起源。在结构因子分析指导下，研究人员沿倒空间M–L方向识别出一支软化声子模（soft phonon mode），该软化效应在L点最强，且此处弹性散射强度随降温增长最快。考虑晶格非谐性（lattice anharmonicity）与电子-声子耦合（electron-phonon coupling, EPC）的第一性原理计算重现了上述观测结果，并确定L点的软模不稳定性是CDW形成的驱动机制。尽管相变为弱一级（weakly first-order），本研究结果明确证明CsV3Sb5中的CDW源于声子软化（phonon softening），阐明了其微观起源并凸显了晶格动力学在Kagome金属中的核心作用。

《Nature Communications》论文解读：超导Kagome金属CsV₃Sb₅中电荷有序的软模起源

研究背景与意义

Kagome晶格因几何阻挫可承载Dirac锥、平带、范霍夫奇点（van Hove singularities, VHSs）及Z₂拓扑不变量，近年发现的AV₃Sb₅（A＝K, Rb, Cs）家族兼具约94 K的电荷密度波（charge-density-wave, CDW）转变温度T_CDW与低温超导（T_c≈2.5 K for CsV₃Sb₅）。然而CsV₃Sb₅中CDW的本质长期存疑：争议包括CDW是2×2×2还是2×2×4超结构、是否源于费米面嵌套（Fermi surface nesting）驱动的Peierls不稳定性，或动量依赖的电子-声子耦合（electron-phonon coupling, EPC）；此前非弹性X射线散射（inelastic X-ray scattering, IXS）与中子散射未见明显声子异常，部分理论指出需考虑离子熵与非谐性方能理解高温对称相稳定性。明确CDW微观机制是理解Kagome金属中交织量子序的关键。本研究发表于《Nature Communications》，通过高精度IXS结合含非谐性的第一性原理计算，首次在实验上直接观测到驱动CDW的软声子模，确立其为EPC主导的声子不稳定性机制。

主要关键技术方法

研究人员以助熔剂法生长Cs:V:Sb＝1:3:5.1的单晶CsV₃Sb₅，经EDS成分分析与磁化率验证T_CDW≈94 K及T_c。先通过密度泛函微扰理论（density functional perturbation theory, DFPT）计算动力结构因子（dynamical structure factor, S(Q,ω)），筛选不稳定声子支具较大散射截面的布里渊区（选定Γ₁₀₃区M点Q＝(0.5,0.5,3)与L点Q＝(0.5,0.5,2.5)，排除先前使用的Γ₄₂₀区因结构因子近零）。在欧洲同步辐射ESRF ID28线站（能量分辨3 meV）与日本RIKEN SPring-8 BL43LXU线站（硅(12 12 12)背散射，能量分辨1.1 meV）进行热漫散射（thermal diffuse scattering, TDS）及变温IXS测量（97–302 K）。第一性原理计算采用随机自洽谐波近似（stochastic self-consistent harmonic approximation, SSCHA）纳入晶格非谐性与电子-声子相互作用，辅以高斯近似势（Gaussian approximation potential, GAP）训练于optB88-vdW泛函DFT数据；DFPT获取电子-声子线宽；混合基赝势法算谐波声子色散与结构因子。IXS能谱用阻尼谐振子（damped harmonic oscillator, DHO）卷积分辨率函数拟合获声子能量E_j、阻尼Γ_j及未阻尼能量E_0j＝√(E_j²＋Γ_j²)。

Results and discussion

谐波声子谱与结构因子计算：DFPT给出高温P6/mmm相谐波声子色散在M与L点附近出现虚频（负频率），沿M–L全段有一支不稳定。计算各布里渊区该支的动力结构因子发现Γ₄₂₀区L点结构因子近乎为零（此前未观测到此模原因），而Γ₁₀₃区M点Q＝(0.5,0.5,3)与L点Q＝(0.5,0.5,2.5)结构因子显著——仅不足10%布里渊区L点结构因子达此强度一半。SSCHA含非谐与EPC的谱权重计算支持此结论。指导了后续IXS实验几何选取。

热漫散射（TDS）确认布里渊区选择：95 K（略高于T_CDW）测(HK 0.5)与(HK 2.5)倒易面，Γ₄₂₀区A–L方向无漫散条纹（对应结构因子可忽略），Γ₁₀₃区L点Q＝(0.5,0.5,2.5)见强CDW前驱体漫散强度，印证该区适宜探测不稳定声子模。

高分辨IXS原始谱与低能谱重分布：在Γ₁₀₃区M点、L点及一低对称点测97–302 K IXS谱。降温至近T_CDW时M与L点低能（＜5 meV）谱权重增强，L点尤为显著；低对称点无此现象（Bose修正后温度无关）。Bose修正沿M–L方向谱在302 K与97 K间显示～11 meV处强度损失、＜5 meV增强，表明声子模软化可直接被观测，无需拟合。

声子软化定量拟合：DHO拟合获各温区声子能量。多数声子温变不显著；M点一支由302 K的～10 meV软化至97 K的～6 meV；L点同支由～10.3 meV软化至102 K的＜3 meV（97 K趋过阻尼难可靠提取）。SSCHA计算重现此温度依赖及M₁⁺与L₂^?不可约表示模软化（谐波计算不稳但被非谐性稳定化），并预言L点两L₂^?模间avoided crossing与谱权重转移。

软化支倒空间弥散与弹性峰增长：沿Q＝(0.5,0.5,L)即M–L线，302 K该支基本无色散～10.3 meV；197 K起全线软化且从M至L单调增强，确认L点为首要晶格不稳定性且软化展宽涵盖半布里渊区；平面内从M向Γ方向软化渐消并于Γ–M中点消失。弹性峰强度I_el(Q)随降温全线增长，L点附近最强，反映本征慢涨落CDW关联发育——除位移型（displacive）软模特征外还含有序-无序（order-disorder）中心峰成分。

机制判定——EPC驱动而非纯嵌套：软化在M–L全线＞4 meV且与计算强EPC区吻合，不符M点VHS嵌套致尖锐Kohn反常特征；压力研究中CDW波矢与T_CDW对压力敏感而VHS与嵌套变化小，支持EPC主导。L点不稳定性属L₂^?对称，M₁⁺/L₂^?模具强EPC而非此对称模几无温变与弱EPC。结论：CsV₃Sb₅CDW由L点L₂^?声子软化驱动，倾向三维2×2×2序，机制为各向异性电子-声子耦合驱动CDW，类2H-NbSe₂等体系，AV₃Sb₅族或共用此机制。

讨论与结论翻译

研究人员实验与理论结果表明CDW由L₂^?声子软化驱动，L点为首要不稳定性并倾向三维2×2×2有序图案。声子自能含电子相空间因子与电子-声子矩阵元；此处报道的软化在倒空间较广——沿M–L全线软化至少4 meV且平面延续至Γ–M半程，匹配计算预测的不稳定支强EPC区域，这与M点VHS嵌套致锐Kohn反常不符而暗示关联（动量依赖EPC）效应，类公认的各向异性EPC驱动CDW材料而非具Kohn反常体系。此结论与静水压研究（CDW波矢及T_CDW对压力敏感而VHS/嵌套几不变）及time-resolved ARPES指EPC为主机制相符。近期KV₃Sb₅亦报软声子关联CDW形成，暗示EPC驱动软声子机制或为AV₃Sb₅共有。此外L点附近软化声子近T_CDW趋过阻尼（线宽＞能量），M点见软支连续展宽（源于强EPC与非谐性）。综上，CsV₃Sb₅中CDW起源于L点声子不稳定性（软模），电子-声子耦合与晶格非谐性共同塑造Kagome金属低能物理，也为其中超导的声子机制提供支持。研究结果确立晶格动力学在几何阻挫Kagome金属交织序中的核心地位。