反铁磁RAg（R = Gd, Tb, Dy）金属间化合物中Ce局域磁矩的温度依赖形成：基于泛函积分方法的中间价态作用

《Solid State Communications》：Temperature-dependent formation of local magnetic moments of Ce in antiferromagnetic RAg ( R = Gd, Tb, and Dy) intermetallics: The role of intermediate valence via functional integral approach

　　

对镧系金属及化合物中局域f电子与巡游价带电子间电荷转移所引发的磁性质的研究，持续吸引着理论和实验物理学家的关注。这一持续兴趣凸显了阐明这些材料中磁现象机制的必要性。反铁磁RAg金属间化合物（其中R为稀土元素，如钆Gd、铽Tb或镝Dy）展现出若干源于其晶体结构和磁相互作用的独特特征。这些化合物通常结晶为CsCl型立方结构，原子占据立方体角顶和中心，这种排列通过间接交换相互作用促进磁有序。在该相互作用中，磁离子通过传导电子而非直接原子轨道重叠发生相互作用。图1展示了RAg化合物的CsCl型立方晶格示意图，稀土原子（R）位于立方体角顶，银（Ag）原子位于立方体中心。R原子的磁矩在面内铁磁（平行）排列，但在面间反铁磁（反平行）排列，形成层状反铁磁有序。该结构以面内铁磁耦合与面间反平行耦合为特征，产生稳定的反铁磁态。低温下，这些材料展现出具有明确奈尔温度（T_N）的反铁磁有序，低于该温度则建立反铁磁序；高于T_N时，化合物呈现顺磁行为（无序磁矩）和负的顺磁居里温度，表明强烈的反铁磁关联。由金属s-p和d带中传导电子介导的鲁德曼-基特尔-卡苏亚-吉田（RKKY）间接交换相互作用主导了该行为，并依赖于电子密度和晶格间距。表1列出了GdAg、TbAg和DyAg主体的磁性质。这些性质使RAg化合物成为研究杂质诱导磁现象的理想体系，尤其对含中间价态元素（如Ce）的系统。其中，含Ce杂质的系统因Ce离子的中间价态行为而引人关注，该行为源于局域4f电子与传导态之间的杂化。Ce能级接近系统费米能级，因此仅需微小能量移动即可改变这些态的相对占据数，从而在局域磁矩与巡游电子行为之间建立平衡。因此，Ce的4f电子预期具有某种程度的巡游性。前期实验研究报道了稀释于RAg金属间化合物中Ce杂质的温度依赖磁超精细场，暗示Ce的中间价态与主体晶格磁有序之间存在关联。然而，能跨不同稀土主体和温度区间捕捉该行为的全面理论框架仍不完善。本研究通过采用基于泛函积分形式与相干势近似（CPA）相结合的理论模型来填补这一空白。第1节概述了理论模型，描述了单个Ce杂质位于反铁磁主体RAg（R = Gd, Tb, Dy）的R位时如何影响磁矩形成及其与磁超精细场的关系。该方法在自洽框架内处理Ce杂质，考虑Ce 4f与5d态之间作为温度函数的中间价态巡游性，并分析了从0 K到各自奈尔温度范围内局域磁矩与磁超精细场的演化。第2节展示了作为温度函数并在中间价态影响下的自洽数值结果，并将计算所得稀释于RAg中Ce杂质的磁超精细场与文献[8]报道的实验数据进行了比较。第3节给出了本研究的结论。

**论文解读文章**

**研究背景**
镧系金属及化合物中，局域4f电子与巡游价带电子间的电荷转移是磁性质的核心来源。反铁磁RAg金属间化合物（R = Gd, Tb, Dy）为研究该问题提供了理想平台：它们具有CsCl型立方晶格结构，在面内铁磁耦合、面间反铁磁耦合，形成稳定的层状反铁磁序，且其磁有序由Ruderman–Kittel–Kasuya–Yosida (RKKY) 间接交换相互作用主导。当引入铈（Ce）杂质时，因其4f能级接近费米能级，易于发生局域与巡游态间的杂化，导致中间价态行为，从而在局域磁矩与巡游电子行为间产生平衡。前期实验已观测到Ce杂质在RAg中的温度依赖磁超精细场，暗示中间价态与主体磁序相关，但缺乏捕捉该行为跨不同稀土主体和温度区间的全面理论框架。因此，研究人员旨在建立一种基于泛函积分形式与相干势近似（CPA）的理论模型，自洽描述Ce杂质在反铁磁主体中的局域磁矩及磁超精细场随温度的变化，并验证其与实验的一致性。

**研究内容与结论**
研究人员通过自洽理论计算，揭示了Ce杂质在反铁磁RAg主体中局域磁矩的温度依赖形成机制，并成功复现了实验观测的磁超精细场趋势。该研究的重要意义在于：首次从理论上系统阐明了中间价态在反铁磁环境中如何调控Ce杂质的局域磁响应，为理解类似稀土杂质体系提供了可推广的理论框架。论文发表在《Solid State Communications》。

