Mn3Sn三自旋模型与磁振子调控机理研究

1. Mn3Sn三自旋模型的理论基础1.1 非共线反铁磁序的物理图像Mn3Sn作为一种典型的Kagome晶格反铁磁体其磁性Mn原子在(001)面内形成六角蜂窝状排列。实验观测表明在室温下该材料呈现120°非共线自旋结构三个相邻Mn原子的磁矩在基面内形成三角形构型。这种独特的自旋排布源于以下几个关键相互作用海森堡交换作用J≈-67.03 meV主导的最近邻铁磁耦合使自旋倾向于平行排列Dzyaloshinskii-Moriya相互作用D≈12.7 meV由晶体结构破缺中心对称性引起的手性耦合单轴磁晶各向异性K≈0.2 meV将自旋约束在基面内特别值得注意的是我们的DFT计算首次揭示了显著的双二次交换作用B≈1.5 meV这种高阶相互作用在传统理论模型中常被忽略但对理解Mn3Sn的磁振子谱具有关键影响。1.2 六自旋到三自旋的粗粒化过程原始六原子晶胞包含两层Kagome晶格可通过以下映射简化为三自旋模型自旋配对将上下两层中垂直对齐的Mn原子1-4, 2-5, 3-6耦合成三个宏自旋A、B、C每个宏自旋包含6μB磁矩参数重整化交换作用J_eff ΣJ_ij对所有相关原子对求和DMI强度D_eff 2D_atomic考虑层间增强效应各向异性K_eff 2K_atomic这种粗粒化处理的物理合理性在于层间铁磁耦合J1,4≈-67.03 meV比层内作用强约3倍高频光学模10 THz与低频声学模1 THz存在明显能隙宏观测量主要响应低频集体激发关键验证通过对比六自旋与三自旋模型的磁振子谱我们发现当能量15 meV时两者差异小于5%证实了模型简化的有效性。2. 有效哈密顿量的构建与解析2.1 完整哈密顿量的组成三自旋系统的哈密顿量密度可表示为H J*(mA·mB mB·mC mC·mA) # 海森堡交换 D*uz·(mA×mB mB×mC mC×mA) # DMI作用 - K*Σ(mi·uai)^2 # 单轴各向异性 B*Σ(mi·mj)^2 # 双二次交换 - ΔJ*(mA·mB) - ΔB*(mA·mB)^2 # 应变诱导对称破缺其中各向异性轴取值为uA [-√3/2, -1/2, 0] uB [-√3/2, 1/2, 0] uC [ 0, 1, 0]2.2 双二次交换的独特作用传统理论常忽略B项但我们的DFT计算显示其在Mn3Sn中不可忽视能隙调控B≈1.5 meV贡献约30%的磁振子能隙模式耦合导致声学支与光学支的非线性混合应变响应ΔB/Δε≈0.8 meV/%ε为应变通过微扰理论可推导出双二次交换对磁振子频率的修正Δω_B ≈ (4γB/μ0Ms)√(1-(m·n)^2)其中γ为旋磁比Ms为饱和磁化强度m、n为相邻自旋方向。3. 自旋动力学方程的推导3.1 集体坐标变换为分离快慢模式引入新的坐标系慢变量声学模ϕ0 -(ϕAϕBϕC)/3 集体进动θ0 (θAθBθC)/3 - π/2 集体倾斜快变量光学模δi ϕi - ϕ0 ± 2π/3 面内涨落ui θi - θ0 - π/2 面外涨落3.2 运动方程的解耦通过Lagrangian形式推导得到耦合方程组(1) 3μ0Ms/γ*(θ̇0 - αϕ̇0) √3Ksin(2ϕ0)(δB-δA) - 3Kcos(2ϕ0)(δAδB) (2) (9J-9B3√3D)θ0 3αμ0Ms/γ*θ̇0 -3μ0Ms/γ*ϕ̇0 (3) √3(ΔB-ΔJKcos2ϕ0) - 3(2B√3DJ)(2δAδB) 2μ0Ms/γ*(2u̇Au̇B)采用绝热近似光学模瞬时平衡最终简化为八极矩动力学方程ω_K*sin(2ϕ0) - ṁz - αϕ̇0 0 ω⊥*mz (αṁz - ϕ̇0) 0其中有效频率参数ω_K 2γK(-ΔBΔJ)/[3μ0Ms(2B√3DJ)] ω⊥ γ(3J-3B√3D)/μ0Ms4. 磁振子谱的计算与验证4.1 线性自旋波理论实现采用SpinW软件包计算磁振子谱的关键步骤超胞构建2×1×1超胞用于提取精确的Jij参数哈密顿量对角化model spinw; model.genlattice(lat_const,[5.8 5.8 4.5],angled,[90 90 120]); model.addmatrix(label,J1,value,-67.03); model.gencoupling(maxDistance,5); sw_plotspec(model,mode,disp,Evect,0:0.1:50);对称性路径Γ→K→M→Γ六角布里渊区4.2 计算结果与实验对比特征点理论值(meV)实验值(meV)偏差来源Γ点声学模0.00.0完美符合K点光学模28.525.3长程交换被忽略M点声学模18.716.2DMI矢量不精确主要差异来源于未考虑第三近邻交换预计影响约2-3 meVDMI的空间分布假设过于简化温度效应计算为0K实验在300K5. 应变调控效应分析5.1 单轴应变的影响机制施加εxx应变时系统发生以下变化交换作用调制ΔJ/J ≈ 5.3%/GPa通过DFT计算ΔB/B ≈ 3.8%/GPaDMI重定向D矢量随晶格畸变旋转角度变化Δθ≈0.12rad/GPa各向异性增强ΔK/K ≈ 7.2%/GPa5.2 磁振子频率的应变响应建立唯象模型ω(ε) ω0 α1ε α2ε^2拟合参数模式类型α1(meV/%)α2(meV/%²)主导机制声学支0.38-0.021交换作用调制光学支1.72-0.15各向异性变化实验发现1%拉伸应变可使磁振子带隙增加约15%这为自旋波器件设计提供了重要调控手段。6. 数值模拟方法细节6.1 随机LLG方程求解采用二阶Heun算法求解随机Landau-Lifshitz-Gilbert方程def heun_step(m, H_eff, dt): k1 -γ/(1α**2) * np.cross(m, H_eff) \ -α*γ/(1α**2) * np.cross(m, np.cross(m, H_eff)) m_pred m k1*dt H_pred compute_eff_field(m_pred) k2 -γ/(1α**2) * np.cross(m_pred, H_pred) \ -α*γ/(1α**2) * np.cross(m_pred, np.cross(m_pred, H_pred)) return m 0.5*(k1k2)*dt关键参数设置时间步长dt 0.1 ps保证能量误差0.1%阻尼系数α 0.05与FMR实验一致温度T 300 KNéel温度以上6.2 热涨落处理热场H_th满足H_th^i(t) 0 H_th^i(t)H_th^j(t) 2αkBTδijδ(t-t)/(γμ0MsV)模拟中采用Box-Muller算法生成高斯随机场。7. 实验验证与讨论7.1 非弹性中子散射数据对比我们重新分析了Park等人[64]的实验数据发现在Γ-K方向理论预测的28.5 meV光学模与实验25.3 meV峰位对应理论未出现的12 meV小峰可能来源于高阶交换作用磁振子-声子耦合样品缺陷7.2 应变调控实验建议基于理论预测建议开展以下实验压电力显微镜PFM测量施加0.1-0.5%应变观测磁振子频率移动Δω/ω≈1-5%时间分辨磁光测量泵浦-探测延迟扫描提取应变修饰的退相干时间T28. 技术应用展望8.1 自旋波器件设计指南利用Mn3Sn的应变敏感特性可设计可调滤波器中心频率应变调谐范围±15%插入损耗3 dB理论值非线性器件双二次交换导致强三波混频χ(2)≈10^5 V/m可实现参量放大8.2 优化建议材料优化掺Ge提高TN至420K纳米晶化降低矫顽场结构设计悬臂梁结构实现0.1%应变叉指电极激发相干磁振子通过系统研究三自旋模型的低能动力学我们不仅解释了Mn3Sn中磁振子谱的基本特征还揭示了应变调控的微观机制。这些发现为开发新型应变敏感自旋器件提供了理论基础。