NSR： 强关联Kagome金属中的交织序

NSR： 强关联Kagome金属中的交织序

 

NATIONAL SCIENCE REVIEW 2026

CsCr₃Sb₅ 强关联Kagome金属中的交织序

Unraveling Intertwined Orders in the Strongly Correlated Kagome Metal CsCr₃Sb₅

导读 导读：通过多探针超快光学技术（TRR/TRB/TRMOKE）和弹性电阻测量，系统揭示了CsCr₃Sb₅中1×4 CDW、电子轨道向列性和三态Potts向列性的共存与交织。CDW能隙极小（Δ₀/k_BT*~1），电子轨道向列性源于多轨道平带简并解除——与铁基超导体惊人相似。该工作为理解强关联Kagome金属中平带驱动的量子序奠定了基础。

一、前言背景

Kagome晶格：拓扑与关联的交汇平台

• Kagome晶格由共角三角形构成，天然包含Dirac点、van Hove奇点和平带——三者在电子结构中同时出现

• 平带极大地压缩带宽（W→0），使电子动能与库仑排斥能比U/W → ∞，从而放大电子关联效应

• AV₃Sb₅（A=K, Rb, Cs）体系：发现超导、CDW、电荷序、向列性，但平带位于费米面以下较深能量处

• 关键局限：平带对输运性质影响有限 → 需要将平带调谐至费米面的新材料平台

CsCr₃Sb₅：首个费米面处内禀平带的Kagome金属

• Cr替代V：将多轨道平带绝热调谐至费米面，同时引入磁性（Cr³⁺, S=3/2）和增强电子关联

• 加压超导：Tc=6.4 K @ 4.2 GPa，首次在Kagome平带体系中实现超导

• 常压相变：T*=55 K处反铁磁+结构相变+1×4周期晶格畸变共存

• 此前争议：STM和ARPES对CDW波矢和对称性的结论不一致，电子向列性的起源未明

• 本文核心：多探针超快光学+弹性电阻，系统揭示CDW、电子轨道向列性和三态Potts向列性的交织

Doniach相图视角下的CsCr₃Sb₅

• 传统Kondo晶格：J_Kρ决定基态——大J_Kρ→Kondo屏蔽→顺磁，小J_Kρ→RKKY→磁有序

• CsCr₃Sb₅的类比：平带增强U/W，类似Kondo晶格中J_Kρ≪1的磁有序区

• 多轨道自由度：Cr³⁺的3d轨道（dxy, dyz, dxz, dz², dx²-y²）在平带中简并 → 向列性解除简并

• 交织序：自旋（反铁磁）+电荷（CDW）+轨道（向列性）+晶格（1×4畸变）四自由度耦合

▸ DFT (PBE):本文使用。能带结构+声子谱，成本~100核时。但低估带隙、关联效应。▸ DFT+U (PBE+U):修正Cr 3d自相互作用误差。U~3-5 eV。成本与PBE相当。▸ HSE06:精确带隙和能带色散。成本~100×PBE。对CsCr₃Sb₅大型超胞不现实。▸ GW:准粒子能带，精确能隙。成本~1000×PBE。仅适合单胞验证。▸ DMFT/GW+DMFT:处理强关联极限。需要Wannier投影+杂质求解器。CsCr₃Sb₅的理想理论工具。

【知识扩展】Rothwarf-Taylor模型：CDW能隙的声子瓶颈探针

• 提出背景：Rothwarf & Taylor (PRL 1967) 为解释超导体中准粒子弛豫的反常慢化而提出

• 物理图像：光激发产生准粒子→准粒子复合发出高频声子（ω>2Δ）→高频声子再次拆散Cooper对→循环

• 关键特征：τ_slow在T_c附近发散（∝Δ(T)⁻¹），这是能隙存在的标志性特征——声子瓶颈效应

• 适用范围：不仅限于超导体，任何有能隙的体系（CDW、SDW、电荷序）均适用

• 本文应用：τ_slow在T*≈55 K附近发散 → CDW能隙打开的直接证据，无需依赖动量空间探测

• 经典参考：Rothwarf & Taylor, PRL 19, 27 (1967); Kabanov et al., PRB 59, 1497 (1999)

• 迁移能力：凡有能隙打开的体系（CDW、超导、Mott绝缘体），RT模型都是超快光谱学的标准分析工具

二、研究方法

多探针超快光学系统：四项技术集成

• 光源：Ti:Sapphire振荡器+OPO，泵浦610 nm/探测780 nm，脉宽250 fs，重复率80 MHz

• EOM调制：473 kHz泵浦调制，锁相检测，信噪比优于10⁻⁶

• TRR模式：探测泵浦诱导的反射率变化ΔR/R，提取CDW能隙和准粒子弛豫