**主要技术方法**
研究人员采用的理论模型基于泛函积分形式（静态近似）与相干势近似（CPA）相结合，用于处理电子关联与杂质效应。模型通过自洽求解含杂质位点的局域格林函数，并利用Friedel电荷屏蔽条件确定杂质能级，同时引入有效电荷差修正以考虑体积效应。关键参数包括d-f有效库仑相互作用U_df、d-f交换耦合参数J_df⁰及其反铁磁修正项δJ_df，以及费米-塞格雷接触耦合参数A(Z_imp)和核心极化常数A_cp^5d、A_cp^4f。所有计算均在0 K至各主体奈尔温度（T_N）范围内进行，温度依赖通过电子占据数和4f能级自洽体现。

**研究结果**
1. **Model and methodology**
研究人员构建了包含sp-d巡游电子（H_cd）、d-f局域相互作用（H_df）及杂质势（V）的哈密顿量，利用泛函积分将d-d电子关联映射为有效单电子哈密顿量中的涨落电荷（ν）和自旋（ξ）场，再通过CPA将系统替换为均匀自能有效介质。Ce杂质替代R位后，通过自洽求解局域格林函数和Friedel条件确定杂质能级，并考虑体积效应修正有效电荷差ΔZ′。局域磁矩由d电子贡献m_l0^γ（通过扰动态密度积分）和f电子贡献m₀^f（通过布里渊函数B_J(J_df)）组成，总磁超精细场B_hf来自6s、5d和4f电子的贡献。

2. **Results and discussions**
- **密度态与Ce价态**：在T=0 K下计算了主体密度态，并引入Ce杂质后修正了局域密度态。通过体积修正有效电荷差，自洽确定杂质能级。Ce的4f占据数（1-δ）随温度从0 K到T_N从0.75降至0.50，相应5d占据数变化（图2）。
- **局域磁矩的温度依赖**：计算了Ce在GdAg、TbAg、DyAg中的局域磁矩贡献（图3-5）。发现反铁磁相邻R位对m₀^f产生额外抑制，不同于典型布里渊函数行为，体现了RKKY相互作用在层状反铁磁相中的调制。d-f交换耦合的修正参数δJ_df/J_df⁰分别取0.75、0.77、0.80以拟合实验趋势。
- **磁超精细场对比**：利用自洽计算得到的磁矩，通过式(8)计算B_hf，并与实验数据对比（图6-8）。结果显示，Ce在GdAg、TbAg、DyAg中的磁超精细场随温度变化的理论曲线与实验点（分别用方形、圆形、三角形标记）吻合良好。各分量（B_hf^c、B_hf^d、B_hf^f）的贡献分析表明，温度依赖主要源于电子占据的重新分布和磁有序的减弱。

**总结与讨论**
研究人员在讨论中指出，尽管自旋-自旋关联及原子间自旋关联对反铁磁态的完整描述可能重要，但本工作采用的局域平均场框架中，主导贡献来自杂质与重稀土传导电子的有效电荷差及Ce杂质4f电荷转移，因此未显式包含这些关联，预期其影响为次要。通过自洽计算，模型成功捕捉了从0 K到T_N范围内局域磁矩和磁超精细场的形成与抑制，仅用少量可调参数即复现了实验趋势。研究结论翻译如下：
**结论**：在本工作中，研究人员采用基于泛函积分形式与CPA的中间价态模型，对反铁磁RAg化合物中Ce杂质局域磁矩的温度依赖形成进行了理论研究。尽管自旋-自旋相互作用和原子间自旋关联可能对反铁磁态的完整微观描述有贡献，但研究人员认为，这些关联对本文讨论的局域磁矩和磁超精细场影响不显著，如同前期工作[10,11]所述。这是因为本计算在局域平均场框架内进行，主导贡献来自杂质与重稀土传导电子间的有效电荷差以及Ce杂质的4f电荷转移，因此未显式包含这些关联。通过自洽计算，模型成功捕捉了从0 K到奈尔温度（T_N）范围内局域磁矩和磁超精细场的出现与抑制，仅用少量可调参数即实现了超精细场和局域磁矩形成的主要趋势。计算得到的Ce杂质在RAg主体（R = Gd, Tb, Dy）中磁超精细场的温度依赖性与实验观测[8]吻合良好，验证了模型的有效性。这些发现凸显了中间价态效应在理解稀释于稀土化合物中Ce杂质磁行为中的重要性。该模型为探索类似系统提供了理论框架，并清晰展示了中间价态如何调控反铁磁环境中Ce杂质的磁响应。特别地，结果强调了电子关联、温度与局域磁矩形成之间的相互作用。