磁振子极化子自旋-晶格耦合

磁振子极化子自旋-晶格耦合

 

CHINESE PHYSICS B 2026 · TOPICAL REVIEW

磁振子极化子自旋-晶格耦合体系中的混合准粒子

Magnon Polarons: Hybrid Quasiparticles in Coupled Spin–Lattice Systems

导读：磁振子极化子（Magnon Polarons）是磁振子与声子通过磁弹耦合杂化形成的混合准粒子，兼具自旋和晶格双重特征。本文系统综述了MPs的物理原理、探测方法（非弹性中子散射、光散射、自旋泵浦、自旋Seebeck/Peltier效应）、材料平台（YIG、稀土铁石榴石、反铁磁体、多铁材料、二维范德华磁体）以及器件应用（SAW器件、毫米级长距离传播）。MPs的形成不仅深刻改变了自旋和热输运的面貌，更为新一代低功耗、相干、多功能信息技术奠定了基础。

一、前言背景

混合准粒子：凝聚态物理的新范式

• 固体中的元激发（电子、光子、磁振子、声子）传统上被独立研究，但它们的相互耦合可产生具有超越个体组分的涌现性质的混合准粒子

• 磁振子（magnon）是磁有序体系中自旋波的量子化激发，具有无电荷、携带角动量、低耗散等优势，被视为超越CMOS计算的有前景信息载体

• 声子（phonon）是晶格振动的量子，主导固体的热学和弹性性质

• 磁振子极化子（Magnon Polarons, MPs）：当磁振子与声子在能量-动量空间中相互靠近时，磁弹耦合（MEC）使两者杂化形成兼具自旋和晶格特征的混合本征态

• 标志性特征：色散关系中的反交叉（anticrossing）行为，伴随磁性与弹性谱权重的重新分配

磁振子极化子的研究意义

• 基础物理：MEC源于自旋-轨道、偶极-偶极和交换相互作用，MPs为研究准粒子杂化、自旋-晶格能量交换和耦合玻色子体系耗散提供理想平台

• 应用价值：为操控自旋、热量和晶格自由度提供新途径，在磁子学、自旋电子学和热电器件中具有重要潜力

• 本综述聚焦：系统总结MPs的实验进展，涵盖物理原理、探测方法、材料平台和新兴功能

• 典型实验手段：非弹性中子散射（INS）、光散射（布里渊/拉曼）、自旋泵浦、自旋Seebeck效应（SSE）、自旋Peltier效应（SPE）

图1 | 磁振子极化子的实验技术、材料和器件概览。涵盖非弹性中子散射、布里渊/拉曼光散射、自旋泵浦、自旋Seebeck效应和声表面波器件等技术手段，以及YIG、稀土铁石榴石、反铁磁体、多铁材料和二维范德华磁体等材料平台。

二、磁弹耦合机制

磁振子与声子的基本描述

• 磁振子哈密顿量：H_mag = Σ_k A_k a_k†a_k + (B_k/2)(a_k†a_{-k}† + a_{-k}a_k)，其中a_k†/a_k为磁振子产生/湮灭算符

• 磁振子色散：ω_k = √(D_ex k² + γμ₀H) · √(D_ex k² + γμ₀H + γμ₀M_s sin²θ_k)，D_ex为交换刚度，γ为旋磁比，M_s为饱和磁化强度

• 声子哈密顿量：H_ph = Σ_{k,λ} ħω_{λk}(c_{λk}†c_{λk} + 1/2)，声学声子色散为ω_{λk} = c_λ|k|，c_λ为对应声子支的声速

• 磁振子色散可通过外磁场H调控（Zeeman效应），而声子色散不受磁场影响——这一差异性是MPs实验探测的关键

磁弹耦合哈密顿量与杂化

• MEC哈密顿量（线性近似）：H_mec ∝ nB_⊥ · Σ k ω_{kλ}^{-1/2} e^{-iφ} a_k(c_{λ,-k}+c_{λk}†) + H.c.，其中B_⊥为磁弹耦合常数，n为自旋数密度

• 总哈密顿量：H = H_mag + H_ph + H_mec，对角化后得到杂化本征模式——磁振子极化子

• 杂化条件：磁振子支与声子支在频率和波矢上同时接近时，MEC使两者杂化，产生反交叉间隙

• 耦合强度：反交叉间隙的大小直接反映了MEC的强度，不同材料体系差异显著

• Kittel先驱性工作：早在1958年就指出MEC可在磁性与弹性激发色散接近时混合两者

磁振子色散关系（含交换、Zeeman和偶极相互作用的完整形式）

线性近似下的磁弹耦合哈密顿量，包含横/纵声子支与磁振子的耦合项

耦合系统的总哈密顿量：对角化后得到磁振子极化子的杂化本征模式

【DFT Tip 1】磁弹耦合常数B_⊥的DFT计算：应力-磁各向异性法

VASP中计算磁弹耦合常数B_⊥的标准方法：① 对晶体施加单轴应变ε=±1%、±2%，② 在每个应变下做SOC自洽计算，提取磁各向异性能MAE(ε)，③ B_⊥=∂(MAE)/∂ε。YIG的B_⊥~10⁶ J/m³量级，2D磁体（如FePS₃）可达~10⁷ J/m³。

关键：应变必须足够小（<2%）以保持线性响应，但也要足够大以克服MAE数值噪声。建议：先测试±0.5%、±1%、±2%三组应变，验证MAE-ε关系的线性度。非线性时需用更小的应变步长。

常见错误：在非SOC下做应变计算——MEC源于SOC，不加SOC得到的B_⊥物理上无意义。另外，ISIF=3（变胞弛豫）时应力张量可直接用于MEC计算，但建议用ISIF=2固定体积做应变扫描以精确控制应变张量。

【DFT Tip 2】磁振子色散的自旋波计算：Phonopy+SpinW vs VASP直接法