• TRB模式：光学桥路+平衡光电二极管检测偏振旋转，灵敏度~μrad，直接探测旋转对称性破缺

• TRMOKE模式：探测磁光克尔效应，外加6 T面外磁场

• 空间分辨率：泵浦光斑3 μm，探测光斑2 μm → 可分辨单个畴（~10 μm）

弹性电阻测量：向列性电子起源的判定

• 改进Montgomery构型：压电堆施加单轴应变，应变片直接监测

• 提取E₂g对称性分量：m_{E₂g} ~ (d(Δρ_xx/ρ₀)/dε_xx - d(Δρ_yy/ρ₀)/dε_xx)

• 判据：电子起源向列性 → E₂g分量在T*附近发散；结构起源 → 台阶式跳跃

• 居里-外斯拟合：T*以上E₂g ~ A/(T-T*) + y₀，提取T*≈48.6 K

第一性原理计算：DFT+声子谱

• VASP, GGA-PBE泛函，平面波截断能300 eV

• CDW相建模：通过交替磁构型诱导人工畸变后完全弛豫

• 声子谱：有限位移法（phonopy），比较正常相、2×2 ISD CDW和1×4 SA-ISD CDW

• 关键预测：仅SA-ISD相无虚频，且声子模式频率和对称性与实验观测一致

TRR/TRB信号的双指数拟合模型。A_fast/A_slow分别为快慢分量振幅，τ_fast/τ_slow为弛豫时间，A₀(t)为长延迟残余信号。卷积高斯函数考虑有限时间分辨率τ_res。

三、实验结果

3.1 基本物性表征

图1 | CsCr₃Sb₅的基本性质。(a) 晶体结构：Cr原子形成Kagome层，Sb原子位于六边形中心，Cs原子层间插层。(b) 电阻率温度依赖：T*=55 K处异常。(c) 磁化率（μ₀H=3 T, H//ab）：T*处转变，Curie-Weiss温度>300 K表明强几何阻挫。(d) 多探针超快光学系统示意图——TRR/TRB/TRMOKE三模式。(e) 温度依赖TRR信号三维图——揭示非单调准粒子弛豫。(f) TRB振幅角分布——低温C₂对称性，高温恢复C₆。

基本物性关键参数

• 晶体生长：自助熔法，典型尺寸0.5×0.5×0.02 mm³

• 电阻率：T*=55 K处出现异常，与相变吻合

• 磁化率：Curie-Weiss温度Θ_CW>300 K >> T*≈55 K → 强几何阻挫

• 反铁磁转变：与结构相变和CDW同时发生，暗示自旋-晶格耦合

• TRB角分布：3 K展示C₂对称性，60 K恢复C₆ → 旋转对称性破缺的直接证据

3.2 1×4 CDW的TRR证据

图2 | 1×4 CDW态。(a) TRR信号及双指数拟合。(b-c) 温度依赖振幅和弛豫时间：慢分量τ_slow在T*附近发散——Rothwarf-Taylor模型的声子瓶颈效应。(d) 残余振荡分量及傅里叶变换：0.66 THz和0.88 THz两个相干声子模式。(e) 探针偏振依赖的FT振幅：展示2重对称性。(f) 温度依赖FT强度图。(g) SA-ISD相示意图：1×4 CDW+2×4 SDW调制，旋转对称性破缺。

CDW能隙的Rothwarf-Taylor分析

• 双指数衰减：快分量τ_fast≈0.3 ps（电子-声子耦合），慢分量τ_slow≈6 ps（3 K）

• τ_slow在T*附近发散 → Rothwarf-Taylor模型的声子瓶颈效应 → 能隙打开的特征标志

• 零温能隙：沿a₁轴Δ₀=6.7±2.0 meV，沿a₂轴Δ₀=3.3±0.7 meV

• Δ₀/k_BT*~1：远小于典型CDW的2-9，也远小于CsV₃Sb₅的3-5

• 极小的CDW能隙解释了ARPES未能检测到能隙的原因——低于标准ARPES能量分辨率

相干声子：CDW结构的指纹

• 两个模式：0.66 THz（~22 cm⁻¹）和0.88 THz（~29 cm⁻¹）

• 偏振依赖：均展示2重旋转对称性 → 与C₂对称的CDW结构一致

• 频率和对称性：仅与DFT预测的1×4 SA-ISD CDW相吻合，排除2×2 ISD和正常相

• 温度依赖：频率软化+线宽展宽（Klemens模型），T*以上消失

• 声子CDW振幅模式：由声学支折叠至Γ点形成，是CDW的直接光谱学证据

沿两个主轴的CDW零温能隙值

CDW能隙与转变温度的比值——异常小，暗示弱耦合CDW

3.3 电子轨道向列性

