ScV6Sn6:对称性选择性 CDW 稳定化
导读:本文利用单轴应变结合高分辨 X 射线衍射,展示了各向异性晶格形变如何通过打破 Kagome 晶格旋转对称性,选择性地稳定和增强 ScV6Sn6 中的电荷密度波序。三态 Potts 模型 + Monte Carlo 模拟揭示了阻挫晶格中序-无序相变的对称性调控机制。
Kagome · CDW · Uniaxial Strain · Potts Model
Symmetry-Selective Stabilization of Charge-Density Wave in ScV6Sn6
一、前言背景
Kagome 金属中的 CDW 与几何阻挫
Kagome 晶格中过渡金属原子的二维三角排列具有几何阻挫特性,天然映射到广义三态 Potts 模型,在临界温度 TC 处发生序-无序相变。
AV3Sb5 家族中 CDW 与超导共存竞争:CsV3Sb5 中 ~1% 压缩应变即可使 CDW 能隙增强近 3 倍,单轴应变作为对称性破缺手段可选择性调控 CDW 态。
ScV6Sn6 是理想平台:HfFe6Ge6 型六方结构(P6/mmm),V 原子构成 Kagome 层,SnT-Sc-SnT 三聚体沿 c 轴排列,几何阻挫+晶格对称性+电子-声子耦合竞争产生多重电荷有序态。
ScV6Sn6 在 TCDW ~ 92 K 发生 sqrt3 x sqrt3 x 3 超结构相变,由 SnT 原子面外集体位移驱动的一级相变,SnT-SnT 键距在相变处发生分裂,同时存在竞争的 q2 = (1/3, 1/3, 1/2) 不稳定性。
与静水压不同(通常压制 CDW),单轴应变引入方向各向异性,打破三重旋转对称性->作为 Ising CDW 序参量的共轭场->择优选择调制方向。
本文核心科学问题
单轴应变如何通过打破 Kagome 晶格旋转对称性(C6->C2),选择性地稳定和增强 q1 CDW 态?
CDW 稳定性增强主要由非谐晶格效应(声子软模稳定化)还是各向异性 Ising 耦合(SnT 三聚体面内有序化)驱动?
三态 Potts 模型中应变依赖的 Ising 耦合能否定量复现实验观测的 TCDW 连续增强?
结合高分辨 X 射线衍射+单轴应变装置+DFT 声子计算+Monte Carlo 模拟的多尺度研究方法。
三态 Potts 模型与序-无序相变的关系
ScV6Sn6 中 SnT-Sc-SnT 三聚体可沿三个等效方向位移(对应 q1 的三个正交畴),天然构成三态变量,映射到三态 Potts 模型。
无应变下三个方向完全简并,体系在 TCDW 处发生一级相变,随机选择一个方向长程有序,形成三个正交 merohderal 畴。
单轴应变打破三个方向简并->引入等效"磁场"倾向 CDW 序参量->降低对称性至正交 Cmmm->TCDW 升高。
该体系为阻挫 Kagome 晶格中序-无序相变的研究提供了理想实验平台,三态 Potts 模型是描述 CDW 相变的最小模型。
二、研究方法
实验手段:单轴应变下的原位 X 射线衍射
高分辨 X 射线衍射(同步辐射 ESRF):BM01 束线(17.5 keV,无应变表征)和 ID15B 束线(30 keV,应变下测量),PILATUS 2M 和 EIGER2 X 9M CdTe 探测器。
单轴应变装置(Razorbill CS200T):沿 [H00] 和 [HH0] 方向压缩,应变幅度 0-1.5%,样品厚度 ~70 um,TCDW+10 K 处施加应变后降温测量。
变温测量:80-200 K 范围追踪 CDW 超晶格峰强度,线性拟合近相变区强度-温度曲线确定 TCDW。
结构精修:CrysAlisPro 数据处理,Jana2006 和 SHELX 精修,ISODISTORT 对称性模式分析,确认六方 P6/mmm->正交 Cmmm 的 GM5+ 不可约表示畸变。
竞争 CDW 波矢:q1 凝聚为长程有序,q2 以漫散射形式存在
TCDW 确定方法:序参量导数的最大值对应的温度
理论方法:DFT 声子谱与 Monte Carlo 模拟
DFT(VASP):PBE-GGA 泛函,六方 16x16x8 k 网格,正交 16x8x8 k 网格,Phonopy 有限位移法计算声子谱(3x2x3 超胞)。
三态 Potts 模型 + 各向异性 Ising 耦合:Ji(d) = Ji0 * exp[lambda*(1 - d/d0)],lambda = 25,J1 = -6.56 K,J2 = +1.74 K,J3 = -1.03 K(未应变参数取自 Alvarado et al.)。
Monte Carlo 模拟:90x90x1 超胞,模拟退火 20 步,每步 50 倍原子数 MC 移动,DISCUS 程序计算漫散射,sigmoid 函数拟合确定 TCDW。
应变编码:通过正交晶格参数变化(a_orth = sqrt3 x b_hex, b_orth = b_hex for [H00];a_orth = sqrt3 b_hex, b_orth = x b_hex for [HH0])将应变映射到最近邻和次近邻耦合参数。
三态 Potts 模型哈密顿量:J1/J2/J3 分别为第一/第二/第三近邻耦合,delta 为 Kronecker delta 函数
应变依赖的 Ising 耦合:J_i^0 为未应变参考值,d 为原子间距,d0 为参考间距,lambda 为衰减常数
CDW 序参量:phi 在三个等价方向中取非零值,sqrt(3) 归一化因子
正交晶格映射与应变参数化
六方->正交对称性降低:P6/mmm(a=b, gamma=120deg)-> Cmmm(ao=sqrt3 ah, bo=bh, co=ch),三个六方方向分裂为两个不等价正交方向。
[H00] 压缩:施加沿 ao 方向(长正交轴),导致 ao 显著压缩、bo 膨胀、c 几乎不变,对应 J1a 增强、J1b 减弱。
[HH0] 压缩:施加沿 bo 方向(短正交轴),与 [H00] 方向正交,产生类似的各向异性晶格响应但方向互换。
DFT 验证:计算不同正交应变下的声子谱,提取虚频软模的稳定化程度,判断声子驱动 vs 序-无序驱动的 CDW 增强机制。
三、实验结果
无应变下的 CDW 相变特征
ScV6Sn6 在 TCDW ~ 95 K 发生一级相变:(i) 晶格参数和 CDW 超晶格峰强度在 TCDW 处突变;(ii) SnT 各向异性位移参数比 U001/U100 随降温显著增大,反映面外位移增强。
SnT-SnT 键距在相变以下发生分裂,表明 SnT-Sc-SnT 三聚体沿 c 轴的集体位移冻结为长程有序的 sqrt3 x sqrt3 x 3 调制结构。
q2 = (1/3, 1/3, 1/2) 漫散射强度随降温增长但不凝聚为长程有序相,q1 和 q2 两个不稳定性在能量上近乎简并,竞争决定了基态选择。
声子谱:无应变晶格中沿 L-A-H 路径存在几乎平坦的虚频支,对应 SnT 原子面外振动,q2 不稳定性与平坦光学声子支软化密切相关。
图 1:(A) 高温 P6/mmm 晶体结构,V 原子(红色)形成 Kagome 晶格,SnH(浅灰)占据蜂窝位,SnT(蓝)和 Sc(紫)形成三聚体;(B) CDW 诱导的 SnT-Sc-SnT 三聚体集体位移;(C) 布里渊区;(D) 晶格参数与 CDW 峰强度随温度演化;(E) SnT 各向异性位移参数比;(F-G) SnT-SnT 和 SnT-Sc 距离的温度依赖性。
单轴应变:对称性破缺驱动的 CDW 增强
应变施加后六方 P6/mmm 对称性降至正交 Cmmm(GM5+ 不可约表示),ao = sqrt3 ah, bo = bh, co = ch,ISODISTORT 分析确认单一对称性允许的正交畸变。
[H00] 压缩(沿长正交轴):ao 轴显著压缩,bo 轴适度膨胀,c 轴几乎不变 -> TCDW 从 92 K 单调提升至 ~105 K(-0.93% 应变)。
[HH0] 压缩(沿短正交轴):产生类似的各向异性晶格响应和 TCDW 增强 -> 两种方向应变效应定性相似,仅存在微弱的面内各向异性。
CDW 超晶格峰沿应变方向和垂直方向温度依赖性几乎相同,证实均匀应变场选择性稳定了三个正交 merohderal 畴中的单一畴。
图 2:(A) 沿 c 方向晶体结构视图,展示三个正交 merohderal 畴及 [H00]/[HH0] 方向应变;(B) 上:应变施加后原子位移示意图,下:单轴应变装置示意图;(C) 特征 Bragg 峰在应变下的位移;(D) 正交晶格参数和体积随 [H00] 应变的演化;(E) HHL 倒空间图;(F) TCDW 随应变的变化。
应变下的声子谱:软模仅中等程度稳定化
无应变:声子谱沿 L-A-H 路径存在近乎平坦的虚频支,由 SnT 原子面外振动主导;正交晶胞中双倍原子数揭示了第二个 SnT 衍生不稳定模式。
[H00] -1.5% 压缩:虚频软模在 L 和 H 点仍不稳定,A 点变为稳定,但 H 点仍为优势 CDW 序矢量(q1 = (0, 1/3, 1/3) 在正交基下)。
[HH0] -1.5% 压缩:声子谱演化与 [H00] 定性相似,未出现强烈的方向依赖性各向异性修改 -> 排除了各向异性电子-声子耦合变化主导 TCDW 增强的情景。
关键结论:声子软模仅中等程度稳定化,不足以解释 ~13 K 的 TCDW 提升,TCDW 增强主要由非声子机制(即 SnT 三聚体面内有序化)驱动。
图 3:原子投影声子色散:(A) 无应变;(B) [H00] 方向 -1.5% 压缩;(C) [HH0] 方向 -1.5% 压缩。Sc/V/SnH/SnT 原子振动模式分别着色,虚频软模以负频率表示。
Monte Carlo 模拟:三态 Potts 模型复现 TCDW 演化
应变依赖的 Ising 耦合参数 Ji(d) 通过指数衰减函数编码,lambda = 25 时模型复现了实验观测的 TCDW 斜率。
[H00] 和 [HH0] 方向压缩的 MC 模拟 TCDW 与实验值高度吻合,验证了应变通过增强面内 Ising 耦合强度来稳定 q1 CDW 的物理图像。
无应变时 MC 模拟显示显著的一级相变特征;压缩应变下三聚体无序逐步抹平锐利转变,序参量向更连续的过渡演化,符合三态 Potts 模型的预期。
Warren-Cowley 短程序参数 alpha 和最终 MC 能量在应变下均展示与 TCDW 一致的演化趋势,证实了序-无序相变本质。
图 4:(A) 正交晶胞中最近邻和次近邻 Ising 耦合 J1a-J3b 的实空间定义;(B) 实验与 MC 模拟的 TCDW 随应变演化对比,[H00] 和 [HH0] 方向模拟结果均成功复现了 CDW 的连续增强。
四、对比分析
[H00] vs [HH0] 应变方向:对称性选择性的面内各向异性
[H00] 压缩:压缩沿长正交轴(ao),TCDW 从 92 K 提升至 ~105 K,ao 压缩显著、bo 膨胀、c 几乎不变,对应压缩沿长轴方向。
[HH0] 压缩:压缩沿短正交轴(bo),TCDW 提升幅度略小于 [H00],但两种方向应变依赖性相似,仅存在微弱的面内各向异性。
两种方向均导致六方->正交结构畸变(Cmmm),但正交轴的长短方向互换 -> CDW 增强程度接近 -> 对称性降低本身而非方向细节是主导因素。
若电子-声子耦合各向异性主导 TCDW 增强,[HH0] 方向应产生更强的效果(因为压缩方向不同),但实验观测不支持 -> 排除了电子-声子矩阵元各向异性为主因。
与静水压和化学掺杂的对比
静水压(Zhang et al.):压制 ScV6Sn6 的 CDW,11 GPa 以下无超导电性 -> 各向同性体积压缩不利于 CDW 稳定。
化学掺杂(Lu/Y 替代 Sc):TCDW 下降,CDW 被压制 -> 化学压力效应与静水压类似,但单轴应变产生相反效果。
CsV3Sb5 中单轴应变:~1% 压缩使 CDW 能隙增强近 3 倍,但 ScV6Sn6 中应变提升 TCDW 而非仅增强能隙 -> 序参量直接被稳定化而非仅增强幅度。
关键区别:ScV6Sn6 中 TCDW 增强是应变驱动的序-无序相变调控,而非简单的能隙调制或 Fermi 面嵌套效应。
与 AV3Sb5 家族 CDW 对称性调控的对比
AV3Sb5:CDW 涉及 V 原子 Kagome 层内 2x2 电荷有序,手性电荷序(星形 vs 三角形图案)已被实验证实,单轴应变可选择性稳定其中之一。
ScV6Sn6:CDW 涉及 SnT 三聚体面外位移而非 Kagome 层内调制,CDW 波矢更长(sqrt3 x sqrt3 x 3),Potts 序参量描述不同畴的选择。
共性:单轴应变作为对称性破缺的共轭场,在两种体系中均能选择性稳定特定 CDW 畴,揭示了 Kagome 体系中应变-电荷有序的普适耦合机制。
五、讨论
CDW 稳定化的微观机制:序-无序相变 vs 声子软模
声子计算表明虚频软模在应变下仅中等程度稳定化,不足以解释 ~13 K 的 TCDW 提升 -> 排除纯粹声子驱动的增强机制。
c 轴晶格参数对应变无响应 -> SnT 面外振动与 SnT pz 轨道相互作用不受面内压缩影响 -> 电子-声子耦合非主要调控因素。
主要驱动力:SnT-Sc-SnT rattling 链在阻挫 Kagome 三角晶格上的面内有序化,通过"无序中产生有序"(order-from-disorder)机制实现。
三态 Potts 模型中应变增强了面内 Ising 耦合参数,等价于将等效 Ising 磁场施加于 CDW 序参量,降低了对称性等价态之间的简并度。
该机制与 FeGe 中观察到的阻挫 CDW 和准长程键取向序一致,表明阻挫 Kagome 晶格中序-无序相变是 CDW 形成的普适特征。
配分函数与自由能:Z 为所有自旋构型的 Boltzmann 求和,F 决定平衡态和相变行为
声子软化机制:chi(q) 为电子极化率,g 为电子-声子耦合,虚频表示 omega_q^2 < 0 的晶格不稳定性
CDW 能隙:TCDW 以下序参量以 sqrt(TCDW-T) 形式增长,典型二级相变行为
对称性降低:P6/mmm->Cmmm 通过 GM5+ 不可约表示,C6 旋转对称性降至 C2
应变能:C_ijkl 为弹性张量,epsilon_ij 为应变张量,面内压缩产生各向异性弹性能
六方->正交畸变对 Kagome 电子结构的深远影响
极小的应变(<1%)即足以打破 P6/mmm 对称性,旋转对称性从 C6 降至 C2,M 点 vHS 和 K 点 Dirac 点的简并被解除。
三个等价方向分裂 -> 择优选择单 q CDW 波矢 -> 不利于手性相但有利于向列相 -> 应变作为向列性的共轭场。
对 (Cs,K,Rb)V3Sb5 体系的影响:其基态在手性和向列序之间争论,单轴应变可能驱动手性->向列相变,或揭示振幅和相位随应变连续演化的新物理。
单轴应变作为一种干净、可控的对称性选择性调控手段,为阻挫晶格中序-无序相变的研究提供了通用实验范式。
三态 Potts 模型的普适性与局限
三态 Potts 模型成功捕捉了 ScV6Sn6 中 CDW 相变的核心特征:一级相变、应变依赖的 TCDW 增强、序-无序转变。
模型局限:(i) 未包含电子-声子耦合的微观细节;(ii) 声子软模效应被简化为有效耦合参数;(iii) 量子涨落效应被忽略。
未来改进:(i) 引入动态声子自由度;(ii) 考虑长程 Coulomb 相互作用;(iii) 构建包含电子自由度的完全量子 Potts 模型。
该模型框架可推广至其他阻挫 Kagome 体系,如 FeGe、Mn3Sn、Co3Sn2S2 等,为对称性破缺驱动的序参量调控提供统一理论描述。
六、总结
核心结论
单轴应变通过打破 Kagome 晶格 C6 旋转对称性(六方 P6/mmm->正交 Cmmm),选择性地稳定和增强 ScV6Sn6 中 q1 = (1/3, 1/3, 1/3) 的 CDW 态。
[H00] 和 [HH0] 方向压缩均将 TCDW 从 ~92 K 提升至 ~105 K(-0.93% 应变),两种方向效应定性相似,仅存在微弱的面内各向异性。
CDW 增强主要由 SnT-Sc-SnT rattling 链的面内有序化(order-from-disorder 机制)驱动,而非声子软模稳定化或电子-声子耦合变化。
三态 Potts 模型 + 应变依赖 Ising 耦合(Ji(d) = Ji0 * exp[lambda*(1 - d/d0)], lambda = 25)在 Monte Carlo 模拟中成功复现 TCDW 连续演化。
本研究建立了单轴应变作为阻挫 Kagome 晶格中对称性选择性序-无序相变调控的通用范式,为多层次 CDW 竞争体系中序参量的解耦提供了新途径。
支撑信息 (Supporting Information)
S1 · 电子能带结构
六方相电子结构
六方 P6/mmm 相电子能带结构显示 Sn pz 轨道特征集中在费米能级附近的 A-L-H 路径,反映了三角 SnT 原子面外轨道在该能量窗口的贡献。
V 原子 Kagome 层贡献的 van Hove 奇点 (vHS) 位于 M 点附近,Dirac 锥位于 K 点,是 CDW 不稳定性起源的电子结构基础。
pz 轨道在费米能级附近的高态密度表明 SnT 原子面外振动与电子态之间存在强耦合,为声子软化和 CDW 形成提供了电子驱动力。
应变与 CDW 重构下的能带演化
[H00] 和 [HH0] 方向 1.5% 压缩应变下的展开能带结构显示,展开谱权重紧密跟随正常态 (NS) 色散,仅在高对称点附近出现轻微展宽,无显著能量位移。
pz 谱权重分布在应变下基本保持不变,表明正交畸变未显著改变低能电子结构,排除了 Fermi 面拓扑变化驱动 TCDW 增强的机制。
sqrt(3) x sqrt(3) x 3 CDW 重构超胞的能带结构显示额外的能带折叠和谱权重再分布,但整体能带色散保持与高对称相一致,证实电子结构的鲁棒性。
图 S1:电子能带结构。(a) 六方相 Sn pz 轨道投影能带;(b-c) [H00] 和 [HH0] 方向 1.5% 压缩应变下的展开能带结构;(d) sqrt(3)xsqrt(3)x3 CDW 重构超胞展开能带,红线为正常态能带参考。
S2 · 声子色散谱
高斯展宽依赖的声子软化
通过增加 Gaussian 展宽宽度模拟电子温度效应,观察到 SnT 软支对应的虚频声子模幅值随展宽增大逐步减小,证实 CDW 不稳定性的电子辅助特征。
即便在最大展宽值下,不稳定支在倒空间广阔区域内依然存在,表明底层晶格不稳定性是鲁棒的,并非纯粹由电子驱动。
[H00] 和 [HH0] 方向压缩的不稳定性热抑制行为定性相同,说明正交畸变的面内取向对热抑制效应无显著方向依赖性。
应变依赖的声子谱演化
无应变参考声子谱显示沿 L-A-H 路径的特征性近乎平坦虚频支,由 SnT 原子面外振动主导。
沿 [H00] 和 [HH0] 方向增加压缩应变,虚频支在布里渊区大部分区域持续存在,证明正交畸变不破坏底层晶格不稳定性。
最大压缩应变下 A 点附近出现部分稳定化,[H00] 方向效应略强于 [HH0],与 [H00] 对应沿较长正交轴的压缩一致。
两种应变方向声子谱定性相似,无强各向异性修改,排除了各向异性电子-声子耦合变化主导 TCDW 增强的情景。
图 S2:声子色散谱。(a-e) [H00] 方向压缩,增加 Gaussian 展宽;(f-j) [HH0] 方向压缩,增加 Gaussian 展宽。颜色编码表示不同原子的振动权重。
S3 · 声子线宽分析
虚频声子模的应变演化
无应变下沿 L-A-H 路径的虚频支几乎完全平坦,对应的声子线宽极窄,表明 SnT 面外振动模式高度局域化,群速度接近零。
施加压缩应变后,虚频支在 L 和 H 点仍保持不稳定,但 A 点在高应变下变为稳定,声子线宽略有展宽,反映应变对晶格动力学的微弱修饰。
正交晶胞中双倍原子数揭示了第二个 SnT 衍生不稳定模式,在 L 和 H 点虚频,A 点稳定,表明 CDW 不稳定性优先选择 H 点作为序矢量。
声子线宽与电子-声子耦合的关系
声子线宽 Gamma_q = 2*pi*omega_q*lambda_q 正比于电子-声子耦合常数 lambda_q,线宽的应变依赖性可直接探测 EPI 矩阵元的变化。
虚频声子模的线宽在 [H00] 和 [HH0] 应变下均未出现显著增大,表明电子-声子耦合矩阵元未因面内压缩而增强。
这一观察排除了各向异性 EPI 变化主导 TCDW 增强的机制,支持序-无序相变(面内 Ising 耦合增强)为主要驱动力。
图 S3:声子色散随应变演化。(a) 无应变参考谱;(b-d) [H00] 方向压缩增加;(e-g) [HH0] 方向压缩增加。虚频支在 L 和 H 点持续不稳定,A 点部分稳定化。
S4 · 电子-声子耦合分析
EPI 矩阵元与能带结构的关系
电子-声子耦合强度由矩阵元 g_{mn,nu}(k,q) = 决定,其中 delta_{nu,q} V 为声子模式 nu 对自洽势的扰动。
SnT pz 轨道与面外振动模的耦合集中在费米能级附近,pz 谱权重在应变下基本不变,表明 EPI 矩阵元对应变不敏感。
展开能带结构显示,正交畸变未在费米能级处引入新的能带交叉或 van Hove 奇点,Fermi 面嵌套条件未因应变显著改变。
EPI 与序-无序机制的区分
若 EPI 变化主导 TCDW 增强,[HH0] 方向压缩应产生比 [H00] 更强的效应(因压缩方向不同导致不同的矩阵元修饰),但实验观测不支持。
c 轴晶格参数对应变无响应,表明 SnT 面外振动与 pz 轨道的相互作用不受面内压缩影响,进一步排除了 EPI 作为主要调控因素。
TCDW 增强的能量尺度 (~13 K ~ 1.1 meV) 远小于典型声子能量 (~10-20 meV),支持 Ising 耦合增强而非电子-声子耦合变化为主导机制。
图 S4:电子-声子耦合证据。(a) 六方相 Sn pz 轨道投影能带;(b-c) [H00] 和 [HH0] 方向 1.5% 压缩应变下的展开能带,pz 谱权重分布基本不变;(d) CDW 重构超胞能带,整体色散保持。
S5 · 各向异性位移参数
SnT 原子 ADP 的温度依赖性
SnT 原子各向异性位移参数 (ADP) 的比值 U001/U100 在 TCDW 以上接近 1,表明面外和面内热振动幅度相当。
TCDW 以下 U001/U100 显著增大,反映 SnT 原子沿 c 轴的大幅面外位移,与 CDW 诱导的三聚体集体位移图像一致。
ADP 的突变行为与晶格参数的一级相变特征同步,证实 CDW 相变是结构序参量驱动的。
ADP 比值与 CDW 序参量的关联
U001/U100 可作为 CDW 序参量的实验探针:比值越大,面外调制幅度越强,长程 CDW 有序度越高。
SnT-SnT 键距在相变以下发生分裂,与 ADP 比值的增大同步,表明 SnT-Sc-SnT 三聚体沿 c 轴集体位移冻结为长程有序调制结构。
应变下 ADP 比值的演化进一步证实了应变对 CDW 序参量的稳定化作用,序参量幅值随应变增大而增强。
图 S5:各向异性位移参数。(A) SnT 原子 ADP 椭球体示意图;(B) U001/U100 比值随温度演化,TCDW 以下显著增大;(C) 不同温度下 SnT 原子位移椭球的取向变化。
S6 · 温度依赖的晶格参数
应变诱导的晶格演化
施加面内单轴应变诱导高温相均匀结构畸变,特征 Bragg 反射发生位移,伪六方设定下 alpha 和 beta 角保持接近 90 度,gamma 角显著偏离 120 度。
正交晶格参数 ao 和 bo 对应变响应强烈:ao 压缩、bo 膨胀,c 轴几乎不变,体积随应变增大而减小。
[H00] 和 [HH0] 方向应变导致的晶格参数演化在实验误差范围内等价,但正交晶胞中方向相反:[H00] 压缩沿长正交轴,[HH0] 压缩沿短正交轴。
零应变参考点的确定
由于样品持续安装在应变装置中冷却,有限残余应变可能来自差热收缩和压电堆残余行程。
零应变条件通过线性外推精修晶格参数 ah 作为正交畸变 ao/sqrt(3)-bo 的函数来估计,与机械断裂样品后直接测量的值一致。
晶格参数的温度依赖性在 80-200 K 范围内平滑,仅 TCDW 处出现突变,证实 CDW 相变的一级特征。
图 S6:应变诱导晶格演化。(A) 伪六方设定下晶胞角度随面内单轴应变的演化;(B) [H00] 和 [HH0] 方向压缩下正交晶格参数对应变的依赖性。
S7 · CDW 调制波矢与结构表征
应变诱导的正交结构
单轴应变选择单一 merohderal 正交畴,其晶格与母体六方基的关系为 [ao, bo, co]^T = [[1,0,0],[1,2,0],[0,0,1]] * [ah, bh, ch]^T。
正交结构保持与母相相同的化学组成:每个化学式含一个 Sc、两个 SnH、一个 SnT,V 原子 Kagome 位置沿 c 方向的分裂可忽略。
ISODISTORT 对称性模式分析确认应变诱导畸变由单一 GM5+ 不可约表示描述,对称性从六方 P6/mmm 降至正交 Cmmm。
CDW 相的结构精修
群-子群分析允许多种正交子群:Cmmm, C222, Cmm2, Amm2,比较精修结果确认 Amm2 空间群给出最佳拟合 (GOF=1.19, RF=6.22%)。
Amm2 模型需要更多独立原子位点和精修参数,但各向异性位移参数的自由精修进一步改善 R 因子,清晰分辨 SnT 和 Sc 原子的大幅位移。
低对称性允许 V 原子 Kagome 位置沿 c 方向分裂,但对应位移在实验误差内可忽略,应变诱导 CDW 相与常压 CDW 结构具有高度相似性。
图 S7:应变诱导正交结构。(A) Cmmm 设定下的正交结构,沿 c 方向视图;(B) 六方原胞以红色轮廓标示;(C) Amm2 空间群精修的 CDW 调制,位移原子着色,非调制位置以黑色阴影显示;(D-E) 高温 Cmmm 相和低温 Amm2 相 Fcalc vs Fobs 对比。
S8 · 超胞构建与结构弛豫
DFT 超胞设置
六方相计算使用 16x16x8 k 网格,正交畸变使用 16x8x8 k 网格,CDW sqrt(3)xsqrt(3)x3 重构使用 9x9x3 k 网格。
所有结构在计算声子谱前进行了晶胞和内部位置的弛豫,CDW 重构结构取自文献报道的原子坐标,未进行内部原子位置弛豫。
谐波声子谱使用 Phonopy 有限位移法计算,正交畸变使用 3x2x3 超胞 (3x3x2 k 网格),展开能带使用 VASPKIT 后处理。
交换耦合参数定义
最近邻和次近邻 Ising 耦合参数在六方晶格上定义:J1 为第一近邻(面内最近邻方向),J2 为第二近邻,J3 为第三近邻。
有效交换哈密顿量:同类近邻对 (UU 或 DD) 在三聚体链中产生 -2J 的能量惩罚,异类对 (UD) 产生 +J 的能量增益。
符号约定:J1 为负值表示最近邻相互作用产生能量惩罚(反铁磁型 Ising 耦合),有利于同类三聚体构型的有序排列。
图 S8:超胞与交换耦合。(A) 六方晶格上交换耦合参数 J1, J2, J3 的实空间定义;(B-C) 有效交换哈密顿量示意图。"UU" 赝原子构型的总惩罚为 -2J1,"UN" 构型的总增益为 +J1。
S9 · 应变依赖的自由能
MC 最终能量与序参量的应变演化
最终 MC 能量在无应变下显示明显的一级相变特征,随应变增大,能量-温度曲线逐渐平滑,反映三聚体无序抹平锐利转变。
Warren-Cowley 短程序参数 alpha 对应变敏感,U-D 对的一阶 alpha 参数在应变增大时单调变化,与 TCDW 演化趋势一致。
积分峰面积 (Area) 在无应变下显示尖锐的 sigmoid 型转变,压缩应变下转变宽度 delta 增大,表明序参量向更连续的过渡演化。
三种探针的一致性
从积分峰面积、最终 MC 能量和 Warren-Cowley 短程序参数三种独立探针提取的 TCDW 在应变依赖性和绝对值上高度一致。
三种探针均显示 [H00] 和 [HH0] 方向压缩产生定性相似的 TCDW 增强,验证了对称性降低而非方向细节是主导因素的结论。
MC 能量和短程序参数的应变演化进一步证实了 CDW 的序-无序本质,与声子软模仅中等程度稳定化的 DFT 结果一致。
图 S9:最终 MC 能量、Warren-Cowley 短程序参数和积分峰面积作为模拟温度和应变的函数。(A) J1a 相互作用的 alpha 参数;(B) 最终 MC 能量;(C) 积分峰面积。
S10 · Monte Carlo 模拟细节
模拟设置与退火方案
MC 模拟在二维 90x90x1 超胞上进行,使用粗粒化方法:Sc-Sn-Sc 三聚体被分配为上 (U)、下 (D) 或中性 (N) 三种构型之一。
模拟退火:20 个温度步,起始温度 Tstart = max(100 K, 3xTfinal),每步尝试 MC 移动次数 = 50 x 总原子数。
每次 MC 移动随机交换两个不同状态的赝原子,接受概率 = min(1, exp(-DeltaE/TMC)),确保细致平衡。
系统假设沿 c 方向完全有序,赝原子标签取决于图 1B 所示构型的底层原子。
漫散射计算与 TCDW 确定
MC 构型扩展为 90x90x3 超胞,使用 DISCUS 程序计算单晶漫散射,NUFFT 傅里叶算法在 181x181x31 体素网格上计算。
5 个独立构型的漫散射取平均,按 2/m Laue 对称性(c 为唯一轴)对称化,以不抑制漫散射信号中的局域正交畸变。
TCDW 通过 sigmoid 函数拟合积分卫星峰面积 Isat(T) 确定:Isat = Ilow + (Ihigh-Ilow)/(1+exp(-(T-TCDW)/delta)) + alpha*(T-150K)。
无应变 + TMC = 0 K 的 hk4.33 层漫散射显示清晰的 sqrt(3)xsqrt(3) 衍射特征,与 Alvarado et al. 结果一致。
图 S10:Monte Carlo 模拟细节。(A) 无应变 TMC=0 K 的 hk4.33 层计算漫散射;(B) h1/3l 层漫散射,l 方向分辨率受限于粗模拟网格;(C) 0% 应变下 sigmoid 函数拟合积分卫星峰面积。
S11 · 序参量演化
CDW 峰强度的温度与应变依赖性
沿应变方向和垂直方向的 CDW 超晶格峰在各应变值下展示几乎相同的温度依赖性,证实均匀应变场选择性稳定了单个正交畴。
CDW 峰强度通过积分卫星反射周围小区域内的散射强度获得,随应变增大,饱和强度总体减小,伴随相变逐渐展宽。
TCDW 通过线性拟合相变附近的积分 CDW 强度-温度曲线确定,拟合斜率随应变增大发生微小变化,反映相变展宽。
Warren-Cowley 短程序参数的应变演化
U-D 对的 Warren-Cowley 短程序参数 alpha 在三种选定应变 (0%, 0.5%, 1%) 下随温度变化,展示从无序到有序的逐渐过渡。
J1a 和 J1b 相互作用的 alpha 参数显示不同的应变依赖性:J1a 沿压缩方向增强,J1b 沿膨胀方向减弱,与各向异性晶格响应一致。
alpha 参数的应变演化幅度与 TCDW 的 ~13 K 提升对应,证实面内 Ising 耦合增强是序参量稳定化的微观起源。
图 S11:序参量演化。(A) 沿应变方向和垂直方向的 CDW 峰温度依赖性;(B) 不同 [H00] 压缩应变下 CDW 峰强度的温度依赖性,虚线为线性拟合;(C) TCDW 和拟合斜率对应变的依赖性。J1a 和 J1b 的 Warren-Cowley 短程序参数 (D) 在三种应变下的温度依赖性。
S12 · 比热数据分析
CDW 相变的热力学特征
CDW 强度-温度曲线的锐利转变对应一级相变的热力学特征:序参量在 TCDW 处不连续跳跃,与比热容的 lambda 型异常一致。
积分 CDW 峰强度可视为 CDW 序参量的实验探针,其温度导数的最大值对应 TCDW,峰宽对应相变宽度。
随应变增大,CDW 强度-温度曲线逐渐展宽,拟合斜率绝对值减小,表明相变从一级向更连续的特征过渡,与三态 Potts 模型预期一致。
应变对比热异常的调制
一级相变在比热容中表现为尖锐的 lambda 峰,峰位对应 TCDW,峰高对应相变潜热。
应变通过增强面内 Ising 耦合强度,一方面提升 TCDW (峰位移动),另一方面通过三聚体无序抹平锐利转变 (峰展宽、峰高降低)。
MC 模拟中 sigmoid 函数拟合的转变宽度 delta 随应变增大,与实验观测的峰展宽趋势一致,支持序-无序相变驱动 TCDW 增强的物理图像。
图 S12:CDW 强度-温度曲线。(A) 不同 [H00] 压缩应变下归一化 CDW 峰强度的温度依赖性;(B) 应变依赖性 TCDW 和拟合斜率,展示相变展宽趋势;(C) 相变附近积分 CDW 强度,虚线为线性拟合。
S13 · 电子结构表征
芯能级光谱与表面终止确认
107 eV 光子测量的芯能级光谱显示清晰的 Sc 3p, V 3s, V 3p 和 Sn 4d 芯能级峰,确认样品的元素组成。
Sn 4d 峰 (~-23.5 eV 和 ~-24.5 eV) 对应自旋轨道分裂的 Sn 4d5/2 和 4d3/2 体芯能级,清晰区分 V3Sn 终止 (红) 和 ScSn2 终止 (蓝) 表面。
Kagome V3Sn 终止仅显示体相 Sn 芯能级峰,ScSn2 终止显示额外的表面相关分量,产生四个明显峰,与先前 ScV6Sn6 的 ARPES 研究一致。
两种表面终止的电子结构差异
V3Sn 终止和 ScSn2 终止的 Fermi 面图谱均显示清晰的 Kagome 晶格六方对称性特征。
仅 V3Sn 终止在低温下显示 M 点附近与 CDW 相变相关的显著电子结构重整化,证实 Kagome 表面在低能电子不稳定性中的关键作用。
两种终止的色散均清晰分辨 Kagome 衍生的 Dirac 锥和 van Hove 奇点,V3Sn 终止显示更强的 CDW 诱导能带重构特征。
图 S13:芯能级光谱。(A) 107 eV 光子全谱,Sc 3p, V 3s, V 3p, Sn 4d 芯能级峰,红/蓝曲线对应 V3Sn/ScSn2 终止;(B) Sn 4d 峰附近放大谱,ScSn2 终止显示四个峰,V3Sn 终止仅显示两个体相峰。
S14 · ARPES 数据
Fermi 面与色散
V3Sn 终止的 Fermi 面在 130 K (TCDW 以上) 和 35 K (TCDW 以下) 均显示六方对称性,低温下出现与 CDW 相变相关的谱权重再分布。
沿 Gamma-M-K-Gamma 方向的色散清晰分辨 Dirac 锥和 van Hove 奇点,V3Sn 终止低温下在 M 点附近显示更强的能带重构。
ScSn2 终止的 Fermi 面和色散特征与 V3Sn 终止定性相似,但 CDW 诱导的低温电子结构重整化明显较弱。
CDW 能隙的打开
沿 Gamma-M 方向,低温 (35 K) 谱在 k|| = 0.55 A^-1 处显示清晰的能隙打开,约 Delta ~ 75 meV。
能量分布曲线 (EDC) 显示两个分裂峰,位于费米能级以下 162 meV 和 235 meV,由 CDW 诱导的能带折叠导致能带反交叉。
能隙形成于 sqrt(3)xsqrt(3) CDW 序相关的布里渊区边界附近,动量依赖的重构与 Kagome 衍生电子态与 CDW 序参量之间的强耦合一致。
图 S14:ARPES 结果。(A-C) V3Sn 终止的 Fermi 面和 Gamma-M-K-Gamma 色散 (130 K/35 K);(D-F) ScSn2 终止的 Fermi 面和色散;(G-H) Gamma-M 方向 CDW 能隙打开,EDC 显示 162 meV 和 235 meV 两个分裂峰。
S15 · 相图细节
TCDW 的应变依赖性:多探针一致性
从积分峰面积 (Area)、最终 MC 能量 (Efinal) 和 Warren-Cowley 短程序参数 (alpha_J1a) 三种独立探针提取的 TCDW 在应变依赖性上高度一致。
[H00] 和 [HH0] 方向压缩的 MC 模拟 TCDW 均显示单调递增,与实验观测的 ~13 K 提升斜率吻合,验证了应变依赖 Ising 耦合模型的定量能力。
三种探针的 TCDW 绝对值在 ~75-95 K 范围内,与实验 TCDW 约 92-105 K 略有差异,源于 MC 模拟中简化模型与真实体系之间的系统偏移。
应变-温度相图与物理机制
应变-温度相图显示 TCDW 随压缩应变单调递增,无饱和迹象,表明在实验可及的应变范围内序参量持续被稳定化。
相图与静水压和化学掺杂形成鲜明对比:后两者压制 CDW,而单轴应变增强 CDW,凸显了对称性破缺(而非体积效应)的关键作用。
三态 Potts 模型框架下,应变通过增强面内 Ising 耦合参数 Ji(d) 降低了三个等价 CDW 畴之间的简并度,等效于向序参量施加共轭场。
该相图为阻挫 Kagome 晶格中序-无序相变的对称性选择性调控提供了定量描述,建立了应变-耦合-序参量-相变温度的完整因果链。
图 S15:拟合 TCDW vs 应变。(A) 积分峰面积提取的 TCDW;(B) 最终 MC 能量提取的 TCDW;(C) J1a 相互作用 Warren-Cowley 短程序参数提取的 TCDW。[H00] 和 [HH0] 方向模拟结果均显示单调递增。
表格
A. Korshunov, C.-Y. Lim, J. Corral-Sertal et al., arXiv:2606.05866 (2026) | ScV6Sn6 · Kagome · CDW · 单轴应变 · Potts 模型