MPs研究需要精确的磁振子色散ω_k——两种方案：① VASP中做非共线磁计算+SpinW后处理（四态法提取交换参数J_ij→SpinW计算自旋波谱），② Phonopy的磁振子模式（有限位移法）。

方案①更灵活：可提取完整的J_ij矩阵，计算任意q点的磁振子能量，且能区分各向同性交换和各向异性项。方案②更直接但仅适用于共线磁结构。

参数建议：超胞至少2×2×2（>100原子），K点≥3×3×3，INCAR中设置LSPIRAL=.TRUE.或手动设置MAGMOM。YIG的交换刚度D_ex~5×10⁻⁴⁰ J·m²，计算误差<10%即可接受。

注意：磁振子-声子交叉点附近的k空间采样必须加密——交叉点附近±0.05 Å⁻¹内至少3个q点，否则反交叉间隙的计算误差可>50%。

【DFT Tip 3】声子谱计算精度：磁振子极化子交叉点的关键

MPs的"触碰"条件μ₀H=c²/(4γD_ex)中声速c是关键参数——声速c来自声子色散在Γ点附近的斜率。VASP中声子计算（IBRION=8+DFPT或Phonopy有限位移法）的精度直接影响H_TA/H_LA的预测。

建议：① 声子谱用Phonopy的有限位移法（比DFPT更稳定），② 超胞≥3×3×3（YIG 160原子），③ 力收敛标准EDIFFG≤1×10⁻⁴ eV/Å（比常规-0.01严格100倍），④ 声速拟合沿Γ→X和Γ→L两个方向取平均。

常见陷阱：PBE泛函系统性低估声速~5-10%（晶格常数偏大→声速偏低），导致H_TA/H_LA预测偏低。如果与实验SSE峰值位置偏差>0.5 T，应怀疑声子精度。修正方案：用实验晶格常数做声子计算（而非DFT弛豫晶格常数），或使用PBEsol泛函。

【知识扩展 1】磁弹耦合：从现象学到第一性原理

【理论解释】磁弹耦合（MEC）描述晶格应变与磁序的相互影响——物理根源是SOC：晶格畸变改变晶体场劈裂→通过SOC调制轨道磁矩→影响自旋取向。MEC有两种描述：① 现象学（磁致伸缩常数λ_s），② 微观（磁弹耦合常数B_⊥）。B_⊥与λ_s的关系：B_⊥∝λ_s·c_ij（c_ij为弹性常数），YIG中λ_111~-2.4×10⁻⁶→B_⊥~10⁶ J/m³。

【方法比较】MEC的第一性原理计算：① 应力-磁各向异性法（直接，适合所有体系），② 磁致伸缩法（需大超胞，适合FM），③ 线性响应DFPT（VASP 6.4+支持，最高效但限于共线磁序）。推荐：共线FM用DFPT，非共线/AFM用应力-磁各向异性法。

【经典参考】Kittel, Phys. Rev. 110, 836 (1958)——MEC混合磁振子与声子的原始理论；Rückriegel et al., Phys. Rev. B 89, 184414 (2014)——磁振子极化子的现代量子理论框架；Streib et al., Phys. Rev. B 99, 184442 (2019)——第一性原理计算MEC常数的方法论。

【迁移能力】MEC的分析框架不仅适用于YIG/ReIGs，也适用于所有磁有序体系——铁磁、反铁磁、亚铁磁、螺旋磁体、斯格明子晶格。关键区别：AFM中MEC常数通常比FM大一个量级（因为AFM晶格畸变直接调制交换相互作用），这解释了为什么AFM中MPs效应更强。

三、实验技术

3.1 非弹性中子散射（INS）

图2 | 磁振子极化与INS探测。(a) FM/AFM/FiM中磁振子手性示意图。(b) YIG晶体结构。(c) 极化中子散射装置示意图。(d) YIG极化磁振子谱。(e-f) YIG中磁振子-声子色散及磁场效应——0 T下杂化散射强度增强，9.1 T下减弱。

INS：磁振子极化子的直接观测

• INS利用中子（无电荷、自旋1/2）与磁性材料的非弹性散射，可同时探测能量和动量转移，直接映射全布里渊区磁振子色散E(q)

• 极化中子散射：区分YIG中α（正极化）和β（负极化）磁振子模式，证实YIG为亚铁磁体

• 反交叉直接证据：磁振子与声子支在交点处杂化而非交叉，伴随谱权重转移和线宽展宽

• Man et al. (2017)发现：YIG中0 T下杂化MP散射强度大于9.1 T，与能隙打开不同，耦合增强杂化散射强度

• 优势：直接探测磁矩，覆盖宽能量-动量范围，可测体单晶；局限：需强中子源和大型单晶样品

3.2 光散射技术

图3 | 非弹性光散射。(a-b) Stokes和anti-Stokes散射示意图，满足能量和动量守恒。(c) 拉曼散射和布里渊散射的光谱范围对比——拉曼探测光学声子（>150 GHz），布里渊探测声学声子和磁振子（500 MHz-1.5 THz）。

布里渊与拉曼光谱：互补的光学探针

• 拉曼光谱：配备滤波器和精细光栅，可探测频率>150 GHz（~5 cm⁻¹）的准粒子，主要用于光学声子和高频磁振子

• 布里渊光谱：多通串联Fabry-Pérot干涉仪，可探测500 MHz-1.5 THz的准粒子，用于声学声子和大多数磁振子

• 磁场依赖光谱：通过调谐外磁场移动磁振子色散，直接观测与声子支的反交叉——MPs形成的明确证据

• 2D材料应用：温度依赖拉曼光谱揭示MnPSe₃中磁振子-声子Fano共振，FePS₃中Zeeman分裂磁振子与声子的强耦合（M_c=87.8 GHz）

3.3 自旋泵浦与自旋Seebeck效应

图4 | 自旋泵浦探测MPs。(a) YIG/Pt自旋泵浦装置。(b) 微波功率依赖P_abs和V_ISHE。(c) 微波参量泵浦示意图。(d) 不同磁场下磁振子与TA/LA声子色散。(e) 归一化自旋泵浦效率κ̄-H曲线——MPs形成导致效率峰。