图3 | 电子轨道向列性。(a) 选定偏振方向的各向异性TRR信号。(b) TRR信号随探针偏振变化。(c) T*以上各向同性TRR。(d-f) TRR振幅角分布：CsCr₃Sb₅与Fe-SCs相似的非对称四分支模式，CsV₃Sb₅几乎圆形。(g) E₂g弹性电阻分量的温度依赖——发散行为确认电子起源。

电子轨道向列性：多轨道平带的指纹

• 3 K TRR：显著的异常各向异性+符号反转+清晰2重对称性

• 80 K TRR：无角度依赖、无符号反转 → 向列性仅存在于T*以下

• 与CsV₃Sb₅对比：后者TRR振幅角分布几乎圆形 → 无电子向列性

• 与Fe-SCs对比：振幅角分布的非对称四分支模式惊人相似 → 多轨道简并解除的普适特征

• 物理机制：向列性解除费米面附近多3d轨道简并，各向异性地扭曲费米面，部分轨道DOS增加、部分减少

弹性电阻：向列性电子起源的判决性证据

• E₂g分量：在T*附近发散 → 电子起源向列性的标志性特征

• 与AV₃Sb₅对比：后者在CDW转变附近为"台阶跳跃" → 结构起源

• T*以上：E₂g遵循居里-外斯行为 y=A/(T-T*)+y₀ → 向列序逐渐建立

• 拟合参数：T*≈48.6 K, y₀=2.6, A=78.8

• 向列转变温度（居里-外斯）与密度波转变温度（Rothwarf-Taylor）一致 → 交织序

E₂g弹性电阻分量的居里-外斯拟合，T*≈48.6 K

E₂g对称性弹性电阻系数的定义

3.4 三态Potts向列性

图4 | 三态Potts向列性与交织序。(a) 3 K时不同探针偏振方向的TRB信号。(b) 三个相邻畴的角分布模式——各旋转120°。(c) 单畴内温度依赖TRB信号。(d-e) TRB振幅和弛豫时间温度依赖。(f) 自旋、轨道、电荷和晶格在平带平台上的相互作用示意图。

三态Potts向列性的直接观测

• TRB直接捕获旋转对称性破缺：振幅随偏振调制，确认C₆→C₂

• 实空间扫描（2 μm光斑）：揭示三种畴，主轴向各旋转120° → Potts型向列性

• 畴尺寸：~10 μm（3 K），比CsV₃Sb₅（~100 μm）小约一个数量级

• 外磁场不敏感：6 T面外磁场对TRB动力学无影响 → 反铁磁序稳健

• τ_slow在T*附近发散：与TRR相同的声子瓶颈效应 → CDW和向列性共享同一序参量

交织序的微观图像

• 四个自由度：自旋（反铁磁）+电荷（CDW）+轨道（向列性）+晶格（1×4畸变）

• 平带平台：多轨道平带在费米面处，强关联放大所有相互作用

• 几何阻挫：Kagome晶格+Curie-Weiss温度>300 K → 强自旋涨落

• 晶格畸变：1×4 CDW是对强自旋涨落的补偿，降低总能量

• 轨道分裂：结构畸变+轨道序的组合效应，进一步放大向列性

• 畴尺寸小：强关联效应限制了畴的生长，导致~10 μm的微畴

【DFT Tip 1】VASP声子计算的ENCUT不要用弛豫值

用phonopy算声子谱时，ENCUT必须比弛豫时高20-30%。本文ENCUT=300 eV对含Cr体系偏低——Cr的3d电子需要更高截断才能准确描述力常数。

建议：对含3d过渡金属的体系，ENCUT≥400 eV，或做ENCUT收敛测试——计算Γ点声子频率随ENCUT的变化，收敛到<1 cm⁻¹。

此外，ADDGRID=.TRUE. 可显著提高力计算的精度，对软模体系尤其重要。

【DFT Tip 2】phonopy超胞尺寸：2×2×1够吗？

本文比较了1×4 SA-ISD和2×2 ISD两种CDW相，用phonopy算声子谱。注意：超胞必须足够大，使力常数在超胞边界处衰减到零。

对于Kagome层状体系，面内至少3×3×1（甚至4×4×1），否则长程力常数截断会导致虚频。

验证方法：增大超胞后声学支在Γ点附近是否仍保持线性色散？若有弯曲，说明超胞不够大。

【DFT Tip 3】CDW声子不稳定性检测：虚频≠数值噪声

声子谱中出现虚频（负频率）是CDW不稳定性的关键判据，但必须区分物理虚频和数值虚频。

物理虚频：量级>1 THz，出现在特定q点，增大超胞不消失。数值虚频：量级<0.5 THz，出现在Γ点附近，增大超胞后消失。

本文S7中正常相和2×2 ISD相的虚频是物理的——它们对应CDW不稳定性，增大超胞后虚频反而加深。

