CsV₃Sb₅：软模驱动的 Kagome 电荷密度波
McGuinness, Henssler, Alkorta et al. | Nat. Commun. (2026) | 高分辨非弹性 X 射线散射（IXS）结合 SSCHA 非谐声子理论和结构因子计算，在 CsV₃Sb₅ 中首次直接观测到 CDW 软模——非谐效应主导，电子-声子耦合辅助，修正了"纯电子驱动 CDW"的主流观点。
Kagome · CDW · 软模 · IXS · 非谐效应 · SSCHA
Soft mode origin of charge ordering in superconducting kagome CsV₃Sb₅
一、研究背景
Kagome 金属 AV₃Sb₅：CDW 与超导共存
AV₃Sb₅（A = K, Cs, Rb）是一类 kagome 晶格超导体，因其 CDW 与超导共存而引发广泛关注。CsV₃Sb₅ 在 TCDW ≈ 94 K 发生 2×2×2 CDW 相变，同时保持超导 Tc ≈ 2.5 K，是研究 CDW-超导相互作用的理想平台。
CDW 起源机制长期争论：电子 Peierls 不稳定性（Fermi 面嵌套驱动）vs 软声子模式驱动（EPC 主导）。两者在物理上密切相关——声子自能同时受电子态密度嵌套和电子-声子矩阵元影响。
此前 Subires et al. (Nat. Commun. 2023) 的高分辨 IXS 研究未发现声子软化，被广泛引用为反对声子驱动 CDW 的证据。但该研究仅测量了特定 Brillouin 区的 L 点，结构因子可能恰好为零。
本文核心突破：通过结构因子计算引导，在 Γ₁₀₃ 等高结构因子 Brillouin 区测量 M-L 方向 IXS，首次直接观测到 CDW 软模。
实验策略：三角验证框架
高分辨 IXS（ESRF ID28 和 SPring-8 BL35XU）：能量分辨 ~1.5-3 meV，系统探测 M-L 方向声子色散和温度演化
漫散射 mapping：同步辐射 X 射线衍射，在 TCDW 上下测量弹性散射强度，追踪 CDW 序参量增长
第一性原理计算：SSCHA 非谐声子谱 + 温度依赖非谐效应 + 电子-声子耦合线宽计算
结构因子分析：S(Q,E) 拟合识别不同声子分支的贡献，定量分离各模式权重，指导 Brillouin 区选择
综合分析：IXS 实验 + 漫散射 + DFT 计算构成的三角验证框架，系统性解决 CDW 起源争议
二、研究方法与理论框架
实验方法：高分辨 IXS 与漫散射
IXS 测量：ESRF ID28 线站（3 meV 分辨率）和 SPring-8 BL35XU（1.5 meV 分辨率），Si 单色器和球面分析器，沿 M-L 方向扫描波矢
漫散射 mapping：Pilatus 2M 探测器，0.25° 角积分，0.5 s 曝光，CrysAlis Pro 软件重构倒空间
温度控制：液氦低温恒温器，97-302 K 范围，覆盖 TCDW 上下
数据分析：阻尼谐振子（DHO）模型拟合，卷积仪器分辨率函数，弹性线扣除 + Bose 因子校正
理论方法：SSCHA 非谐声子计算
SSCHA（随机自洽简谐近似）：超越传统 DFT 简谐声子，考虑温度依赖的非谐效应和量子涨落
DFT 参数：Quantum ESPRESSO，PBEsol 泛函，超胞法计算力常数，SSCHA 最小化自由能
非谐谱函数：从 SSCHA 声子自能计算温度依赖的谱函数 A(q,ω)，包含声子频移和展宽
EPC 计算：Wannier 插值获取精细 k/q 网格上的电子-声子矩阵元，计算动量依赖的 EPC 线宽
二·附、核心理论公式
散射截面与动态结构因子 S(Q,E) 的关系。IXS 测量的是 d²σ/dΩdE，正比于 S(Q,E)。
涨落-耗散定理：动态结构因子 S(Q,E) 通过 Bose 因子与动态磁化率虚部 χ″(Q,E) 关联。温度依赖来源于 Bose 因子和 χ″ 本身的非谐效应。
阻尼谐振子（DHO）模型：χ″(Q,E) 表示为多个 DHO 的叠加，每个 DHO 由振幅 Aⱼ、线宽 Γⱼ 和本征能量 Eⱼ 描述。
Bose 因子：n(ω,T) = 1/(exp(ħω/kBT) − 1)。低温下高频声子被冻结，低频声子被热激发，是区分平凡温度效应与真正声子软化的关键。
软模温度依赖：ω₀²(T) = ω₀²(0) + αT。非谐声子理论预言软模能量平方与温度呈线性关系，α < 0 表示软化。
电子-声子耦合（EPC）线宽：Γ_ep = 2πω λ_{Qν}。EPC 贡献的声子展宽正比于声子频率和耦合常数。
EPC 耦合常数：λ = N_F⟨I²⟩/(M⟨ω²⟩)。由费米面态密度 N_F、电声子矩阵元 I 和原子质量 M 共同决定。
无阻尼声子频率：E₀ = √(E² + Γ²)。实验拟合得到的是有阻尼的峰值能量 E，需通过线宽 Γ 修正得到无阻尼频率 E₀。
SSCHA 自由能：F_SSCHA = F₀ + ⟨V − V₀⟩₀。通过最小化自由能得到温度依赖的有效简谐力常数，自然包含非谐效应。
三、主要实验结果
Fig 1 核心发现：计算声子强度与 Brillouin 区选择
计算揭示 CsV₃Sb₅ 沿 M-L 方向存在一支低频光学声子模式（L₂⁻ 分支），在 L 点频率最低，主要涉及 Cs 原子沿 c 轴的振动
结构因子 S(Q,E) 在不同 Brillouin 区间差异极大：<10% 的 Brillouin 区在 L 点有可观的结构因子，解释了此前 Subires et al. 未观测到软化的原因
Γ₁₀₃ 区（Q=(0.5, 0.5, 2.5) 附近）L 点结构因子最强，是本实验的关键测量窗口
非谐计算（SSCHA）显示声子频率在 TCDW 附近显著降低，而简谐计算无法捕捉这一软化行为
图 1 | 计算声子强度。(a,c) 简谐和非谐计算在 Γ₄₂₀ 和 Γ₁₀₃ 区的动态结构因子。(b,d) 对应的声子色散。红色箭头标记 L₂⁻ 软模。非谐计算（c,d）显示 M-L 方向显著软化。
Fig 2 核心发现：TCDW 之上的热漫散射
95 K（略高于 TCDW ≈ 94 K）的热漫散射 mapping 揭示 M-L 方向存在显著的漫散射强度
Γ₁₀₃ 区 L 点附近漫散射最强，与结构因子计算完全一致，直接验证了 Brillouin 区选择策略
漫散射强度沿 M-L 线具有明显的波矢选择性，L 点附近最强，M 点较弱
弹性散射强度在 TCDW 以下急剧增长，而漫散射强度在 TCDW 以上已显示前驱特征
图 2 | 95 K 热漫散射（TCDW 之上）。(a-c) Γ₄₂₀、Γ₁₀₃ 和 (H,K,0) 面的漫散射 mapping。(d) 沿 M-L 方向的强度分布，L 点附近最强。
Fig 3 核心发现：IXS 原始谱与 Bose 校正
原始 IXS 谱（未归一化）在 302 K 至 97 K 之间系统测量，M 点和 L 点各温度下的能谱叠加对比
Bose 因子校正后，302 K 和 97 K 的谱权重差异显著减小——但 L 点仍存在超出 Bose 因子的额外谱权重转移
M 点：谱权重温度依赖基本由 Bose 因子解释，声子软化幅度较小（~6 meV → ~4 meV）
L 点：谱权重从 ~8 meV 向低能（~3 meV）显著转移，超出 Bose 因子预期，直接证明声子软化
拟合采用 DHO 模型卷积仪器分辨率，分离弹性线、低频软模和高频声子分支的贡献
图 3 | 原始 IXS 谱。(a-c) M 点 97-302 K 非归一化 IXS 谱。(d-f) L 点 97-302 K 非归一化 IXS 谱。(g-i) Bose 因子校正后的谱，L 点仍显示低能谱权重转移。
Fig 4 核心发现：拟合声子能量与线宽的温度依赖
M 点：6 支声子模式拟合，其中一支在 ~10 meV 处从 302 K 的 10 meV 软化至 97 K 的 6 meV
L 点：软模能量从 302 K 的 ~8 meV 单调降至 97 K 的 ~3 meV，接近完全软化
声子线宽 Γ 在 TCDW 以上基本恒定为 ~2 meV，说明 EPC 贡献有限，非谐效应主导
阻尼比 Γ/E₀ 在 TCDW 接近时因 E₀ 减小而增大，但 EPC 线宽无明显异常增强
综合判断：CDW 软模主要由晶格非谐效应驱动，EPC 提供辅助但非主导作用
图 4 | 拟合 IXS 谱与声子能量。(a-f) M 点和 L 点各温度下 DHO 拟合结果。(g-l) 拟合声子能量和线宽的温度依赖。
Fig 5 核心发现：相变级数、软化空间尺度与弹性峰增长
软模能量 ω₀ 与温度的对数-对数图显示 ω₀ ∝ (T − TCDW)^(1/2) 的标度关系，符合二级相变（mean-field）的经典软模行为
沿 M→Γ 方向测量软化效应的空间尺度：远离 M 点向 Γ 移动时，软化逐渐减弱，说明 CDW 不稳定性局域在 M-L 线附近
弹性散射强度 I_elastic(T) 在 TCDW 以下急剧增长，零能弹性峰强度与 CDW 序参量 |Δ|² 成正比
弹性强度在 TCDW 附近显示微小的滞后回线（~1-2 K），暗示 CDW 相变具有弱的一级相变特征
图 5 | 相变级数与软化空间尺度。(a) 软模能量 ω₀ 的对数-对数标度。(b) 沿 M→Γ 方向软化效应的空间衰减。(c) 弹性强度 I_elastic 的温度依赖，展示 CDW 序参量增长。
四、讨论
软模起源：非谐效应 vs 电子-声子耦合
ω₀²(T) ∝ T 的线性关系是经典非谐声子软化的标志性行为，与 PbTe、SnTe 等铁电软模体系类似
EPC 线宽 Γ_ep 在 TCDW 附近无明显增强，说明电子关联非 CDW 主要驱动力——但这不意味着 EPC 不重要
更准确的物理图像：非谐晶格不稳定性是 CDW 的"引擎"，EPC 通过重正化声子频率辅助软化，两者协同作用
这一结论与纯电子驱动（Fermi 面嵌套）的 CDW 图像形成对比，但与近期对 AV₃Sb₅ 体系的理论理解趋于一致
与 KV₃Sb₅ 的对比：软模普适性
近期 KV₃Sb₅ 中也报道了类似的软声子模式（Phys. Rev. Lett. 2024），暗示声子软化是 AV₃Sb₅ 体系的普适特征
CsV₃Sb₅ 的软模行为与 KV₃Sb₅ 定性相似但细节不同：CsV₃Sb₅ 的软化幅度更大，可能与 Cs 原子更大的离子半径和更强的层间耦合有关
CsV₃Sb₅ 的 CDW 相图更复杂（2×2×4 相、向列性、手性输运），软模差异可能为理解这些额外序的起源提供线索
量子非谐效应与软模临界行为
SSCHA 计算表明量子非谐效应（零点振动）在 CsV₃Sb₅ 中起关键作用：即使在 0 K，量子涨落也能显著重正化声子频率
L₂⁻ 分支在 SSCHA 框架下在 TCDW 以下变为动力学不稳定（ω² < 0），对应 CDW 序参量的凝聚
97 K 时软模能量降至 ~3 meV，接近仪器分辨率极限——声称"过阻尼"（overdamped）需谨慎，实际可能是仪器分辨率限制
弹性峰的增长在 TCDW 以上已显示前驱特征，可能源于 TCDW 的微小空间分布（不均匀性）或临界涨落
四·附、讨论中的关键公式
阻尼比 Γ/E₀：当 Γ/E₀ > 1 时声子过阻尼（overdamped）。CsV₃Sb₅ 中软模在 97 K 时 Γ/E₀ 接近 1，但尚未达到过阻尼条件。
弹性散射强度与 CDW 序参量：I_elastic(T) ∝ |Δ(T)|² ∝ (TCDW − T)^(2β)，β 为临界指数（mean-field 预期 β=1/2）。
非谐谱函数：A(q,ω) = −(1/π) Im G(q,ω)，其中 G(q,ω) 为 SSCHA 声子 Green 函数。谱函数峰位对应重正化声子频率，峰宽对应声子寿命。
温度依赖的谱权重：S(Q,ω,T) = χ″(Q,ω) [1+n(ω,T)]。低温下 [1+n(ω,T)] 因子压低高频声子强度，增强低频弹性散射。
声子结构因子：S(Q,E) = Σⱼ |Q·eⱼ|²/ωⱼ · Fⱼ(E)。散射截面正比于 |Q·eⱼ|²，即波矢与声子本征矢的内积平方，解释了不同 Brillouin 区间结构因子的巨大差异。
五、结论与意义
核心结论
CsV₃Sb₅ 中 CDW 由 M-L 方向的软声子模式驱动——非谐晶格不稳定性为主，电子-声子耦合为辅
L 点软模能量从 300 K 的 ~8 meV 单调降至 95 K 的 ~3 meV，ω₀²(T) 与温度呈线性关系，符合非谐软模理论
EPC 线宽在 TCDW 附近无明显增强，说明电子关联非 CDW 主要驱动力，修正了"纯电子驱动 CDW"的主流观点
结构因子计算揭示此前实验未观测到软化的原因：<10% 的 Brillouin 区在 L 点有可观结构因子
实验-计算-理论三角验证框架为研究 kagome 体系 CDW 机制提供了系统性方法论，可推广至其他量子材料体系
六、支撑信息 (Supplementary Information)
S1 · 非谐谱函数
图 S1 | 非谐谱函数。Q=(3.5, 2.5, 0) 处 240 K、180 K、120 K 和 100 K 的非谐谱函数 A(q,ω)。峰位随温度降低向低能移动，展示声子软化的理论预期。
S1 要点
SSCHA 非谐谱函数在 240 K 至 100 K 间显示清晰的峰位下移，与实验 IXS 谱的软化趋势一致
谱函数峰宽随温度降低而增大，反映非谐散射增强导致的声子寿命缩短
S2 · L₂⁻ 模式避免交叉与杂化
图 S2 | L₂⁻ 模式的避免交叉与杂化。(a) 温度依赖的声子频率演化，展示两支 L₂⁻ 模式在 ~5 meV 处的避免交叉。(b) 模式杂化导致的谱权重转移。
S2 要点
在 L 点附近，两支 L₂⁻ 模式在 ~5 meV 处发生避免交叉，频谱权重在模式间转移
避免交叉的缺失是 M 点软化行为与 L 点不同的主要原因之一
S3 · 计算结构因子
图 S3 | 不稳定声子分支在不同 M 点的计算结构因子。(H,K,0)、(H,K,1)、(H,K,2) 和 (H,K,3) 面的代表性 M 点结构因子。
S3 要点
结构因子在不同 Brillouin 区差异极大，解释了为什么只有特定 Brillouin 区能观测到软模
8 个面板对应 8 个代表性 M 点，直观展示结构因子对 Brillouin 区选择的敏感性
S4 · 弹性强度温度依赖
图 S4 | CDW 转变温度附近弹性散射强度的温度依赖。归一化弹性强度在 TCDW 以下急剧增长，显示微小的滞后回线。
S4 要点
弹性强度在 TCDW 附近显示 ~1-2 K 的滞后，暗示 CDW 相变具有弱的一级相变特征
TCDW 以上弹性强度非零，可能源于临界涨落或样品中 TCDW 的微小空间分布
S5 · 声子能量温度依赖（Γ-M 中点）
图 S5 | Q=(0.75, 0.25, 3)（Γ-M 中点）处六支拟合声子模式的能量温度依赖，97-302 K。
S5 要点
远离 M 点向 Γ 移动时，声子软化幅度显著减小，说明 CDW 不稳定性局域在 M-L 线附近
六支模式中仅有一支显示明显软化，其他模式能量基本不变
S6 · M 和 L 点实验谱对比
图 S6 | M 点和 L 点各温度下实验 IXS 谱与仅含 Bose 因子温度依赖的模拟谱对比。偏离表明存在超出 Bose 因子的本征软化。
S6 要点
若声子仅受平凡 Bose 因子温度依赖影响，模拟谱与实验谱应一致
L 点实验谱在低能区的额外谱权重直接证明本征声子软化的存在
S7 · 拟合参数温度依赖
图 S7 | 软模（Phonon B）的拟合能量 (a)、宽度 (b) 和阻尼比 (c) 的温度依赖。
S7 要点
162 K 处 L 点宽度数据点偏离趋势，原因是软模峰值与 ~5 meV 处的强声子模式重叠，增大了拟合不确定性
阻尼比 Γ/E₀ 在 TCDW 接近时因 E₀ 减小而增大，但 Γ 本身无异常增强
S8 · 面内漫散射 mapping
图 S8 | 95 K 时沿 Q=(0, 0, L) 方向的面内漫散射 mapping。展示 CDW 关联的二维特征。
S8 要点
漫散射强度沿 L 方向（出平面方向）变化，L=0.5 整数附近最强
面内漫散射分布呈各向异性，反映 CDW 关联的二维特征和 kagome 晶格的对称性
S9 · ID28 低分辨 IXS 谱
图 S9 | ESRF ID28 线站 3 meV 分辨率的 IXS 谱。(a) 300 K 原始数据。(b) 基于 300 K 数据模拟的 94 K 谱。(c) 实测 94 K 谱。
S9 要点
ID28 的低分辨率（3 meV）测量覆盖了更宽的动量空间范围（Γ-M-L-A 路径）
低分辨率下软模表现为宽化的谱权重包络，与高分辨率 DHO 拟合结果一致
谱权重从高能向低能转移的趋势在不同分辨率下均可复现，增强了结果的可靠性
S10 · EPC 动量依赖
图 S10 | M-L 方向电子-声子耦合线宽的动量依赖。SSCHA 100 K 非谐声子色散（黑色虚线）与 EPC 贡献线宽（彩色标度）对比。
S10 要点
EPC 线宽在 M-L 路径上具有显著的动量依赖，L 点附近 EPC 贡献最大
但 EPC 线宽的绝对值（~0.5-1 meV）远小于实验观测到的总展宽（~2 meV），证实非谐贡献主导
七、表格
表格
未找到表格文件。
八、同行评审 (Peer Review)
Version 0 — 初始投稿审稿意见
审稿人 #1 (正面)
总体评价：认为本文报道了重要的实验数据，对 kagome 材料和 CDW 领域具有重要参考价值。
建议 1：更深入地探讨 Fermi 面嵌套 vs 电子-声子耦合作为 CDW 起源的二分法是否有意义。
建议 2：面向一般读者解释 IXS 强度的两个决定因素（矩阵元效应和声子谱函数），以及本文如何克服矩阵元效应。
建议 3：Fig 4 中的误差棒可能误导读者——误差棒仅代表频移而非展宽的不确定度，应同时展示展宽的误差棒。
审稿人 #2 (拒稿)
总体评价：认为软模证据不充分，无法支持 Nature Communications 发表。
主要问题 1：软模的识别完全依赖拟合，在 127 K 和 102 K（M 点）及 197 K 和 127 K（L 点）的原始谱中软模几乎无法与背景区分，质疑其为拟合伪影。
主要问题 2：应提供与 Subires et al. 研究的直接对比，在相同条件下测量两个 Brillouin 区，验证结构因子差异是否确实导致软模可探测性差异。
主要问题 3：EPC 线宽趋势不一致——L 点线宽在低温下保持极小且误差小，与预期的大展宽矛盾，质疑拟合可靠性。
次要问题：声称"<10% 的 BZ 有可观结构因子"缺乏数据支持；SI 引用存在多处 "??" 缺失。
审稿人 #3 (正面)
总体评价：认为软模发现意义重大，鉴于此前多个理论工作试图解释软模缺失，本文为 CDW 电子-声子机制提供了关键证据，适合 Nature Communications。
主要问题：软模在原始数据中不够明显，虽经 Fig 3 详细分析后变得可信，但对拟合参数的可靠性存疑——误差棒是否可靠？若释放更多固定参数，误差棒如何变化？建议在 SI 中展示所有模型参数的温度依赖。
次要问题 1：量子非谐效应为何起关键作用？
次要问题 2：Fig 1 中非谐和 EPC 导致的声子展宽在彩色图中难以辨认，这些展宽对理解本文有多重要？与实验拟合阻尼值如何比较？
次要问题 3：97 K 数据声称"过阻尼"可能过于强烈——声子能量从 3 meV 突降至 0 meV 仅 5 K 内，应审慎评估或软化措辞。
次要问题 4：理论与实验的连接需更好解释——L=0.5 处理论预期两个模式，实验仅观测到一个，对应哪个模式？另一模式为何未观测到？
作者回复与修改
主要修改内容
添加 Bose 因子校正后的原始谱图（Fig 3 新版），使声子软化在原始数据中直接可见，回应审稿人 #2 和 #3 的核心关切。
补充 SI 中所有拟合参数（能量、线宽、振幅）的温度依赖，明确标注固定参数和自由参数，回应审稿人 #3 的拟合可靠性问题。
增加与 Subires et al. 的详细对比讨论，展示结构因子计算如何解释此前实验的"阴性"结果。
修正所有 SI 引用 "??" 错误，补充结构因子统计分布数据。
软化"过阻尼"措辞，明确说明 97 K 时 Γ/E₀ 接近 1 但尚未达到过阻尼条件。
增加量子非谐效应的物理解释，补充 SSCHA 计算中零点振动对声子频率重正化的定量分析。
Version 1 — 修订稿审稿意见
审稿人 #1 (接受)
作者已成功解决我提出的问题，推荐在 Nature Communications 发表。
审稿人 #2 (接受)
在回复信和修订稿中，作者详细回应并澄清了我提出的所有疑虑。特别是 Bose 因子校正后的谱图（Fig 3）现在直接展示了声子模式演化。因此，我支持在 Nature Communications 发表。
审稿人 #3 (接受，附次要建议)
作者做出了显著努力改进稿件，Bose 校正原始数据的加入使得声子温度依赖变化变得明确。鉴于此前大量工作将 CDW 非常规性质归因于 van Hove 奇点嵌套，本文提出的电子-声子机制具有重要意义，适合 Nature Communications。
剩余次要问题 1：误差棒是否在所有自由参数均允许变化时确定？有时人们错误地在固定参数下计算误差棒，导致误差棒过小。
剩余次要问题 2：弹性峰增长是本征的还是源于 TCDW 分布？即使极小部分样品有稍高的 TCDW，也可能在平均 TCDW 之上产生大的弹性信号。如果是本征的，则同时存在有序-无序和声子软化过程，与 Ref. 45 一致。
剩余次要问题 3：KV₃Sb₅ 中近期也报道了软声子，与 CsV₃Sb₅ 是否类似？若不同，差异是否与 CsV₃Sb₅ 更复杂的相图（2×2×4 相、向列性、手性输运）有关？
剩余次要问题 4：SI 中仍有多处 "??" 引用。
CsV₃Sb₅ · Kagome 金属 · CDW 软模 · 非弹性 X 射线散射 · SSCHA 非谐声子 · 电子-声子耦合