自旋泵浦：电学探测MPs的接口

• 机制：YIG中磁化进动泵浦纯自旋流进入Pt层，通过逆自旋霍尔效应（ISHE）转换为可测电压V_ISHE

• 参量激发：微波光子分裂为频率f_p/2、波矢相反的磁振子对，通过磁场调谐磁振子色散，使激发磁振子波矢接近磁振子-声子交叉点

• 关键发现：归一化自旋泵浦效率κ̄在特定磁场（~78 mT）出现峰值，对应磁振子与TA声子耦合形成MPs的条件

• 优势：直接敏感于自旋动力学，高分辨率探测界面杂化态

图5 | 自旋Seebeck效应示意图。(a) 纵向SSE构型。(b) 磁振子流J_m。(c) 逆自旋霍尔效应。(d) 非局域SSE构型。

图6 | YIG中SSE磁振子极化子异常。(a) LSSE装置与磁振子/声子传播。(b) 50 K LSSE电压H依赖。(c) 20 K nlSSE电压H依赖。(d) 不同磁场下的磁振子-声子色散关系。(e) 散射率比η决定峰/谷形态。(f-g) 温度驱动的峰-谷转变。

SSE：磁振子极化子最灵敏的电学探针

• LSSE和nlSSE中均观测到MPs异常：在临界场H_TA和H_LA处，SSE电压出现尖锐峰/谷结构

• "触碰"条件：当磁振子支与声子支相切（而非交叉）时，磁振子-声子耦合在动量空间中覆盖最广范围，MP效应最强

• 触碰场公式：μ₀H_TA(LA) = c²_T(L)/(4γD_ex)，只依赖于声速和交换刚度

• Boltzmann输运理论（Flebus et al.）：峰/谷形态取决于磁性与非磁性杂质散射率之比η=τ⁻¹_mag/τ⁻¹_ph

• η>1→峰（MP寿命>纯磁振子），η<1→谷（MP寿命<纯磁振子），η=1→异常消失

• 温度驱动峰-谷转变（Shi et al. 2021）：高温下四磁振子散射率急剧上升，η从<1变为>1

图7 | 自旋Peltier效应（SPE）。(a) BiGa:LuIG/Pt SPE装置。(b) 方波电流输入与热电偶电压输出。(c) H依赖ΔV_TC——H_TA=0.16 T和H_LA=1.0 T处出现MP峰。

SPE：SSE的Onsager逆过程

• SPE中自旋流驱动热流：Pt层中电荷流→自旋流→界面非平衡磁振子布居→热流→温度变化

• BiGa:LuIG中MP异常发生在更低磁场（H_TA=0.16 T, H_LA=1.0 T），归因于掺杂效应改变磁振子色散

• SPE与SSE的一致性验证了MPs作为自旋-热-晶格耦合媒介的普适性

【DFT Tip 4】INS谱模拟：SpinW vs McPhase vs UppASD

INS实验数据需要理论模拟来验证MPs的形成——三种主流工具：① SpinW（MATLAB，线性自旋波理论，适合共线磁序），② McPhase（C++，支持平均场+RPA，适合稀土磁体），③ UppASD（原子尺度Landau-Lifshitz-Gilbert动力学，适合有限温度效应）。

对于MPs的INS模拟：SpinW最直接——输入从DFT提取的J_ij矩阵→计算S(q,E)→与实验对比反交叉位置和间隙。UppASD的优势在于可以模拟温度效应：模拟磁振子-声子耦合随温度升高如何减弱。

关键参数：S(q,E)的分辨率取决于超胞大小——要模拟q~0.1 Å⁻¹附近的MPs，超胞需≥30×30×1（~900原子），计算成本~1000核时/温度点。对于YIG（80原子原胞），建议先用SpinW做线性近似，再用UppASD验证温度效应。

【DFT Tip 5】SSE异常的Boltzmann输运模拟：η参数的第一性原理计算

SSE中MP异常的峰/谷由η=τ⁻¹_mag/τ⁻¹_ph决定——理论上可以从第一性原理计算：① τ⁻¹_mag来自磁振子-杂质散射和磁振子-磁振子散射，② τ⁻¹_ph来自声子-杂质散射和声子-声子散射。

τ⁻¹_mag的计算：用DFT提取的J_ij→自旋波色散→态密度→Fermi黄金规则计算散射率。磁振子-杂质散射率∝(ΔJ/J)²·DOS_mag(ω)，其中ΔJ是杂质处交换常数的变化。

τ⁻¹_ph的计算：Phonopy+Phono3py（三阶力常数）→声子寿命。YIG中声子寿命通常~10⁻¹⁰-10⁻⁹ s（高纯晶体），磁振子寿命~10⁻⁸-10⁻⁷ s→η<1→SSE谷。

实践建议：如果无法完整计算τ⁻¹_mag和τ⁻¹_ph，可以用实验SSE峰/谷方向反推η——这是更可靠的验证方式。BFO/LSMO中η≈0.3-0.5（谷），Cr₂O₃中η≈1.5-2.0（峰）。

【DFT Tip 6】布里渊光谱的DFT辅助分析：磁振子频率的温度依赖性

布里渊光谱探测的磁振子频率f_m通常在1-100 GHz——与DFT的零温结果有系统性偏差。温度效应：① 交换刚度D_ex(T)∝M(T)²（低温），② 磁振子-磁振子相互作用导致频率重正化。

DFT修正方法：① 计算T=0 K的D_ex和声速c，② 用实验M(T)曲线标度D_ex(T)≈D_ex(0)·[M(T)/M(0)]²，③ 代入触碰场公式H_TA(T)=c²/(4γD_ex(T))。

YIG中D_ex(T)从0 K到300 K降低~15%，导致H_TA从2.51 T降至2.14 T——与实验趋势一致。GdIG中磁补偿点附近M(T)→0，D_ex→0，H_TA→∞——解释了为什么GdIG在T_comp附近无MP信号。