【DFT Tip 4】Cr³⁺体系的DFT+U：不要盲信默认U值

PBE严重低估Cr 3d电子的在位库仑排斥，导致d带过宽、磁性描述不准。DFT+U可以修正，但U值选择是关键。

对Cr³⁺（3d³），常用U=3-5 eV, J=0.87 eV。但U值强烈依赖于：赝势类型（PAW vs US）、氧化态、配位环境。

建议：做U值扫描（U=2,3,4,5,6 eV），比较能带、磁矩、CDW能隙对U的敏感性。本文用PBE（U=0）已能定性描述CDW，但定量能隙可能偏低。

【DFT Tip 5】Kagome体系的SOC计算：层状材料的方向依赖

Kagome层状材料中SOC可能有显著的面内/面外各向异性。本文未包含SOC，但若研究磁各向异性或拓扑性质，SOC是必需的。

VASP中SOC计算注意事项：① 使用LSORBIT=.TRUE.；② 必须用非共线磁结构（LNONCOLLINEAR=.TRUE.）；③ SAXIS用于设定磁矩方向。

对Cr体系，SOC耦合常数~50 meV，可能不足以改变基态磁序，但足以影响Dirac点能隙和Berry曲率。

【DFT Tip 6】磁构型建模：从ISD到SA-ISD

本文通过交替磁构型（↑↓↑↓）诱导1×4晶格畸变来构建SA-ISD相。这种"人工畸变+弛豫"策略是处理CDW的标准方法。

注意事项：① 初始畸变幅度需合理（~0.05-0.1 Å），过大可能弛豫到错误极小值，过小弛豫回正常相；② 弛豫时用IBRION=2, ISIF=3（同时优化晶格参数）；③ 比较不同磁构型的能量——最低能量构型才是基态。

【科研经验】为什么ARPES看不到meV级CDW能隙？

问题：本文CDW能隙Δ₀~3-7 meV，但ARPES报告"几乎不变的能带结构"——为什么？

原因：① 标准ARPES能量分辨率~10-15 meV（激光ARPES好一些，~5 meV），本文Δ₀已接近或低于分辨率极限；② ARPES光斑>30 μm，而畴尺寸~10 μm——畴平均效应抹平了各向异性信号；③ 光电子逃逸深度限制表面灵敏度，表面可能重构。

解决方案：① 使用激光ARPES（6 eV）提高分辨率至<2 meV；② 微区ARPES（<5 μm光斑）分辨单个畴；③ 结合超快光谱学（如本文）——不依赖动量分辨率，声子瓶颈效应是能隙的普适探针。

建议：当ARPES和输运/光谱学结果矛盾时，不要轻易否定任何一方。要考虑分辨率、空间平均、表面效应等系统性因素。

【科研经验】畴平均效应：如何避免对称性破缺被"平均掉"？

问题：CsV₃Sb₅中STM看到了C₂对称性破缺，但ARPES和输运测到的是C₆——为什么？

原因：三种120°旋转畴（Potts型）统计等概率分布时，宏观测量平均恢复C₆对称性。探测尺度>畴尺寸时，对称性破缺被"平均掉"。

解决方案：① 提高空间分辨率——使探测光斑<畴尺寸（本文2 μm光斑<10 μm畴，成功）；② 施加外场（单轴应变、磁场）打破畴的简并分布，使某一畴取向占优；③ 用局域探针（STM、微区光学）而非宏观探针。

建议：发现对称性破缺信号时，先问：我的探测尺度是否小于畴尺寸？如果否，很可能被平均掉了。

【知识扩展】三态Potts模型：从Ising到Potts的向列性

• Ising向列性（C₄→C₂）：见于Fe-SCs，由x和y方向的Ising型对称性破缺描述。序参量是标量，两种畴（相差90°）。

• Potts向列性（C₆→C₂）：本文CsCr₃Sb₅的特征。三态Potts模型描述三个等价方向上的对称性破缺，三种畴（相差120°）。

• 物理差异：Ising向列性用Z₂序参量，Potts向列性用Z₃序参量——后者更丰富：畴壁结构、临界行为、外场响应均不同。

• 经典参考：Potts, Proc. Camb. Phil. Soc. 48, 106 (1952); Wu, Rev. Mod. Phys. 54, 235 (1982) — Potts模型综述。

• 本文贡献：首次在Kagome体系中通过超快光学+弹性电阻同时证实CDW和Potts向列性的共存，为Z₃序参量理论提供实验基准。

四、对比分析

CsCr₃Sb₅ vs CsV₃Sb₅：Cr替代的根本性变革

• 磁性：CsCr₃Sb₅含Cr³⁺(S=3/2)反铁磁序，CsV₃Sb₅非磁性V⁴⁺(S=0)

