Dy6(Fe,Mn)Bi2化合物磁性的第一性原理(First-Principle)计算研究

《Crystals》：First-Principle Calculations of Magnetic Properties of Dy6(Fe,Mn)Bi2 Compounds

　　

研究人员在密度泛函理论(Density Functional Theory, DFT)框架下研究了Dy6MnBi2和Dy6FeBi2金属间化合物的磁性。该类材料因潜在的磁制冷(Magnetocaloric Effect, MCE)应用及显著的磁热效应受到广泛关注。本研究计算了状态方程(Equation of State, EOS)、电子态密度(Electronic Density of States, DOS)及磁矩，并将计算结果与已有实验数据进行对比。计算所得物理量均与实验结果吻合良好，再次佐证了DFT作为探索三元稀土金属间化合物磁行为理论框架的可靠性。

Dy₆(Fe,Mn)Bi₂化合物磁性的第一性原理计算研究——发表于《Crystals》

一、研究背景与目的

三元稀土金属间化合物R₆TX₂（R为稀土元素，T为过渡金属，X = Bi, Te或Sb）因非常规磁行为及磁热效应(Magnetocaloric Effect, MCE)受到关注，部分成员具双磁相变及较大MCE。Dy₆(Fe,Mn)Bi₂系列中顺磁(PM)–铁磁(FM)转变温度约129–370 K，低温存在自旋重取向转变。实验中Fe被Mn替代后居里温度(Curie Temperature, T_C)升高，而按平均场理论磁离子间距增大应使T_C降低，此反常现象难以用常规理论解释，轨道杂化假说缺乏完备微观理论支撑。解析法困难，密度泛函理论(DFT)可定量预测结构、电子及磁性。研究人员此前已成功将DFT应用于Ho₆MnBi₂和Ho₆FeBi₂，本文扩展至Dy体系以评估DFT在不同稀土离子下的稳健性，并探究稀土替代对MCE性能的调控作用。本研究限制于共线自旋组态（铁磁相），未涉及低温非共线相。

二、主要关键技术方法

研究人员采用Quantum ESPRESSO软件包进行第一性原理计算。Dy使用投影缀加平面波(Projector-Augmented-Wave, PAW)赝势，Bi使用GBRV超软(Ultrosoft, US)赝势，Fe使用PSlibrary PAW赝势，Mn使用GBRV US赝势。交换关联能采用广义梯度近似(Generalized Gradient Approximation, GGA)—Perdew–Burke–Ernzerhof(PBE)及PBEsol泛函，最终选用与实验更符合的PBEsol。考虑自旋极化与非自旋极化、共线铁磁(FM)与反铁磁(AFM)组态。平面波截断能120 Ry，电荷密度截断能1440 Ry，布里渊区积分用Γ中心6×6×12 k点网格，Marzari–Vanderbilt冷 smear宽度0.075 Ry。晶胞及原子位置用BFGS准牛顿算法弛豫。状态方程(EOS)通过各向同性缩放体积(0.96–1.04V_eq)并在定容下弛豫内坐标获得，拟合Murnaghan EOS。总及分波态密度(PDOS)用四面体方法计算。

三、研究结果

4.1. Geometry Optimization and EOS（几何优化与状态方程）

研究人员以Herrero等及Morozkin等的实验结构为初猜进行全弛豫，优化晶格参数与原子位置。自旋极化计算的平衡体积与实验值吻合远优于非自旋极化情形，自旋极化态能量更低，证实磁性基态特征。Murnaghan EOS拟合精度高（偏差<10??%），平衡体积与全弛豫差异<0.1%。沿EOS曲线，a轴随体积线性增加；c/a比在平衡体积附近出现斜率变化（交叉现象），反映各向异性压缩率及可能的磁弹效应，Dy₆FeBi₂更显著（与Fe–Bi共价性强有关）。六方对称性得以保持，无结构相变。与Ho类似物相比，Dy体系c/a异常更明显，表明更强的磁弹耦合，可能影响层间Dy₁–Dy₂磁交换。

4.2. Geometry Optimization and Magnetic Moments（几何优化与磁矩）

在平衡体积附近，两个不等价Dy位（Dy1占3f位，Dy2占3g位）磁矩约4.6–4.7 μ_B（Dy1略大于Dy2），随体积增大轻微减小，归因于晶格膨胀改变传导电子与Dy 5d轨道重叠进而影响4f–5d交换极化，非直接3d–4f杂化所致。过渡金属(T=M n,Fe)位DFT给出非零磁矩（Mn > 3 μ_B，Fe ≈ 1.5–2.5 μ_B）且随体积增大而增大，但实验中过渡金属磁矩淬灭(quenched)，此为标准共线GGA忽略自旋–轨道耦合(Spin–Orbit Coupling, SOC)、高估Stoner交换劈裂及未考虑动态自旋涨落的已知局限。Ho类似物呈相同趋势。间接RE₁–RE₂交换由过渡金属与Bi介导，决定高有序温度。

4.3. Density of States（态密度分析）

总态密度自旋向上与向下通道明显不对称，确认自旋极化。两不等价Dy位多数自旋同向（铁磁排列），过渡金属与Dy铁磁耦合。分波态密度显示Dy 4f态为窄峰位于约?4 eV（高度局域化），过渡金属3d带分布于?3 eV至+2 eV，二者无显著能级交叠，排除强3d–4f杂化。Dy₆MnBi₂费米面(E_F)处态密度略高且不对称，表明更多巡游载流子，可能增强Dy离子间间接交换(RKKY-like mediation via itinerant electrons/Bi 6p)，解释Mn替代Fe后T_C升高的实验现象。Ho系列中RE–TM呈反铁磁耦合而Dy系列呈铁磁耦合，但因TM磁矩在DFT中为计算假象，真实物理由传导电子媒介的RE₁–RE₂间接交换主导。

四、结论（Discussion & Conclusions总结翻译）

研究人员通过DFT计算了Dy₆MnBi₂和Dy₆FeBi₂的平衡晶体结构、状态方程、电子态密度及磁矩并与实验对比。DFT对该体系给出一致可靠描述，优化结构再现Fe/Mn替代趋势，正确捕捉自旋极化态能量稳定性。与Ho类似物相比，Dy化合物平衡体积更大（符合镧系收缩），磁态稳定化更显著（尤Dy₆MnBi₂），c/a各向异性响应增强使Mn/Fe体系区分更清晰。这表明Ho替换为Dy改变层间与层内交换相互作用平衡，可能影响磁热性能。本研究强化了DFT对稀土基金属间化合物的预测能力，为未来系统探究稀土替代在调控磁热性质中的作用奠定基础。