常见错误：直接拿DFT的D_ex(0)代入触碰场公式与室温实验对比——误差可达20-30%。必须做温度修正。

四、自旋与热输运中的磁振子极化子

SSE中MP异常的微观机制

• 热梯度激发宽谱磁振子，当磁振子能量与声子色散相切时，MPs形成→SSE电压出现共振特征

• 触碰vs反交叉：触碰条件下（H=H_TA, H_LA）磁振子与声子在宽动量范围内耦合，效应最强

• 峰/谷转变：η=τ⁻¹_mag/τ⁻¹_ph控制MPs相对纯磁振子的弛豫时间——这是区分薄膜磁性和声学质量的重要判据

• 四磁振子散射：高温下显著增强，是SSE温度依赖中MP信号变化的关键因素

• 非局域SSE：通过改变注入-探测距离，可精确测定自旋载流子的输运性质

磁振子-声子"触碰"临界磁场：由声速和交换刚度唯一确定

磁性与非磁性杂质散射率之比：决定SSE异常的峰/谷形态

反交叉能级：描述FePS₃等2D材料中磁振子-声子杂化的能级分裂，M_c为耦合强度（87.8 GHz）

【DFT Tip 7】非磁性离子替代的DFT建模：虚晶近似 vs 超胞法

材料工程中Bi替代Y（Bi:YIG）或Ga替代Fe的DFT建模有两种策略：① 虚晶近似（VCA）——VASP中无原生支持，需用特殊赝势或混合原子法，精度有限，② 超胞法——建2×2×2超胞（160原子），手动替换特定位置的原子，是最可靠的方法。

超胞法的关键：① 必须测试多种替代构型（如Bi替代四面体d位vs八面体a位），② 用SQS（特殊准随机结构）模拟无序替代，ATAT或mcsqs工具可生成，③ 弛豫时允许晶格常数和原子位置同时变化（ISIF=3）。

Bi替代对声子色散的影响：Bi原子质量大→TA/LA声速降低10-20%→H_TA降低。但DFT通常低估质量效应对声速的影响（因为PBE低估晶格常数），建议用实验晶格常数约束弛豫。Ga替代对磁振子色散的影响更复杂——需同时计算交换参数J_ij的变化和磁矩的重分布。

【DFT Tip 8】拓扑磁振子极化子的计算：Berry曲率+Chern数

FMO（Fe₂Mo₃O₈）中拓扑MPs的计算需要：① 构建磁振子-声子耦合的BdG哈密顿量（~4N×4N矩阵，N为原子数），② 对角化得到杂化能带，③ 计算Berry曲率Ω_n(k)和Chern数C_n。

工具链：VASP（磁构型+声子）→输出J_ij和力常数矩阵→Python脚本构建BdG哈密顿量→BerryCurvature（Wannier90风格的k空间积分）。

关键参数：① k空间网格≥200×200（六方晶格）以保证Chern数收敛到整数，② DM相互作用的符号和大小对拓扑性质有决定性影响——FMO中DM诱导的MEC是拓扑MPs的必要条件（无DM时C=0）。

常见错误：只计算磁振子带的Berry曲率而忽略磁振子-声子杂化带——拓扑性质来自杂化能带的能带反演，纯磁振子带可能拓扑平庸。另外，时间反演对称性破缺是Berry曲率非零的前提——确认体系有净磁矩或外加磁场。

【科研经验 1】SSE中的MP异常 vs 磁化翻转：如何区分？

问题：SSE信号在特定磁场下的峰/谷结构不一定来自MPs——磁化翻转、磁畴壁运动、自旋翻转转变均可产生类似特征。许多初学者将所有SSE异常归因于MPs，导致误判。

原因：MPs的SSE异常和磁化翻转的SSE异常都表现为H依赖的尖锐结构，且两者都出现在磁转变附近。TbIG中97.37 K的额外符号反转起初被怀疑是MPs，但最终确认是DUMS磁结构引起的磁化翻转。

解决方案：① 同时测量M-H曲线——如果SSE异常与M-H的不可逆行为（磁滞）重合，大概率是磁化翻转而非MPs，② 做角度依赖测量——MPs的触碰场H_TA/H_LA应随磁场方向连续变化，而磁化翻转是离散的，③ 做温度依赖——MPs异常随温度平滑演化，磁化翻转在T_comp附近突变，④ 理论验证——用DFT计算D_ex和c后预测H_TA/H_LA，与实验对比。

建议：任何声称"SSE中观测到MPs"的工作，必须提供H_TA/H_LA与理论预测的定量对比。仅凭SSE峰形不足以证明MPs。

五、材料体系

5.1 YIG及其非磁性离子替代

图8 | 非磁性离子替代调控MPs。(a) YIG和Bi:YIG声子色散。(b) 触碰条件下的磁振子-声子色散。(c) Bi_xY_{3-x}IG/Pt的LSSE系数S。(d-e) 非磁性替代示意与磁性补偿效应。(f-i) (Lu₂Bi)Fe₄GaO₁₂中MPs——H_TA=0.42 T，峰强比YIG增强~700%。

Bi替代：调控声子色散

• Bi替代Y位（Bi:YIG）：Bi原子质量（209）远大于Y（89），降低TA和LA声速

• 效果：触碰场H_TA和H_LA随Bi含量x增加而降低（x=0→0.5→0.9：H_TA从2.51 T降至2.06 T）

• MP峰强度增强：归因于Bi替代增强的磁弹耦合

• Ga共掺杂(Lu₂Bi)Fe₄GaO₁₂：Ga优先占据d位，引入磁性补偿效应，磁振子色散从抛物线型（类FM）向线性（类AFM）转变

磁性补偿：最大化MP效应

• 在补偿态附近，磁振子群速度接近声速→磁振子与声子色散重叠最大化→耦合显著增强

• (Lu₂Bi)Fe₄GaO₁₂中H_TA=0.42 T、H_LA=1.86 T（300 K），远低于纯YIG