• 平带位置：CsCr₃Sb₅平带在费米面，CsV₃Sb₅平带在费米面以下~0.5 eV

• CDW能隙：CsCr₃Sb₅ Δ₀~meV, Δ₀/k_BT*~1；CsV₃Sb₅ Δ₀~40 meV, Δ₀/k_BT*~3-5

• 向列性：CsCr₃Sb₅有显著电子轨道向列性，CsV₃Sb₅无可检测的电子向列性

• 畴尺寸：CsCr₃Sb₅~10 μm，CsV₃Sb₅~100 μm

• 核心差异：Cr替代V将平带调至费米面+引入磁性+增强关联 → 全新交织序物理

CsCr₃Sb₅ vs Fe-SCs：向列性的跨体系类比

• 相似性：TRR符号反转+异常各向异性，弹性电阻发散，多轨道简并解除驱动

• 差异性：CsCr₃Sb₅为Potts型（C₆→C₂），Fe-SCs为Ising型（C₄→C₂）

• 耦合方式：CsCr₃Sb₅向列性与CDW、反铁磁交织；Fe-SCs向列性与条纹磁序耦合

• 普适性：多轨道平带中轨道简并解除是向列性的普适驱动力，不受晶格对称性约束

与ARPES结果的表观矛盾与解释

• ARPES：报告几乎不变的能带结构和C₆对称化等能轮廓 → 与CDW和对称性破缺矛盾

• 原因1：CDW能隙极小（~meV），低于标准ARPES能量分辨率（~10 meV）

• 原因2：典型ARPES光斑>30 μm，平均了不同畴（~10 μm）的信号 → 畴平均掩盖旋转对称性破缺

• 超快光学的优势：高空间分辨率（2 μm）直接探测单个畴内动力学，避免畴平均

五、讨论

平带物理的新范式

• 传统平带研究：魔角石墨烯等扭曲体系，平带由moire势产生

• CsCr₃Sb₅的新范式：内禀平带+Kagome几何+磁性+强关联 → 交织序

• 多轨道自由度：5个3d轨道在平带中近似简并 → 向列性作为轨道序

• 超快动力学：飞秒时间分辨率直接探测序参量的弛豫和相干激发

• 理论意义：为理解多轨道平带中的关联效应和交织序提供实验基准

超快光谱学的方法论价值

• 声子瓶颈效应：CDW能隙存在的普适探针，不依赖动量分辨率

• 相干声子：CDW结构的高精度指纹，可区分不同CDW构型

• 偏振依赖TRR/TRB：旋转对称性破缺的直接探针，可分辨单个畴

• 弹性电阻：区分电子起源和结构起源向列性的判决性方法

• 多探针协同：TRR（电荷）+TRB（旋转对称性）+TRMOKE（磁性）+弹性电阻（电子-晶格耦合）

未解决的问题与展望

• 精确转变温度：多个序参量在窄温度窗口内交织，需要更高精度实验区分各自T_c

• 动量空间探测：微区ARPES或STM直接解析不同畴的电子结构，验证轨道向列性微观机制

• 加压下的超快动力学：Tc=6.4 K超导态与向列性/CDW的竞争关系

• 理论发展：包含多轨道平带+几何阻挫+强关联的统一微观理论

• 材料探索：寻找其他费米面处内禀平带的Kagome材料

六、总结

核心发现

(1) 1×4 CDW的强光谱学证据：声子瓶颈效应+相干声子频率和对称性与DFT一致

(2) 电子轨道向列性：多轨道平带简并解除驱动，与Fe-SCs类似，弹性电阻发散确认电子起源

(3) 三态Potts向列性：TRB实空间扫描直接观测到三种120°旋转畴

(4) 交织序物理图像：自旋+电荷+轨道+晶格四自由度在平带平台上耦合

(5) 方法论贡献：多探针超快光学+弹性电阻是研究交织序的强大工具组合

参考文献

[1] Liu L, Li Y, Tan H, et al. Natl. Sci. Rev. 13, nwag044 (2026) — 本工作

[2] Ortiz BR, et al. Phys. Rev. Mater. 3, 094407 (2019) — AV₃Sb₅的发现

[3] Jiang K, et al. Natl. Sci. Rev. 10, nwac199 (2022) — AV₃Sb₅综述

[4] Xiang L, et al. Nature 632, 63 (2024) — CsCr₃Sb₅加压超导

[5] Xu Y, et al. Nat. Commun. 16, 1312 (2025) — CsCr₃Sb₅磁性与结构

[6] Tan H, Yan B. Phys. Rev. Lett. 134, 126501 (2025) — CsCr₃Sb₅ DFT

[7] Wu Q, et al. Nat. Commun. 14, 526 (2023) — CsV₃Sb₅超快动力学