• 关键成果：MP峰强度比YIG增强约700%，归因于磁振子寿命降低使η≥1，杂化后磁振子寿命大幅增强

5.2 补偿型稀土铁石榴石（ReIGs）

图9 | YIG和GdIG的磁性与SSE。(a) YIG M-T曲线。(b) YIG S-T曲线。(c) YIG自旋波谱。(d) GdIG自旋波谱。(e) α（右手）和β（左手）磁振子。(f) GdIG三套磁性子晶格。(g) GdIG磁化强度和SSE信号温度依赖。

GdIG：选择性磁振子激发的平台

• GdIG含两个手性磁振子支（α和β），可通过温度选择性激发——Gd子晶格对温度敏感

• 补偿温度T_comp：Gd和Fe子晶格磁矩反平行，T>T_comp时SSE信号由β支主导，T

• 进一步降温：Gd主导的α支贡献超过β支→第二次SSE符号反转

图10 | GdIG中温度驱动的MP异常演化。(a) 不同温度下SSC-H曲线。(b) 120 K和150 K自旋波谱。(c) α模磁振子与TA声子群速度关系。(d) 双峰形成的色散关系。(e) GdIG/Pt实验SSC-H曲线——单峰到双峰转变。

GdIG中单峰→双峰MP异常

• 理论预测（Shen 2019）：随温度升高，α模磁振子群速度v_m-k曲线逐渐抬升

• T<~60 K：v_m远小于v_t，无MP异常

• 60-108 K：v_m≈v_t→单峰MP异常，108 K附近强度最大

• T>108 K：v_m与v_t有两个交点→双峰结构，间距随温度增大

• 不同于YIG中TA/LA双峰，GdIG双峰均来自α磁振子与TA声子的耦合——温度驱动的群速度变化是根本原因

5.2b TbIG：非共线磁结构与巨MP异常

图11 | TbIG中的LSSE与磁结构。(a) TbIG/Pt样品中LSSE测量示意图，H沿[100]轴扫描。(b) 不同温度下SSC-H曲线——97.37 K时出现±2 T处额外符号反转和±6 T处磁滞异常。(c) 体相TbIG在4.2 K的M-H曲线。(d) 模拟的[100]方向总磁化强度M[100]。(e) DUMS（双伞形磁结构）在不同磁场H₁、H₂、H₃下的示意图。

TbIG：非共线磁结构中的巨MP异常

• TbIG具有强各向异性和大磁致伸缩常数（Tb离子非零轨道角动量），暗示潜在的大MP诱导SSE异常

• Li et al. LSSE测量（H // [100]）：207 K和157 K时SSC表现为常规台阶状符号翻转

• 97.37 K的异常：±2 T处额外符号反转 + ±6 T处磁滞状异常，SSC放大一个量级以上

• 物理根源：TbIG在低温下呈现"双伞形磁结构"（DUMS）——Tb子晶格的非共线排列导致复杂的磁化翻转过程

• 模拟表明：H扫描过程中，Fe和Tb子晶格磁矩的相对取向发生多次重构，每次重构伴随MP耦合条件的变化

图12 | TbIG中SSC补偿与MP增强。(a) SSC-T和M-T曲线。(b) T_comp以下SSC-H曲线。(c) MP导致的SSC增强计算值随H变化。(d) 体相TbIG低频磁振子谱（<0.5 THz），∆T=-9.5 K，三个典型磁场H₁、H₂、H₃下的磁振子色散。(e) T_comp以上SSC-H曲线和Fe子晶格磁矩投影。(f) H-∆T相图——SSC和θ_Fe（M_Fe与H夹角）的精细温度映射。

TbIG中MP增强的物理机制

• 补偿温度T_comp以下：SSC在特定磁场出现MP增强峰，强度与GdIG相当

• 磁振子谱计算显示：低频磁振子支与TA/LA声子在多个磁场下满足"触碰"条件

• H-∆T相图揭示：Fe子晶格磁矩取向θ_Fe随H和∆T的系统演化——磁矩翻转过程直接调制MP耦合强度

• 非共线磁结构中的MPs：不同于YIG和GdIG的共线磁结构，TbIG的伞形磁结构为MPs提供了更丰富的调控维度

• 意义：非共线ReIGs系列（TbIG、DyIG、HoIG等）为探索磁织构与MPs相互作用提供了独特平台

5.3 反铁磁绝缘体

图13 | Cr₂O₃中的反铁磁MPs。(a) Cr₂O₃晶体结构。(b) 单轴反铁磁体中MPs实空间示意图。(c-d) 磁振子模式。(e-f) SF转变前后的磁振子-声子色散。(g-h) Cr₂O₃/Ta和Cr₂O₃/Pt的SSE信号。

Cr₂O₃：反铁磁MPs的范式体系

• Cr₂O₃为单轴AFM绝缘体，易轴沿c轴，AFM共振频率~0.165 THz，自旋翻转场H_SF=6.0 T

• 两套简并磁振子模式（α和β），H沿c轴可解除简并——β支下移与LA/TA声子相交形成MPs

• SSE异常均为峰→AFM磁振子散射强于声子（η>1）

• H_SF处SSE突变+符号反转：磁振子模式从β切换到m-3，自旋极化反转

• 理论拓展（Liu & Shen 2022）：不同各向异性场和交换场下，MP异常可出现1-3个峰

5.4 多铁材料中的拓扑磁振子极化子

图14 | Fe₂Mo₃O₈中的拓扑MPs。(a) FMO晶体和磁结构。(b) 磁振子-声子反交叉示意图。(c-d) INS测量的MP激发谱。(e) DM相互作用诱导的磁振子-声子耦合。(f) 无/有DM耦合的MP能带计算。

FMO：DM相互作用驱动的新型MPs

• FMO（Fe₂Mo₃O₈）：极性六方结构（P6₃mc），Fe²⁺蜂窝子晶格，反演和时间反演对称性均破缺

• 强单轴（Ising型）各向异性，面内DM相互作用，巨线性磁电效应

• INS全景映射（Bao et al. 2023）：高能磁振子（~11 meV, 14 meV）与声子杂化形成MPs

• 关键：无DM耦合时计算无法重现反交叉；加入DM诱导的磁振子-声子耦合后完美复现实验

• 拓扑MPs：能带反演+磁振子-声子成分混合→拓扑非平庸态，为低损耗磁子学逻辑器件开辟新途径

5.5 二维范德华磁体

图15 | 二维AFM中的MPs。(a) MnPSe₃温度依赖拉曼光谱——Fano共振。(b) FePS₃磁场依赖拉曼光谱——Zeeman分裂和反交叉。(c) 圆偏振度。(d-e) FePSe₃中手性选择性杂化。

2D磁体：手性选择性与维度调控

• MPX₃（M=Fe, Mn, Ni; X=S, Se）系列：2D极限下表面自旋完全暴露+面外声子振幅增强→强磁振子-声子耦合

• MnPSe₃：温度降至T_N以下，E_g声子与双磁振子连续谱杂化，出现非对称Fano线型→磁振子退相干途径的确立

• FePS₃：磁振子Zeeman分裂为M↑和M↓，与声子反交叉，耦合强度M_c=87.8 GHz——比FM体系大2个量级

• FePSe₃手性选择性：非对称螺旋轴对称性限制声子P±仅与同角动量方向的磁振子Mα/β耦合——角动量匹配条件

• 应用前景：纠缠磁振子-声子-光子腔体系，角动量编码量子信息

【DFT Tip 9】2D vdW磁体中MPs的DFT计算：维度效应与基底选择

2D磁体（MPX₃系列）的MPs计算面临独特挑战：① 面外声子振幅在单层极限下增强→MEC增大，② 基底效应（如SiO₂/Si或hBN）显著改变声子色散。

VASP设置：① 真空层≥15 Å（避免层间镜像相互作用），② 偶极修正LDIPOL=.TRUE.+IDIPOL=3，③ vdW修正DFT-D3（IVDW=11）——层间距影响磁交换J_ij。

基底效应建模：建议测试两种极限——① 悬空单层（无基底），② 刚性基底（固定底层原子）。FePS₃计算表明：悬空vs基底的声子频率差可达10-20 GHz，足以影响MPs的交叉位置。

关键：2D磁体的J_ij通常比体材料小（维度效应+量子限域），但MEC常数B_⊥反而更大（面外声子振幅增强）。FePS₃的M_c=87.8 GHz（比YIG大2个量级）即源于此——不要用体材料的MEC常数估算2D体系的MPs耦合强度。

【DFT Tip 10】DM相互作用的VASP计算：手性MPs的关键输入

FMO中DM相互作用诱导的磁振子-声子耦合是拓扑MPs的核心——VASP中计算DM相互作用：① 打开LSORBIT=.TRUE.（SOC是DM的物理来源），② 对多个自旋构型做非共线计算（四态法），③ 从能量差中提取D_ij矢量。

参数：自旋构型需包含FM和至少三种不同螺旋构型——D_ij的方向和大小都敏感于自旋取向。建议用2×2×1超胞（≥48原子），K点≥4×4×2。

输出分析：D_ij的符号决定磁振子-声子杂化的手性——D_ij>0时左旋MPs能量更低，D_ij<0时右旋更低。FMO中D_ij~0.5 meV（面内），与交换J~3 meV竞争→能带反演→拓扑非平庸。

常见错误：① 在非SOC下提取DM相互作用——DM严格由SOC导致，非SOC结果物理上无意义，② 超胞太小时周期性镜像DM相互作用干扰→用足够大的超胞（DM相互作用随距离衰减快，但周期性边界条件可能引入伪影）。

【科研经验 2】INS中MPs的交叉 vs 反交叉：分辨率陷阱

问题：INS实验中色散交叉（crossing）和反交叉（anticrossing）的区分取决于能量分辨率——分辨率不足时，小的反交叉间隙（<1 meV）可能被误判为交叉，导致MPs耦合强度被低估或漏检。

原因：INS的能量分辨率ΔE由中子能量、单色器/分析器类型和飞行时间决定。YIG中磁振子-TA声子反交叉间隙~0.2 meV，如果ΔE>0.2 meV则无法分辨。Man et al. (2017)使用高分辨率三轴谱仪（ΔE~0.1 meV）才成功观测到反交叉。

解决方案：① 引用INS谱仪的具体分辨率参数（ΔE和ΔQ），② 在反交叉附近做更精细的能量扫描（步长≤0.05 meV），③ 做磁场依赖——如果间隙随H变化则来自MPs，不随H变化则可能是无关的能带交叉，④ 理论辅助——先用DFT+SpinW计算预期的反交叉间隙，判断是否在INS分辨率范围内。

建议：反交叉间隙<仪器分辨率时，不能排除MPs的存在——可用光散射（布里渊/拉曼，分辨率~0.1 GHz≈4×10⁻⁴ meV）或SSE（通过输运异常间接探测）交叉验证。

【科研经验 3】2D磁体中MPs的基底效应：从悬空到异质结

问题：2D磁体（如FePS₃、CrI₃）的MPs测量通常在SiO₂/Si或hBN基底上进行——但基底如何影响声子色散和磁弹耦合？许多研究忽略基底效应，导致DFT预测与实验MPs交叉位置偏差>10 GHz。

原因：基底引入两个效应：① 机械耦合——基底刚度改变2D磁体的面外声子频率（悬空→低，基底→高），② 静电掺杂——基底与2D磁体之间的电荷转移改变磁交换J_ij。FePS₃/SiO₂中声子频率比悬空高~15 GHz，MPs交叉点偏移~0.05 T。

解决方案：① DFT计算中明确包含基底（至少2-3层基底原子），② 对比悬空和基底的MPs预测——如果两者差异<10%，基底效应可忽略，③ 实验上：用hBN封装（最接近悬空）、SiO₂/Si（常规）、金属基底（强耦合）三种基底做对比测量。