[8] Chu JH, et al. Science 329, 824 (2010) — Fe-SCs向列性

[9] Rothwarf A, Taylor M. Phys. Rev. Lett. 19, 27 (1967) — RT模型

[10] Klemens PG. Phys. Rev. 148, 845 (1966) — 声子非谐衰减

【DFT Tip 7】Kagome体系的K点采样：Γ中心 vs Monkhorst-Pack

六方晶系（Kagome属于此）的布里渊区采样必须用Γ-centered网格（VASP默认），不能切换到Monkhorst-Pack。

原因：Γ-centered保证高对称点（Γ, K, M）被精确采样——这对Dirac点和van Hove奇点的描述至关重要。

K点密度：对于弛豫，至少9×9×5；对于声子超胞（如2×2×1），可以降低到3×3×5（超胞倒空间缩小）。

常见错误：对超胞仍用与单胞相同的K点密度 → 计算量无谓增加10倍。

【DFT Tip 8】ISMEAR选择：金属体系别用-5

CsCr₃Sb₅是金属，ISMEAR=1（Methfessel-Paxton）是最佳选择。SIGMA建议0.05-0.1 eV。

ISMEAR=-5（四面体法）给出精确DOS但需要极密K点（>20×20×10），且对金属体系稳定性差。

ISMEAR=0（Gaussian）简单但展宽过大时可能模糊费米面附近的精细结构。

声子计算：力和能量对SIGMA敏感，建议SIGMA≤0.05 eV，且始终用ISMEAR=1。

【DFT Tip 9】弛豫精度：EDIFF和力收敛对软模体系

对于有声子软模的体系（如CsCr₃Sb₅的CDW不稳定性），弛豫收敛标准必须更严格。

EDIFF=1E-6 eV（默认1E-4不够），EDIFFG=-0.001 eV/Å（力收敛阈值），否则残余应力可能导致虚假虚频。

验证：弛豫结束后检查离子受力是否全部<0.001 eV/Å，应力张量对角项是否<1 kBar。

更保险的做法：用IBRION=6（DFPT）代替IBRION=8（finite differences），对软模体系更稳定。

【DFT Tip 10】声子频率与实验对比：DFT的系统性偏差

DFT声子频率通常与实验有3-5%的系统偏差（PBE通常偏软，LDA偏硬）。本文DFT预测0.66/0.88 THz，应与Raman实验（S12, 120/138 cm⁻¹≈3.6/4.1 THz）区分——这是不同模式。

对比策略：① 比较色散形状和对称性，而非绝对频率；② 做频率缩放：实验/DFT比值通常0.95-1.05；③ 检查模式的原子位移矢量是否一致——这是确认模式对应关系的最可靠方法。

常见错误：直接比较DFT声子频率和实验值，发现差5%就认为"DFT不准"——实际上PBE的系统偏差在正常范围内。

反常霍尔进阶计算（按需）DFT+U / HSE06关联修正EPC 电声耦合超导Tc估计Monte Carlo磁相变模拟DMFT强关联极限▸ 关键参数: ENCUT≥400 eV, KPOINTS Γ-centered 9×9×5, ISMEAR=1, SIGMA=0.05, EDIFF=1E-6, EDIFFG=-0.001▸ 声子计算: 超胞≥3×3×1, ADDGRID=.TRUE., 区分物理虚频与数值虚频▸ CDW建模: 人工畸变+弛豫, 比较多种磁构型能量, 声子谱验证稳定性▸ 进阶分析: DFT+U/HSE06修正关联, Wannier90提取TB+Berry曲率, EPC估算超导Tc

【知识扩展】Klemens声子非谐衰减模型

• 提出背景：Klemens (PRB 1966) 为解释晶体中声子寿命的温度依赖而提出。核心假设：光学声子通过三声子过程衰变为两个声学声子（能量和动量守恒）。

• 频率温度依赖：ω(T)=ω₀ - A[2/(exp(ħω₀/2k_BT)-1)]，T升高→声子占据数增加→非谐耦合增强→频率软化。

• 线宽温度依赖：Γ(T)=Γ₀[1+2/(exp(ħω₀/2k_BT)-1)]，T*以上Phonon寿命急剧缩短→线宽展宽→模式消失。

• 本文应用（S6）：0.88 THz声子的频率软化和线宽展宽完美符合Klemens模型→确认该模式是CDW相的特征声子，而非缺陷或杂质模式。

• 迁移能力：Klemens模型是分析相干声子温度依赖的标准工具，适用于任何有声子非谐衰减的体系（CDW、铁电、热电材料）。

 

总览：从DFT声子到强关联Kagome的完整理论图景

本文DFT计算聚焦于声子谱验证CDW结构，已经做得很好——但以本文的丰富物理（平带、磁性、向列性、CDW、超导），还有大量值得深入的计算方向。以下推荐8项后续计算，按优先级排序。