建议：任何2D磁体MPs的DFT计算，必须报告基底模型（悬空/hBN/SiO₂）。如果基底效应>10%，需用包含基底的超胞做完整计算——超胞≥3×3×1 2D磁体+3层基底，~200原子，~500核时。

六、器件应用与展望

6.1 声表面波器件

图16 | SAW器件中磁振子-声子耦合。(a) Ni/LiNbO₃中SAW激发FMR。(b) 不同频率SAW传输的H依赖。(c) LiNbO₃/Co/Pt中SAW自旋泵浦。(d) SAW功率差和ΔV_DC的H依赖。

SAW：集成化磁声耦合平台

• SAW频率（MHz-GHz）与FMR天然兼容，叉指换能器（IDT）可高效激发和探测

• Weiler et al. (2011)：GHz SAW在Ni/LiNbO₃中实现全弹性FMR激发和探测——无需外加射频磁场

• 声子驱动自旋泵浦：SAW脉冲提供相干弹性波，在Co/Pt中产生自旋流，通过ISHE检测

• 应用：磁声传感器、声波辅助磁记录、可调谐滤波器、雷达和通信核心组件

6.2 长距离磁振子极化子传播

图17 | BFO/LSMO中长距离MP传播。(a-b) BFO/LSMO界面和磁振子-声子耦合示意图。(c-d) 自旋波传输谱和衰减长度。(e-f) 反交叉测量与模拟。(g-n) BTO/LSMO中Rayleigh/Love SAW与自旋波杂化。

毫米级磁振子传播：MPs的里程碑突破

• BFO/LSMO异质结：室温下磁振子衰减长度达1 mm——比裸LSMO金属薄膜大2个量级以上

• 两种模式：Mode Y（常规Damon-Eshbach模，衰减~5 μm，群速~0.8 km/s），Mode X（MP杂化模，群速~2.5 km/s，毫米级衰减长度）

• BTO/LSMO：Rayleigh和Love SAW与自旋波相干杂化，传播距离>80 μm，波长可通过H强度和取向精确调控

• 耦合强度：有效层间MEC约2×10⁵ J/m³，理论模拟定量复现反交叉行为

• 意义：克服金属磁子体系本征阻尼限制，为低损耗、高速混合磁子器件提供可行路径

6.3 未来展望

磁振子极化子研究的五大前沿方向

(1) 交变磁体（Altermagnets）：动量依赖的高各向异性磁振子带+手性分裂→新型MPs平台，有望实现高频、长寿命信息载体

(2) 超快非平衡动力学：飞秒光脉冲和太赫兹脉冲操控MPs→超快信息处理，实时追踪MP形成与退相干

(3) 拓扑磁振子极化子：非共线磁织构（斯格明子、自旋螺旋）提供对称性破缺环境→拓扑保护MP边态，实现无耗散手性边电流

(4) 量子信息处理：长寿命声子作为磁振子态的量子存储器，微波-光学转换，拓扑MP边态作为量子信息通道

(5) 强非线性效应：声子介导的磁振子非线性→宏观量子态生成，量子计量学应用

 

Further 1 | 磁振子-声子耦合常数B_⊥的第一性原理高通量计算

为什么值得算：综述覆盖了YIG、ReIGs、Cr₂O₃、FMO、MPX₃等多个体系，但各体系的MEC常数B_⊥仅来自实验拟合或单独计算，缺乏系统性对比。高通量计算B_⊥可以揭示MEC的化学趋势——哪些元素/结构产生强MEC？

能回答的问题：为什么AFM中MPs效应比FM强？B_⊥与SOC常数的标度关系？2D vs 3D磁体的MEC差异有何普适规律？

适合体系：Materials Project中所有磁有序绝缘体（~2000种），筛选D_ex>0.1 meV、c<5000 m/s的候选体系。

输入：VASP SOC计算（自动应变+MAE），弹性常数，Phonopy声子谱。每个体系~6个应变×2个方向=12个SOC计算。

 

Further 2 | 有限温度磁振子极化子：从零温DFT到Monte Carlo+LLG

为什么值得算：综述中所有DFT计算都是T=0 K，但SSE实验在10-300 K进行。非零温下磁振子-磁振子散射（四磁振子过程）和声子-声子散射显著改变η参数——温度驱动的峰-谷转变（Shi et al. 2021）需要第一性原理定量解释。

能回答的问题：给定材料，η(T)如何演化？η=1的临界温度T_c是多少？T_c与材料参数（D_ex、c、缺陷浓度）的关系？

适合体系：YIG、GdIG、Cr₂O₃（已有SSE实验数据可验证）。

输入：DFT→J_ij+力常数→UppASD（LLG模拟）+Phono3py（声子寿命）。

 

Further 3 | 手性声子与磁振子的角动量转移：非弹性X射线散射+DFT

为什么值得算：FePSe₃的实验揭示了手性选择性的MPs杂化——声子P±仅与同角动量方向的磁振子Mα/β耦合。但这一机制是否普遍？手性声子是否存在于更多材料中？

能回答的问题：哪些材料具有手性声子？手性声子-磁振子的角动量转移效率？能否通过手性选择性实现"磁振子二极管"？

适合体系：所有非中心对称磁体（FMO、CrI₃、CrSBr、Cu₂OSeO₃、TbIG）。

输入：DFT（声子本征矢分析手性）+非弹性X射线散射（IXS，比INS更高能量分辨率）。

 

Further 4 | 反铁磁体MPs的SSE：Cr₂O₃之外还有哪些候选？

为什么值得算：Cr₂O₃是目前唯一系统研究SSE中反铁磁MPs的体系。但AFM家族庞大——NiO、MnF₂、Fe₂O₃、MnPS₃、Cu₃TeO₆等都可能展现MPs。

能回答的问题：AFM中MPs的SSE异常是否普遍？AFM磁振子手性（R/L模式）如何影响SSE符号？Néel温度vs MP耦合强度有相关性吗？