Further 1 | DFT+U / HSE06 能带与PDOS

为什么值得算：PBE低估Cr 3d关联，d带过宽，平带色散可能不准。DFT+U/HSE06能给出更准确的轨道分辨能带结构。

能回答的问题：多轨道平带在费米面的精确能量分布？各轨道对DOS的贡献？向列性解除简并后哪些轨道位移最大？

适合体系：所有含3d过渡金属的Kagome材料。

输入：PBE弛豫后的CONTCAR→POSCAR，设置LDAU=.TRUE., LDAUTYPE=2, U=3-5 eV on Cr d。

 

Further 2 | Wannier90 紧束缚模型 + 费米面分析

为什么值得算：从DFT能带投影到Wannier函数，可以构建最小TB模型，提取跃迁参数和轨道成分——这是理解平带形成机制和向列性的理论起点。

能回答的问题：平带来自哪些轨道？轨道间跃迁积分是多少？向列性如何改变跃迁参数？费米面的嵌套矢量是什么？

适合体系：所有Kagome平带体系。

输入：DFT能带（PBE或PBE+U），Cr 3d轨道投影，Wannier90输入文件。

 

Further 3 | 电子-声子耦合（EPC）+ 超导Tc估算

为什么值得算：本文提到加压Tc=6.4 K，但EPC机制和Tc尚未被DFT研究。EPC可以回答超导配对机制是否为声子介导。

能回答的问题：电声耦合常数λ是多少？哪些声子模式贡献最大？McMillan公式估算Tc与实验是否一致？

适合体系：Kagome超导体（CsCr₃Sb₅, AV₃Sb₅）。

输入：DFT声子谱+Wannier TB模型，EPW代码或Quantum ESPRESSO的ph.x。

 

Further 4 | Berry曲率 + 反常霍尔电导（AHC）

为什么值得算：Kagome晶格天然支持非零Berry曲率，Cr磁矩提供时间反演对称性破缺→预期有大的AHC。

能回答的问题：CsCr₃Sb₅的AHC有多大？Berry曲率在动量空间如何分布？CDW/向列性对AHC的调制？

适合体系：磁性Kagome材料。

输入：Wannier90 TB模型，或VASP直接计算（LORBIT=11, ICHARG=11）。

 

Further 5 | Monte Carlo磁相变模拟

为什么值得算：本文Curie-Weiss温度>300 K但T*=55 K——强几何阻挫。Monte Carlo可以模拟Kagome自旋的有限温度行为。

能回答的问题：磁有序温度为何远低于Θ_CW？阻挫参数是多少？磁序与CDW/向列性的耦合？

适合体系：几何阻挫磁体。

输入：DFT提取的交换耦合参数J₁, J₂, J₃（通过总能差方法或TB2J）。

 

Further 6 | NEB计算畴壁迁移

为什么值得算：Potts向列性有三种畴（120°旋转），畴壁的能量和迁移势垒决定了畴的稳定性和外场响应。

能回答的问题：两种畴之间的畴壁能是多少？畴壁迁移的能垒？外场（应变、电场）能否驱动畴壁运动？

适合体系：多畴体系（向列性、铁电、铁磁）。

输入：SA-ISD超胞，两种畴构型作为初末态，5-7个中间图像。

 

Further 7 | Bader电荷 + ELF分析

为什么值得算：Cr替代V如何改变电荷分布？Kagome层的化学键特征？电荷转移与CDW的关系？

能回答的问题：Cr和V的Bader电荷差？电子局域化函数（ELF）在Kagome层的分布？

适合体系：任何含电荷转移的体系。

输入：PBE静态计算的CHGCAR。

 

Further 8 | DMFT：强关联Kagome平带的终极理论工具

为什么值得算：当U/W→∞时，DFT+U和HSE06的微扰处理失效。DMFT精确处理局域关联，是平带强关联体系的"金标准"。

能回答的问题：平带中的准粒子有效质量？Mott物理是否出现？轨道选择性的关联效应？光谱函数与ARPES对比？

适合体系：U/W比值大的强关联体系（包括CsCr₃Sb₅平带）。

输入：Wannier90 TB模型→DMFT杂质求解器（如TRIQS, w2dynamics）。

 

七、支撑信息

支撑信息总览

以下13个Supplementary Figures涵盖：实验系统详解、TRR/TRB数据拟合、磁场依赖、相干声子分析、第一性原理计算、弹性电阻测量等。每张图配详细解读。

S1 | 多探针超快光学系统

S1 | 多探针超快光学系统

实验室自建多探针超快光学系统示意图及Jones矩阵理论描述。系统集成低温（3 K）、强磁场（6 T）、高时空分辨率（250 fs, 2 μm）和光子能量可调性。TRR/TRB/TRMOKE三模式协同工作。

支撑信息图 S1

S2 | TRR数据双指数拟合细节

S2 | TRR数据双指数拟合细节

探针偏振和温度依赖TRR数据的完整拟合过程。展示沿a₁和a₂轴的双指数拟合结果（黑色拟合线），验证模型适用性。提取的振幅和弛豫时间汇总于正文图2b-c和图S4c-d。

支撑信息图 S2

S3 | 零场与6 T磁场下的TRR/TRB对比

S3 | 零场与6 T磁场下的TRR/TRB对比

零场和6 T面外磁场下的TRR和TRB数据对比。所有数据几乎完全相同，证实反铁磁相互作用稳健，不受面外磁场影响。TRB在磁场下等效于TRMOKE测量。

支撑信息图 S3

S4 | 沿主轴a₂的温度依赖TRR数据

S4 | 沿主轴a₂的温度依赖TRR数据

探针偏振沿主轴a₂的TRR三维图、选定温度TRR信号和双指数拟合参数。与沿a₁轴（正文图2a-c）类似，展示清晰的相变和Rothwarf-Taylor行为，CDW能隙Δ₀=3.3±0.7 meV。

支撑信息图 S4

S5 | 沿主轴a₁的相干声子振荡

S5 | 沿主轴a₁的相干声子振荡

残余振荡分量（蓝色）和傅里叶变换（绿色）在3 K沿a₁轴。两个声子模式：0.66 THz和0.88 THz，温度依赖FT强度图。T*以上模式逐渐消失，与a₂轴（正文图2d-f）一致。

支撑信息图 S5

S6 | 0.88 THz相干声子的非谐衰减

S6 | 0.88 THz相干声子的非谐衰减

温度依赖频率和半高宽（FWHM），使用Klemens模型拟合。频率软化+FWHM展宽→非谐声子相互作用。零温频率ω₀和线宽Γ₀由拟合提取。

支撑信息图 S6

S7 | 三种相的计算声子谱

S7 | 三种相的计算声子谱

正常相、2×2 ISD CDW和1×4 SA-ISD CDW的声子谱对比。正常相和2×2 ISD相存在虚频→结构不稳定。仅SA-ISD相无虚频，声学支折叠至Γ点形成光学支，与实验观测的0.66/0.88 THz模式一致。

支撑信息图 S7

S8 | 轨道投影态密度与轨道简并解除

S8 | 轨道投影态密度与轨道简并解除

SA-ISD CDW相中Cr1和Cr2（相反自旋）的轨道投影DOS。dyz（蓝色）和dxz（黄色）轨道简并解除→轨道序。示意图总结超快动力学全过程：光激发→轨道简并解除→各向异性CDW能隙形成→弛豫。

支撑信息图 S8

S9 | 探针波长依赖的TRR和TRB

S9 | 探针波长依赖的TRR和TRB

不同探针波长（740-860 nm）下的TRR和TRB数据。振幅随波长增加略有减小，弛豫时间不变→底层物理过程对波长不敏感，主要反映费米面附近DOS。

支撑信息图 S9

S10 | TRB数据双指数拟合细节

S10 | TRB数据双指数拟合细节

温度依赖TRB和三个畴中偏振依赖TRB的双指数拟合。拟合模型（公式1）很好地描述所有数据。提取的振幅和弛豫时间汇总于正文图4d-e。

支撑信息图 S10

S11 | TRB振幅的温度依赖与临界指数

S11 | TRB振幅的温度依赖与临界指数

双指数拟合提取的TRB快慢分量振幅的温度演化。使用Landau相变理论的临界指数模型拟合：O(T)=O₀(1-T/T_c)^λ（T

支撑信息图 S11

S12 | 不同偏振几何下的Raman光谱

S12 | 不同偏振几何下的Raman光谱

CsCr₃Sb₅单晶在XX、XY、LR、LL四种偏振几何下的Raman光谱（6 K, 532 nm激发）。两个声子峰：120 cm⁻¹和138 cm⁻¹，主要来自Sb原子振动。Raman选择规则确认Kagome晶体结构。

支撑信息图 S12

S13 | 改进Montgomery构型的弹性电阻测量

S13 | 改进Montgomery构型的弹性电阻测量

压电装置和Montgomery方法示意图。不同温度下应变依赖的相对电阻变化。E₂g分量在T*以上遵循居里-外斯行为：y=A/(T-T*)+y₀，T*≈48.6 K。向列涨落在T*附近峰值→电子起源。

支撑信息图 S13

L. Liu, Y. Li, H. Tan, et al. | National Science Review 13, nwag044 (2026) | Kagome金属 · 强关联 · 交织序 · 超快光谱