适合体系：NiO（T_N=523 K，高稳定性）、MnF₂（单轴AFM，类似Cr₂O₃）、FePS₃（2D AFM，已有拉曼MPs证据）。

输入：DFT→J_ij→自旋波谱→预测H_TA/H_LA→设计SSE实验。

 

Further 5 | 磁振子极化子介导的界面热导：NEGF+DFT

为什么值得算：MPs改变了磁振子和声子的输运性质——在界面处，MPs的形成是否增强或抑制界面热导？这对自旋热电器件的热管理至关重要。

能回答的问题：YIG/Pt界面热导中MPs的贡献比例？MPs对界面Kapitza热阻的影响？能否通过外磁场调控界面热导？

适合体系：YIG/Pt、YIG/Au、BFO/LSMO、Cr₂O₃/Pt——所有磁性绝缘体/金属界面。

输入：DFT（界面结构+力常数）→NEGF（非平衡格林函数）界面热导计算。

 

Further 6 | 磁振子极化子的非线性效应：Kerr频率梳+四波混频

为什么值得算：综述提及"强非线性效应"作为未来方向之一，但未深入。磁振子Kerr非线性（源自磁各向异性和偶极相互作用）与声子非线性结合，可产生磁振子极化子频率梳。

能回答的问题：MPs的Kerr系数n_2是否比纯磁振子大？能否实现MPs孤子？磁振子极化子频率梳的阈值功率？

适合体系：YIG（超低阻尼，~10⁻⁵）、GdIG（补偿点附近非线性增强）、FMO（拓扑保护边态+非线性）。

输入：DFT（磁各向异性常数）→LLG模拟（非线性自旋波动力学）→耦合声子模式。

 

Further 7 | Altermagnets中的磁振子极化子：CrSb、MnTe、RuO₂的预测

为什么值得算：交变磁体（Altermagnets）是2023年才发现的新磁性相——动量空间中的d波/g波自旋劈裂，无净磁矩但有自旋流。交变磁体的磁振子色散具有独特的各向异性（节点线/节点面），与声子耦合将产生全新的MPs拓扑结构。

能回答的问题：交变磁体中MPs的拓扑性质（Chern数、边态手性）？节点线MPs的输运特征？与常规FM/AFM的MPs有何本质区别？

适合体系：CrSb（T_N=703 K，高稳定性）、MnTe（T_N=307 K）、RuO₂（T_N>300 K）、V₂Se₂O（2D altermagnet）。

输入：DFT（交变磁序+SOC）→J_ij矩阵（含各向异性）→构建BdG哈密顿量→Berry曲率+Chern数。

 

【知识扩展 2】自旋Seebeck效应：从自旋卡路里电子学到磁振子极化子

【理论解释】SSE（自旋Seebeck效应）是由温度梯度驱动的自旋流产生现象——2008年Uchida等人在YIG/Pt中首次发现。物理过程：热梯度→磁振子布居非平衡→界面自旋流→Pt层ISHE→电压。SSE不依赖净电荷流，只依赖自旋自由度，是自旋卡路里电子学（Spin Caloritronics）的核心效应。

【方法比较】SSE的探测方式：① 纵向SSE（LSSE，最简单，热梯度垂直于界面），② 横向SSE（TSSE，需要面内磁化），③ 非局域SSE（nlSSE，分离注入和探测端，可测自旋扩散长度）。nlSSE是最干净的MPs探针——排除Seebeck效应和异常Nernst效应的干扰。

【经典参考】Uchida et al., Nature 455, 778 (2008)——SSE首次发现；Bauer, Saitoh, van Wees, Nat. Mater. 11, 391 (2012)——自旋卡路里电子学综述；Kikkawa et al., PRL 117, 207203 (2016)——SSE中MPs的首次观测；Flebus et al., PRB 95, 144420 (2017)——MPs输运的Boltzmann理论。

【迁移能力】SSE中MPs的分析框架（η参数→峰/谷）适用于所有磁绝缘体/金属双层体系——不仅限于YIG/Pt。关键迁移条件：① 磁绝缘体需有低阻尼（α<10⁻³），② 界面需有强自旋-轨道耦合（Pt, Ta, W比Au, Cu好），③ 温度梯度需足够大（>0.1 K/μm）以产生可测信号。

参考文献

主要参考文献

[1] Li Y, Gu Q, Wang M, et al. Chin. Phys. B 35, 047507 (2026) — 本综述

[2] Kikkawa T, Shen K, Flebus B, et al. Phys. Rev. Lett. 117, 207203 (2016) — YIG中LSSE磁振子极化子首次发现

[3] Cornelissen LJ, et al. Phys. Rev. B 96, 104441 (2017) — 非局域SSE探测MPs

[4] Flebus B, et al. Phys. Rev. B 95, 144420 (2017) — Boltzmann输运理论，η参数

[5] Shi Z, et al. Phys. Rev. Lett. 127, 277203 (2021) — 温度驱动峰-谷转变

[6] Bao S, et al. Nat. Commun. 14, 6093 (2023) — FMO拓扑MPs

[7] Zhang J, et al. Nat. Commun. 12, 7258 (2021) — BFO/LSMO毫米级MP传播

[8] Li J, et al. Phys. Rev. Lett. 125, 217201 (2020) — Cr₂O₃反铁磁MPs

[9] Shen K, Phys. Rev. B 99, 024417 (2019) — GdIG单峰→双峰理论

[10] Liu S, et al. Phys. Rev. Lett. 127, 097401 (2021) — FePS₃强磁振子-声子耦合

[11] Cui J, et al. Nat. Commun. 14, 3396 (2023) — FePSe₃手性选择性杂化

[12] Ramos R, et al. Nat. Commun. 10, 5162 (2019) — 磁性补偿增强MPs 700%

Y. Li, Q. Gu, M. Wang, et al. | Chinese Physics B 35, 047507 (2026) | 磁振子极化子 · 磁弹耦合 · 自旋Seebeck效应 · 混合准粒子